[Dart-dev] DART/branches Revision: 12705

[dart at ucar.edu]

nancy at ucar.edu
2018-07-03 09:22:25 -0600 (Tue, 03 Jul 2018)
595
start of a converter for NCEP obs distributed
in NetCDF format instead of prepbufr format.
this converter isn't finished because we discovered
the files were missing some of the data we needed
to determine exactly which obs types were there.
also the expected errors were missing.  these were
going to be added to the output files by the people
writing the converters, so this program may be the
way forward instead of continuing to try to parse
prepbufr format.  see info on the NCAR RDA for
dataste 337.0 - look for info on the subsetting
option and asking for netcdf output instead of bufr.




Added: DART/branches/rma_trunk/observations/obs_converters/NCEP/netcdf/README
===================================================================
--- DART/branches/rma_trunk/observations/obs_converters/NCEP/netcdf/README	                        (rev 0)
+++ DART/branches/rma_trunk/observations/obs_converters/NCEP/netcdf/README	2018-07-03 15:22:25 UTC (rev 12705)
@@ -0,0 +1,23 @@
+# DART software - Copyright UCAR. This open source software is provided
+# by UCAR, "as is", without charge, subject to all terms of use at
+# http://www.image.ucar.edu/DAReS/DART/DART_download
+#
+# DART $Id$
+
+This is the start of a converter from the new NetCDF format
+files with NCEP obs in them, directly to DART obs_seq files.
+
+THIS CONVERTER IS NOT FINISHED YET.  anyone interested in
+doing the rest, let dart at ucar.edu know.  we'd love to use 
+it instead of parsing prepbufr files.
+
+Thanks to Soyoung Ha for the initial work on this converter.
+We stopped when some of the needed information was not in
+the netcdf file, but later versions of the converter are
+expected to contain what's needed.
+
+# <next few lines under version control, do not edit>
+# $URL$
+# $Revision$
+# $Date$
+


Property changes on: DART/branches/rma_trunk/observations/obs_converters/NCEP/netcdf/README
___________________________________________________________________
Added: svn:eol-style
## -0,0 +1 ##
+native
\ No newline at end of property
Added: svn:keywords
## -0,0 +1 ##
+Date Rev Author HeadURL Id
\ No newline at end of property
Added: svn:mime-type
## -0,0 +1 ##
+text/plain
\ No newline at end of property
Added: DART/branches/rma_trunk/observations/obs_converters/NCEP/netcdf/convert_pb_netcdf.f90
===================================================================
--- DART/branches/rma_trunk/observations/obs_converters/NCEP/netcdf/convert_pb_netcdf.f90	                        (rev 0)
+++ DART/branches/rma_trunk/observations/obs_converters/NCEP/netcdf/convert_pb_netcdf.f90	2018-07-03 15:22:25 UTC (rev 12705)
@@ -0,0 +1,453 @@
+! DART software - Copyright UCAR. This open source software is provided
+! by UCAR, "as is", without charge, subject to all terms of use at
+! http://www.image.ucar.edu/DAReS/DART/DART_download
+!
+! $Id$
+
+program convert_pb_netcdf
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!
+!   convert_pb_netcdf - program that reads a prep_bufr file that has
+!                       been converted to netcdf format by MET version 5.
+!                       (Model Evaluation Tools) Developmental Testbed Center
+!
+! adapted from MADIS converters by Soyoung Ha and Nancy Collins,  3/27/2018
+!
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+use         types_mod, only : r8, missing_r8
+use     utilities_mod, only : nc_check, initialize_utilities, finalize_utilities
+use  time_manager_mod, only : time_type, set_calendar_type, set_date, set_time, &
+                              increment_time, get_time, operator(-), GREGORIAN, &
+                              set_time_missing, print_date
+use      location_mod, only : VERTISSURFACE, VERTISPRESSURE
+use  obs_sequence_mod, only : obs_sequence_type, obs_type, read_obs_seq, &
+                              static_init_obs_sequence, init_obs, write_obs_seq, & 
+                              init_obs_sequence, get_num_obs, & 
+                              set_copy_meta_data, set_qc_meta_data
+use           meteor_mod, only : sat_vapor_pressure, specific_humidity, & 
+                                 wind_dirspd_to_uv, pres_alt_to_pres, &
+                                 temp_and_dewpoint_to_rh
+use          obs_err_mod  ! all
+use dewpoint_obs_err_mod, only : dewpt_error_from_rh_and_temp, &
+                                 rh_error_from_dewpt_and_temp
+use         obs_kind_mod  ! FIXME, all for now
+use             sort_mod,  only : index_sort
+use obs_def_altimeter_mod, only : compute_altimeter
+use     obs_utilities_mod, only : getvar_real_2d, get_or_fill_QC, add_obs_to_seq, &
+                                  create_3d_obs, getvar_int, getdimlen, set_missing_name, &
+                                  getvar_char
+
+use netcdf
+
+implicit none
+
+character(len=128), parameter :: pb_netcdf_file = 'prep_bufr_input.nc'
+character(len=128), parameter :: pb_out_file    = 'obs_seq.prep_bufr'
+
+! the following logical parameters control which water-vapor variables appear in the output file,
+! whether to use the NCEP error or Lin and Hubbard (2004) moisture error model, and if the
+! input file has data quality control fields, whether to use or ignore them.
+logical, parameter :: LH_err                    = .false.