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

[Jean-Michel Campin]

Hi Dimitris,

On Mon, May 05, 2008 at 02:45:50PM -0700, Dimitris Menemenlis wrote:
> 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.

I don't think you need to output also SI_Fract when you select
SI_Tice1 ; but you do need to provide enough (memory) space:
> >(PID.TID 0000.0001) *** ERROR *** DIAGNOSTICS_SET_POINTERS: Not  
> >enough space for all active diagnostics (from data.diagnostics)

Jean-Michel

> 
> 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
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel