<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<title></title>
</head>
<body>
<div name="messageBodySection">
<div dir="auto">Dear all,<br />
<br />
I’m running a simulation for a freshwater-discharged coastal channel. The model is run with multi processors and it could restart from the pickup files and end normally without reporting any error. However, the output results (T, S, U) become different from those at the same timestep in a non-restart run, after initial conditions being loaded correctly from the pickup. That is to say, the 1st-iteration output in the restart is identical to the same-realtime output in a non-restart run but consequent results begin to differ. A plot showing the divergence at one grid point is attached here:<br />
<br />
<br />
<img height="auto" src="cid:237103AF5139482C98D1A2A37923D40B" /><br />
<br />
The data file I used is as follow:<br />
<br />
&PARM01<br />
tRef=90*0.,<br />
sRef=90*34.,<br />
viscAz=0.00001,<br />
viscC2smag=2.2,<br />
diffKzT=0.00001,<br />
diffKhT=20,<br />
diffKzS=0.00001,<br />
diffKhS=20,<br />
bottomDragQuadratic=2.5E-3,<br />
implicitViscosity=.TRUE.,<br />
implicitDiffusion=.TRUE.,<br />
no_slip_sides=.FALSE.,<br />
no_slip_bottom=.TRUE.,<br />
rigidLid=.FALSE.,<br />
implicitFreeSurface=.TRUE.,<br />
hFacMin=0.05,<br />
readBinaryPrec=64,<br />
staggerTimeStep=.TRUE.,<br />
tempAdvScheme=33,<br />
saltAdvScheme=33,<br />
eosType='JMD95Z',<br />
f0=1.2E-4,<br />
beta=0.E-11,<br />
selectAddFluid=1.,<br />
useSingleCpuIO=.TRUE.,<br />
&<br />
<br />
# Elliptic solver parameters<br />
&PARM02<br />
cg2dMaxIters=300,<br />
cg2dTargetResidual=1.E-13,<br />
cg3dMaxIters=20,<br />
cg3dTargetResidual=1.E-8,<br />
&<br />
<br />
# Time stepping parameters<br />
&PARM03<br />
nIter0=0,<br />
endTime=20000.,<br />
deltaT=2.,<br />
abEps=0.01,<br />
pChkptFreq=1000.,<br />
chkptFreq=0.,<br />
pickupSuff='',<br />
dumpFreq=1000.,<br />
monitorSelect=1,<br />
monitorFreq=1000.,<br />
periodicExternalForcing=.FALSE.,<br />
&<br />
<br />
# Gridding parameters<br />
&PARM04<br />
usingCartesianGrid=.TRUE.,<br />
delXfile='delx.bin',<br />
delYfile='dely.bin',<br />
delZ= 20*2,30*4,40*6.,<br />
&<br />
<br />
# Input datasets<br />
&PARM05<br />
bathyFile='bathymetry.bin',<br />
hydrogThetaFile='tempini.bin',<br />
hydrogSaltFile ='saltini.bin',<br />
&<br />
<br />
I only modified nIter0 every time before restart, is there anything else needs attention? Thanks a lot for your time!</div>
</div>
<div name="messageSignatureSection"><br />
<div class="matchFont">Best, Weiyang</div>
</div>
</body>
</html>