Changeset 3692 for palm

Timestamp:

Jan 23, 2019 2:45:49 PM (6 years ago)

Author:

suehring

Message:

Add check for soil moisture higher than its saturation value

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Revise check for soil moisture higher than its saturation value
+! 
+! 3685 2019-01-21 01:02:11Z knoop
 ! Some interface calls moved to module_interface + cleanup
 ! 
@@ -1653 +1656 @@
           CALL message( 'lsm_check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
        ENDIF
-    ENDIF    
-    
-    
-!
-!-- Check for proper setting of soil moisture, must not be larger than its 
-!-- saturation value. 
-    DO  k = nzb_soil, nzt_soil
-       IF ( soil_moisture(k) > saturation_moisture )  THEN 
-          message_string = 'soil_moisture must not exceed its saturation' //    &
-                            ' value'
-          CALL message( 'lsm_check_parameters', 'PA0458', 1, 2, 0, 6, 0 )
-       ENDIF
-    ENDDO
- 
+    ENDIF     
 !
 !-- Calculate grid spacings. Temperature and moisture are defined at
@@ -4441 +4431 @@
           ENDIF
 !
+!--       After soil moisture and temperature are finally initialized, check 
+!--       if soil moisture is higher than its saturation value. Else, this 
+!--       will produce floating point errors in the soil model parametrization.
+          DO  m = 1, surf_lsm_h%ns
+             IF ( surf_lsm_h%vegetation_surface(m)  .OR.                       &
+                  surf_lsm_h%pavement_surface(m) )  THEN
+                DO  k = nzb_soil, nzt_soil
+                   IF ( m_soil_h%var_2d(k,m) >= surf_lsm_h%m_sat(k,m) )  THEN
+                      WRITE( message_string, * ) 'soil moisture is higher '//  &
+                            'than its saturation value at (k,j,i) ', k,        &
+                            surf_lsm_h%i(m), surf_lsm_h%j(m)
+                      CALL message( 'lsm_init', 'PA0458', 2, 2, 0, 6, 0 )
+                   ENDIF               
+                ENDDO
+             ENDIF
+          ENDDO
+          DO  l = 0, 3
+             DO  m = 1, surf_lsm_v(l)%ns
+                IF ( surf_lsm_v(l)%vegetation_surface(m)  .OR.                 &
+                     surf_lsm_v(l)%pavement_surface(m) )  THEN
+                   DO  k = nzb_soil, nzt_soil
+                      IF ( m_soil_v(l)%var_2d(k,m) >= surf_lsm_v(l)%m_sat(k,m) )  &
+                      THEN
+                         WRITE( message_string, * )                               &
+                            'soil moisture is higher than ' //                 &
+                            'its saturation value at (k,j,i) ', k,             &
+                            surf_lsm_v(l)%i(m), surf_lsm_v(l)%j(m)
+                         CALL message( 'lsm_init', 'PA0458', 2, 2, 0, 6, 0 )                   
+                      ENDIF
+                   ENDDO
+                ENDIF
+             ENDDO
+          ENDDO
+
+!
 !--       Further initialization
           DO  m = 1, surf_lsm_h%ns