[MITgcm-support] Diagnostics Package
Kirk Bryan
kbryan at Princeton.EDU
Fri Apr 25 21:03:14 EDT 2014
I am using the mitgcm to test stochastic backscatter in ocean mesoscale eddy simulations.
The model does a great job. Presently my set up works fine for 6 diagnostic variables. I want to extend to 12 diagnostic
quantities. I doubled all dimensions in
DIAGNOSTICS_SIZE.h
and recompiled and then specified the new set of 12 diagnostics in data.diagnostics. The two
files are shown below. It seems to hang up in "diagnostics_set_pointers.f" after doing only 6
of the diagnostics. All of the diagnostics are in the "available_diagnostics.log". Can you see what
I might be doing wrong?
Kirk
files
Comparison of the old and new data.diagnostics file
energy_apr14/T8_160_Tz2_Ty4LF_20_energy3> diff data.diagnostics data.diagnostics.old
16,22c16,17
< # fields(1,1) = 'UVEL ','VVEL ','WVEL ',
< # 'THETA ','momKE ','PHIHYD ',
< fields(1,1) = 'UVEL ','VVEL ','Um_dPHdx',
< 'Um_Advec','Um_Diss','Um_Ext ',
< 'Vm_Ext' ,'Vm_dPHdy','Vm_Advec',
< 'Vm_Diss','AB_gU ','AB_gV' ,
<
---
> fields(1,1) = 'UVEL ','VVEL ','WVEL ',
> 'THETA ','momKE ','PHIHYD ',
Comparison of the new to old DIAGNOSTICS_SIZE.h
T8_160_Tz2_Ty4LF_20_energy3/code> diff DIAGNOSTICS_SIZE.h DIAGNOSTICS_SIZE.old
23,27c23,27
< PARAMETER( ndiagMax = 10000 )
< PARAMETER( numlists = 40, numperlist = 200, numLevels=2*Nr )
< PARAMETER( numdiags = 6000*Nr )
< PARAMETER( nRegions = 0 , sizRegMsk = 1 , nStats = 8 )
< PARAMETER( diagSt_size = 100*Nr )
---
> PARAMETER( ndiagMax = 5000 )
> PARAMETER( numlists = 20, numperlist = 100, numLevels=2*Nr )
> PARAMETER( numdiags = 3000*Nr )
> PARAMETER( nRegions = 0 , sizRegMsk = 1 , nStats = 4 )
> PARAMETER( diagSt_size = 50*Nr )
33d32
Results of the new run with 12 diagnostics
more STDERR.0005
(PID.TID 0005.0001) *** ERROR *** DIAGNOSTICS_SET_POINTERS: is not a Diagnostic
more STDOUT.0001
(PID.TID 0001.0001)
(PID.TID 0001.0001) // ======================================================
(PID.TID 0001.0001) // MITgcm UV
(PID.TID 0001.0001) // =========
(PID.TID 0001.0001) // ======================================================
(gap)
(PID.TID 0001.0001) -----------------------------------------------------
(PID.TID 0001.0001) DIAGNOSTICS_READPARMS: active diagnostics summary:
(PID.TID 0001.0001) -----------------------------------------------------
(PID.TID 0001.0001) Creating Output Stream: main1
(PID.TID 0001.0001) Output Frequency: 172800.000000 ; Phase: 0.000000
(PID.TID 0001.0001) Averaging Freq.: 172800.000000 , Phase: 0.000000 , Cycle: 1
(PID.TID 0001.0001) missing value: -9.990000000000E+02
(PID.TID 0001.0001) Levels: will be set later
(PID.TID 0001.0001) Fields: UVEL VVEL Um_dPHdx Um_Advec Um_Diss Um_Ext Vm_Ext Vm_dPHdy Vm_Advec
(PID.TID 0001.0001) Fields: Vm_Diss AB_gU AB_gV
(PID.TID 0001.0001) -----------------------------------------------------
(PID.TID 0001.0001) DIAGNOSTICS_READPARMS: statistics diags. summary:
(PID.TID 0001.0001) -----------------------------------------------------
(PID.TID 0001.0001)
(PID.TID 0001.0001) SET_PARMS: done
(gap)
=================================
(PID.TID 0001.0001) tempAdvScheme = /* Temp. Horiz.Advection scheme selector */
(PID.TID 0001.0001) 7
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) tempVertAdvScheme = /* Temp. Vert. Advection scheme selector */
(PID.TID 0001.0001) 7
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) tempMultiDimAdvec = /* use Muti-Dim Advec method for Temp */
(PID.TID 0001.0001) T
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) tempSOM_Advection = /* use 2nd Order Moment Advection for Temp */
(PID.TID 0001.0001) F
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) AdamsBashforthGt = /* apply Adams-Bashforth extrapolation on Gt */
(PID.TID 0001.0001) F
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) AdamsBashforth_T = /* apply Adams-Bashforth extrapolation on Temp */
(PID.TID 0001.0001) F
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) saltAdvScheme = /* Salt. Horiz.advection scheme selector */
(PID.TID 0001.0001) 7
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) saltVertAdvScheme = /* Salt. Vert. Advection scheme selector */
(PID.TID 0001.0001) 7
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) saltMultiDimAdvec = /* use Muti-Dim Advec method for Salt */
(PID.TID 0001.0001) T
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) saltSOM_Advection = /* use 2nd Order Moment Advection for Salt */
(PID.TID 0001.0001) F
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) AdamsBashforthGs = /* apply Adams-Bashforth extrapolation on Gs */
(PID.TID 0001.0001) F
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) AdamsBashforth_S = /* apply Adams-Bashforth extrapolation on Salt */
(PID.TID 0001.0001) F
(PID.TID 0001.0001) ;
(PID.TID 0001.0001) // ===================================
(PID.TID 0001.0001) ------------------------------------------------------------
(PID.TID 0001.0001) DIAGNOSTICS_SET_LEVELS: done
(PID.TID 0001.0001) Total Nb of available Diagnostics: ndiagt= 202
(PID.TID 0001.0001) SETDIAG: Allocate 8 x 1 Levels for Diagnostic # 30 UVEL
(PID.TID 0001.0001) SETDIAG: Allocate 8 x 1 Levels for Diagnostic # 31 VVEL
(PID.TID 0001.0001) SETDIAG: Allocate 8 x 1 Levels for Diagnostic # 165 Um_dPHdx
(PID.TID 0001.0001) SETDIAG: Allocate 8 x 1 Levels for Diagnostic # 161 Um_Advec
(PID.TID 0001.0001) SETDIAG: Allocate 8 x 1 Levels for Diagnostic # 159 Um_Diss
(PID.TID 0001.0001) SETDIAG: Allocate 8 x 1 Levels for Diagnostic # 167 Um_Ext
more T8_160_Tz2_Ty4LF_20.e1202380
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS
[NID 01154] 2014-04-25 13:59:44 Apid 33404492: initiated application termination
======== BEGIN EPILOGUE ========
Moab JobID: gaea.9651226
RM JobID: 1202380.c1-sys0.ncrc.gov
State: Running
Name: T8_160_Tz2_Ty4LF_20
User: Kirk.Bryan
--More--(64%)
Hangs up after only 6 diagnostics
Kirk Bryan
Sayre Hall, AOS Program
Princeton University
Princeton, NJ 08544
tel 609-258-3688
fax 609-258-2850
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20140426/6b75c460/attachment.htm>
More information about the MITgcm-support
mailing list