Changeset 1824

Timestamp:

Apr 7, 2016 9:55:29 AM (9 years ago)

Author:

gronemeier

Message:

check if roughness_length < dz/2

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-!
+! Check if roughness_length < dz/2. Moved location_message(finished) to the end.
 ! 
 ! Former revisions:
@@ -4070 +4070 @@
     ENDIF
 
-    CALL location_message( 'finished', .TRUE. )
-
 !
 !-- Prevent empty time records in volume, cross-section and masked data in case of 
@@ -4095 +4093 @@
     ENDIF
 
+!
+!-- Check roughness length, which has to be smaller than dz/2
+    IF ( ( constant_flux_layer .OR.  &
+           INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) &
+         .AND. roughness_length <= 0.5 * dz )  THEN
+       message_string = 'roughness_length must be smaller than dz/2'
+       CALL message( 'check_parameters', 'PA0424', 1, 2, 0, 6, 0 )
+    ENDIF
+
+    CALL location_message( 'finished', .TRUE. )
 !
 !-- Check &userpar parameters