/[MITgcm]/MITgcm_contrib/gael/profilesMatlabProcessing/ecco_v4/insitu_diags.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/profilesMatlabProcessing/ecco_v4/insitu_diags.m

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

revision 1.3 by gforget, Sun Dec 30 00:35:53 2012 UTC revision 1.6 by gforget, Tue May 10 21:38:51 2016 UTC
# Line 5  function []=insitu_diags(dirMat,doComp,d Line 5  function []=insitu_diags(dirMat,doComp,d
5  %            doComp states whether to compute (1) or display (0)  %            doComp states whether to compute (1) or display (0)
6  %            dirTex is the directory where tex and fig files will be created  %            dirTex is the directory where tex and fig files will be created
7  %            nameTex is the tex file name (default : 'myPlots')  %            nameTex is the tex file name (default : 'myPlots')
8    %
9    %notes : MITprof files will be used from the myenv.profiles directory
10    %            or [dirMat '/profiles/output/'] is myenv.profiles is missing
11    
12    gcmfaces_global;
13    
14  dirMat=[dirMat '/'];  dirMat=[dirMat '/'];
15    
16  dirData=[dirMat '/profiles/output/'];  if isfield(myenv,'profiles');
17  while ~isdir(dirData);    dirData=myenv.profiles;
18    else;
19      dirData=[dirMat '/profiles/output/'];
20    end;
21    while ~isdir(dirData)&doComp;
22    fprintf(['directory : ' dirData '\n'])    fprintf(['directory : ' dirData '\n'])
23    dirData=input(['does not exist. Specify directory of nc file : \n']);    dirData=input(['does not exist. Specify directory of nc file : \n']);
24  end;  end;
25    
26    listData=dir([dirData '*.nc']);
27    
28  dirMat={'dirMat',dirMat};  dirMat={'dirMat',dirMat};
29  dirData={'dirData',dirData};  dirData={'dirData',dirData};
30    listData={'listData',listData(:).name};
31    
32  if isempty(who('dirTex'));  if isempty(who('dirTex'));
33    addToTex={'addToTex',0}; dirTex={'dirTex',''}; nameTex={'nameTex',''};    addToTex={'addToTex',0}; dirTex={'dirTex',''}; nameTex={'nameTex',''};
# Line 28  else; Line 40  else;
40  end;  end;
41    
42  if doComp;  if doComp;
43  insitu_misfit(1,dirData,dirMat);  fprintf('starting insitu_misfit\n'); clock
44  insitu_cost(1,dirData,dirMat);  insitu_misfit(1,dirData,dirMat,listData);
45    fprintf('starting insitu_cost\n'); clock
46    insitu_cost(1,dirData,dirMat,listData);
47    fprintf('done with insitu_diags\n'); clock
48  else;  else;
49  insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);  insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
50  insitu_cost(0,dirMat,addToTex,dirTex,nameTex);  insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22