Changeset 2729 for palm/trunk/SOURCE

Timestamp:

Jan 9, 2018 11:22:28 AM (7 years ago)

Author:

maronga

Message:

input of deep soil temperature separated from prognostic soil temperature

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Separated deep soil temperature from soil_temperature array
+! 
+! 2724 2018-01-05 12:12:38Z maronga
 ! Added security check for insufficient soil_temperature values
 ! 
@@ -434 +437 @@
                 canopy_resistance_coefficient = 9999999.9_wp, & !< NAMELIST g_d
                 c_surface = 9999999.9_wp,               & !< Surface (skin) heat capacity (J/m2/K)
+                deep_soil_temperature =  9999999.9_wp,  & !< Deep soil temperature (bottom boundary condition)
                 drho_l_lv,                              & !< (rho_l * l_v)**-1
                 exn,                                    & !< value of the Exner function
@@ -1495 +1499 @@
     ENDIF
 
+
+!!! these checks are not needed for water surfaces??
+
 !
 !-- Determine number of soil layers to be used and check whether an appropriate
@@ -1514 +1521 @@
 !
 !-- Check whether valid soil temperatures are prescribed
-    IF ( ANY( soil_temperature(nzb_soil:nzt_soil+1) == 9999999.9_wp ) )  THEN
+    IF ( COUNT( soil_temperature /= 9999999.9_wp ) /= nzs )  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, ')'
+                                   soil_temperature /= 9999999.9_wp ), ')'
           CALL message( 'check_parameters', 'PA0471', 1, 2, 0, 6, 0 )
+    ENDIF
+
+    IF ( deep_soil_temperature == 9999999.9_wp ) THEN
+          message_string = 'deep_soil_temperature is not set but must be'//    &
+                           '/= 9999999.9'
+          CALL message( 'check_parameters', 'PA0472', 1, 2, 0, 6, 0 )
     ENDIF
 
@@ -3847 +3860 @@
                    m_soil_h%var_2d(k,m) = soil_moisture(k)
                 ENDDO
-                t_soil_h%var_2d(nzt_soil+1,m) = soil_temperature(nzt_soil+1)
+                t_soil_h%var_2d(nzt_soil+1,m) = deep_soil_temperature
              ENDIF
           ENDDO
@@ -3858 +3871 @@
                       m_soil_v(l)%var_2d(k,m) = soil_moisture(k)
                    ENDDO
-                   t_soil_v(l)%var_2d(nzt_soil+1,m) = soil_temperature(nzt_soil+1)
+                   t_soil_v(l)%var_2d(nzt_soil+1,m) = deep_soil_temperature
                 ENDIF
              ENDDO
@@ -4489 +4502 @@
                                   constant_roughness,                          &
                                   conserve_water_content,                      &
+                                  deep_soil_temperature,                       &
                                   dz_soil,                                     &
                                   f_shortwave_incoming, field_capacity,        &