[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