[MITgcm-devel] diagnostics for SI_Tice1, SI_Tice2, SI_Qice1, and SI_Qice2

[Dimitris Menemenlis]

A little more info towards solving the puzzler. SI_Fract is defined as  
the diagmate of SI_Tice1:

         diagName  = 'SI_Tice1'
         diagTitle = 'Sea-Ice Temperature, 1srt layer (area weighted)'
         diagUnits = 'degC            '
         diagCode  = 'SM  C   M1      '
         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
      I     diagName, diagCode, diagUnits, diagTitle, numFract, myThid )

but why does pkg/diagnostics request to write the two outputs in the  
same file.  The behavior below
remains even when diagnostic output for SI_Fract is specified in a  
separate output file.

D.

Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Institute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  cell: 818-625-6498;  fax: 818-393-6720

On May 5, 2008, at 2:15 PM, Dimitris Menemenlis wrote:

> There is some strange behavior in requesting diagnostics for  
> SI_Tice1, SI_Tice2, SI_Qice1, and SI_Qice2.
> For example if I put the following in the data.diagnostics file:
>
> &diagnostics_list
>  frequency(1) = 2635200.,
>   fields(1,1) = 'SI_Tice1',
>   filename(1) = 'SI_Tice1',
> &
>
> The model stops and I get the following error message:
>
> (PID.TID 0000.0001) DIAGNOSTICS_SET_LEVELS: done
> (PID.TID 0000.0001)  Total Nb of available Diagnostics: ndiagt=   273
> (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic  
> #   205 SI_Tice1
> (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Count.Diag  
> #   201 SI_Fract
> DIAGNOSTICS_SET_POINTERS: Not enough space for all active  
> diagnostics (from data.diagnostics)
> DIAGNOSTICS_SET_POINTERS: numperlist=     1 but needs at least     2
> STOP ABNORMAL END: S/R DIAGNOSTICS_SET_POINTERS statement executed
>
> (PID.TID 0000.0001) *** ERROR *** DIAGNOSTICS_SET_POINTERS: Not  
> enough space for all active diagnostics (from data.diagnostics)
> (PID.TID 0000.0001) *** ERROR *** DIAGNOSTICS_SET_POINTERS:  
> numperlist=     1 but needs at least     2
>
> It's as if I was also requesting diagnostics for SI_Fract in the sam  
> file, but I am not.  Same thing happens for
> SI_Tice2, SI_Qice1, and SI_Qice2.
>
> Is this normal-known behavior or is this a bug?
>
> Dimitris
>
> Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
> Jet Propulsion Lab, California Institute of Technology
> MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
> tel: 818-354-1656;  cell: 818-625-6498;  fax: 818-393-6720
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel