<div dir="ltr">Thanks very much Martin. In case anyone is interested (or has further enlightening thoughts), the problem is in <br>gad_som_advect.F, where:<br><br> noFlowAcrossSurf = rigidLid .OR. nonlinFreeSurf.GE.1<br> & .OR. select_rStar.NE.0<div><br></div><div>yields <span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">noFlowAcrossSurf=.TRUE. even though <span style="text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">rigidLid=.FALSE., <span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">nonlinFreeSurf=0 and <span style="text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">select_rStar=0 in my setup. I have solved the problem for now by setting <span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">noFlowAcrossSurf=.FALSE. in my code directory, but obviously it would be nice if it could be fixed properly. I don't know why it's not working as expected though.</span></span></span></span></span></div><div><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Thanks, </span></span></span></span></span></div><div><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Spencer</span></span></span></span></span></div><div><p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures">
<br></span></p><p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><br></span></p><p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><br></span></p>
</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 24, 2018 at 3:59 AM, Martin Losch <span dir="ltr"><<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Spencer,<br>
<br>
in doc/tag-index I found that there were changes to the advection packages in checkpoint64e. Maybe you can diff the files (e.g. like this: <<a href="https://github.com/MITgcm/MITgcm/commit/d077f4a7023faea230432dc180400e4970aff0ff#diff-5cf10394edefefcf714f07bf62cf6d40" rel="noreferrer" target="_blank">https://github.com/MITgcm/<wbr>MITgcm/commit/<wbr>d077f4a7023faea230432dc180400e<wbr>4970aff0ff#diff-<wbr>5cf10394edefefcf714f07bf62cf6d<wbr>40</a>> or using other git capabilities) and figure out which change causes the in instability in your setup. <br>
<span class="HOEnZb"><font color="#888888"><br>
Martin<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
> On 21. Aug 2018, at 21:23, Spencer Jones <<a href="mailto:spencerj@ldeo.columbia.edu">spencerj@ldeo.columbia.edu</a>> wrote:<br>
> <br>
> Hi everyone, <br>
> I'm trying to update my code tree before I start a new project, but I'm finding that SOM is unstable in my setup from checkpoint 64e onwards. With the flux limiter off (advection scheme 80), the temperature becomes large and negative. With the flux limiter on (advection scheme 81), the temperature remains positive, but I get an unphysical-looking standing wave in the zonal direction at mid-latitudes. <br>
> <br>
> I have tried changing the viscosity, diffusivity and various other things but to no avail. Does anyone know what's going on? My data file is pasted below. <br>
> Thanks very much, <br>
> Spencer <br>
> <br>
> # ====================<br>
> # | Model parameters |<br>
> # ====================<br>
> #<br>
> # Continuous equation parameters<br>
> &PARM01<br>
> tRef = 32*0.,<br>
> sRef = 32*35.,<br>
> no_slip_bottom=.TRUE.,<br>
> no_slip_sides=.TRUE.,<br>
> viscAr=1.E-4,<br>
> viscAh=40.E3,<br>
> # diffKrT=1.E-5,<br>
> diffKrNrT=0.0083, 0.0065, 0.0038,0.0014, 0.00039,0.000092,<br>
> 0.00003, 0.00002,0.00002,0.00002,0.<wbr>00002,0.00002,<br>
> 0.00002,0.00002,0.00002,0.<wbr>00002,0.00002,0.00002,0.00002,<wbr>0.00002,<br>
> 0.00002,0.00002,0.00002,0.<wbr>00002,0.00002,0.00002,0.00002,<wbr>0.00002,<br>
> 0.00002,0.00002,0.00002,0.<wbr>00002,<br>
> # diffKrS=2.E-5,<br>
> diffKrNrS=0.0083, 0.0065, 0.0038,0.0014, 0.00039,0.000092,<br>
> 0.00003, 0.00002,0.00002,0.00002,0.<wbr>00002,0.00002,<br>
> 0.00002,0.00002,0.00002,0.<wbr>00002,0.00002,0.00002,0.00002,<wbr>0.00002,<br>
> 0.00002,0.00002,0.00002,0.<wbr>00002,0.00002,0.00002,0.00002,<wbr>0.00002,<br>
> 0.00002,0.00002,0.00002,0.<wbr>00002,<br>
> rhonil=1035.,<br>
> rhoConstFresh=1000.,<br>
> eosType = 'LINEAR',<br>
> ivdc_kappa=10.,<br>
> implicitDiffusion=.TRUE., <br>
> bottomDragLinear=1.E-3,<br>
> exactConserv=.TRUE.,<br>
> useRealFreshWaterFlux=.FALSE.,<br>
> useSingleCpuIO=.TRUE.,<br>
> # turn on looped cells<br>
> hFacMin=.1,<br>
> hFacMindr=50.,<br>
> # set precision of data files<br>
> readBinaryPrec=32,<br>
> tempAdvScheme=81,<br>
> saltAdvScheme=81,<br>
> staggerTimeStep=.TRUE.,<br>
> debugLevel=-1, <br>
> &<br>
> <br>
> # Elliptic solver parameters<br>
> &PARM02<br>
> cg2dMaxIters=500,<br>
> cg2dTargetResidual=1.E-13,<br>
> printResidualFreq=37606140,<br>
> &<br>
> <br>
> # Time stepping parameters<br>
> &PARM03<br>
> nIter0=13000000,<br>
> nTimeSteps = 4800000,<br>
> # 100 years of integration will yield a reasonable flow field<br>
> # startTime = 0.,<br>
> # endTime = 3110400000.,<br>
> deltaTmom = 900.,<br>
> deltaTtracer= 10800.,<br>
> deltaTClock = 10800.,<br>
> deltaTfreesurf= 5400.,<br>
> abEps = 0.1,<br>
> pChkptFreq= 4320000000.,<br>
> chkptFreq= 4320000000.,<br>
> dumpFreq= 43200000000000.,<br>
> monitorFreq=86400000.,<br>
> # 10 days restoring timescale for temperature<br>
> tauThetaClimRelax= 864000.,<br>
> # 6 months restoring timescale for salinity<br>
> # tauSaltClimRelax = 15552000.,<br>
> # periodicExternalForcing=.<wbr>FALSE.,<br>
> # externForcingPeriod=2592000.,<br>
> # externForcingCycle=31104000.,<br>
> pickupStrictlyMatch=.FALSE., <br>
> &<br>
> <br>
> # Gridding parameters<br>
> &PARM04<br>
> usingSphericalPolarGrid=.TRUE.<wbr>,<br>
> delR= 32*125.,<br>
> ygOrigin=-70.,<br>
> dySpacing=0.9722,<br>
> dxSpacing=0.9722,<br>
> &<br>
> <br>
> # Input datasets<br>
> &PARM05<br>
> bathyFile= 'bathy.bin',<br>
> zonalWindFile= 'taux.bin',<br>
> meridWindFile= 'tauy.bin',<br>
> thetaClimFile= 'sst2.bin',<br>
> # saltClimFile= 'lev_sss.bin',<br>
> # surfQFile= 'ncep_qnet.bin',<br>
> # the_run_name= 'global_oce_latlon',<br>
> # fresh water flux is turned on, comment next line to it turn off<br>
> # (maybe better with surface salinity restoring)<br>
> EmPmRFile= 'sflux.bin',<br>
> &<br>
> <br>
> <br>
</div></div><div class="HOEnZb"><div class="h5">> ______________________________<wbr>_________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
<br>
______________________________<wbr>_________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
</div></div></blockquote></div><br></div>