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
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