[MITgcm-support] NaNs in output

Yuan Lian lian at ashimaresearch.com
Mon Oct 1 02:54:25 EDT 2012


I think the model runs global simulation in cube sphere grid, so no delX and  
delY are defined. 

It looks like the model runs unstable. What are the CFL values before NaNs  
appear? Dose the temperature field look  suspicious too?

Yuan

-----Original message-----
From: Martin Losch <Martin.Losch at awi.de>
To: mitgcm-support at mitgcm.org
Sent: Sun, Sep 30, 2012 23:40:02 PDT
Subject: Re: [MITgcm-support] NaNs in output

The model becomes numerically unstable (and this can happen after some time,  
if the forcing is variable) and you need to investigate your input  
parameters (time step, forcing etc.).

I am surprised it runs at all. What is your grid resolution? There is no  
dx/ySpacing or delX/delY in your PARM04 namelist.
It is also strange that your cg2d: line has zeros for the rhsMax.

Martin

On Oct 1, 2012, at 8:01 AM, Fengchao Yao wrote:

> Dear all,
> 
> The MITgcm I am running is constantly giving NaNs after a couple years of  
stable running.  Below is the information from the STDOUT file where the  
NaNs occur. Anybody has some idea on how to handle this? Thanks for your  
help.
> 
> Fengchao
> 
> STDOUT:
>  
----------------------------------------------------------------------------- 
---------------------------------
> cg2d: Sum(rhs),rhsMax =   7.95491814529655E+03  0.00000000000000E+00
> cg2d: Sum(rhs),rhsMax =   7.95490803828502E+03  0.00000000000000E+00
> SURF_ADJUSTMENT: Iter=   3348930 Nb_pts,Vol=         1  7.60072406E+06
> cg2d: Sum(rhs),rhsMax =   7.95489793127349E+03  0.00000000000000E+00
> SURF_ADJUSTMENT: Iter=   3348931 Nb_pts,Vol=         5  1.96015022E+08
> cg2d: Sum(rhs),rhsMax =  -1.77088743107612E+21  0.00000000000000E+00
> SURF_ADJUSTMENT: Iter=   3348932 Nb_pts,Vol=    128878  6.73412758E+43
> cg2d: Sum(rhs),rhsMax =                    NAN  0.00000000000000E+00
> cg2d: Sum(rhs),rhsMax =                    NAN                   NAN
> cg2d: Sum(rhs),rhsMax =                    NAN                   NAN
> cg2d: Sum(rhs),rhsMax =                    NAN                   NAN
> cg2d: Sum(rhs),rhsMax =                    NAN                   NAN
> cg2d: Sum(rhs),rhsMax =                    NAN                   NAN
> cg2d: Sum(rhs),rhsMax =                    NAN                   NAN
> cg2d: Sum(rhs),rhsMax =                    NAN                   NAN
> 
> 
> data file
>  
----------------------------------------------------------------------------- 
--------------------------------
> &PARM01
> tRef               =  
3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13.,
>                      12.,11.,2*9.,8.,7.,6.,
> sRef               = 25*38,
> no_slip_sides=.FALSE.,
> no_slip_bottom=.TRUE.,
> 
> viscAh = 3.0e1,
> 
> viscA4 = 0e7,
> viscAhGrid = 0.00,
> 
> viscA4Grid = 0.03,
> 
> viscAr = 7.E-4,
> 
> diffK4T=0e7,
> diffK4S=0e7,
> diffKrT=1.E-5,
> diffKrS=1.E-5,
> diffKrBL79surf=0.1E-4,
> diffKrBL79deep=1.0E-4,
> bottomDragQuadratic= 0.002,
> implicitDiffusion=.TRUE.,
> implicitViscosity=.TRUE.,
> gravity=9.81,
> rhonil=1029.,
> rhoConstFresh=1000.,
> eosType='JMD95Z',
> staggerTimeStep=.TRUE.,
> vectorInvariantMomentum=.TRUE.,
> useEnergyConservingCoriolis=.true.,
> implicitFreeSurface=.TRUE.,
> exactConserv=.TRUE.,
> useSingleCpuIO=.TRUE.,
> tempAdvScheme=30,
> saltAdvScheme=30,
> nonlinFreeSurf=1,
> useRealFreshWaterFlux=.FALSE.,
> allowFreezing=.false.,
> hFacMin=.3,
> hFacMinDr=30.,
> readBinaryPrec=64,
> /
> &PARM02
> cg2dMaxIters=5000,
> cg2dTargetResWunit=1.E-11,
> /
> &PARM03
> nIter0=_stTIME_,
> nTimeSteps=_nsteps_,
> deltaT = _delt_,
> pickupStrictlyMatch=.FALSE.,
> cAdjFreq=0.,
> abEps = 0.01,
> forcing_In_AB=.FALSE.,
> pChkptFreq=31104000.0,
> chkptFreq=31104000.0,
> taveFreq=259200.0,
> periodicExternalForcing=.FALSE.
> #dumpFreq=3600.,
> monitorFreq=2635200000.,
> # tauThetaClimRelax= 2592000.,
> # tauSaltClimRelax= 7776000.,
> /
> &PARM04
> usingCurvilinearGrid=.TRUE.,
> delR= 10.2652, 11.8262, 13.6398, 15.7495, 18.2060, 21.0677, 24.4002,
> 28.2749, 32.7643, 37.9348, 43.8340, 50.4735, 57.8103, 65.7382,
> 74.1055, 82.7805, 91.7766,101.4306,112.6418,127.2439,148.7415,
> 183.8161,244.2570,340.5531,442.0510,
> /
> &PARM05
> bathyFile='redsea_bathy_fixed.bin',
> # bathyFile='red2_bathyed-10.bin',
> hydrogThetaFile='lev01_temp.bin',
> hydrogSaltFile='lev01_salt.bin',
> zonalWindFile ='uf_198002.data',
> meridWindFile='vf_198002.data',
> surfQFile='qf_198002.data',
> EmPmRFile='ef_198002.data',
> surfQswFile='qsf_198002.data',
> /
> 
> 
> This message and its contents including attachments are intended solely  
for the original recipient. If you are not the intended recipient or have  
received this message in error, please notify me immediately and delete this  
message from your computer system. Any unauthorized use or distribution is  
prohibited. Please consider the environment before printing this email.
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support


_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org
http://mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120930/00e8fa27/attachment-0001.htm>


More information about the MITgcm-support mailing list