AIRS Retrieval V5 decoder

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
  • No labels