[MITgcm-support] diagnostics DFrE_TH

[jdeshayes at whoi.edu]

Hi,
I am trying to use the diagnostics package to measure on-line the fluxes of
temperature. The simulation works fine, except that I get the following warning
message :

- WARNING - from DIAGNOSTICS_OUT at iter=    414778
- WARNING -   diag.#  88 : DFrE_TH  (#  1 ) in outp.Stream: vdiffth
- WARNING -   has not been filled (ndiag= 0 )
WARNING DIAGNOSTICS_OUT  => write ZEROS instead

for each time step which I asked for the diagnostics !
and when checking the output, the data corresponding to DFrE_TH is not only
zeros but have numbers which do not look unrealistic. so what does the warning
message actually mean ?

Besides, as far as I understood, oceanic convection is represented by a local
increase of the vertical diffusivity parameter, hence a local increase of the
vertical diffusive flux of temperature. How can I track these events ? by
comparing DFrE_TH (calculated with "normal" vertical diffusivity) and DFrI_TH
(with locally increased vertical diffusivity to reproduce convection) ?

following is my packages.conf file, my data file, and my data.diagnostics file.

Thank you in advance for any help !
Julie


************************** packages.conf ******************************
# $Header: /u/gcmpack/MITgcm/verification/lab_sea/code/packages.conf,v 1.6
2004/10/07 19:31:14 edhill Exp $
# $Name:  $

cd_code
obcs
mdsio
mom_fluxform
mom_vecinv
monitor
rw
sbo
shap_filt
timeave
generic_advdiff
diagnostics


************************** data file ******************************
# Model parameters
# Continuous equation parameters
 &PARM01
 tRef=  5.35, 5.20, 5.05, 4.9, 4.75, 4.6, 4.45, 4.3, 4.15, 4.0,
 sRef=10*10.,
 viscAz=2.E-4,
 viscAh=200.,
 no_slip_sides=.TRUE.,
 no_slip_bottom=.TRUE.,
 diffKhT=120.,
 diffKzT=1.E-5,
 f0=1.0e-4,
 beta=0.E-11,
 tAlpha=2.E-4,
 sBeta =0.,
#gravity=9.81,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 eosType='LINEAR',
 tempAdvScheme=4,
 saltAdvScheme=4,
 readBinaryPrec=64,
 hFacMin=0.2,
 bottomDragLinear=1.e-6,
 implicitDiffusion=.TRUE.,
 useCDscheme=.FALSE.,
 useNHMTerms=.FALSE.,
 useSingleCpuIO=.TRUE.,
 saltStepping=.FALSE.,
 ivdc_kappa=1000.0,
 &
# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-11,
 &
# Time stepping parameters
 &PARM03
startTime=622080000,
 endTime=622252800,
 deltaT=1500.0,
 abEps=0.1,
 pChkptFreq=0.0,
 chkptFreq=86400,
 dumpFreq=86400,
 cAdjFreq=0,
# tbudgetstart=3600.0,
 &
# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.FALSE.,
 dXspacing=10.e+3,
 dYspacing=10.e+3,
 delZ=10*100.,
 &
 &PARM05
 bathyFile='topog_10km_obcs.bin',
 hydrogThetaFile='',
 hydrogSaltFile=,
 zonalWindFile=,
 meridWindFile=,
 uVelInitFile='',
 vVelInitFile='',
 surfQfile='Qnet_10km_obcs.bin',
 &
# time-dependent wind forcing
 &PARM06
 HfluxP=360.,
 HfluxW=60.,
 qmax=1.2,
 relaxtime=10.,
 xr0=375.e+3,xr1=750.e+3,yr0=0.0,yr1=100.e+3,
 &


************************** data.diagnostics ******************************
# 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
  frequency(1)=86400.,
   fields(1,1)='DFrE_TH ','DFrI_TH ',
   filename(1)='vdiffth',
 &
 &diag_statis_parms
 &





----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.