<div dir="ltr">Martin,<div>I tried a few more things.</div><div>1) When I use the EXF package, but turn off the OBCS package (and comment out relevant parts in the data.exf file), then the solver seems to be fine.</div><div>2) When I use both the EXF and OBCS package, and also turn on the SEAICE package, the solver also seems to be fine, but I get:</div><div>** WARNING ** SEAICE_LSR (it= 0, 1) did not converge :<br>(PID.TID 0808.0001) nIt,wFU,dU= 1500 0.950 2.126087E-09 ; nIt,wFV,dV= 1500 0.950 1.679611E-09<br>(PID.TID 0808.0001) ** WARNING ** SEAICE_LSR (it= 0, 2) did not converge :<br>(PID.TID 0808.0001) nIt,wFU,dU= 1500 0.950 3.301846E-09 ; nIt,wFV,dV= 1500 0.950 5.170461E-09<br>(PID.TID 0808.0001) ** WARNING ** SEAICE_LSR (it= 1, 1) did not converge :<br>(PID.TID 0808.0001) nIt,wFU,dU= 1500 0.950 4.122604E-09 ; nIt,wFV,dV= 1500 0.950 1.121440E-08<br>(PID.TID 0808.0001) ** WARNING ** SEAICE_LSR (it= 1, 2) did not converge :<br>(PID.TID 0808.0001) nIt,wFU,dU= 1500 0.950 4.035429E-09 ; nIt,wFV,dV= 1500 0.950 9.547553E-09</div><div>Also, when I run this for longer, I get a very unrealistic ocean, with deep mixed layers and surface cooling (the restoring seems to try and compensate this, but not enough) where they shouldn't be (also far away from any ice) etc, so while the solver seems to be fine in these runs, there's still something very wrong. </div><div>So it seems to be that there is something wrong with how I use the OBCS package, but since my OBCS input files work fine without the exf package, I assume those files are fine.</div><div>Any ideas?</div><div>Andreas</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jan 23, 2021 at 12:52 PM Andreas Klocker <<a href="mailto:aklocker42@gmail.com">aklocker42@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Martin,<div><br></div><div>- # useExfCheckRange = .TRUE.,</div><div><br></div><div>I don't get any complaints here.</div><div><br></div><div>- 32 bit</div><div><br></div><div>Yes, I use 32 bit everywhere.</div><div><br></div><div>- readStressOnAgrid</div><div><br></div><div>Correct, that was a mistake. But it didn't solve my problem.</div><div><br></div><div>This should be so simple, hence my confusion...</div><div><br></div><div>Andreas</div><div><br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jan 23, 2021 at 12:26 PM Martin Losch <<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Andreas,<br>
<br>
when you uncomment<br>
# useExfCheckRange = .TRUE.,<br>
does the model complain about weird forcing values?<br>
<br>
Wild guess: you use<br>
> exf_iprec = 32,<br>
> exf_iprec_obcs = 32,<br>
are your input fields really 32 bit (and you used readBinaryPrec=32 in data&PARM01)?<br>
<br>
Then, if you use the same files, readStressOnAgrid should be false, shouldn’t it? without exf the model expects the forcing to be on the appropriate grid points, ie. U and V points for stress.<br>
<br>
Not important but I think you don’t need to set ustressRepCycle etc to the same value if you set repeatPeriod=31104000.<br>
<br>
Martin<br>
<br>
> On 22. Jan 2021, at 08:02, Andreas Klocker <<a href="mailto:aklocker42@gmail.com" target="_blank">aklocker42@gmail.com</a>> wrote:<br>
> <br>
> Hi all,<br>
> <br>
> I'm having some issues with the EXF package which I just can't understand. Maybe one of you can help.<br>
> <br>
> I am running a regional model configuration using the OBCS package. In my first experiment I set<br>
> <br>
> periodicExternalForcing=.TRUE.,<br>
> externForcingPeriod=2592000.,<br>
> externForcingCycle=31104000.,<br>
> tauThetaClimRelax = 864000.,<br>
> tauSaltClimRelax = 2592000.,<br>
> <br>
> in the data file, i.e. I restore surface T/S with above restoring time scale/forcing period/forcing cycle. This works fine. <br>
> <br>
> Then I wanted to add sea ice to my open boundary conditions which, if I understand correctly, means that I have to use the EXF package to set the restoring time scale/forcing period/forcing cycle. To move forward in baby steps, I didn't use sea ice yet, but just tried to let the EXF package do exactly the same as the restoring and obcs did before I switched to the EXF package. See below for my data.exf file. But when I try and run this, the solver immediately writes NaNs, despite me using the same initial conditions, obcs files, and restoring files as before. I would have thought that both ways of using OBCS and restoring should give the same results, but one runs fine while the other one crashes immediately without me being able to find the issue.<br>
> <br>
> Please let me know if anyone is seeing something I'm doing wrong here.<br>
> <br>
> Thanks in advance!<br>
> <br>
> Andreas Klocker<br>
> <br>
> <br>
> #<br>
> # *********************<br>
> # External Forcing Data<br>
> # *********************<br>
> &EXF_NML_01<br>
> #<br>
> # useExfCheckRange = .TRUE.,<br>
> # repeatPeriod = 2592000.0,<br>
> exf_iprec = 32,<br>
> exf_iprec_obcs = 32,<br>
> useAtmWind=.FALSE.,<br>
> readStressOnAgrid=.TRUE.,<br>
> #<br>
> &<br>
> <br>
> # *********************<br>
> &EXF_NML_02<br>
> #<br>
> ustressstartdate1 = 19790101,<br>
> ustressstartdate2 = 000000,<br>
> ustressperiod = 2592000.0,<br>
> ustressRepCycle = 31104000.,<br>
> #<br>
> vstressstartdate1 = 19790101,<br>
> vstressstartdate2 = 000000,<br>
> vstressperiod = 2592000.0,<br>
> vstressRepCycle = 31104000.,<br>
> #<br>
> climsststartdate1 = 19790101,<br>
> climsststartdate2 = 000000,<br>
> climsstperiod = 2592000.0,<br>
> climsstTauRelax = 2592000.0,<br>
> climsstRepCycle = 31104000.,<br>
> #<br>
> climsssstartdate1 = 19790101,<br>
> climsssstartdate2 = 00000,<br>
> climsssperiod = 2592000.0,<br>
> climsssTauRelax = 864000.0,<br>
> climsssRepCycle = 31104000.,<br>
> #<br>
> ustressfile = 'TauX_BSOSE ',<br>
> vstressfile = 'TauY_BSOSE ',<br>
> climsstfile = 'SST_BSOSE.res ',<br>
> climsssfile = 'SSS_BSOSE.res',<br>
> #<br>
> &<br>
> <br>
> # *********************<br>
> &EXF_NML_03<br>
> &<br>
> <br>
> # *********************<br>
> # old open64 compiler (4.2.1) cannot skip this namelist to read in the next one;<br>
> # comment out this namelist (not read).<br>
> #&EXF_NML_04<br>
> #&<br>
> <br>
> # *********************<br>
> &EXF_NML_OBCS<br>
> <br>
> # useOBCSYearlyFields = .TRUE.,<br>
> <br>
> obcsSstartdate1 = 19790101,<br>
> obcsSstartdate2 = 000000,<br>
> obcsSperiod = 2592000.0,<br>
> obcsSRepCycle = 31104000.,<br>
> #<br>
> obcsNstartdate1 = 19790101,<br>
> obcsNstartdate2 = 000000,<br>
> obcsNperiod = 2592000.0,<br>
> obcsNRepCycle = 31104000.,<br>
> #<br>
> obcsWstartdate1 = 19790101,<br>
> obcsWstartdate2 = 000000,<br>
> obcsWperiod = 2592000.0,<br>
> obcsWRepCycle = 31104000.,<br>
> #<br>
> obcsEstartdate1 = 19790101,<br>
> obcsEstartdate2 = 000000,<br>
> obcsEperiod = 2592000.0,<br>
> obcsERepCycle = 31104000.,<br>
> #<br>
> &<br>
> <br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>
</blockquote></div>