<div dir="ltr">That's a good point Oliver which I haven't thought about! I'll try to fix these timing differences and hopefully it'll work...<div>Andreas</div></div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jan 23, 2021 at 3:06 PM Oliver Jahn <<a href="mailto:jahn@mit.edu" target="_blank">jahn@mit.edu</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>
I think you will have to place the fields in the middle of the month as<br>
non-exf forcing assumes time averages,<br>
<br>
obcsSstartdate1 = 19790116,<br>
<br>
etc. Alternatively, you can just set<br>
<br>
obcsSperiod = -12,<br>
<br>
and not set startdate1/2 or RepCycle. This means monthly climatology<br>
(even with Gregorian calendar). This applies to all fields, of course,<br>
but I imagine for obcs a mismatch with the initial condition could make<br>
the model blow up early on. Not sure this is the only problem, but it<br>
might help.<br>
<br>
Best,<br>
Oliver<br>
<br>
<br>
On 2021-01-22 22:44, Andreas Klocker wrote:<br>
> Martin,<br>
> I tried a few more things.<br>
> 1) When I use the EXF package, but turn off the OBCS package (and<br>
> comment out relevant parts in the data.exf file), then the solver seems<br>
> to be fine.<br>
> 2) When I use both the EXF and OBCS package, and also turn on the SEAICE<br>
> package, the solver also seems to be fine, but I get:<br>
> ** 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=<br>
> 1500 0.950 1.679611E-09<br>
> (PID.TID 0808.0001) ** WARNING ** SEAICE_LSR (it= 0, 2) did<br>
> not converge :<br>
> (PID.TID 0808.0001) nIt,wFU,dU= 1500 0.950 3.301846E-09 ; nIt,wFV,dV=<br>
> 1500 0.950 5.170461E-09<br>
> (PID.TID 0808.0001) ** WARNING ** SEAICE_LSR (it= 1, 1) did<br>
> not converge :<br>
> (PID.TID 0808.0001) nIt,wFU,dU= 1500 0.950 4.122604E-09 ; nIt,wFV,dV=<br>
> 1500 0.950 1.121440E-08<br>
> (PID.TID 0808.0001) ** WARNING ** SEAICE_LSR (it= 1, 2) did<br>
> not converge :<br>
> (PID.TID 0808.0001) nIt,wFU,dU= 1500 0.950 4.035429E-09 ; nIt,wFV,dV=<br>
> 1500 0.950 9.547553E-09<br>
> Also, when I run this for longer, I get a very unrealistic ocean, with<br>
> deep mixed layers and surface cooling (the restoring seems to try and<br>
> compensate this, but not enough) where they shouldn't be (also far away<br>
> from any ice) etc, so while the solver seems to be fine in these runs,<br>
> there's still something very wrong. <br>
> So it seems to be that there is something wrong with how I use the OBCS<br>
> package, but since my OBCS input files work fine without the exf<br>
> package, I assume those files are fine.<br>
> Any ideas?<br>
> Andreas<br>
> <br>
> <br>
> <br>
> On Sat, Jan 23, 2021 at 12:52 PM Andreas Klocker <<a href="mailto:aklocker42@gmail.com" target="_blank">aklocker42@gmail.com</a><br>
> <mailto:<a href="mailto:aklocker42@gmail.com" target="_blank">aklocker42@gmail.com</a>>> wrote:<br>
> <br>
> Hi Martin,<br>
> <br>
> - # useExfCheckRange = .TRUE.,<br>
> <br>
> I don't get any complaints here.<br>
> <br>
> - 32 bit<br>
> <br>
> Yes, I use 32 bit everywhere.<br>
> <br>
> - readStressOnAgrid<br>
> <br>
> Correct, that was a mistake. But it didn't solve my problem.<br>
> <br>
> This should be so simple, hence my confusion...<br>
> <br>
> Andreas<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> 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><br>
> <mailto:<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>>> wrote:<br>
> <br>
> 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<br>
> readBinaryPrec=32 in data&PARM01)?<br>
> <br>
> Then, if you use the same files, readStressOnAgrid should be<br>
> false, shouldn’t it? without exf the model expects the forcing<br>
> 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<br>
> etc to the same value if you set repeatPeriod=31104000.<br>
> <br>
> Martin<br>
> <br>
> > On 22. Jan 2021, at 08:02, Andreas Klocker<br>
> <<a href="mailto:aklocker42@gmail.com" target="_blank">aklocker42@gmail.com</a> <mailto:<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<br>
> understand. Maybe one of you can help.<br>
> ><br>
> > I am running a regional model configuration using the OBCS<br>
> 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<br>
> restoring time scale/forcing period/forcing cycle. This works fine.<br>
> ><br>
> > Then I wanted to add sea ice to my open boundary conditions<br>
> which, if I understand correctly, means that I have to use the<br>
> EXF package to set the restoring time scale/forcing<br>
> period/forcing cycle. To move forward in baby steps, I didn't<br>
> use sea ice yet, but just tried to let the EXF package do<br>
> exactly the same as the restoring and obcs did before I switched<br>
> to the EXF package. See below for my data.exf file. But when I<br>
> try and run this, the solver immediately writes NaNs, despite me<br>
> using the same initial conditions, obcs files, and restoring<br>
> files as before. I would have thought that both ways of using<br>
> OBCS and restoring should give the same results, but one runs<br>
> fine while the other one crashes immediately without me being<br>
> able to find the issue.<br>
> ><br>
> > Please let me know if anyone is seeing something I'm doing<br>
> 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<br>
> 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> <mailto:<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> <mailto:<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>
> _______________________________________________<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>