/[MITgcm]/MITgcm_contrib/jmc_script/grph_StD.m
ViewVC logotype

Diff of /MITgcm_contrib/jmc_script/grph_StD.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.5 by jmc, Mon Aug 25 21:35:12 2014 UTC revision 1.6 by jmc, Thu Mar 5 20:57:26 2015 UTC
# Line 7  Line 7 
7  % $Header$  % $Header$
8  % $Name$  % $Name$
9    
   
10  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;
11  nItMx=2400*ones(1,Nexp);  nItMx=2400*ones(1,Nexp);
12  namLg=namA ; namLg=strrep(namLg,'_','\_');  namLg=namA ; namLg=strrep(namLg,'_','\_');
13    undef=123456.7;
14  %-----------  %-----------
15  %- test if the variable krd is define :  %- test if the variable krd is define :
16  if size(who('krd'),1) > 0,  if size(who('krd'),1) > 0,
# Line 32  if krd > 0, Line 32  if krd > 0,
32    [ntA(n),rList,tim,vv1,listV] = ...    [ntA(n),rList,tim,vv1,listV] = ...
33      read_StD(prefix,namA(n),listV);      read_StD(prefix,namA(n),listV);
34    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
35      vv1(find(vv1==undef))=NaN;
36  %- set global dims: & load vvA --> vvB  %- set global dims: & load vvA --> vvB
37    nbV=size(listV,2);    nbV=size(listV,2);
38    nrec=nIt; n3d=nk; nReg=nRg;    nrec=nIt; n3d=nk; nReg=nRg;
# Line 42  if krd > 0, Line 43  if krd > 0,
43    [ntA(n),rList,tim,vv1,listV] = ...    [ntA(n),rList,tim,vv1,listV] = ...
44      read_StD(prefix,namA(n),listV);      read_StD(prefix,namA(n),listV);
45    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
46      vv1(find(vv1==undef))=NaN;
47    if (nrec < nIt),    if (nrec < nIt),
48      fprintf('\n');      fprintf('\n');
49      error([' Nb of records=',int2str(nIt),' exceeds nrec=',int2str(nrec)]);      error([' Nb of records=',int2str(nIt),' exceeds nrec=',int2str(nrec)]);
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22