170 likes | 319 Views
Level 2 Ocean Salinity Processor. 8 - 9 July 2009. Paul Spurgeon Mike McCulloch. New processor 03_07. Operational & prototype still in sync Removed old MIRE Many changes & bug fixes... To be used in DPGS tests next week Updated documentation: OPSRD, PPSRD, ATBD, IODD, DPM.
E N D
Level 2 Ocean Salinity Processor 8 - 9 July 2009 Paul Spurgeon Mike McCulloch
New processor 03_07 • Operational & prototype still in sync • Removed old MIRE • Many changes & bug fixes... • To be used in DPGS tests next week • Updated documentation: OPSRD, PPSRD, ATBD, IODD, DPM...
New default values • -999 • New versions of SDV & smos-box know this value – consistent with L2 SM • May need to adjust MATLAB scripts?
No poor quality grid points retrieved in 03_06 • Classified as poor quality if • fg_ctrl_many_outliers | fg_ctrl_sunglint | fg_ctrl_moonglint | fg_ctrl_gal_noise | fg_ctrl_gal_noise_pol | fg_ctrl_num_meas_low | fg_sc_rain | fg_sc_TEC_gradient
Poor quality failed if • fg_ctrl_range | fg_ctrl_sigma | fg_ctrl_chi2 | fg_ctrl_chi2_P | fg_OoR_Gam1_dim1_RA | fg_OoR_Gam1_dim2_Dec | fg_OoR_Sunglint_dim2_ThetaSun | fg_OoR_Sunglint_dim3_Phi | fg_OoR_Sunglint_dim4_Theta | fg_OoR_Sunglint_dim5_WS | fg_OoR_Foam_dim1_WS | fg_OoR_Foam_dim2_TseaAir | fg_OoR_Foam_dim3_SSS | fg_OoR_Foam_dim4_SST | fg_OoR_Foam_dim5_Theta | fg_OoR_Gam2_Dec | fg_OoR_Gam2_Ra | fg_OoR_Gam2_WSn | fg_OoR_Gam2_Theta | fg_OoR_Gam2_Psi
Poor quality failed ignores • fg_OoR_Rough1_dim1_SST • fg_OoR_Rough1_dim2_SSS • fg_OoR_Rough1_dim3_U10 • fg_OoR_Rough1_dim4_Theta • fg_OoR_Rough2_dim5_SST • fg_OoR_Rough2_dim4_SSS • fg_OoR_Rough2_dim3_Theta • fg_OoR_Rough2_dim1_Ust • fg_OoR_Rough2_dim2_Omega
Detecting sea-ice • SPH Sea_Ice is count of fg_sc_suspect_ice • fg_sc_suspect_ice if (Dg_Suspect_ice)/nMeas > Tg_suspect_ice //default 50% • Dg_Suspect_ice = count(fm_suspect_ice) • fm_suspect_ice if (tb > (tbFwdModel + Tm_DT_ice ) //default 20°
Detecting sea-ice old MIRE, Dg_Suspect_ice
Other problems fixed • Cleaned up and corrected GUI breakpoint list • Corrected tests for missing ECMWF data: • Sea_Ice_Cover needed by all config • Drag_Coefficient_With_Waves needed by roughness model 1 • Minor changes to AUX_CNFL2OS for DPGS • removed L2_Creator • added List_of_Hosts: new definition of HW_Identifier
Other problems fixed • Now all %: • Tg_gal_noise_max, Tg_gal_noise_pol_max, Tg_moonglint_max, Tg_sun_glint_max • Flag Tg_sunglint_max ,Tg_suspect_ice, Tg_lat_ice_Acard can now be modified in GUI • Fg_ctrl_ECMWF = 0 means missing ECMWF data (corrected in ATBD)
Other problems fixed • Tg_sunglint_max%, Tg_moonglint_max%, Tg_gal_noise_max% - all 10% (in TGRD) • Dg_Suspect_ice, Dg_galactic_Noise_Error, Dg_galactic_Noise_Pol, Dg_quality_Acard now written to UDP
Sun glint • Stopped fml1c_sun_glint_area setting fm_l1c_sun (otherwise measurement is discarded)
Moon glint Scenario 20
Missing ECMWF Wind_10m_Wave_Model
Tools Processors Readers filter csv csv L1C ECMWF L2OS UDP L2OS DAP csv csv csv csv csv merge csv csv count csv stats csv XML
Tools sample script • extract='ID,Latitude,Longitude,Sea_Surface_Temperature,Drag_Coefficient_With_Waves,Wind_10m_Wave_Model' • test='Drag_Coefficient_With_Waves < 0 & Sea_Surface_Temperature > 274 & Latitude > -40 & Latitude < -30' • readECMWF -v3 -file"$ecmwf" -extract"$extract" | filter -test"$test" > e.txt • extract='ID,SST,Control_Flags_1.Fg_ctrl_selected,Science_Flags_1.Fg_sc_suspect_ice' • test='Control_Flags_1.Fg_ctrl_selected = 1' • readL2OSUDP -file"$udp" -extract"$extract" | filter -test"$test" > u.txt • merge -files"e.txt,u.txt" -field"ID" > m.txt • echo “Grid point count =" `count < m.txt`
Tools sample script result • Item 177475 • ID 2254687 • Latitude -25.139 • Longitude 17.736 • Sea_Surface_Temperature 298.158 • Drag_Coefficient_With_Waves -99998 • Wind_10m_Wave_Model -99998 • SST 25.008 • Control_Flags_1.Fg_ctrl_selected 1 • Science_Flags_1.Fg_sc_suspect_ice 0