| 1 |
function []=insitu_misfit(doComp,varargin); |
| 2 |
%object: computes or displays cost function statistics |
| 3 |
%input: doComp states whether to compute & save (1) or load & display (0) |
| 4 |
%optional must take the following form {'name',param1,param2,...} Active ones are |
| 5 |
% {'dirData',dirData} is the data directory name ('mat/profiles/output/' by default) |
| 6 |
% {'listData',data1,data2,...} where data1 e.g. is 'argo_indian.nc' |
| 7 |
% {'listVar',var1,var2} where var1 and var2 are 'T' and 'S' |
| 8 |
% {'years',year0,year1} to limit computation to [year0 year1] interval |
| 9 |
% {'dirMat',dirMat} where dirMat is the output directory ('mat/' by default) |
| 10 |
% {'suffMat',suffMat} is the mat file suffix ('all' by default) |
| 11 |
% {'dirTex',dirTex} where dirTex is the tex directory ('tex/' by default) |
| 12 |
% {'nameTex',nameTex} is the tex file name ('myPlots' by default) |
| 13 |
% {'addToTex',addToTex} states whether (1) or not (0 default) to |
| 14 |
% augment an ongoing tex file (see write2tex.m, dirTex) |
| 15 |
% |
| 16 |
%example: insitu_misfit('./',1,{'listData','argo_in*'}); |
| 17 |
|
| 18 |
vecDepth=[[75 125];[250 350];[900 1100]]; |
| 19 |
vecCaxis=[[3 3 1];[0.5 0.5 0.2]]; |
| 20 |
%if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end; |
| 21 |
|
| 22 |
gcmfaces_global; |
| 23 |
global myparms; |
| 24 |
|
| 25 |
listVar={'T','S'}; |
| 26 |
%shorter test case: listData={'argo_in*'}; listVar={'T'}; |
| 27 |
dirData='mat/profiles/output/'; |
| 28 |
dirMat='mat/'; suffMat='_all'; |
| 29 |
dirTex='tex/'; addToTex=0; nameTex='myPlots'; |
| 30 |
%set more optional paramaters to user defined values |
| 31 |
for ii=1:nargin-1; |
| 32 |
if strcmp(varargin{ii}{1},'listData'); listData={varargin{ii}{2:end}}; |
| 33 |
elseif strcmp(varargin{ii}{1},'listVar'); listVar={varargin{ii}{2:end}}; |
| 34 |
elseif strcmp(varargin{ii}{1},'years'); year0=varargin{ii}{2}; year1=varargin{ii}{3}; |
| 35 |
elseif strcmp(varargin{ii}{1},'dirMat'); dirMat=varargin{ii}{2}; |
| 36 |
elseif strcmp(varargin{ii}{1},'suffMat'); suffMat=['_' varargin{ii}{2}]; |
| 37 |
elseif strcmp(varargin{ii}{1},'dirTex'); dirTex=varargin{ii}{2}; |
| 38 |
elseif strcmp(varargin{ii}{1},'addToTex'); addToTex=varargin{ii}{2}; |
| 39 |
elseif strcmp(varargin{ii}{1},'nameTex'); nameTex=varargin{ii}{2}; |
| 40 |
elseif strcmp(varargin{ii}{1},'dirData'); dirData=varargin{ii}{2}; |
| 41 |
else; |
| 42 |
warning('inputCheck:insitu_misfit',... |
| 43 |
['unknown option ''' varargin{ii}{1} ''' was ignored']); |
| 44 |
end; |
| 45 |
end; |
| 46 |
|
| 47 |
if isempty(myparms)|isempty(mygrid); |
| 48 |
load([dirMat 'diags_grid_parms.mat']); |
| 49 |
end; |
| 50 |
|
| 51 |
if isempty(whos('year0')); |
| 52 |
year0=myparms.yearInAve(1); year1=myparms.yearInAve(2); |
| 53 |
end; |
| 54 |
|
| 55 |
if isempty(whos('listData')); |
| 56 |
listData=dir([dirMat 'profiles/output/*.nc']); |
| 57 |
listData={listData(:).name}; |
| 58 |
for ff=1:length(listData); listData{ff}=[listData{ff}(1:end-3) '*']; end; |
| 59 |
end; |
| 60 |
|
| 61 |
if doComp; |
| 62 |
|
| 63 |
%initialize delaunay triangulation: |
| 64 |
gcmfaces_bindata; |
| 65 |
%time limits: |
| 66 |
date0=datenum(year0,1,1); date1=datenum(year1,12,31); |
| 67 |
%number of levels: |
| 68 |
nk=size(vecDepth,1); |
| 69 |
|
| 70 |
for vv=1:length(listVar); |
| 71 |
varCur=listVar{vv}; |
| 72 |
%get the misfits (here not normalized): |
| 73 |
[MITprof]=MITprof_stats_load(dirData,listData,varCur,1); |
| 74 |
%mask out values that are not in year range: |
| 75 |
ii=find(MITprof.prof_date<date0|MITprof.prof_date>date1); |
| 76 |
MITprof.prof(ii,:)=NaN; |
| 77 |
%misfit maps: |
| 78 |
for kkk=1:size(vecDepth,1); |
| 79 |
kk=find(MITprof.prof_depth>=vecDepth(kkk,1)&MITprof.prof_depth<=vecDepth(kkk,2)); |
| 80 |
lon=MITprof.prof_lon; lat=MITprof.prof_lat; obs=nanmean(MITprof.prof(:,kk),2); |
| 81 |
ii=find(~isnan(obs)); lon=lon(ii); lat=lat(ii); obs=obs(ii); |
| 82 |
misfit_map=gcmfaces_bindata(lon,lat,obs); |
| 83 |
eval(['misfit_map_' varCur num2str(kkk) '=misfit_map;']); |
| 84 |
end; |
| 85 |
end; |
| 86 |
|
| 87 |
if ~isdir([dirMat 'cost/']); mkdir([dirMat 'cost/']); end; |
| 88 |
eval(['save ' dirMat 'cost/insitu_misfit' suffMat '.mat misfit_map_* vecDepth year* listData listVar;']); |
| 89 |
|
| 90 |
else; |
| 91 |
|
| 92 |
if isdir([dirMat 'cost/']); dirMat=[dirMat 'cost/']; end; |
| 93 |
|
| 94 |
eval(['load ' dirMat 'insitu_misfit' suffMat '.mat;']); |
| 95 |
|
| 96 |
nk=size(vecDepth,1); |
| 97 |
figureL; |
| 98 |
for vv=1:length(listVar); |
| 99 |
varCur=listVar{vv}; |
| 100 |
for kkk=1:size(vecDepth,1); |
| 101 |
eval(['misfit_map=misfit_map_' varCur num2str(kkk) ';']); |
| 102 |
dep=mean(vecDepth(kkk,:)); |
| 103 |
|
| 104 |
if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end; |
| 105 |
subplot(nk,length(listVar),vv+(kkk-1)*length(listVar)); |
| 106 |
m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc},{'myShading','flat'}); |
| 107 |
title(sprintf('%s misfits at %dm (%d to %d)',varCur,dep,year0,year1)); |
| 108 |
end; |
| 109 |
end; |
| 110 |
|
| 111 |
myCaption={'Time mean misfit (model-data) for in situ profiles, ',... |
| 112 |
'at various depths (rows), for T (left; in K) and S (right; in psu).'}; |
| 113 |
if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end; |
| 114 |
|
| 115 |
end; |
| 116 |
|
| 117 |
|
| 118 |
|
Parent Directory
| 
Revision Log
| 
Revision Graph