[MITgcm-support] NaN output data

[yulinghui]

Dear all:
       I apply the OBCS (OBWU=U0*sin(wt)) to generate the internal waves. If my compute zone is extend to 3000km and dx(min)=1000m, the numerical model could continue running without NaN.  But  if my compute zone is 1500km and dx(min)=200m, the numerical result would have NaN in the output data after about 3 tide cycle.  I want to use the high resolution to study the internal solitary wave, how should i do to improve it except reducing the time step(now is 10s).
Best regards
   




input file data:
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRefFile  = 'MD.tRef', 
 sRef= 200*35.,
 viscAz=1.E-3,
 viscAh=1.E-2,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.FALSE.,
 diffKhT=1.E-2
 diffKzT=1.E-3,
 f0=1.07E-4,
 beta=0.E-11,
 eosType='LINEAR',
 tAlpha=2.E-4,
 sBeta =0.E-4,
 gravity=9.81,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.
 nonHydrostatic=.TRUE.,
 hFacMin=0.2,
 implicSurfPress=0.5,
 implicDiv2DFlow=0.5,
 nonlinFreeSurf=0,
 hFacInf=0.2,
 hFacSup=1.8,
 saltStepping=.FALSE.,
 staggerTimeStep=.TRUE.
 tempAdvScheme=77, 
 tempVertAdvScheme=77, 
#- not safe to use globalFiles in multi-processors runs
#globalFiles=.TRUE.,
 readBinaryPrec=64,
 writeBinaryPrec=64,
 writeStatePrec=64,
 &


# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-13,
 cg3dMaxIters=400,
 cg3dTargetResidual=1.E-13,
 &


# Time stepping parameters
 &PARM03
 nIter0=0,
 nTimeSteps=44600,
 deltaT=10,
 abEps=0.1,
 pChkptFreq=446000.,
 chkptFreq=0.,
 dumpFreq=11150
 monitorFreq=2500.,
 monitorSelect=200,
 &


# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 delXfile='delXvar',
 delY=1.E3,
 delRfile='delZvar',
 &


# Input datasets
 &PARM05
 hydrogThetaFile='T.init',
 bathyFile='topog',
 &


input data data.obcs
# Open-boundaries
 &OBCS_PARM01
 OB_Ieast=-1,
 OB_Iwest=1,
 useOrlanskiEast=.TRUE.,
# useOBCSbalance=.TRUE.,
# OBCS_balanceFacE= 1.,
# OBCS_monitorFreq=4000.,
 &


# Orlanski parameters
 &OBCS_PARM02
 Cmax=0.45,
 cVelTimeScale=50.,
 &
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20131209/af0857a4/attachment.htm>