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

[Dimitris Menemenlis]

Jean-Michel, the problem is not memory space.  The problem is numperlist
For some reason, these four diagnostics, on their own, request  
numperlist=2.
It's a bit strange behavior.  First, I don't understand why SI_Tice1  
should be mated
with SI_Fract.  I thought that mating was only needed for vector  
quantities to take
care of edges properly.  Second, even mated quantities, e.g., UVELMASS  
and
VVELMASS or SIuice and SIvice, do not require numperlist=2, unless of  
course
I try to save two fields in the same file.

This is not a high priority item.  I just thought I would report it so  
we have a record
and maybe there is a simple explanation that I am not seeing and also  
because I
did end up wasting a few hours tracking it down so others don't  
stumble on same issue.

How is the downslope parameterization coming along?  Will you have any  
time
to work on it in next few weeks?

Cheers, 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

On May 5, 2008, at 3:05 PM, Jean-Michel Campin wrote:

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