[Dart-dev] DART/branches Revision: 11450

dart at ucar.edu dart at ucar.edu
Sat Apr 8 08:45:20 MDT 2017


mizzi at ucar.edu
2017-04-08 08:45:20 -0600 (Sat, 08 Apr 2017)
30
Update for PANDA surface obs




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-04-08 14:43:56 UTC (rev 11449)
+++ DART/branches/mizzi/models/wrf_chem/WRF_DART_utilities/wrf_dart_obs_preprocess.f90	2017-04-08 14:45:20 UTC (rev 11450)
@@ -42,8 +42,8 @@
                              MOPITT_CO_RETRIEVAL,  &
                              IASI_CO_RETRIEVAL, &
                              IASI_O3_RETRIEVAL, &
-                             AIRNOW_CO, &
-                             AIRNOW_O3
+                             AIRNOW_CO, AIRNOW_O3, &
+                             PANDA_CO, PANDA_O3, PANDA_PM25
 ! APM/JB ---
 use time_manager_mod, only : time_type, set_calendar_type, GREGORIAN, set_time
 use        model_mod, only : static_init_model
@@ -73,7 +73,10 @@
                       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'
+                      airnow_o3_extra    = 'obs_seq.airnow_o3',  &
+                      panda_co_extra     = 'obs_seq.panda_co',   &
+                      panda_o3_extra     = 'obs_seq.panda_o3',   &
+                      panda_pm25_extra   = 'obs_seq.panda_pm25'
 character(len=80)  :: name, sgday, sgsec
 ! APM/JB ---
 integer            :: max_num_obs              = 600000   ! Largest number of obs in one sequence
@@ -137,6 +140,21 @@
 real(r8)           :: airnow_o3_horiz_int         = 100.0_r8   ! horizontal interval for super-ob
 real(r8)           :: airnow_o3_pres_int          = 2500.0_r8  ! pressure interval for super-ob
 logical            :: overwrite_ncep_airnow_o3_qc = .false.    ! true to overwrite NCEP QC (see instructions)
+!  PANDA CO specific parameters
+logical            :: superob_panda_co            = .false.    ! super-ob sat wind data
+real(r8)           :: panda_co_horiz_int          = 100.0_r8   ! horizontal interval for super-ob
+real(r8)           :: panda_co_pres_int           = 2500.0_r8  ! pressure interval for super-ob
+logical            :: overwrite_ncep_panda_co_qc  = .false.    ! true to overwrite NCEP QC (see instructions)
+!  PANDA O3 specific parameters
+logical            :: superob_panda_o3            = .false.    ! super-ob sat wind data
+real(r8)           :: panda_o3_horiz_int          = 100.0_r8   ! horizontal interval for super-ob
+real(r8)           :: panda_o3_pres_int           = 2500.0_r8  ! pressure interval for super-ob
+logical            :: overwrite_ncep_panda_o3_qc  = .false.    ! true to overwrite NCEP QC (see instructions)
+!  PANDA PM25 specific parameters
+logical            :: superob_panda_pm25           = .false.    ! super-ob sat wind data
+real(r8)           :: panda_pm25_horiz_int         = 100.0_r8   ! horizontal interval for super-ob
+real(r8)           :: panda_pm25_pres_int          = 2500.0_r8  ! pressure interval for super-ob
+logical            :: overwrite_ncep_panda_pm25_qc = .false.    ! true to overwrite NCEP QC (see instructions)
 ! APM/JB ---
 
 !  surface obs. specific parameters
@@ -157,7 +175,10 @@
          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
+         superob_airnow_o3, airnow_o3_pres_int, airnow_o3_extra, airnow_o3_horiz_int, &
+         superob_panda_co, panda_co_pres_int, panda_co_extra, panda_co_horiz_int, &
+         superob_panda_o3, panda_o3_pres_int, panda_o3_extra, panda_o3_horiz_int, &
+         superob_panda_pm25, panda_pm25_pres_int, panda_pm25_extra, panda_pm25_horiz_int
 ! APM/JB ---
 
 ! ----------------------------------------------------------------------
@@ -177,7 +198,8 @@
                            seq_prof, seq_tc, seq_gpsro, seq_other, &
 ! APM/JB +++
                            seq_modis_aod, seq_mopitt_co, seq_iasi_co, seq_iasi_o3, &
-                           seq_airnow_co, seq_airnow_o3
+                           seq_airnow_co, seq_airnow_o3, &
+                           seq_panda_co, seq_panda_o3, seq_panda_pm25
 ! APM/JB ---
 
 type(time_type)         :: anal_time
@@ -257,6 +279,9 @@
 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)
+call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_panda_co)
+call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_panda_o3)
+call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_panda_pm25)
 ! APM/JB ---
 call create_new_obs_seq(num_copies, num_qc, max_obs_seq, seq_other)
 
@@ -267,7 +292,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_co, seq_airnow_o3, seq_other)
+seq_airnow_co, seq_airnow_o3, seq_panda_co, seq_panda_o3, seq_panda_pm25, seq_other)
 ! APM/JB ---
 
 !  add supplimental rawinsonde observations from file
@@ -356,6 +381,21 @@
 AIRNOW_O3, 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(panda_co_extra, seq_panda_co, max_obs_seq, &
+PANDA_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(panda_o3_extra, seq_panda_o3, max_obs_seq, &


More information about the Dart-dev mailing list