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

Diff of /MITgcm_contrib/jmc_script/read_StD.m

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

revision 1.1 by jmc, Thu Mar 1 04:35:35 2007 UTC revision 1.2 by jmc, Thu May 29 23:47:13 2008 UTC
# Line 67  end Line 67  end
67  fclose(fid);  fclose(fid);
68  if rf > 0,  if rf > 0,
69  %- rename fields (consistent with script "extract_StD"):  %- rename fields (consistent with script "extract_StD"):
70     listV=strrep(listV,'ETAN', 'Eta'); listV=strrep(listV,'ETANSQ','Et2');     for j=1:rf,
71     listV=strrep(listV,'THETA' ,'T');  listV=strrep(listV,'SALT'  ,'S');       var1=char(listV(j));
72     listV=strrep(listV,'UVEL'  ,'U');  listV=strrep(listV,'VVEL'  ,'V');       switch var1
73     listV=strrep(listV,'WVEL'  ,'W');  listV=strrep(listV,'PHIHYD','Phi');       case 'ETAN'    ,  var2='Eta';
74     listV=strrep(listV,'UVELSQ','U2'); listV=strrep(listV,'VVELSQ','V2');       case 'ETANSQ'  ,  var2='Et2';
75         case 'THETA'   ,  var2='T';
76         case 'SALT'    ,  var2='S';
77         case 'UVEL'    ,  var2='U';
78         case 'VVEL'    ,  var2='V';
79         case 'WVEL'    ,  var2='W';
80         case 'PHIHYD'  ,  var2='Phi';
81         case 'UVELSQ'  ,  var2='U2';
82         case 'VVELSQ'  ,  var2='V2';
83         case 'WVELSQ'  ,  var2='W2';
84         case 'THETASQ' ,  var2='T2';
85         otherwise     var2=var1;
86         end
87         listV(j)=cellstr(var2);
88         if strcmp(var1,var2), fprintf(' %s\n',var2);
89              else fprintf(' %s --> %s\n',var1,var2); end
90       end
91    %  listV
92  end  end
93  if flag ~= 0 | k > 0, frq,rList,kList, end  if flag ~= 0 | k > 0, frq,rList,kList, end
94  if flag ~= 0, error(['not normal end after reading ',int2str(l),' lines']); end  if flag ~= 0, error(['not normal end after reading ',int2str(l),' lines']); end
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22