[Dart-dev] DART/branches Revision: 11465

[dart at ucar.edu]

thoar at ucar.edu
2017-04-12 18:50:37 -0600 (Wed, 12 Apr 2017)
215
model_mod_140716.f90
   correct the serious mistake of 
       soilcolumn%ssm = h2osoi_liq(6) * 0.001 / ((LEVGRND(6)+LEVGRND(7))/2)
   to
       soilcolumn%ssm = h2osoi_liq(6) * 0.001 / ((LEVGRND(1)+LEVGRND(2))/2)




Modified: DART/branches/Tb/models/clm/model_mod.f90
===================================================================
--- DART/branches/Tb/models/clm/model_mod.f90	2017-04-13 00:49:45 UTC (rev 11464)
+++ DART/branches/Tb/models/clm/model_mod.f90	2017-04-13 00:50:37 UTC (rev 11465)
@@ -4576,9 +4576,9 @@
 ! area-weight the average
 obs_val = sum(tb * weights) / sum(weights)
 
-if (obs_val > 350.0_r8 .or. obs_val < 200.0_r8 ) then
-   obs_val=MISSING_R8
-endif  
+!if (obs_val > 350.0_r8 .or. obs_val < 200.0_r8 ) then
+!   obs_val=MISSING_R8
+!endif  
 
 !==================Long
 ! if (loc_lon > 330_r8 .or. loc_lon < 180_r8) then
@@ -4838,7 +4838,7 @@
 ! Fill the output array ... finally
 soilcolumn%nlayers = nlayers
 ! Currently, only consider the first layer of topsoil, i.e., the 6th layer of levtot
-soilcolumn%ssm = h2osoi_liq(6) * 0.001 / ((LEVGRND(6)+LEVGRND(7))/2) ! convert unit from kg/m2 to m3/m3
+soilcolumn%ssm = h2osoi_liq(6) * 0.001 / ((LEVGRND(1)+LEVGRND(2))/2) ! convert unit from kg/m2 to m3/m3 !!! mistakely using the 6th and 7th layer before, corected
 soilcolumn%stg = t_soisno(6) 
 
 deallocate(h2osoi_liq, h2osoi_ice, t_soisno)