Changeset 1856 for palm/trunk/SOURCE

Timestamp:

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

Author:

maronga

Message:

bugfix in land surface / radiation model

Location:

palm/trunk/SOURCE

Files:

• land_surface_model_mod.f90 (modified) (3 diffs)
• radiation_model_mod.f90 (modified) (2 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.

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: allocation of rad_lw_out for radiation_scheme = 'clear-sky'
 ! 
 ! Former revisions:
@@ -806 +806 @@
              ALLOCATE ( rad_lw_in(0:0,nysg:nyng,nxlg:nxrg) )
           ENDIF
+          IF ( .NOT. ALLOCATED ( rad_lw_out ) )  THEN
+             ALLOCATE ( rad_lw_out(0:0,nysg:nyng,nxlg:nxrg) )
+          ENDIF
 
           IF ( .NOT. ALLOCATED ( rad_lw_in_av ) )  THEN