<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">The difference between the “truly” incompressible and the semi-compressible Boussinesq equations is all in the equation of state. In the “truly” impressible equations, density does not depend on pressure at all while the semi-compressible equations allow the density to depend on the hydrostatic pressure due to the reference density. This is the standard set used in oceanography and most people just call these the Boussinesq or incompressible equations since the compressibility doesn’t affect dynamics, just thermodynamics.<div class=""><br class=""></div><div class="">The ISOMIP reference experiment uses the JMD95Z equation of state, which depends on the above hydrostatic pressure. So the ISOMIP experiment would be using the semi-compressible equations.</div><div class=""><br class=""></div><div class="">Christopher</div><div class=""><br class=""><div class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">----------------------------------------------------------------------------<br class="">Christopher L. Pitt Wolfe<br class="">Associate Professor (Physical Oceanography)<br class="">School of Marine and Atmospheric Sciences<br class="">Stony Brook University<br class=""><a href="mailto:christopher.wolfe@stonybrook.edu" class="">christopher.wolfe@stonybrook.edu</a>             631-632-3152<br class="">----------------------------------------------------------------------------<br class=""><br class=""></div></div></div>
</div>
<div><br class=""><blockquote type="cite" class=""><div class="">On Nov 4, 2021, at 6:11 AM, Martin Losch <<a href="mailto:Martin.Losch@awi.de" class="">Martin.Losch@awi.de</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Arran,<br class=""><br class="">the isomip verification experiment uses hydrostatic (because nonHydrostatic=.FALSE. in data, defining ALLOW_NONHYDROSTATIC only compiles the code) and incompressible equations.<br class=""><br class="">In general, I have always thought the the non-hydrostatic code always uses incompressible equations.<br class=""><br class="">I am not entirely sure about the documentation, but 1.5.1.4 (where you find eq 1.98-1.103) also contains the statement that “it is the hydrostatic () form of these equations that are used throughout the ocean modelling community …”<br class="">I agree that this is a little misleading. But note, that 1.100 (non-divergent velocity field) implies incompressibility, so I am not quite sure what the section heading means here, maybe someone else can explain.<br class=""><br class="">Martin<br class=""><br class=""><blockquote type="cite" class="">On 3. Nov 2021, at 21:11, Arran Whiteford <<a href="mailto:arran.whiteford@vuw.ac.nz" class="">arran.whiteford@vuw.ac.nz</a>> wrote:<br class=""><br class="">Hullo,<br class=""><br class="">I'm trying to figure out what versions of the Boussinesq equations the ISOMIP tutorial example is solving <br class=""><a href="https://github.com/MITgcm/MITgcm/tree/master/verification/isomip" class="">https://github.com/MITgcm/MITgcm/tree/master/verification/isomip</a><br class=""><br class="">It includes non hydrostatic code.<br class="">#define ALLOW_NONHYDROSTATIC https://github.com/MITgcm/MITgcm/blob/master/verification/isomip/code/CPP_OPTIONS.h<br class=""><br class="">In the MITgcm documents it explains:<br class="">1.3.4.4.2. Ocean<br class="">Non-hydrostatic forms of the incompressible Boussinesq equations in z- coordinates are supported - see eqs. (1.98) to (1.103).<br class=""><br class="">BUT eqns (1.98) to (1.103) are not the incompressible eqns, they are the semi-compressible.<br class="">https://mitgcm.readthedocs.io/en/latest/overview/hydrostatic.html<br class=""><br class="">In the non-hydrostatic ocean model all terms in equations Eqs.(1.29 →<br class="">1.31) are retained.<br class=""><br class="">Which is true? Is the nonhydrostatic code running incompressible or semi-compressible Boussinesq equations?<br class=""><br class="">Thanks,<br class="">Arran<br class=""><br class="">_______________________________________________<br class="">MITgcm-support mailing list<br class="">MITgcm-support@mitgcm.org<br class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br class=""></blockquote><br class="">_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br class=""></div></div></blockquote></div><br class=""></div></body></html>