[Dart-dev] [5604] DART/branches/development/models/wrf: add wind profiler obs as an option for the wrf/dart obs
nancy at ucar.edu
nancy at ucar.edu
Thu Mar 15 14:51:19 MDT 2012
Revision: 5604
Author: nancy
Date: 2012-03-15 14:51:19 -0600 (Thu, 15 Mar 2012)
Log Message:
-----------
add wind profiler obs as an option for the wrf/dart obs
preprocess program. note that this update requires you
to add the obs_def_gts_mod to your input_files list in
the &preprocess namelist. i've updated the default one
in the work directory to include gts. also updated the
.nml file to reflect the current namelist and defaults.
Modified Paths:
--------------
DART/branches/development/models/wrf/WRF_DART_utilities/wrf_dart_obs_preprocess.f90
DART/branches/development/models/wrf/WRF_DART_utilities/wrf_dart_obs_preprocess.nml
DART/branches/development/models/wrf/work/input.nml
-------------- next part --------------
Modified: DART/branches/development/models/wrf/WRF_DART_utilities/wrf_dart_obs_preprocess.f90
===================================================================
--- DART/branches/development/models/wrf/WRF_DART_utilities/wrf_dart_obs_preprocess.f90 2012-03-15 20:45:30 UTC (rev 5603)
+++ DART/branches/development/models/wrf/WRF_DART_utilities/wrf_dart_obs_preprocess.f90 2012-03-15 20:51:19 UTC (rev 5604)
@@ -40,7 +40,7 @@
use obs_kind_mod, only : RADIOSONDE_U_WIND_COMPONENT, ACARS_U_WIND_COMPONENT, &
MARINE_SFC_U_WIND_COMPONENT, LAND_SFC_U_WIND_COMPONENT, &
METAR_U_10_METER_WIND, GPSRO_REFRACTIVITY, &
- SAT_U_WIND_COMPONENT, VORTEX_LAT
+ SAT_U_WIND_COMPONENT, PROFILER_U_WIND_COMPONENT, VORTEX_LAT
use time_manager_mod, only : time_type, set_calendar_type, GREGORIAN, set_time
use model_mod, only : static_init_model
use netcdf
@@ -60,6 +60,7 @@
metar_extra = 'obs_seq.metar', &
marine_sfc_extra = 'obs_seq.marine', &
sat_wind_extra = 'obs_seq.satwnd', &
+ profiler_extra = 'obs_seq.profiler', &
gpsro_extra = 'obs_seq.gpsro', &
trop_cyclone_extra = 'obs_seq.tc'
integer :: max_num_obs = 100000 ! Largest number of obs in one sequence
@@ -100,7 +101,7 @@
sfc_elevation_check, overwrite_ncep_sfc_qc, overwrite_ncep_satwnd_qc, &
aircraft_pres_int, sat_wind_pres_int, sfc_elevation_tol, &
obs_pressure_top, obs_height_top, obs_boundary, sonde_extra, metar_extra, &
- acars_extra, land_sfc_extra, marine_sfc_extra, sat_wind_extra, &
+ acars_extra, land_sfc_extra, marine_sfc_extra, sat_wind_extra, profiler_extra, &
trop_cyclone_extra, gpsro_extra, tc_sonde_radii, increase_bdy_error, &
maxobsfac, obsdistbdy, sat_wind_horiz_int, aircraft_horiz_int, &
overwrite_obs_time
@@ -120,7 +121,7 @@
logical :: file_exist, pre_I_format
type(obs_sequence_type) :: seq_all, seq_rawin, seq_sfc, seq_acars, seq_satwnd, &
- seq_tc, seq_gpsro, seq_other
+ seq_prof, seq_tc, seq_gpsro, seq_other
type(time_type) :: anal_time
@@ -177,6 +178,7 @@
call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_sfc)
call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_acars)
call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_satwnd)
+call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_prof)
call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_gpsro)
call create_new_obs_seq(num_copies, num_qc, 100, seq_tc)
call create_new_obs_seq(num_copies, num_qc, max_num_obs, seq_other)
@@ -224,6 +226,12 @@
obs_pressure_top, obs_height_top, sfc_elevation_check, sfc_elevation_tol, &
overwrite_obs_time, anal_time)
+! add supplimental profiler observations from file
+call add_supplimental_obs(profiler_extra, seq_prof, max_obs_seq, &
+PROFILER_U_WIND_COMPONENT, nx, ny, obs_boundary, include_sig_data, &
+obs_pressure_top, obs_height_top, sfc_elevation_check, sfc_elevation_tol, &
+overwrite_obs_time, anal_time)
+
! add supplimental GPSRO observations from file
call add_supplimental_obs(gpsro_extra, seq_gpsro, max_obs_seq, &
GPSRO_REFRACTIVITY, nx, ny, obs_boundary, include_sig_data, &
@@ -250,8 +258,8 @@
max_obs_seq = get_num_obs(seq_tc) + get_num_obs(seq_rawin) + &
get_num_obs(seq_sfc) + get_num_obs(seq_acars) + &
- get_num_obs(seq_satwnd) + get_num_obs(seq_gpsro) + &
- get_num_obs(seq_other)
+ get_num_obs(seq_satwnd) + get_num_obs(seq_prof) + &
+ get_num_obs(seq_gpsro) + get_num_obs(seq_other)
call create_new_obs_seq(num_copies, num_qc, max_obs_seq, seq_all)
@@ -273,6 +281,9 @@
call build_master_sequence(seq_satwnd, seq_all)
call destroy_obs_sequence(seq_satwnd)
+call build_master_sequence(seq_prof, seq_all)
+call destroy_obs_sequence(seq_prof)
+
call build_master_sequence(seq_other, seq_all)
call destroy_obs_sequence(seq_other)
Modified: DART/branches/development/models/wrf/WRF_DART_utilities/wrf_dart_obs_preprocess.nml
===================================================================
--- DART/branches/development/models/wrf/WRF_DART_utilities/wrf_dart_obs_preprocess.nml 2012-03-15 20:45:30 UTC (rev 5603)
+++ DART/branches/development/models/wrf/WRF_DART_utilities/wrf_dart_obs_preprocess.nml 2012-03-15 20:51:19 UTC (rev 5604)
@@ -1,18 +1,42 @@
&wrf_obs_preproc_nml
- obs_boundary = 0.0,
- increase_bdy_error = .false.,
- maxobsfac = 2.5,
- obsdistbdy = 15.0,
- sfc_elevation_check = .false.,
- sfc_elevation_tol = 300.0,
- obs_pressure_top = 0.0,
- obs_height_top = 200000000000000.0,
- include_sig_data = .true.,
- tc_sonde_radii = -1.0,
- superob_aircraft = .false.,
- aircraft_horiz_int = 36.0,
- aircraft_pres_int = 2500.0,
- superob_sat_winds = .false.,
- sat_wind_horiz_int = 100.0,
- sat_wind_pres_int = 2500.0,
-/
+
+ file_name_input = 'obs_seq.old'
+ file_name_output = 'obs_seq.new'
+
+ sonde_extra = 'obs_seq.rawin'
+ land_sfc_extra = 'obs_seq.land_sfc'
+ metar_extra = 'obs_seq.metar'
+ marine_sfc_extra = 'obs_seq.marine'
+ sat_wind_extra = 'obs_seq.satwnd'
+ profiler_extra = 'obs_seq.profiler'
+ gpsro_extra = 'obs_seq.gpsro'
+ trop_cyclone_extra = 'obs_seq.tc'
+
+ max_num_obs = 100000
+ overwrite_obs_time = .false.
+
+ obs_boundary = 0.0_r8
+ increase_bdy_error = .false.
+ maxobsfac = 2.5_r8
+ obsdistbdy = 15.0_r8
+
+ sfc_elevation_check = .false.
+ sfc_elevation_tol = 300.0_r8
+ obs_pressure_top = 0.0_r8
+ obs_height_top = 2.0e10_r8
+
+ include_sig_data = .true.
+ tc_sonde_radii = -1.0_r8
+
+ superob_aircraft = .false.
+ aircraft_horiz_int = 36.0_r8
+ aircraft_pres_int = 2500.0_r8
+
+ superob_sat_winds = .false.
+ sat_wind_horiz_int = 100.0_r8
+ sat_wind_pres_int = 2500.0_r8
+
+ overwrite_ncep_satwnd_qc = .false.
+ overwrite_ncep_sfc_qc = .false.
+
+/
Modified: DART/branches/development/models/wrf/work/input.nml
===================================================================
--- DART/branches/development/models/wrf/work/input.nml 2012-03-15 20:45:30 UTC (rev 5603)
+++ DART/branches/development/models/wrf/work/input.nml 2012-03-15 20:51:19 UTC (rev 5604)
@@ -114,6 +114,7 @@
'../../../obs_def/obs_def_altimeter_mod.f90',
'../../../obs_def/obs_def_gps_mod.f90',
'../../../obs_def/obs_def_vortex_mod.f90',
+ '../../../obs_def/obs_def_gts_mod.f90',
/
&obs_kind_nml
More information about the Dart-dev
mailing list