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

Dimitris Menemenlis menemenlis at sbcglobal.net
Mon May 5 18:55:16 EDT 2008 

Hey Andrea!  It's so nice to hear from you.
Thanks for solving this puzzler (pronounced as in NPR's car talk).
Hope all is well.  Cheers, Dimitris

Dimitris Menemenlis <menemenlis at sbcglobal.net>
5056 Oakwood Ave, La Canada, CA 91011-2450
tel/fax: 818-790-6735;  cell: 818-625-6498

On May 5, 2008, at 3:49 PM, Andrea Molod wrote:

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

More information about the MITgcm-devel mailing list