[Dart-dev] DART/trunk Revision: 11958

[dart at ucar.edu]

thoar at ucar.edu
2017-09-27 16:15:16 -0600 (Wed, 27 Sep 2017)
412
Documenting the alternatives to the kuo observation error variance
specification in convert_cosmic_gps_cdf.

Updating convert_cosmic_ionosphere to be consistent with the rma_trunk
so the TIEGCM users on the trunk/classic will be able to use the only
observation converter we have for electron density observations.
The documentation for the data sources, etc., are contained in an
update and expanded gps.html
 



Modified: DART/trunk/models/tiegcm/work/input.nml
===================================================================
--- DART/trunk/models/tiegcm/work/input.nml	2017-09-27 21:27:46 UTC (rev 11957)
+++ DART/trunk/models/tiegcm/work/input.nml	2017-09-27 22:15:16 UTC (rev 11958)
@@ -119,7 +119,7 @@
 #              'TI',    'KIND_TEMPERATURE_ION',           'NA',      'NA',      'restart',   'NO_COPY_BACK',
 #              'TE',    'KIND_TEMPERATURE_ELECTRON',      'NA',      'NA',      'restart',   'NO_COPY_BACK',
 #              'ZG',    'KIND_GEOMETRIC_HEIGHT',          'NA',      'NA',      'secondary', 'NO_COPY_BACK',
-#              'VTEC',  'KIND_VERTICAL_TEC',              'NA',      'NA',      'calculate'  'NO_COPY_BACK' 
+#              'VTEC',  'KIND_VERTICAL_TEC',              'NA',      'NA',      'calculate'  'NO_COPY_BACK'
 
 &model_nml
    debug                       = 1
@@ -138,10 +138,9 @@
                'O2',    'KIND_MOLEC_OXYGEN_MIXING_RATIO', '0.00001', '0.99999', 'restart',    'NO_COPY_BACK',
                'TN',    'KIND_TEMPERATURE',               '0.0',     '6000.0',  'restart',    'NO_COPY_BACK',
                'ZG',    'KIND_GEOMETRIC_HEIGHT',          'NA',      'NA',      'secondary',  'NO_COPY_BACK',
-               'VTEC',  'KIND_VERTICAL_TEC',              'NA',      'NA',      'calculate',  'NO_COPY_BACK' 
+               'VTEC',  'KIND_VERTICAL_TEC',              'NA',      'NA',      'calculate',  'NO_COPY_BACK'
    /
 
-
 &model_to_dart_nml
    file_out      = 'dart_ics'
    /
@@ -168,8 +167,8 @@
    /
 
 &obs_kind_nml
-   assimilate_these_obs_types = 'CHAMP_DENSITY', 'GPS_VTEC_EXTRAP', 'GPS_PROFILE'
-   evaluate_these_obs_types   = 'GND_GPS_VTEC'
+   assimilate_these_obs_types = 'GPS_PROFILE', 'COSMIC_ELECTRON_DENSITY'
+   evaluate_these_obs_types   = 'GND_GPS_VTEC', 'CHAMP_DENSITY', 'GPS_VTEC_EXTRAP'
    /
 
 &assim_model_nml
@@ -237,11 +236,11 @@
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final'
    obs_sequence_list = ''
-   first_bin_center =  2005, 9,  9, 0, 0, 0 
-   last_bin_center  =  2005, 9, 10, 0, 0, 0 
-   bin_separation   =     0, 0,  0, 1, 0, 0 
-   bin_width        =     0, 0,  0, 1, 0, 0 
-   time_to_skip     =     0, 0,  0, 1, 0, 0 
+   first_bin_center =  2005, 9,  9, 0, 0, 0
+   last_bin_center  =  2005, 9, 10, 0, 0, 0
+   bin_separation   =     0, 0,  0, 1, 0, 0
+   bin_width        =     0, 0,  0, 1, 0, 0
+   time_to_skip     =     0, 0,  0, 1, 0, 0
    max_num_bins     = 1000
    trusted_obs      = 'null'
    Nregions   = 4

Modified: DART/trunk/obs_def/obs_def_upper_atm_mod.f90
===================================================================
--- DART/trunk/obs_def/obs_def_upper_atm_mod.f90	2017-09-27 21:27:46 UTC (rev 11957)
+++ DART/trunk/obs_def/obs_def_upper_atm_mod.f90	2017-09-27 22:15:16 UTC (rev 11958)
@@ -46,6 +46,7 @@
 ! SAT_F107,                        KIND_1D_PARAMETER,               COMMON_CODE
 ! SAT_RHO,                         KIND_DENSITY
 ! GPS_PROFILE,                     KIND_ELECTRON_DENSITY,           COMMON_CODE
+! COSMIC_ELECTRON_DENSITY,         KIND_ELECTRON_DENSITY,           COMMON_CODE
 ! GND_GPS_VTEC,		           KIND_GND_GPS_VTEC
 ! CHAMP_DENSITY,                   KIND_DENSITY
 ! MIDAS_TEC,                       KIND_VERTICAL_TEC

Modified: DART/trunk/observations/gps/convert_cosmic_gps_cdf.f90
===================================================================
--- DART/trunk/observations/gps/convert_cosmic_gps_cdf.f90	2017-09-27 21:27:46 UTC (rev 11957)
+++ DART/trunk/observations/gps/convert_cosmic_gps_cdf.f90	2017-09-27 22:15:16 UTC (rev 11958)
@@ -84,7 +84,7 @@
 real(r8) :: ray_ds                 = 5000.0_r8    ! delta stepsize (m) along ray, nonlocal op
 real(r8) :: ray_htop               = 15000.0_r8 ! max height (m) for nonlocal op
 character(len=256) :: gpsro_netcdf_file     = 'cosmic_gps_input.nc'
-character(len=256) :: gpsro_netcdf_filelist = 'cosmic_gps_input_list'
+character(len=256) :: gpsro_netcdf_filelist = ''
 character(len=256) :: gpsro_out_file        = 'obs_seq.gpsro'
 
 namelist /convert_cosmic_gps_nml/ obs_levels, local_operator, ray_ds,   &
@@ -931,7 +931,7 @@
 ! routines below here were lifted from soyoung ha's version of
 ! the ncep bufr format converter for gps obs, including lidia's error spec
 
-! soyoung says this error spec from bill is only 1 of 5 possible profiles
+! soyoung says the error spec from Bill Kuo is only 1 of 5 possible profiles
 ! and most appropriate for the CONUS domain.  if used for global obs, it
 ! should be updated to include the other functions.  nsc jan 2016
 

Modified: DART/trunk/observations/gps/convert_cosmic_gps_cdf.nml
===================================================================
--- DART/trunk/observations/gps/convert_cosmic_gps_cdf.nml	2017-09-27 21:27:46 UTC (rev 11957)
+++ DART/trunk/observations/gps/convert_cosmic_gps_cdf.nml	2017-09-27 22:15:16 UTC (rev 11958)
@@ -1,10 +1,11 @@
 &convert_cosmic_gps_nml
-  obs_levels(nmaxlevels) = -1.0
+  gpsro_netcdf_file      = 'cosmic_gps_input.nc'
+  gpsro_netcdf_filelist  = ''
+  gpsro_out_file         = 'obs_seq.gpsro'
   local_operator         = .true.