[Dart-dev] DART/branches Revision: 12251

[dart at ucar.edu]

thoar at ucar.edu
2017-12-29 13:56:34 -0700 (Fri, 29 Dec 2017)
314
Confirmed the assimilation of Temperature observations. Had to open 
the time window to 240 seconds in the navdas converter to get some 
temperature obs. Increased total number of observations assimilated 
to about 2800. Took almost 16 minutes end-to-end 
(debugging on - 4 members, 16 tasks on 1 cheyenne node).




Modified: DART/branches/coamps/models/coamps_nest/externals/obs_def/obs_def_navdas_mod.f90
===================================================================
--- DART/branches/coamps/models/coamps_nest/externals/obs_def/obs_def_navdas_mod.f90	2017-12-29 18:54:36 UTC (rev 12250)
+++ DART/branches/coamps/models/coamps_nest/externals/obs_def/obs_def_navdas_mod.f90	2017-12-29 20:56:34 UTC (rev 12251)
@@ -255,7 +255,7 @@
     ! Pressure - for now, assume that we are working in pressure
     ! coordinates only.  Once we add other capabilities, change this
     ! to call model_interpolate for QTY_PRESSURE
-    if (.not. is_vertical(location, 'pressure')) then
+    if (.not. is_vertical(location, 'PRESSURE')) then
        t(:) = missing_r8
        istatus(:) = 11
        return

Modified: DART/branches/coamps/models/coamps_nest/work/input.nml
===================================================================
--- DART/branches/coamps/models/coamps_nest/work/input.nml	2017-12-29 18:54:36 UTC (rev 12250)
+++ DART/branches/coamps/models/coamps_nest/work/input.nml	2017-12-29 20:56:34 UTC (rev 12251)
@@ -56,7 +56,7 @@
    async                        = 4
    adv_ens_command              = "./advance_wrapper.csh"
 
-   obs_sequence_in_name         = "obs_seq.out"
+   obs_sequence_in_name         = "obs_seq.processed"
    obs_sequence_out_name        = "obs_seq.final"
    num_output_obs_members       = 40
    init_time_days               = -1
@@ -142,7 +142,7 @@
 &obs_sequence_tool_nml
    num_input_files   = 1
    filename_seq      = 'obs_seq.out'
-   filename_out      = 'obs_seq.subset'
+   filename_out      = 'obs_seq.processed'
    first_obs_days    = -1
    first_obs_seconds = -1
    last_obs_days     = -1
@@ -158,7 +158,7 @@
                        'U_WSAT_WIND' , 'V_WSAT_WIND',
                        'U_ASCAT_WIND', 'V_ASCAT_WIND',
                        'U_TERRA_WIND', 'V_TERRA_WIND',
-                       'T_RAOB'
+                       'T_RAOB', 'T_AMDAR', 'T_MDCRS'
    keep_types        = .true.
    print_only        = .false.
    min_lat           = -90.0
@@ -178,11 +178,6 @@
   /
 
 
-#  evaluate_these_obs_types   = 'EVAL_POT_TEMPERATURE',
-#                               'RADIOSONDE_U_WIND_COMPONENT',
-#                               'RADIOSONDE_LOG_SPECIFIC_HUMIDITY', 
-#                               'P_SFC_LAND' 
-
 &obs_kind_nml
    assimilate_these_obs_types = 'U_RAOB'      , 'V_RAOB', 'Z_RAOB' , 'T_RAOB' ,
                                 'U_AIREP'     , 'V_AIREP',
@@ -363,7 +358,7 @@
    innov_file_name = 'innov.out'
    ngt_file_name   = 'ngt.out'
    obs_seq_in_name = 'obs_seq.out'
-   obs_window      = 60
+   obs_window      = 240
    verbose         = .true.
   /