[MITgcm-support] output frequencies NOT consistent

Jean-Michel Campin jmc at ocean.mit.edu
Wed Sep 21 14:45:06 EDT 2005


Hi Paola,

The default behavior is to write a pickup & to dump the state
variables at the end of the run, whatever the frequencies are.
This is the reason why you get
> the instantaneous fields and pickup.chkptA files output at the correct time
at iteration 650000 (= the end of this run, niter0+nTimeSteps=650000)

For the time-average and pickup.000NNNN, they are written according to
taveFreq & pChkptFreq , each time the model time is a multiple of
those frequencies (in your case, both are 250.E6 (s)).
Starting at t=nIter0*deltaT = 400.E6 (s), I expect to get those file
written at t=500.E6 (s), 750.E6 (s) ... and so on.
and this is what the model did.

I hope this helps to clarify things.

Jean-Michel

On Wed, Sep 21, 2005 at 05:04:07AM -0700, Paola Cessi wrote:

> 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
> 

> # ====================
> # | 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=,
>  &

> # 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.,
>  &
> 

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




More information about the MITgcm-support mailing list