[Dart-dev] [3337] DART/trunk/models/wrf: Update WRF model_mod to only print a single copy of the domain details

nancy at subversion.ucar.edu nancy at subversion.ucar.edu
Thu May 15 10:28:11 MDT 2008


An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20080515/df487958/attachment.html
-------------- next part --------------
Modified: DART/trunk/models/wrf/WRF_DART_utilities/dart_tf_wrf.f90
===================================================================
--- DART/trunk/models/wrf/WRF_DART_utilities/dart_tf_wrf.f90	2008-05-14 22:26:24 UTC (rev 3336)
+++ DART/trunk/models/wrf/WRF_DART_utilities/dart_tf_wrf.f90	2008-05-15 16:28:11 UTC (rev 3337)
@@ -53,12 +53,17 @@
 real (kind=r8) :: center_search_half_length = 500000.0_r8
 integer :: center_spline_grid_scale = 10
 integer :: vert_localization_coord  =  3  ! 1,2,3 == level,pressure,height
+!nc -- we are adding these to the model.nml until they appear in the NetCDF files
+logical :: polar = .false.
+logical :: periodic_x = .false.
 
+
 namelist /model_nml/ output_state_vector, num_moist_vars, &
                      num_domains, calendar_type, surf_obs, soil_data, h_diab, &
                      adv_mod_command, assimilation_period_seconds, &
-                     vert_localization_coord, allow_obs_below_vol, &
-                     center_search_half_length, center_spline_grid_scale
+                     allow_obs_below_vol, vert_localization_coord, &
+                     center_search_half_length, center_spline_grid_scale, &
+                     polar, periodic_x
 
 !-------------------------------------------------------------
 

Modified: DART/trunk/models/wrf/WRF_DART_utilities/ensemble_init.f90
===================================================================
--- DART/trunk/models/wrf/WRF_DART_utilities/ensemble_init.f90	2008-05-14 22:26:24 UTC (rev 3336)
+++ DART/trunk/models/wrf/WRF_DART_utilities/ensemble_init.f90	2008-05-15 16:28:11 UTC (rev 3337)
@@ -42,6 +42,12 @@
                              wrf_dealloc, wrfbdy_dealloc, &
                              wrf_io, wrfbdy_io, &
                              set_wrf_date
+use      location_mod, only : location_type, get_location, set_location, & 
+                              vert_is_surface, vert_is_level, vert_is_pressure, vert_is_height, &
+                              VERTISUNDEF, VERTISSURFACE, VERTISLEVEL, VERTISPRESSURE, &
+                              VERTISHEIGHT
+                             
+
 use    utilities_mod, only : get_unit, file_exist, open_file, &
                              close_file, error_handler, E_ERR, E_MSG, initialize_utilities, &
                              register_module, logfileunit, timestamp, &
@@ -61,24 +67,36 @@
 ! Model namelist parameters with default values.
 !-----------------------------------------------------------------------
 
-logical :: output_state_vector  = .false.  ! state vs. prognostic format
+logical :: output_state_vector  = .false.     ! output prognostic variables
 integer :: num_moist_vars       = 3
 integer :: num_domains          = 1
 integer :: calendar_type        = GREGORIAN
-logical :: surf_obs             = .true.
-logical :: h_diab               = .false.
-logical :: allow_obs_below_vol  = .false.
 integer :: assimilation_period_seconds = 21600
+logical :: surf_obs             = .true.      
+logical :: soil_data            = .true.      
+logical :: h_diab               = .false.     
 character(len = 72) :: adv_mod_command = './wrf.exe'
-integer :: center_search_size       = 25
-integer :: center_spline_grid_scale = 10
-integer :: vert_localization_coord  =  3  ! 1,2,3 == level,pressure,height
+real (kind=r8) :: center_search_half_length = 500000.0_r8
+integer :: center_search_half_size            
+integer :: center_spline_grid_scale = 10      
+integer :: vert_localization_coord = VERTISHEIGHT
+! Allow (or not) observations above the surface but below the lowest
+! sigma level.
+logical :: allow_obs_below_vol = .false.
+! Max height a surface obs can be away from the actual model surface
+! and still be accepted (in meters)
+!real (kind=r8) :: max_surface_delta = 500.0
+!nc -- we are adding these to the model.nml until they appear in the NetCDF files
+logical :: polar = .false.  
+logical :: periodic_x = .false.
 
+
 namelist /model_nml/ output_state_vector, num_moist_vars, &
-                     num_domains, calendar_type, surf_obs, h_diab, &
+                     num_domains, calendar_type, surf_obs, soil_data, h_diab, &
                      adv_mod_command, assimilation_period_seconds, &
-                     vert_localization_coord, allow_obs_below_vol, &
-                     center_search_size, center_spline_grid_scale
+                     allow_obs_below_vol, vert_localization_coord, &
+                     center_search_half_length, center_spline_grid_scale, &
+                     polar, periodic_x
 
 !-----------------------------------------------------------------------
 

Modified: DART/trunk/models/wrf/model_mod.f90
===================================================================
--- DART/trunk/models/wrf/model_mod.f90	2008-05-14 22:26:24 UTC (rev 3336)
+++ DART/trunk/models/wrf/model_mod.f90	2008-05-15 16:28:11 UTC (rev 3337)
@@ -38,7 +38,8 @@
 
 use     utilities_mod, only : file_exist, open_file, close_file, &
                               register_module, error_handler, E_ERR, E_WARN, &
-                              E_MSG, logfileunit, find_namelist_in_file, check_namelist_read
+                              E_MSG, logfileunit, do_output, &
+                              find_namelist_in_file, check_namelist_read
 
 use      obs_kind_mod, only : KIND_U_WIND_COMPONENT, KIND_V_WIND_COMPONENT, &
                               KIND_SURFACE_PRESSURE, KIND_TEMPERATURE, &
@@ -301,11 +302,14 @@
 
 do id=1,num_domains
 
-   write( idom , '(I1)') id
+   ! only print this once, no matter how many parallel tasks are running
+   if (do_output()) then
+      write( idom , '(I1)') id
 
-   write(*,*) '******************'
-   write(*,*) '**  DOMAIN # ',idom,'  **'
-   write(*,*) '******************'
+      write(*,*) '******************'
+      write(*,*) '**  DOMAIN # ',idom,'  **'
+      write(*,*) '******************'
+   endif
 
    if(file_exist('wrfinput_d0'//idom)) then
 

Modified: DART/trunk/models/wrf/work/path_names_ensemble_init
===================================================================
--- DART/trunk/models/wrf/work/path_names_ensemble_init	2008-05-14 22:26:24 UTC (rev 3336)
+++ DART/trunk/models/wrf/work/path_names_ensemble_init	2008-05-15 16:28:11 UTC (rev 3337)
@@ -3,3 +3,6 @@
 common/types_mod.f90
 time_manager/time_manager_mod.f90
 utilities/utilities_mod.f90
+location/threed_sphere/location_mod.f90
+random_nr/random_nr_mod.f90
+random_seq/random_seq_mod.f90


More information about the Dart-dev mailing list