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

Diff of /MITgcm_contrib/jmc_script/plot_StD.m

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

revision 1.3 by jmc, Tue Sep 30 22:10:26 2014 UTC revision 1.4 by jmc, Thu Mar 5 20:57:26 2015 UTC
# Line 11  Line 11 
11  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;
12  %nItMx=60*ones(1,Nexp);  %nItMx=60*ones(1,Nexp);
13  namLg=namA ; namLg=strrep(namLg,'_','\_');  namLg=namA ; namLg=strrep(namLg,'_','\_');
14    undef=123456.7;
15  %-----------  %-----------
16  %- test if the variable krd is define :  %- test if the variable krd is define :
17  if size(who('krd'),1) > 0,  if size(who('krd'),1) > 0,
# Line 32  if krd > 0, Line 33  if krd > 0,
33    [ntA(n),rList,tim,vv1,listV,kList] = ...    [ntA(n),rList,tim,vv1,listV,kList] = ...
34      read_StD(prefix,namA(n,:),listV);      read_StD(prefix,namA(n,:),listV);
35    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
36      vv1(find(vv1==undef))=NaN;
37  %- set global dims: & load vvA --> vvB  %- set global dims: & load vvA --> vvB
38    nbV=size(listV,2);    nbV=size(listV,2);
39    nrec=nIt; n3d=nk; nReg=nRg;    nrec=nIt; n3d=nk; nReg=nRg;
# Line 42  if krd > 0, Line 44  if krd > 0,
44    [ntA(n),rList,tim,vv1,listV] = ...    [ntA(n),rList,tim,vv1,listV] = ...
45      read_StD(prefix,namA(n,:),listV);      read_StD(prefix,namA(n,:),listV);
46    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
47      vv1(find(vv1==undef))=NaN;
48    if (nrec < nIt),    if (nrec < nIt),
49      fprintf('\n');      fprintf('\n');
50      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.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22