OBSPROC
- OBSPROC runs on bluefire, using the following wrapper script: /mmm/users/auligne/ampsrt/da_run_suite_wrapper_amps1.ksh. The specific namelist parameters are:
export NL_QC_TEST_VERT_CONSISTENCY=false
export NL_QC_TEST_CONVECTIVE_ADJ=false
export MAX_OB_RANGE=3
- OBSPROC code is located on bluefire:/ptmp/auligne/WRFDA/var/obsproc/src. To compile "make clean", then "make", then link the executable to /ptmp/auligne/WRFDA/var/obsproc/3dvar_obs.exe
- The code has been modified to read U & V and write into T_error and Q_error (module_decoded.F90 and module_err_afwa.F90)
- Polarday/night is stored in T Quality Flag
- numCloud is stored in Z Quality Flag
WRF-Var
Code located on bluefire:/ptmp/auligne/WRFDA/. Modifications to:
- Avoid conversion from RH to Q in da_obs_proc_station.inc (since we now read Q directly) and conversion of RH_error to Q_error in da_fill_obs_structures.inc
- Write Q in g/kg (instead of kg/kg) in gts_omb_oma_01 (da_write_obs.inc)
- Calculate surface type for each AIRSRET observation (call da_detsurtyp.inc in da_get_innov_vector_airsr.inc). Put surface type information in T Quality Control Flag according to the formula:
QI = 10*ISFLG + DAYNIGHT, with DAYNIGHT=0 for night and 4 for day and ISFLG is the surface flag:
- 0 sea
- 1 sea ice
- 2 land
- 3 snow
- 4 mixed predominately sea
- 5 mixed predominately sea ice
- 6 mixed predominately land
- 7 mixed predominately snow
NCL plotting scripts
Modified to read and compute AIRSRET, located on bluefire:/blhome/auligne/monitoring:
- plot_ob_ascii_loc.ncl represents the map of observation locations from ob.ascii (output from OBSPROC)
- plot_gts_omb_oma.ncl shows maps and scatterplots of OMB and OMA for various vertical levels