[Dart-dev] DART/branches Revision: 11158

[dart at ucar.edu]

mizzi at ucar.edu
2017-02-27 15:12:17 -0700 (Mon, 27 Feb 2017)
33
Update WRF-Chem/DART repository




Modified: DART/branches/mizzi/models/wrf_chem/WRF_DART_utilities/wrf_dart_obs_preprocess.f90
===================================================================
--- DART/branches/mizzi/models/wrf_chem/WRF_DART_utilities/wrf_dart_obs_preprocess.f90	2017-02-27 20:31:07 UTC (rev 11157)
+++ DART/branches/mizzi/models/wrf_chem/WRF_DART_utilities/wrf_dart_obs_preprocess.f90	2017-02-27 22:12:17 UTC (rev 11158)
@@ -42,6 +42,7 @@
                              MOPITT_CO_RETRIEVAL,  &
                              IASI_CO_RETRIEVAL, &
                              IASI_O3_RETRIEVAL, &
+                             AIRNOW_CO, &
                              AIRNOW_O3
 ! APM/JB ---
 use time_manager_mod, only : time_type, set_calendar_type, GREGORIAN, set_time
@@ -71,6 +72,7 @@
                       mopitt_co_extra    = 'obs_seq.mopitt_co',  &
                       iasi_co_extra      = 'obs_seq.iasi_co',    &
                       iasi_o3_extra      = 'obs_seq.iasi_o3',    &
+                      airnow_co_extra    = 'obs_seq.airnow_co',  &
                       airnow_o3_extra    = 'obs_seq.airnow_o3'
 character(len=80)  :: name, sgday, sgsec
 ! APM/JB ---
@@ -125,6 +127,11 @@
 real(r8)           :: iasi_o3_horiz_int           = 100.0_r8   ! horizontal interval for super-ob
 real(r8)           :: iasi_o3_pres_int            = 2500.0_r8  ! pressure interval for super-ob
 logical            :: overwrite_ncep_iasi_o3_qc   = .false.    ! true to overwrite NCEP QC (see instructions)
+!  AIRNOW CO specific parameters
+logical            :: superob_airnow_co           = .false.    ! super-ob sat wind data
+real(r8)           :: airnow_co_horiz_int         = 100.0_r8   ! horizontal interval for super-ob
+real(r8)           :: airnow_co_pres_int          = 2500.0_r8  ! pressure interval for super-ob
+logical            :: overwrite_ncep_airnow_co_qc = .false.    ! true to overwrite NCEP QC (see instructions)
 !  AIRNOW O3 specific parameters
 logical            :: superob_airnow_o3           = .false.    ! super-ob sat wind data
 real(r8)           :: airnow_o3_horiz_int         = 100.0_r8   ! horizontal interval for super-ob
@@ -149,6 +156,7 @@
          superob_mopitt_co, mopitt_co_pres_int, mopitt_co_extra, mopitt_co_horiz_int, &
          superob_iasi_co, iasi_co_pres_int, iasi_co_extra, iasi_co_horiz_int, &
          superob_iasi_o3, iasi_o3_pres_int, iasi_o3_extra, iasi_o3_horiz_int, &
+         superob_airnow_co, airnow_co_pres_int, airnow_co_extra, airnow_co_horiz_int, &
          superob_airnow_o3, airnow_o3_pres_int, airnow_o3_extra, airnow_o3_horiz_int
 ! APM/JB ---
 
@@ -169,7 +177,7 @@
                            seq_prof, seq_tc, seq_gpsro, seq_other, &
 ! APM/JB +++
                            seq_modis_aod, seq_mopitt_co, seq_iasi_co, seq_iasi_o3, &
-                           seq_airnow_o3
+                           seq_airnow_co, seq_airnow_o3
 ! APM/JB ---
 
 type(time_type)         :: anal_time
@@ -248,6 +256,7 @@
 call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_mopitt_co)
 call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_iasi_co)
 call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_iasi_o3)
+call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_airnow_co)
 call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_airnow_o3)
 ! APM/JB ---
 call create_new_obs_seq(num_copies, num_qc, max_obs_seq, seq_other)
@@ -259,7 +268,7 @@
 sfc_elevation_tol, overwrite_ncep_sfc_qc, overwrite_ncep_satwnd_qc, &
 overwrite_obs_time, anal_time, seq_rawin, seq_sfc, seq_acars, seq_satwnd, & 
 seq_tc, seq_gpsro, seq_modis_aod, seq_mopitt_co, seq_iasi_co, seq_iasi_o3, &
-seq_airnow_o3, seq_other)
+seq_airnow_co, seq_airnow_o3, seq_other)
 ! APM/JB ---
 
 !  add supplimental rawinsonde observations from file
@@ -339,6 +348,11 @@
 obs_pressure_top, obs_height_top, sfc_elevation_check, sfc_elevation_tol, &
 overwrite_obs_time, anal_time)
 !
+call add_supplimental_obs(airnow_co_extra, seq_airnow_co, max_obs_seq, &
+AIRNOW_CO, nx, ny, obs_boundary, include_sig_data, &
+obs_pressure_top, obs_height_top, sfc_elevation_check, sfc_elevation_tol, &
+overwrite_obs_time, anal_time)
+!
 call add_supplimental_obs(airnow_o3_extra, seq_airnow_o3, max_obs_seq, &
 AIRNOW_O3, nx, ny, obs_boundary, include_sig_data, &
 obs_pressure_top, obs_height_top, sfc_elevation_check, sfc_elevation_tol, &
@@ -370,8 +384,12 @@
 !  super-ob IASI_O3 data
 ! if ( superob_iasi_o3 ) call superob_iasi_o3_data(seq_iasi_o3, anal_time, &
 !                                     iasi_co_horiz_int, iasi_co_pres_int)
+!  super-ob AIRNOW_CO data
+! if ( superob_airnow_co ) call superob_airnow_co_data(seq_airnow_co, anal_time, &
+!                                     airnow_co_horiz_int, airnow_co_pres_int)
 !  super-ob AIRNOW_O3 data
 ! if ( superob_airnow_o3 ) call superob_airnow_o3_data(seq_airnow_o3, anal_time, &
+!                                     airnow_o3_horiz_int, airnow_o3_pres_int)
 ! APM/JB ---
 
 max_obs_seq = get_num_obs(seq_tc)     + get_num_obs(seq_rawin) + &
@@ -381,7 +399,7 @@
 ! APM/JB +++
               get_num_obs(seq_modis_aod) + get_num_obs(seq_mopitt_co) + &
               get_num_obs(seq_iasi_co) + get_num_obs(seq_iasi_o3) + &
-              get_num_obs(seq_airnow_o3) 
+              get_num_obs(seq_airnow_co) + get_num_obs(seq_airnow_o3) 
 ! APM/JB ---
 
 call create_new_obs_seq(num_copies, num_qc, max_obs_seq, seq_all)
@@ -423,6 +441,9 @@