| 1 |
gforget |
1.1 |
|
| 2 |
|
|
function []=profiles_prep_test_atlas(thresholdVAL); |
| 3 |
|
|
|
| 4 |
|
|
global rep_out filename_out; |
| 5 |
|
|
global t_std s_std t_test s_test t_w s_w t_equi s_equi; |
| 6 |
|
|
global ymd hms pnum lon lat direc dmod ilon ilat imonth; |
| 7 |
|
|
global fill_value_output initChkFile; |
| 8 |
|
|
|
| 9 |
|
|
global dmod_ref vec_lon vec_lat tmask3D; |
| 10 |
|
|
global T_climatology S_climatology T_climatology2 S_climatology2; |
| 11 |
|
|
global T_climatology3 S_climatology3 T_climatology4 S_climatology4; |
| 12 |
|
|
global T_weights S_weights; |
| 13 |
|
|
|
| 14 |
|
|
|
| 15 |
|
|
%cost for T : |
| 16 |
|
|
|
| 17 |
|
|
t_equi=squeeze(T_climatology(imonth,ilon,ilat,:))'; t_equi = interp1q(dmod_ref',t_equi',dmod')'; |
| 18 |
|
|
t_equi2=squeeze(T_climatology2(imonth,ilon,ilat,:))'; t_equi2 = interp1q(dmod_ref',t_equi2',dmod')'; |
| 19 |
|
|
t_equi3=squeeze(T_climatology3(imonth,ilon,ilat,:))'; t_equi3 = interp1q(dmod_ref',t_equi3',dmod')'; |
| 20 |
|
|
t_equi4=squeeze(T_climatology4(imonth,ilon,ilat,:))'; t_equi4 = interp1q(dmod_ref',t_equi4',dmod')'; |
| 21 |
|
|
|
| 22 |
|
|
t_cost=t_w.*((t_std-t_equi).^2); |
| 23 |
|
|
t_cost2=t_w.*((t_std-t_equi2).^2); |
| 24 |
|
|
t_cost3=t_w.*((t_std-t_equi3).^2); |
| 25 |
|
|
t_cost4=t_w.*((t_std-t_equi4).^2); |
| 26 |
|
|
|
| 27 |
|
|
%cost for S : |
| 28 |
|
|
|
| 29 |
|
|
s_equi=squeeze(S_climatology(imonth,ilon,ilat,:))'; s_equi = interp1q(dmod_ref',s_equi',dmod')'; |
| 30 |
|
|
s_equi2=squeeze(S_climatology2(imonth,ilon,ilat,:))'; s_equi2 = interp1q(dmod_ref',s_equi2',dmod')'; |
| 31 |
|
|
s_equi3=squeeze(S_climatology3(imonth,ilon,ilat,:))'; s_equi3 = interp1q(dmod_ref',s_equi3',dmod')'; |
| 32 |
|
|
s_equi4=squeeze(S_climatology4(imonth,ilon,ilat,:))'; s_equi4 = interp1q(dmod_ref',s_equi4',dmod')'; |
| 33 |
|
|
|
| 34 |
|
|
s_cost=s_w.*((s_std-s_equi).^2); |
| 35 |
|
|
s_cost2=s_w.*((s_std-s_equi2).^2); |
| 36 |
|
|
s_cost3=s_w.*((s_std-s_equi3).^2); |
| 37 |
|
|
s_cost4=s_w.*((s_std-s_equi4).^2); |
| 38 |
|
|
|
| 39 |
|
|
|
| 40 |
|
|
%remarks: 1) the climatology fields are assumed to be masked with NaN |
| 41 |
|
|
% and the next two lines does not do anything for cost=NaN |
| 42 |
|
|
% logic: we cannot test the data, this is no reason to reject observations |
| 43 |
|
|
% 2) the test must fail for all atlases for rejection |
| 44 |
|
|
% logic: an individual atlas may be wrong in some places, this is no reason to reject observations |
| 45 |
|
|
|
| 46 |
|
|
t_test(find( (t_test==0)&(t_std~=fill_value_output)&(t_cost>thresholdVAL)&(t_cost2>thresholdVAL)&(t_cost3>thresholdVAL)&(t_cost4>thresholdVAL) ))=5; |
| 47 |
|
|
s_test(find( (s_test==0)&(s_std~=fill_value_output)&(s_cost>thresholdVAL)&(s_cost2>thresholdVAL)&(s_cost3>thresholdVAL)&(s_cost4>thresholdVAL) ))=5; |
| 48 |
|
|
|
| 49 |
|
|
|
Parent Directory
| 
Revision Log
| 
Revision Graph