[MITgcm-support] rdnctiles

Martin Losch Martin.Losch at awi.de
Thu May 10 09:28:14 EDT 2007


Hi all,

sorry, partially my mistake, which 'bytile' it seems to work half the  
way:
s=rdnctiles('diags2D*','ETAN',[],'bytile',20);
works
but
tlev.tdname='T';
tlev.tvname='T';
tlev.tvals=3079296000;
s=rdnctiles('diags2D*','ETAN',tlev,'bytile',20);
does not:
> Total files found:  18
>   diags2D.0000000000.t001.nc
>   diags2D.0000000000.t002.nc
>   diags2D.0000000000.t003.nc
>   diags2D.0000000000.t004.nc
>   diags2D.0000000000.t005.nc
>   diags2D.0000000000.t006.nc
>   diags2D.0000000000.t007.nc
>   diags2D.0000000000.t008.nc
>   diags2D.0000000000.t009.nc
>   diags2D.0000000000.t010.nc
>   diags2D.0000000000.t011.nc
>   diags2D.0000000000.t012.nc
>   diags2D.0000000000.t013.nc
>   diags2D.0000000000.t014.nc
>   diags2D.0000000000.t015.nc
>   diags2D.0000000000.t016.nc
>   diags2D.0000000000.t017.nc
>   diags2D.0000000000.t018.nc
> ??? Undefined function or method 'eq' for input arguments of type  
> 'struct'.
>
> Error in ==> ismember at 78
>           found = find(a(i)==s(:));  % FIND returns indices for LOC.
>
> Error in ==> intersect at 86
>     [tf,pos] = ismember(b,a);
>
> Error in ==> rdnctiles at 116
>         v = intersect(vit.vars.(allvars{i}),tlev.tvals);
>

Ed, you can find the files in
eddy.csail.mit.edu:~mlosch/scratch/rdnctiles

Martin



More information about the MITgcm-support mailing list