<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="">Hi Rita,<div class=""><br class=""></div><div class="">if you do this, you should keep in mind that before the sea ice model is called, Qnet = Qnet_ocean = lwdown + remaining terms.</div><div class="">After the sea ice model has been called Qnet = Qnet_ocean*(1-AREA) + Qnet_ice*AREA (Area = seaice fractional area \elem [0,1]), and lwdown is not used anymore.</div><div class="">If you don’t care about the sea ice model, you could add any heat anomaly dQ to Qnet **after** calling the sea ice model (in do_oceanic_phys), but this will likely change the surface temperature so that is inconsistent with the sea ice model.</div><div class="">If you add your anomaly before calling the sea ice model (ie. somewhere to the fields of the exf pkg), this would be similar as adding these heat anomalies like this: Qnet + dQ*(1-AREA), after the seaice, i.e only over the open ocean part (except that you dQ affects the sea ice model, too). So I’d add dQ*(1-AREA) to Qnet after calling the sea ice model to be on the safe side.</div><div class=""><br class=""></div><div class="">Martin</div><div class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 8. Apr 2022, at 12:37, Margarita Markina <<a href="mailto:margarita.markina@earth.ox.ac.uk" class="">margarita.markina@earth.ox.ac.uk</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta charset="UTF-8" class=""><div style="font-style: normal; font-variant-caps: normal; font-weight: 400; 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; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);" class="">Thanks again Dimitris and Martin - all suggestions are very useful!</div><div style="font-style: normal; font-variant-caps: normal; font-weight: 400; 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; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: 400; 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; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);" class="">Since my aim is to investigate oceanic response to perturbations in <span style="background-color: rgb(255, 255, 255); display: inline !important;" class="">heat flux that we want to apply in certain locations - we were</span><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""> thinking that another way to do it (without turning off ice model so we have all ice feedbacks) - would be to stick with using bulk formulae - but add desired localized anomalies to downward longwave radiation flux? </span><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">since the Qnet is anyways a sum of DLW, DSW and latent and sensible heat fluxes. Do you think it could be a viable option?</span></div><div style="font-style: normal; font-variant-caps: normal; font-weight: 400; 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; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);" class=""><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></span></div><div style="font-style: normal; font-variant-caps: normal; font-weight: 400; 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; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: 400; 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; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);" class="">Thanks again, </div><div style="font-style: normal; font-variant-caps: normal; font-weight: 400; 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; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);" class="">Rita</div><div id="appendonsend" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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;" class=""></div><hr tabindex="-1" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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; display: inline-block; width: 830.046875px;" class=""><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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; float: none; display: inline !important;" class=""></span><div id="divRplyFwdMsg" dir="ltr" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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;" class=""><font face="Calibri, sans-serif" style="font-size: 11pt;" class=""><b class="">From:</b><span class="Apple-converted-space"> </span>MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" class="">mitgcm-support-bounces@mitgcm.org</a>> on behalf of Menemenlis, Dimitris (US 329B) <<a href="mailto:dimitris.menemenlis@jpl.nasa.gov" class="">dimitris.menemenlis@jpl.nasa.gov</a>><br class=""><b class="">Sent:</b><span class="Apple-converted-space"> </span>Thursday, April 7, 2022 6:46 PM<br class=""><b class="">To:</b><span class="Apple-converted-space"> </span>MITgcm Support <<a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a>><br class=""><b class="">Subject:</b><span class="Apple-converted-space"> </span>Re: [MITgcm-support] [EXTERNAL] forcing model with hflux and sflux</font><div class=""> </div></div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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; line-break: after-white-space;">Hi Rita, although hflux and sflux are not compatible with pkg/seaice or pkg/thsice, there is a super-simple MITgcm sea ice model that is:<div class=""><br class=""></div><div class="">You enable it with option “allowFreezing=.TRUE.” under namelist PARM01 in runtime parameter file “data"</div><div class=""><br class=""></div><div class=""><div class="">C     allowFreezing  :: Allows surface water to freeze and form ice</div><div class=""><br class=""></div><div class=""><div class="">C     adjustColdSST_diag :: diagnostic field for how much too cold (below</div><div class="">C              Tfreezing) SST has been adjusted (with allowFreezing=T).</div><div class="">C              > 0 for increase of SST (up to Tfreezing).</div><div class="">C              Units are r_unit.K/s (=Kelvin.m/s if r=z).</div><div class="">C        Note: 1) allowFreezing option is a crude hack to fix too cold SST that</div><div class="">C              results from missing seaice component. It should never be used</div><div class="">C              with any seaice component, neither current seaice pkg (pkg/seaice</div><div class="">C              or pkg/thsice) nor a seaice component from atmos model when</div><div class="">C              coupled to it.</div><div class=""><br class=""></div></div><div class=""><br class=""><blockquote type="cite" class=""><div class="">On Apr 5, 2022, at 3:17 AM, Martin Losch <<a href="mailto:Martin.Losch@awi.de" class="">Martin.Losch@awi.de</a>> wrote:</div><br class="x_Apple-interchange-newline"><div class=""><div class="" style="word-wrap: break-word; line-break: after-white-space;">Hi Rita,<div class=""><br class=""></div><div class="">the code cannot apply hflux/sflux with the sea ice model turned on, because the sea ice model computes these fluxes where there is ice. Frankly, what do you need the sea ice model for, if you want to prescribe the fluxes anyway? In this sense I’d support Dimitris’ suggestion. The fluxes that you prescribe should know about the sea ice distribution. It’s like prescribing the sea ice cover. If you want to have the sea ice modify these fluxes, more work is needed (and it will be complicated to get it “right”, i.e. consistent).<div class=""><br class=""></div><div class="">You’ll have to find a solution where you allow the model to compute the fluxes (from atmospheric fields that it reads) **and** read them. Depending what you want to do exactly I would have a look at pkg/exf and pkg/seaice and comment out the stop statements that bother you. Then you’d have to prescribe both the fluxes and the atmospheric fields and modify the code so that the fluxes are not overwritten by the ones that the exf-package computes. The sea ice model will then use the atmospheric fluxes (which you have prescribed) and recompute them in grid cells with sea ice. That could work, but I have never tried it.</div><div class=""><br class=""></div><div class="">Martin</div><div class=""><div class=""><br class=""><blockquote type="cite" class=""><div class="">On 5. Apr 2022, at 12:03, Margarita Markina <<a href="mailto:margarita.markina@earth.ox.ac.uk" class="">margarita.markina@earth.ox.ac.uk</a>> wrote:</div><br class="x_Apple-interchange-newline"><div class=""><div class="" style="font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);">It is subpolar latitudes, so I need the sea ice unfortunately<span class="x_Apple-converted-space"> </span><span id="x_🙁" class="">🙁</span></div><div class="" style="font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);"><br class=""></div><div class="" style="font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);">is there a way to force the model with hflux/sflux when sea ice module is on?</div><div class="" style="font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);"><br class=""></div><div class="" style="font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);">Or alternatively (if the former is impossible) - is there a way to add some fixed values (in certain locations) to heat fluxes computed within the model? </div><div class="" style="font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);"><br class=""></div><div class="" style="font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);"><br class=""></div><div class="" style="font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);">All the best, </div><div class="" style="font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);">Rita</div><div id="x_appendonsend" class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"></div><hr tabindex="-1" class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; display: inline-block; width: 830.046875px;"><span class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; float: none; display: inline !important;"></span><div id="x_divRplyFwdMsg" dir="ltr" class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><font face="Calibri, sans-serif" class="" style="font-size: 11pt;"><b class="">From:</b><span class="x_Apple-converted-space"> </span>MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" class="">mitgcm-support-bounces@mitgcm.org</a>> on behalf of Menemenlis, Dimitris (US 329B) <<a href="mailto:dimitris.menemenlis@jpl.nasa.gov" class="">dimitris.menemenlis@jpl.nasa.gov</a>><br class=""><b class="">Sent:</b><span class="x_Apple-converted-space"> </span>Tuesday, April 5, 2022 10:49 AM<br class=""><b class="">To:</b><span class="x_Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a><span class="Apple-converted-space"> </span><<a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a>><br class=""><b class="">Subject:</b><span class="x_Apple-converted-space"> </span>Re: [MITgcm-support] [EXTERNAL] forcing model with hflux and sflux</font><div class=""> </div></div><div class="x_BodyFragment" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><font size="2" class=""><span class="" style="font-size: 11pt;"><div class="x_PlainText">Turn off the sea ice module ;-)<br class=""><br class="">> On Apr 5, 2022, at 2:43 AM, Margarita Markina <<a href="mailto:margarita.markina@earth.ox.ac.uk" class="">margarita.markina@earth.ox.ac.uk</a>> wrote:<br class="">><span class="x_Apple-converted-space"> </span><br class="">> So, my question is -  how do I set my model so it is forced by hflux and sflux without getting errors from sea ice module?<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=""><a href="https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!eTeWcOxOTAji-EIjGhQHvYQArFlltbqlfjJG7zzui2pWTp-kGH19A5VFBt_9BuT7KbhV_vmX4bw$" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class=""></div></span></font></div><span class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; float: none; display: inline !important;">_______________________________________________</span><br class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><span class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; float: none; display: inline !important;">MITgcm-support mailing list</span><br class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><a href="mailto:MITgcm-support@mitgcm.org" class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;">MITgcm-support@mitgcm.org</a><br class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><a href="https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!eTeWcOxOTAji-EIjGhQHvYQArFlltbqlfjJG7zzui2pWTp-kGH19A5VFBt_9BuT7KbhV_vmX4bw$" class="" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a></div></blockquote></div><br class=""></div></div></div>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class=""><a href="https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!eTeWcOxOTAji-EIjGhQHvYQArFlltbqlfjJG7zzui2pWTp-kGH19A5VFBt_9BuT7KbhV_vmX4bw$" class="">https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!eTeWcOxOTAji-EIjGhQHvYQArFlltbqlfjJG7zzui2pWTp-kGH19A5VFBt_9BuT7KbhV_vmX4bw$</a><span class="Apple-converted-space"> </span><br class=""></div></blockquote></div><br class=""></div></div><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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; float: none; display: inline !important;" class="">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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;" class=""><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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; float: none; display: inline !important;" class="">MITgcm-support mailing list</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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;" class=""><a href="mailto:MITgcm-support@mitgcm.org" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">MITgcm-support@mitgcm.org</a><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; 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;" class=""><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a></div></blockquote></div><br class=""></div></body></html>