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]);