[Dart-dev] DART/branches Revision: 10926

[dart at ucar.edu]

thoar at ucar.edu
2017-01-20 11:43:15 -0700 (Fri, 20 Jan 2017)
490
This test is bitwise identical.

This implements trunk  wrf/model_mod.f90 changes in r8397 
to the wrf_chem/model_mod.f90

Paraphrasing from r8397:
"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: DART/branches/mizzi/models/wrf_chem/model_mod.f90
===================================================================
--- DART/branches/mizzi/models/wrf_chem/model_mod.f90	2017-01-20 18:08:01 UTC (rev 10925)
+++ DART/branches/mizzi/models/wrf_chem/model_mod.f90	2017-01-20 18:43:15 UTC (rev 10926)
@@ -5451,7 +5451,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))
@@ -8610,7 +8610,6 @@
    ! surface field, otherwise you do the full-up 3d interpolation.
    if ( wrf_state_variables(1, i) == 'PSFC' ) then
       in_state_vector(KIND_SURFACE_PRESSURE) = .true.
-      in_state_vector(KIND_LANDMASK) = .true.
    endif
 
 enddo
@@ -8679,10 +8678,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.