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

Andrea Molod molod at ocean.mit.edu
Mon May 5 18:49:19 EDT 2008 

hey dimitris,

i couldn't resist putting in my 2 cents here. you say its not a memory 
issue, but indeed it is. and this 'mate' is different from a vector mate.
this 'mate', the fraction, is an averaging weight for the sea ice temp.
so when you ask for sea ice temp it MUST also allocate the space for
the averaging weight which needs to be computed at the same frequency
at which you have asked for temp. whether or not the fraction appears
in the output file, the space for it is needed and must be allocated.

andrea

On Mon, 5 MayK 2008, Dimitris Menemenlis wrote:

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

-- 
----------------------------------------------------------------
Andrea M. Molod, PhD.
Department of Earth, Atmospheric & Planetary Sciences
MIT
Tel: (617) 253-5458, Email: amolod at mit.edu
WWW:http://paoc.mit.edu/paoc/people/person.asp?position=Postdoc&who=molod
----------------------------------------------------------------

More information about the MITgcm-devel mailing list