profilesMatlabProcessing/profiles_devel/MITprof_bin.m



doPart1=0;
doDisplay=1;


if doPart1;
    
    %load paths and grid:
    %--------------------
    startup_gcmfaces_and_mitprof;
    
    %load data set:
    %--------------
    
    dirData='./MITprof/profiles_devel/';
    fileData='seals_latest_MITprof.nc';
    %fileData='argo_atlantic_MITprof_latest.nc';
    %fileData='argo_indian_MITprof_latest.nc';
    MITprof=MITprof_load([dirData fileData]);
    
    %assign nameData:
    %----------------
    
    ii=strfind(fileData,'_MITprof_latest.nc');
    if ~isempty(ii);
        nameData=fileData(1:ii-1);
    else;
        ii=strfind(fileData,'_MITprof.nc');
        if ~isempty(ii);
            nameData=fileData(1:ii-1);
        else;
            ii=strfind(fileData,'.nc');
            nameData=fileData(1:ii-1);
        end;
    end;
    
    eval(['MITprof_' nameData '=MITprof;']);
    
    
end;

kk=25;
lon=MITprof.prof_lon; lat=MITprof.prof_lat;
obs=MITprof.prof_T(:,kk); estim=MITprof.prof_Testim(:,kk);

ii=find(~isnan(obs.*estim));
lon=lon(ii); lat=lat(ii);
obs=obs(ii); estim=estim(ii);

gcmfaces_bindata;
tmp_point=gcmfaces_bindata(lon,lat);
[tmp_sum,tmp_nb]=gcmfaces_bindata(lon,lat,obs);
[tmp_sum2,tmp_nb]=gcmfaces_bindata(lon,lat,obs.^2);
%[tmp_]=gcmfaces_bindata(MITprofur.prof_lon,MITprofur.prof_lat);

tmp_var=(tmp_sum2-tmp_sum.^2./tmp_nb)./tmp_nb; tmp_var(find(tmp_nb==0|tmp_var<1e-3))=NaN;
%tmp1=convert2array(tmp_var); tmp2=log10(tmp1); tmp_var_log=convert2array(tmp2,tmp_var);
tmp_var_log=log10(tmp_var); 
figure; m_map_gcmfaces(tmp_var_log,0,[-2 1]);


1	gforget	1.1
2
3			doPart1=0;
4			doDisplay=1;
5
6
7			if doPart1;
8
9			%load paths and grid:
10			%--------------------
11			startup_gcmfaces_and_mitprof;
12
13			%load data set:
14			%--------------
15
16			dirData='./MITprof/profiles_devel/';
17			fileData='seals_latest_MITprof.nc';
18			%fileData='argo_atlantic_MITprof_latest.nc';
19			%fileData='argo_indian_MITprof_latest.nc';
20			MITprof=MITprof_load([dirData fileData]);
21
22			%assign nameData:
23			%----------------
24
25			ii=strfind(fileData,'_MITprof_latest.nc');
26			if ~isempty(ii);
27			nameData=fileData(1:ii-1);
28			else;
29			ii=strfind(fileData,'_MITprof.nc');
30			if ~isempty(ii);
31			nameData=fileData(1:ii-1);
32			else;
33			ii=strfind(fileData,'.nc');
34			nameData=fileData(1:ii-1);
35			end;
36			end;
37
38			eval(['MITprof_' nameData '=MITprof;']);
39
40
41			end;
42
43			kk=25;
44			lon=MITprof.prof_lon; lat=MITprof.prof_lat;
45			obs=MITprof.prof_T(:,kk); estim=MITprof.prof_Testim(:,kk);
46
47			ii=find(~isnan(obs.*estim));
48			lon=lon(ii); lat=lat(ii);
49			obs=obs(ii); estim=estim(ii);
50
51			gcmfaces_bindata;
52			tmp_point=gcmfaces_bindata(lon,lat);
53			[tmp_sum,tmp_nb]=gcmfaces_bindata(lon,lat,obs);
54			[tmp_sum2,tmp_nb]=gcmfaces_bindata(lon,lat,obs.^2);
55			%[tmp_]=gcmfaces_bindata(MITprofur.prof_lon,MITprofur.prof_lat);
56
57			tmp_var=(tmp_sum2-tmp_sum.^2./tmp_nb)./tmp_nb; tmp_var(find(tmp_nb==0\|tmp_var<1e-3))=NaN;
58			%tmp1=convert2array(tmp_var); tmp2=log10(tmp1); tmp_var_log=convert2array(tmp2,tmp_var);
59			tmp_var_log=log10(tmp_var);
60			figure; m_map_gcmfaces(tmp_var_log,0,[-2 1]);
61
62