[MITgcm-support] rdmds.m: how to read a given field in a multi-field file?
Menemenlis, Dimitris (3248)
Dimitris.Menemenlis at jpl.nasa.gov
Sun Mar 13 17:43:14 EDT 2011
Alastair, in example below you expect
(8*Nr + 3) records
since there are 8 3D field and 3 2D fields.
Nr is the vertical dimension.
I assume that Nr=5 in your example below?
On Mar 13, 2011, at 1:37 PM, Alastair McKinstry wrote:
> I think I gave a bad example, as the above has both 3 records and 3 fields.
> A better example is:
>
> octave:1> [AA,ITER,MM]=rdmds('pickup.ckptA');
> octave:2> MM
> MM = nDims = [ 2 ]; dimList = [ 192 32 ]; dataprec = [ 'float64' ]; nrecords = [ 43 ]; timeStepNumber = [ 144 ]; timeInterval = [ 8.640000000000E+04 ]; nFlds = [ 11 ]; fldList = { 'Uvel ' 'Vvel ' 'Theta ' 'Salt ' 'GuNm1 ' 'GvNm1 ' 'GtNm1 ' 'GsNm1 ' 'EtaN ' 'dEtaHdt ' 'EtaH ' };
> octave:3> ITER
> ITER = -1
> octave:4> size(AA)
> ans =
>
> 192 32 43
>
>
> This has 43 records but 11 fields. I think the squeeze example above works to pick out the records, but not the fields.
More information about the MITgcm-support
mailing list