[MITgcm-support] Can not run a regional model for long time
Matthew Mazloff
mmazloff at ucsd.edu
Sun Sep 18 12:16:57 EDT 2022
You need to find out why it is becoming unstable. Have a look at the solution and see if you can see where and why it is blowing up.
Matt
> On Sep 18, 2022, at 6:08 AM, Sabri MUTLU (MAM-ÇTÜE) <sabri.mutlu at tubitak.gov.tr> wrote:
>
> Hi all,
>
> As a newbie in this field, I need help for the following problem, which I have been trying to deal with for a long time.
>
> I built a regional model for my thesis study. When I try to run it for 12 or more days, it gives "NaN" values after 8 or 9 days later in the model.
> Could you please examine the following model configuration and give me an idea for possible solution?
>
> (I had to choose Orlanski method since don't have any boundary condition to feed the model.)
>
> Thank you in advance.
> Sabri
>
> &PARM01
> viscAh=2.E-2,
> viscAz=1.5E-5,
> no_slip_sides=.TRUE.,
> no_slip_bottom=.TRUE.,
> diffKhT=2.E-2,
> diffKhS=2.E-2,
> diffKzT=2.E-5,
> diffKzS=2.E-5,
> f0=0.948E-4,
> implicitFreeSurface=.TRUE.,
> exactConserv=.TRUE.,
> implicitDiffusion=.TRUE.,
> implicitViscosity=.TRUE.,
> eosType = 'JMD95Z',
> tRef=22.2803,22.3383,22.3646,22.4220,22.5469,22.5081,21.0406,19.6075,18.3619,16.7902,15.8310,15.4016,15.3625,15.3935,15.4334,15.4465,15.4508,15.4710,15.5226,
> 15.5428,15.5485,15.5589,15.5829,15.5860,15.5859,15.5767,15.5572,15.5477,15.5357,15.5068,15.4782,15.4640,15.4337,15.3816,15.3480,15.2998,15.2423,15.2127,
> 15.1271,15.0639,15.0086,14.9722,14.9229,14.8937,14.8470,14.8073,14.7838,14.7677,14.7392,14.7210,14.7037,14.7036,14.7035,14.7034,14.7033,
> sRef=24.3775,24.4042,24.4730,24.5782,24.7042,24.8407,25.8279,26.4026,28.0163,30.1542,31.7531,33.9951,37.0132,37.5629,37.8303,37.9592,38.0011,38.1433,38.2744,
> 38.3184,38.3548,38.4269,38.5535,38.6043,38.6738,38.7411,38.7479,38.7496,38.7510,38.7512,38.7469,38.7529,38.7569,38.7644,38.7666,38.7683,38.7698,38.7696,
> 38.7635,38.7563,38.7513,38.7476,38.7427,38.7399,38.7354,38.7326,38.7305,38.7280,38.7254,38.7236,38.7220,38.7220,38.7220,38.7220,38.7220,
> staggerTimeStep = .TRUE.,
> tempAdvScheme = 33,
> saltAdvScheme = 33,
> useSingleCpuIO=.TRUE.,
> bottomDragQuadratic = 0.02,
> select_rStar = 0,
> &
>
> &PARM02
> cg2dMaxIters=300,
> cg2dTargetResidual=1.E-13,
> cg3dMaxIters=20,
> cg3dTargetResidual=1.E-8,
> &
> ---------
> &PARM03
> startTime = 0.,
> endTime = 1036800.,
> deltaT=5,
> dumpFreq=10800.,
> monitorFreq=3600.,
> monitorSelect=2,
> periodicExternalForcing=.TRUE.,
> externForcingPeriod=10800.,
> externForcingCycle=1036800.,
> &
>
> &PARM04
> usingSphericalPolarGrid=.TRUE.,
> delz = 32*2,22*8,10,
> delX = 528*0.0013,
> delY = 164*0.001,
> ygOrigin=40.656,
> xgOrigin=29.2612,
> &
>
> &PARM05
> bathyFile='bathy_v6_528-164.bin',
> zonalWindFile='windU_2309-0510.bin',
> meridWindFile='windV_2309-0510.bin',
> &
>
> # Open-boundaries
> &OBCS_PARM01
> OB_Iwest=164*1,
> useOrlanskiWest=.TRUE.,
> useOBCSbalance=.TRUE.,
> useOBCSsponge=.FALSE.,
> OBCS_balanceFacW = -1.,
> &
>
> # Packages
> &PACKAGES
> useOBCS=.TRUE.,
> useKPP =.TRUE.,
> &
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!FESXL3f13-HEtS0Gpg-VY86MgAy-PM4Ai4Q1CWC9SWho-hAR-5s2KHK44htyZaLDvWLlzjo3dQZ2mwDVIV5cfG84t0U2dw$
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20220918/7ddab79e/attachment.html>
More information about the MITgcm-support
mailing list