Changeset 1856 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Apr 13, 2016 12:56:17 PM (9 years ago)

Author:

maronga

Message:

bugfix in land surface / radiation model

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: for water surfaces, the initial water surface temperature is set equal
+! to the intital skin temperature. Moreover, the minimum value of r_a is now
+! 1.0 to avoid too large fluxes at the first model time step
 ! 
 ! Former revisions:
@@ -1074 +1076 @@
 
 !
-!--          Make sure that the resistance does not drop to zero
-             IF ( ABS(r_a(j,i)) < 1.0E-10_wp )  r_a(j,i) = 1.0E-10_wp
+!--          Make sure that the resistance does not drop to zero for neutral 
+!--          stratification
+             IF ( ABS(r_a(j,i)) < 1.0_wp )  r_a(j,i) = 1.0_wp
 
 !
@@ -1885 +1888 @@
              IF ( veg_type_2d(j,i) == 14  .OR.  veg_type_2d(j,i) == 15 )  THEN
                 water_surface(j,i) = .TRUE.
+                t_soil(:,j,i) = t_surface(j,i)
              ELSEIF ( veg_type_2d(j,i) == 20 )  THEN
                 pave_surface(j,i) = .TRUE.