[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