[Dart-dev] [3875] DART/trunk/models/wrf: I'd overlooked adding: default_state_variables = .true.
nancy at ucar.edu
nancy at ucar.edu
Fri May 15 08:23:37 MDT 2009
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20090515/6d11d722/attachment.html
-------------- next part --------------
Modified: DART/trunk/models/wrf/model_mod.f90
===================================================================
--- DART/trunk/models/wrf/model_mod.f90 2009-05-14 04:01:18 UTC (rev 3874)
+++ DART/trunk/models/wrf/model_mod.f90 2009-05-15 14:23:35 UTC (rev 3875)
@@ -145,7 +145,7 @@
! center_spline_scale: coarse grid to spline interp. fine grid ratio
!-----------------------------------------------------------------------
-logical :: output_state_vector = .false. ! output prognostic variables
+logical :: output_state_vector = .false. ! output prognostic variables
logical :: default_state_variables = .true. ! use default state list?
character(len=129) :: wrf_state_variables(num_state_table_columns,max_state_variables) = 'NULL'
character(len=129) :: wrf_state_bounds(num_bounds_table_columns,max_state_variables) = 'NULL'
@@ -308,6 +308,13 @@
if ( default_state_variables ) then
wrf_state_variables = 'NULL'
call fill_default_state_table(wrf_state_variables)
+ call error_handler(E_MSG, 'static_init_model:', &
+ 'Using predefined wrf variable list for dart state vector.')
+ call error_handler(E_MSG, 'static_init_model:', &
+ 'Set "default_state_variables" to .false. in the namelist')
+ call error_handler(E_MSG, 'static_init_model:', &
+ 'to use the "wrf_state_variables" list instead.')
+
endif
if ( debug ) then
@@ -483,6 +490,9 @@
if ( debug ) then
print*,'variable size ',trim(wrf_state_variables(1,my_index)),' ',wrf%dom(id)%var_size(:,ind)
endif
+
+ write(errstring, '(A,I4,2A)') 'state vector array ', ind, ' is ', trim(wrf_state_variables(1,my_index))
+ call error_handler(E_MSG, 'static_init_model: ', errstring)
enddo
! close data file, we have all we need
@@ -5544,7 +5554,7 @@
'static_init_model', 'get_att DY')
call nc_check( nf90_get_att(ncid, nf90_global, 'DT', wrf%dom(id)%dt), &
'static_init_model', 'get_att DT')
- if (do_output()) print*,'dt from wrfinput file is: ', wrf%dom(id)%dt
+ if (do_output()) print*,'dt from wrfinput_d0X file is: ', wrf%dom(id)%dt
if(debug) write(*,*) ' dx, dy, dt are ',wrf%dom(id)%dx, &
wrf%dom(id)%dy, wrf%dom(id)%dt
Modified: DART/trunk/models/wrf/work/input.nml
===================================================================
--- DART/trunk/models/wrf/work/input.nml 2009-05-14 04:01:18 UTC (rev 3874)
+++ DART/trunk/models/wrf/work/input.nml 2009-05-15 14:23:35 UTC (rev 3875)
@@ -133,22 +133,33 @@
# 2 = pressure
# 3 = height
+# set default_state_variables to .false. to use the explicit list.
+# otherwise it uses a hardcoded default list: U, V, W, PH, T, MU only.
+# other possible common fields are:
+# wrf_state_variables = 'U','KIND_U_WIND_COMPONENT','TYPE_U','UPDATE','999',
+# 'V','KIND_V_WIND_COMPONENT','TYPE_V','UPDATE','999',
+# 'W','KIND_VERTICAL_VELOCITY','TYPE_W','UPDATE','999',
+# 'T','KIND_POTENTIAL_TEMPERATURE','TYPE_T','UPDATE','999',
+# 'PH','KIND_GEOPOTENTIAL_HEIGHT','TYPE_GZ','UPDATE','999',
+# 'MU','KIND_PRESSURE','TYPE_MU','UPDATE','999',
+# 'QVAPOR','KIND_VAPOR_MIXING_RATIO','TYPE_QV','UPDATE','999',
+# 'QCLOUD','KIND_CLOUD_LIQUID_WATER','TYPE_QC','UPDATE','999',
+# 'QRAIN','KIND_RAINWATER_MIXING_RATIO','TYPE_QR','UPDATE','999',
+# 'U10','KIND_U_WIND_COMPONENT','TYPE_U10','UPDATE','999',
+# 'V10','KIND_V_WIND_COMPONENT','TYPE_V10','UPDATE','999',
+# 'T2','KIND_TEMPERATURE','TYPE_T2','UPDATE','999',
+# 'TH2','KIND_POTENTIAL_TEMPERATURE','TYPE_TH2','UPDATE','999',
+# 'Q2','KIND_SPECIFIC_HUMIDITY','TYPE_Q2','UPDATE','999',
+# 'PSFC','KIND_PRESSURE','TYPE_PS','UPDATE','999',
+
&model_nml
- wrf_state_variables = 'U','KIND_U_WIND_COMPONENT','TYPE_U','UPDATE','999',
- 'V','KIND_V_WIND_COMPONENT','TYPE_V','UPDATE','999',
- 'W','KIND_VERTICAL_VELOCITY','TYPE_W','UPDATE','999',
- 'PH','KIND_GEOPOTENTIAL_HEIGHT','TYPE_GZ','UPDATE','999',
- 'T','KIND_POTENTIAL_TEMPERATURE','TYPE_T','UPDATE','999',
- 'MU','KIND_PRESSURE','TYPE_MU','UPDATE','999',
- 'QVAPOR','KIND_VAPOR_MIXING_RATIO','TYPE_QV','UPDATE','999',
- 'QCLOUD','KIND_CLOUD_LIQUID_WATER','TYPE_QC','UPDATE','999',
- 'QRAIN','KIND_RAINWATER_MIXING_RATIO','TYPE_QR','UPDATE','999',
- 'U10','KIND_U_WIND_COMPONENT','TYPE_U10','UPDATE','999',
- 'V10','KIND_V_WIND_COMPONENT','TYPE_V10','UPDATE','999',
- 'T2','KIND_TEMPERATURE','TYPE_T2','UPDATE','999',
- 'TH2','KIND_POTENTIAL_TEMPERATURE','TYPE_TH2','UPDATE','999',
- 'Q2','KIND_SPECIFIC_HUMIDITY','TYPE_Q2','UPDATE','999',
- 'PSFC','KIND_PRESSURE','TYPE_PS','UPDATE','999',
+ default_state_variables = .true.,
+ wrf_state_variables = 'U', 'KIND_U_WIND_COMPONENT', 'TYPE_U', 'UPDATE','999',
+ 'V', 'KIND_V_WIND_COMPONENT', 'TYPE_V', 'UPDATE','999',
+ 'W', 'KIND_VERTICAL_VELOCITY', 'TYPE_W', 'UPDATE','999',
+ 'PH','KIND_GEOPOTENTIAL_HEIGHT', 'TYPE_GZ','UPDATE','999',
+ 'T', 'KIND_POTENTIAL_TEMPERATURE','TYPE_T', 'UPDATE','999',
+ 'MU','KIND_PRESSURE', 'TYPE_MU','UPDATE','999',
wrf_state_bounds = 'QVAPOR','0.0','NULL','CLAMP',
'QRAIN','0.0','NULL','CLAMP',
'QCLOUD','0.0','NULL','CLAMP',
@@ -157,8 +168,6 @@
num_domains = 1,
calendar_type = 3,
surf_obs = .true.,
- soil_data = .true.,
- h_diab = .false.,
assimilation_period_seconds = 21600,
adv_mod_command = './wrf.exe',
vert_localization_coord = 3,
More information about the Dart-dev
mailing list