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

Timestamp:

Apr 30, 2019 9:48:33 AM (5 years ago)

Author:

suehring

Message:

Move checks for correct dimensions in static input file; improve checks concerning buildings; check whether at least one surface type is set at a natural-type surface element

File:

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

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Check that at least one surface type is set at surface element.
+! 
+! 3933 2019-04-25 12:33:20Z kanani
 ! Remove unused subroutine and allocation of pt_2m, this is done in surface_mod
 ! now (surfaces%pt_2m)
@@ -2277 +2280 @@
        ENDIF
 
-
-!        pt(k+k_off,j+j_off,i+i_off) = surf_t_surface_p%var_1d(m) / exner(nzb)  !is actually no air temperature
        surf%pt_surface(m)          = surf_t_surface_p%var_1d(m) / exner(nzb)
 
@@ -2881 +2882 @@
                 IF ( water_type_f%var(j,i)      /= water_type_f%fill )         &
                    surf_lsm_h%water_surface(m) = .TRUE.
+!
+!--             Check if at least one type is set.
+                IF ( .NOT. surf_lsm_h%vegetation_surface(m)  .AND.             &
+                     .NOT. surf_lsm_h%pavement_surface(m)    .AND.             &
+                     .NOT. surf_lsm_h%water_surface(m) )  THEN
+                   WRITE( message_string, * ) 'Horizontal surface element ' // &
+                                       ' at i, j = ',  i, j,                   &
+                                       ' is neither a vegetation, ' //         &
+                                       'pavement, nor a water surface.'
+                   CALL message( 'land_surface_model_mod', 'PA0619',          &
+                                  2, 2, myid, 6, 0 )
+                ENDIF
+                
              ENDDO
 !
@@ -2959 +2973 @@
                       surf_lsm_v(l)%water_surface(m) = .TRUE.
 !
-!--                Check if surface element has the flag %building_covered_all,
-!--                this case, set vegetation_type appropriate
+!--                Check if at least one type is set.
+                   IF ( .NOT. surf_lsm_v(l)%vegetation_surface(m)  .AND.       &
+                        .NOT. surf_lsm_v(l)%pavement_surface(m)    .AND.       &
+                        .NOT. surf_lsm_v(l)%water_surface(m) )  THEN
+                      WRITE( message_string, * ) 'Vertical surface element ' //&
+                                       ' at i, j = ',  i, j,                   &
+                                       ' is neither a vegetation, ' //         &
+                                       'pavement, nor a water surface.'
+                      CALL message( 'land_surface_model_mod', 'PA0619',        &
+                                     2, 2, myid, 6, 0 )
+                   ENDIF
                 ENDDO
              ENDDO