TLS

GCEL · Mar 29, 2019 · d654107 · d654107
1 parent 72fa5f7
commit d654107
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/ACM_GPP_ET.f90 b/src/ACM_GPP_ET.f90
@@ -1944,7 +1944,10 @@ subroutine calculate_update_soil_water(ET_leaf,ET_soil,rainfall_in,corrected_ET)
        ! pass information to waterloss variable and zero watergain
        ! convert kg.m-2 (or mm) -> Mg.m-2 (or m)
        waterloss = 0d0 ; watergain = 0d0
-       waterloss(1:nos_root_layers) = evaporation_losses(1:nos_root_layers)*1d-3
+       ! adjust water already committed to evaporation
+       ! convert kg.m-2 (or mm) -> Mg.m-2 (or m)
+       soil_waterfrac(1:nos_root_layers) = soil_waterfrac(1:nos_root_layers) &
+                                         + ((-evaporation_losses(1:nos_root_layers)*1d-3) / layer_thickness(1:nos_root_layers))
 
        !!!!!!!!!!
        ! Gravitational drainage