[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