[MITgcm-support] 2d and 3d-diagnostics in one netcdf file?

Martin Losch Martin.Losch at awi.de
Tue Jan 23 12:37:46 EST 2007


Hi Andrea,

thanks, it's amazing what the package can do: "multi-year january  
averages"! Excellent.
I'll put ETAN into a separate file rather than filling it with zeros.

Martin

On 23 Jan 2007, at 18:22, Andrea Molod wrote:

> hi martin,
>
> On Tue, 23 Jan 2007, Martin Losch wrote:
>
>> Hi all,
>>
>> I am embarrassed to ask because it may have been discussed already  
>> and I just cannot find it:
>> 1. Time-averaging with the diagnostics package: the documentation  
>> says that frequency > 0 produces averages, frequency < 0 produces  
>> snapshots. What does averagingFeq do?
>
> martin i would have to look into this a bit, but averagingFeq
> has to do with the periodic averaging option. for example, multi- 
> year january average from a multi-year run, or monthly mean
> diurnal cycles. i don't remember exactly how this works, ie,
> exactly which freq is specified in averagingFeq.
>
>>
>> 2. Is it possible in the diagnostics package to have a 2D and 3D  
>> field in the same output stream, e.g.,
>>
>> frequency(5) = 43200.,
>> filename(5)  = 'state_diags',
>> fields(1,5)  = 'ETAN','THETA','SALT','UVEL','VVEL','WVEL',
>>
>> ? With these parameters I seem to be getting the level1 of ETAN  
>> (which is correct), but then I also get only the first levels for  
>> THETA, SALT, etc. here's the header of the netcdf file, note, that  
>> both Zmd000001 and Zld000001 = 1. Is that to be expected?
>
>
> yeah, martin. when you don't specify the levels in data.diagnostics
> it uses the minimum of the number of levels of all the diagnostics
> in the output stream (see diagnostics_set_pointers) - if you specify
> the whole list of levels in data.diagnostics it will write ETAN and
> then a bunch of zeros. worth a try?
>
> andrea
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list