[Dart-dev] DART/branches Revision: 11467

[dart at ucar.edu]

thoar at ucar.edu
2017-04-12 18:57:26 -0600 (Wed, 12 Apr 2017)
269
model_mod_140822.f90
   correct the serious mistake!
   comment "where((data_2d_array > 1.0_r8)) data_2d_array = org_array"
   when updating the first and second layer's soil moisutre!!!
   now all the increaments to SM are accepted.
   using Qh RTM instead of wg RTM




Modified: DART/branches/Tb/models/clm/model_mod.f90
===================================================================
--- DART/branches/Tb/models/clm/model_mod.f90	2017-04-13 00:51:42 UTC (rev 11466)
+++ DART/branches/Tb/models/clm/model_mod.f90	2017-04-13 00:57:26 UTC (rev 11467)
@@ -2811,7 +2811,7 @@
 
 !      where(isnan(data_2d_array)) data_2d_array = org_array
 
-      where((data_2d_array > 1.0_r8)) data_2d_array = org_array 
+!      where((data_2d_array > 1.0_r8)) data_2d_array = org_array 
       !===========================================================Long
       where((data_2d_array < 0.0_r8)) data_2d_array = org_array
    elseif (trim(progvar(ivar)%varname) == 'H2OSOI_ICE') then
@@ -4357,7 +4357,7 @@
 ! be part of the DART state vector. They are currently directly harvested from the CLM
 ! restart file. As such, the posteriors are not informative.
 
-use   radiative_transfer_mod, only : forward_wg
+use   radiative_transfer_mod, only : forward_Qh
 
 type(time_type),        intent(in)  :: state_time      ! valid time of DART state
 type(location_type),    intent(in)  :: location        ! observation location
@@ -4535,7 +4535,7 @@
    if (landcovercode >= 0 ) then
       ! the tb_out array contains the calculated brightness temperature outputs
       ! at each polarization (rows) and frequency (columns).
-      call forward_wg(N_FREQ, freq, tetad, aux_ins, tb_out)
+      call forward_Qh(N_FREQ, freq, tetad, aux_ins, tb_out)
    else
       ! call to alternative radiative transfer model goes here.
    endif