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

Timestamp:

Jan 5, 2018 12:12:38 PM (6 years ago)

Author:

maronga

Message:

some changes in spinup mechanism and additional check in land surface model

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added security check for insufficient soil_temperature values
+! 
+! 2723 2018-01-05 09:27:03Z maronga
 ! Bugfix for spinups (end_time was increased twice in case of LSM + USM runs)
 ! 
@@ -478 +481 @@
 
                                             
-    REAL(wp), DIMENSION(0:20)  ::  root_fraction = 9999999.9_wp, & !< (NAMELIST) distribution of root surface area to the individual soil layers
-                                   soil_moisture = 0.0_wp,       & !< NAMELIST soil moisture content (m3/m3)
-                                   soil_temperature = 300.0_wp,  & !< NAMELIST soil temperature (K) +1
-                                   dz_soil  = 9999999.9_wp,      & !< (NAMELIST) soil layer depths (spacing)
+    REAL(wp), DIMENSION(0:20)  ::  root_fraction = 9999999.9_wp,     & !< (NAMELIST) distribution of root surface area to the individual soil layers
+                                   soil_moisture = 0.0_wp,           & !< NAMELIST soil moisture content (m3/m3)
+                                   soil_temperature = 9999999.9_wp,  & !< NAMELIST soil temperature (K) +1
+                                   dz_soil  = 9999999.9_wp,          & !< (NAMELIST) soil layer depths (spacing)
                                    zs_layer = 9999999.9_wp         !< soil layer depths (edge)
                                  
@@ -1509 +1512 @@
     nzs = nzt_soil + 1
 
-
+!
+!-- Check whether valid soil temperatures are prescribed
+    IF ( ANY( soil_temperature(nzb_soil:nzt_soil+1) == 9999999.9_wp ) )  THEN
+       WRITE( message_string, * ) 'number of soil layers (', nzs, ') does not',&
+                                  ' match to the number of layers specified',  &
+                                  ' in soil_temperature (', COUNT(             &
+                                   soil_temperature /= 9999999.9_wp )-1, ')'
+          CALL message( 'check_parameters', 'PA0471', 1, 2, 0, 6, 0 )
+    ENDIF
+
+!
+!-- Check whether the sum of all root fractions equals one
     IF ( vegetation_type == 0 )  THEN
        IF ( SUM( root_fraction(nzb_soil:nzt_soil) ) /= 1.0_wp )  THEN