[MITgcm-support] rdmnc.m

[David Wang]

Hello all,

I've got trouble loading multi-tile two-dimensional mnc diagnostics 
(i.e., surface quantities like ETAN, TAUX) output into matlab. rdmnc.m 
works just fine when I read in one single iteration, but failed when I 
tried to load all the iterations available. What puzzles me is that 
rdmnc.m works well for 3-dimensional diagnostics (THETA, UVEL etc) in 
both cases. Did I get anything wrong? I'm attaching my data.diagnostics 
and data.mnc namelists for reference. The code I use is 
checkpoint57w_post. matlab is R14 (SP3) on linux FC4.

Thanks, D.W.

 >> S = rdmnc('diag2d.*.nc')   % not work for multiple iter.
??? Error using ==> eval
Subscripted assignment dimension mismatch.

Error in ==> rdmnc>rdmnc_local at 194
        eval(['S.(cvar)(',Sstr(1:end-1),')=tmpdata;'])

Error in ==> rdmnc at 98
    S=rdmnc_local(nc,varlist,iters,S);

 >> S = rdmnc('diag2d.*.nc',80)   % works for one iter.

S =

              attributes: [1x1 struct]
    iters_read_from_file: 80
                       T: 864000
                    iter: 80
             diag_levels: 1
                     Xp1: [161x1 double]
                       Y: [100x1 double]
                    TAUX: [161x100 double]
                       X: [160x1 double]
                     Yp1: [101x1 double]
                    TAUY: [160x101 double]
                   TFLUX: [160x100 double]

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20051126/8042d439/attachment.htm>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: data.diagnostics
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20051126/8042d439/attachment.el>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: data.mnc
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20051126/8042d439/attachment-0001.el>