[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