function []=profiles_prep_test_profiler2(); global rep_out filename_out; global t_std s_std t_test s_test t_w s_w t_equi s_equi; global ymd hms pnum pnum_txt lon lat direc dmod ilon ilat imonth; global fill_value_output initChkFile; %general remarks : % -most of the decisions were taken by looking at deep values of S, % expected not to vary strongly, but some suggestivity remains. % -the point in having these tests is that they will not be affected % by evolutions of the climatology and weights (compare with test_cost). % %additional remarks : % -4900221 : first profiles might be ok. % -the biases identified could be double checked. % -51886 : deep S is strange and detected by the cost. % -4900388 : not fully understood. % -15855 : this number is associated with profiles all over the globe. % -39037 : not fully understood. % -2900335 : irregular last profiles. %for the test of the date : tmp_year=floor(ymd*1e-4); tmp_month=floor((ymd-tmp_year*1e4)*1e-2); tmp_day=ymd-tmp_year*1e4-tmp_month*1e2; tmp_jul=jul_0h([tmp_year tmp_month tmp_day 0 0 0]); % pesky profiler 15855 - this shows up multiple times in PACIFIC files % lat=-5.8940 lon=6.978 ydm=19970820 (location is ATLANTIC) tmp0=jul_0h([1997 8 20 0 0 0]); if (pnum==15855 )&(tmp_jul==tmp0);s_test(:)=3;t_test(:)=3;end; %exemple of a dying profiler : tmp0=jul_0h([2002 4 1 0 0 0]); if (pnum==29011)&(tmp_jul>tmp0); s_test(:)=3; end; %exemple of a strong drift : if (pnum==5900198); s_test(:)=3; end; %a profiler that says "30" : if (pnum==21851); s_test(:)=3; end; %two profilers that say "33.5" : if (pnum==2900236); s_test(:)=3; end; if (pnum==2900237); s_test(:)=3; end; %instrument flagged after a couple of profiles. if (pnum==5900853); s_test(:)=3; end; if (pnum==4900221); s_test(:)=3; end; %exemple of a profiler that stays silent for some time... and then provides crap. tmp0=jul_0h([2003 1 1 0 0 0]); if (pnum==1900141)&(tmp_jul>tmp0); s_test(:)=3; end; tmp0=jul_0h([2004 10 1 0 0 0]); if (pnum==1900141)&(tmp_jul>tmp0); t_test(:)=3; end; %exemple of a profiler whose last unflagged profiles are boubtful : tmp0=jul_0h([2003 12 1 0 0 0]); if (pnum==5900382)&(tmp_jul>tmp0); s_test(:)=3; t_test(:)=3; end; %profilers that drift : tmp0=jul_0h([2001 06 1 0 0 0]); if (pnum==29036)&(tmp_jul>tmp0); s_test(:)=3; end; tmp0=jul_0h([2004 6 16 0 0 0]); if (pnum==4900388)&(tmp_jul>tmp0); s_test(:)=3; t_test(:)=3;end; %tmp0=jul_0h([2004 2 1 0 0 0]); %if (pnum==1900067)&(tmp_jul>tmp0); s_test(:)=3; end; tmp0=jul_0h([2002 5 1 0 0 0]); if (pnum==3900137)&(tmp_jul>tmp0); s_test(:)=3; end; tmp0=jul_0h([2002 8 1 0 0 0]); if (pnum==3900132)&(tmp_jul>tmp0); s_test(:)=3; end; %tmp0=jul_0h([]); %if (pnum==)&(tmp_jul>tmp0); s_test(:)=3; end; %profilers that stay silent for some time... and then provide crap. %tmp0=jul_0h([2003 11 1 0 0 0]); %if (pnum==3900116)&(tmp_jul>tmp0); s_test(:)=3; t_test(:)=3; end; tmp0=jul_0h([2004 6 1 0 0 0]); if (pnum==5900311)&(tmp_jul>tmp0); s_test(:)=3; t_test(:)=3; end; tmp0=jul_0h([2004 9 1 0 0 0]); if (pnum==1900306)&(tmp_jul>tmp0); s_test(:)=3; t_test(:)=3; end; %tmp0=jul_0h([2004 2 1 0 0 0]); %if (pnum==69027)&(tmp_jul>tmp0); s_test(:)=3; end; tmp0=jul_0h([2001 11 27 0 0 0]); if (pnum==6900197)&(tmp_jul>tmp0); s_test(:)=3; end; %t_test(:)=3; end; %profilers whose last unflagged profiles are boubtful : tmp0=jul_0h([2004 11 1 0 0 0]); if (pnum==5900391)&(tmp_jul>tmp0); s_test(:)=3; t_test(:)=3; end; tmp0=jul_0h([2004 8 10 0 0 0]); if (pnum==4900391)&(tmp_jul>tmp0); s_test(:)=3; t_test(:)=3; end; tmp0=jul_0h([2005 3 1 0 0 0]); if (pnum==1900203)&(tmp_jul>tmp0); s_test(:)=3; end; tmp0=jul_0h([2004 12 15 0 0 0]); if (pnum==3900149)&(tmp_jul>tmp0); s_test(:)=3; end; tmp0=jul_0h([2004 8 1 0 0 0]); if (pnum==1900278)&(tmp_jul>tmp0); s_test(:)=3; t_test(:)=3; end; %profilers judged biased : typically S very different from climatology at 1500-2000m if (pnum==1900425); s_test(:)=3; end; if (pnum==1900424); s_test(:)=3; end; if (pnum==1900447); s_test(:)=3; end; tmp0=jul_0h([2003 6 1 0 0 0]); if (pnum==5900180)&(tmp_jul>tmp0); s_test(:)=3; end; %if (pnum==69027); s_test(:)=3; end; %...and more anecdotic : if (pnum==1900208); s_test(:)=3; end; if (pnum==1900498); s_test(:)=3; end; if (pnum==1900139); s_test(:)=3; end; %particular case : profiler with way too large S at depth... % and the opposite around 300m (!?) if (pnum==5900319); s_test(:)=3; end;