[MITgcm-support] output frequencies NOT consistent
Paola Cessi
pcessi at ucsd.edu
Wed Sep 21 08:04:07 EDT 2005
I downloaded the latest CVS version on 9/19/2005, but I cannot get pickup
files and time averages to output at the final times of my runs. They
output at an earlier time. In the example here, the instantaneous fields
and pickup.chkptA files output at the correct time
(niter0+dumpfreq=650000), but the pickup.000NNNN and the time-averaged
files output at 500000.
What am I doing wrong?
Attached are the data and data.diagnostics that I used. They have:
nIter0=400000,
nTimeSteps=250000,
deltaT=1000,
pChkptFreq=250.E6,
chkptFreq=250.E6,
dumpFreq=250.E6,
monitorFreq=250.E6,
taveFreq=250.E6,
and
&diagnostics_list
frequency(1) = 250.E6,
fields(1,1) = 'UVEL ','VVEL ','WVEL ',
'THETA ','THETASQ ',
filename(1) = 'dyn',
frequency(2) = 250.E6,
fields(1,2) = 'UVELTH ','VVELTH ','WVELTH ',
'DIFx_TH ','DIFy_TH ',
# 'DFrI_TH ',
filename(2) = 'thetaflx',
Thanks,
paola
-------------------------------------------------------------------------
Paola Cessi Tel: +1 858 534 0622
Scripps Institution of Oceanography Fax: +1 858 534 8045
UCSD-0213 e-mail: pcessi at ucsd.edu
La Jolla, CA 92093-0213
USA Web: http://www-pord.ucsd.edu/~pcessi
On sabbatical leave until November 30th, 2005 at:
Courant Institute of Mathematical Sciences Tel: +1 212 998 3154
New York University Fax: +1 212 995 4121
251 Mercer Street
New York, NY, 10012-1185
-------------- next part --------------
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
tRef=-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,-9.8,
sRef=30*35.,
viscAh=30.E0,
viscAz=2.E-3,
no_slip_sides=.TRUE.,
no_slip_bottom=.FALSE.,
viscA4=1.E9,
diffK4T=1.E9,
diffKhT=8.E-5,
diffKzT=8.E-5,
diffKhS=1.E0,
diffKzS=1.E0,
f0=1.E-4,
beta=1.6E-11,
tAlpha=2.E-4,
sBeta =0.,
gravity=10.,
rhoConst=1000.,
rhoNil=1000.,
heatCapacity_Cp=3900.,
rigidLid=.TRUE.,
implicitFreeSurface=.FALSE.,
saltAdvection=.FALSE.,
saltForcing=.FALSE.,
saltStepping=.FALSE.,
eosType='LINEAR',
nonHydrostatic=.FALSE.,
momAdvection=.TRUE.,
implicitViscosity=.TRUE.,
implicitDiffusion=.FALSE.,
readBinaryPrec=64,
tempAdvScheme=77,
staggerTimeStep=.TRUE.,
bottomDragLinear=4.4E-4,
debugLevel=-1,
useJamartWetPoints=.TRUE.,
# hFacMin=0.1,
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=40,
cg2dTargetResidual=1.E-7,
cg3dMaxIters=40,
cg3dTargetResidual=2.E-9,
&
# Time stepping parameters
&PARM03
nIter0=400000,
nTimeSteps=250000,
deltaT=1000,
abEps=0.1,
pChkptFreq=250.E6,
chkptFreq=250.E6,
dumpFreq=250.E6,
monitorFreq=250.E6,
tauThetaClimRelax=764400.,
cAdjFreq=0,
taveFreq=250.E6,
&
# Gridding parameters
&PARM04
usingCartesianGrid=.TRUE.,
usingSphericalPolarGrid=.FALSE.,
dXspacing=10.4166666E3,
dYspacing=10.4166666E3,
#
delZ=7.84,13.00,18.16,23.26,28.30,33.26,38.13,42.90,47.54,52.06,56.43,60.65,64.70,68.58,73.27,75.76,79.04,82.11,84.95,87.56,89.82,91.95,93.81,95.43,96.78,97.85,98.67,99.21,99.49,99.49,
# delZ=1.7155,2.8127,4.5955,7.4654,12.0173,19.0699,29.6094,44.5400,64.1643,87.5647,112.4353,135.8357,155.4600,170.3906,180.9301,187.9827,192.5346,195.4045,197.1873,198.2845,
# delZ=26.0217, 30.8931, 36.4851, 42.8330, 49.9480, 57.8101, 66.3624, 75.5081, 85.1115, 95.0042, 104.9958, 114.8885, 124.4919, 133.6376, 142.1899, 150.0520, 157.1670, 163.5149, 169.1069, 173.9783,
&
# Input datasets
&PARM05
thetaClimFile='TempSurf.semi.sin_y',
bathyFile='topo_semi.box',
hydrogThetaFile=,
hydrogSaltFile=,
zonalWindFile='windx.semi.sin_y_clean',
meridWindFile=,
&
-------------- next part --------------
# Diagnostic Package Choices
#-----------------
# for each output-stream:
# filename(n) : prefix of the output file name (only 8.c long) for outp.stream n
# frequency(n):< 0 : write snap-shot output every |frequency| seconds
# > 0 : write time-average output every frequency seconds
# timePhase(n) : write at time = timePhase + multiple of |frequency|
# levels(:,n) : list of levels to write to file (Notes: declared as REAL)
# when this entry is missing, select all common levels of this list
# fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics" file
# for the list of all available diag. in this particular config)
#-----------------
&diagnostics_list
diag_mnc=.FALSE.,
frequency(1) = 250.E6,
fields(1,1) = 'UVEL ','VVEL ','WVEL ',
'THETA ','THETASQ ',
filename(1) = 'dyn',
frequency(2) = 250.E6,
fields(1,2) = 'UVELTH ','VVELTH ','WVELTH ',
# 'ADVx_TH ',
#'ADVy_TH ',
#'ADVr_TH ',
'DIFx_TH ','DIFy_TH ',
# 'DFrI_TH ',
filename(2) = 'thetaflx',
&
# Parameter for Diagnostics of per level statistics:
#-----------------
# for each output-stream:
# stat_fname(n) : prefix of the output file name (only 8.c long) for outp.stream n
# stat_freq(n):< 0 : write snap-shot output every |stat_freq| seconds
# > 0 : write time-average output every stat_freq seconds
# stat_phase(n) : write at time = stat_phase + multiple of |stat_freq|
# stat_region(:,n) : list of "regions" (default: 1 region only=global)
# stat_fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics.log"
# file for the list of all available diag. in this particular config)
#-----------------
&DIAG_STATIS_PARMS
# an example just to check the agreement with MONITOR output:
# stat_fields(1,1)= ','UVEL ','VVEL ','WVEL ', 'THETA ',
# stat_fname(1)= 'dynStDiag',
## stat_freq(1)= -864000.,
# stat_freq(1)= 10000.,
# stat_phase(1)= 10000.,
&
More information about the MITgcm-support
mailing list