[Dart-dev] [8397] DART/trunk/models/wrf/model_mod.f90: This fixes JIRA ticket DARTSUP-247 and addresses DARTSUP-245

[nancy at ucar.edu]

Revision: 8397
Author:   thoar
Date:     2015-08-18 13:39:55 -0600 (Tue, 18 Aug 2015)
Log Message:
-----------
This fixes JIRA ticket DARTSUP-247 and addresses DARTSUP-245
We believe it fixes the original issue in DARTSUP-245, but are
awaiting confirmation that it fixes the subsequent issues in the ticket.

DARTSUP-245 pertains to the ability to interpolate the land mask.

DARTSUP-247 is the shape of the 2D variables in the DART diagnostic files.

Modified Paths:
--------------
    DART/trunk/models/wrf/model_mod.f90

-------------- next part --------------
Modified: DART/trunk/models/wrf/model_mod.f90
===================================================================
--- DART/trunk/models/wrf/model_mod.f90	2015-08-18 17:29:01 UTC (rev 8396)
+++ DART/trunk/models/wrf/model_mod.f90	2015-08-18 19:39:55 UTC (rev 8397)
@@ -4914,7 +4914,7 @@
                          wrf%dom(id)%var_size(2,ind)/)
 
          allocate ( temp2d(dimsizes_2D(1),dimsizes_2D(2)) )
-         temp2d  = reshape(statevec(i:j), (/ dimsizes_3D(1),dimsizes_3D(2) /) )
+         temp2d  = reshape(statevec(i:j), (/ dimsizes_2D(1),dimsizes_2D(2) /) )
          call nc_check(nf90_put_var( ncFileID, VarID, temp2d, &
                                   start=(/ 1, 1, copyindex, timeindex /) ), &
                     'nc_write_model_vars','put_var '//trim(varname))
@@ -7870,10 +7870,11 @@
                        source, revision, revdate)
 endif
 
-! surface elevation is read in outside the state vector mechanism,
+! surface elevation and land mask are read outside the state vector mechanism,
 ! directly from the wrfinput template file, and does not vary from
 ! one ensemble member to another.
 in_state_vector(KIND_SURFACE_ELEVATION) = .true.
+in_state_vector(KIND_LANDMASK) = .true.
 
 ! there is no field that directly maps to the vortex measurements.
 ! if you have all the fields it needs, allow them.
@@ -7902,6 +7903,8 @@
 in_state_vector(KIND_RADAR_REFLECTIVITY) = .true.
 in_state_vector(KIND_POWER_WEIGHTED_FALL_SPEED) = .true.
 
+
+
 ! FIXME:  i was going to suggest nuking this routine all together because it makes
 ! the default behavior be to exit with an error when requesting to interpolate an
 ! unknown kind, instead of simply returning with a failed interpolation and letting