function []=insitu_cost(dirData,doComp,varargin); %[]=insitu_cost(dirData,listData,listVar,varargin); %object: computes or displays cost function statistics %input: dirData is the data main directory name % doComp states whether to compute & save (1) or load & display (0) %optional must take the following form {'name',param1,param2,...} Active ones are % {'subdirData',subdirData} is the data subdirectory name ('profiles/output/' by default) % {'listData',data1,data2,...} where data1 e.g. is 'argo_indian.nc' % {'listVar',var1,var2} where var1 and var2 are 'T' and 'S' % {'years',year0,year1} to limit computation to [year0 year1] interval % {'subdirMat',subdirMat} where subdirMat is the output subdirectory ('mat/' by default) % {'suffMat',suffMat} is the mat file suffix ('all' by default) % {'subdirTex',subdirTex} where subdirTex is the tex directory ('tex/' by default) % {'addToTex',addToTex} states whether (1) or not (0 default) to % augment an ongoing tex file (see write2tex.m, subdirTex) % %example: insitu_cost('./',1,{'listData','argo_in*'}); %load paths and grid: %-------------------- global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end; %set more optional paramaters to default values year0=1992; year1=2008; %listData={'WOD09_XBT*','WOD09_CTD*','argo_in*','argo_at*','argo_pa*','seals*'}; %listVar={'T','S'}; listData={'argo_in*'}; listVar={'T'}; subdirData='/profiles/output/'; subdirMat='mat/'; suffMat='_all'; subdirTex='tex/'; addToTex=0; %set more optional paramaters to user defined values for ii=1:nargin-2; if strcmp(varargin{ii}{1},'listData'); listData={varargin{ii}{2:end}}; elseif strcmp(varargin{ii}{1},'listVar'); listVar={varargin{ii}{2:end}}; elseif strcmp(varargin{ii}{1},'years'); year0=varargin{ii}{2}; year1=varargin{ii}{3}; elseif strcmp(varargin{ii}{1},'subdirMat'); subdirMat=varargin{ii}{2}; elseif strcmp(varargin{ii}{1},'suffMat'); suffMat=['_' varargin{ii}{2}]; elseif strcmp(varargin{ii}{1},'subdirTex'); subdirTex=varargin{ii}{2}; elseif strcmp(varargin{ii}{1},'addToTex'); addToTex=varargin{ii}{2}; elseif strcmp(varargin{ii}{1},'subdirData'); subdirData=varargin{ii}{2}; else; warning('inputCheck:insitu_cost',... ['unknown option ''' varargin{ii}{1} ''' was ignored']); end; end; if doComp; %time limits: date0=datenum(year0,1,1); date1=datenum(year1,12,31); for vv=1:length(listVar); varCur=listVar{vv}; %get the cost square root: [MITprof]=MITprof_stats_load([dirData '/' subdirData],listData,varCur); %mask out values that are not in year range: ii=find(MITprof.prof_datedate1); MITprof.prof(ii,:)=NaN; %mean and median: costCur=[nanmean(MITprof.prof(:).^2) nanmedian(MITprof.prof(:).^2)]; fprintf('mean cost for %s: %0.3g (median: %0.3g) \n',varCur,costCur); eval(['costAve' varCur '.mean=costCur(1);']); eval(['costAve' varCur '.median=costCur(2);']); %depth/time mean cost: tmp_date=(MITprof.prof_date-date0)/365; % [x,y,z,n]=MITprof_stats(MITprof.prof_depth,MITprof.prof_depth,... [x,y,z,n]=MITprof_stats(MITprof.prof_depth,[0:200:1000 1500:500:6000],... tmp_date,[0:1/4:1+(year1-year0)],'mean',MITprof.prof.^2); z(n==0)=NaN; n(n==0)=NaN; eval(['depthTimeCost' varCur '.x=x;']); eval(['depthTimeCost' varCur '.y=y;']); eval(['depthTimeCost' varCur '.z=z;']); eval(['depthTimeCost' varCur '.n=n;']); %latitudinal distribution: [x,y,z,n]=MITprof_stats(MITprof.prof_lat,[-90:5:90],MITprof.prof,[-5:0.25:5]); z(n==0)=NaN; n(n==0)=NaN; eval(['misfitDistrib' varCur '.x=x;']); eval(['misfitDistrib' varCur '.y=y;']); eval(['misfitDistrib' varCur '.z=z;']); eval(['misfitDistrib' varCur '.n=n;']); end; eval(['save ' dirData '/' subdirMat 'insitu_cost' suffMat '.mat costAve* misfitDistrib* depthTimeCost* year* listData listVar;']); else;%display result eval(['load ' dirData '/' subdirMat 'insitu_cost' suffMat '.mat;']); figureL; for vv=1:length(listVar); varCur=listVar{vv}; eval(['costCur=[costAve' varCur '.mean costAve' varCur '.median];']); eval(['x=depthTimeCost' varCur '.x;']); eval(['y=depthTimeCost' varCur '.y;']); eval(['z=depthTimeCost' varCur '.z;']); eval(['n=depthTimeCost' varCur '.n;']); z(n<1e2)=NaN; subplot(2,length(listVar),vv); depthStretchPlot('pcolor',{y,x,z},[0:200:1000 1500:500:6000],[0 1000 4000]); caxis([0 8]); shading flat; colorbar; ylabel('depth (in m)'); xlabel(['date (in years since ' num2str(year0) ')']); title(sprintf('mean cost for %s: %0.3g (median: %0.3g)',varCur,costCur)); %latitudinal distribution: eval(['x=misfitDistrib' varCur '.x;']); eval(['y=misfitDistrib' varCur '.y;']); eval(['z=misfitDistrib' varCur '.z;']); eval(['n=misfitDistrib' varCur '.n;']); tmp1=(nansum(n,2)>1e4)*ones(1,size(n,2)); z(tmp1==0)=NaN; subplot(2,length(listVar),vv+length(listVar)); contourf(x,y,z,[0:0.05:0.5]); grid on; caxis([0 0.5]); colorbar; ylabel('normalized misfit'); xlabel('latitude'); title('pdf'); end; myCaption={'in situ profiles cost function'}; if addToTex; write2tex([dirData '/' subdirTex '/myPlots.tex'],2,myCaption,gcf); end; end;