Changeset 240 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Feb 18, 2009 5:50:38 PM (15 years ago)

Author:

letzel

Message:

• External pressure gradient (check_parameters, init_3d_model, header, modules, parin, prognostic_equations)
• New topography case 'single_street_canyon'

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -5 +5 @@
 ! -----------------
 ! Bugfix: pressure included for profile output
+! Check pressure gradient conditions
 !
 ! Former revisions:
@@ -2816 +2817 @@
 
 !
+!-- Check pressure gradient conditions
+    IF ( dp_external .AND. conserve_volume_flow )  THEN
+       WRITE( message_string, * )  'Both dp_external and conserve_volume_flow', &
+            ' are .TRUE. but one of them must be .FALSE.'
+       CALL message( 'check_parameters', 'PA0150', 1, 2, 0, 6, 0 )
+    ENDIF
+    IF ( dp_external )  THEN
+       IF ( dp_level_b < zu(nzb) .OR. dp_level_b > zu(nzt) )  THEN
+          WRITE( message_string, * )  'dp_level_b = ', dp_level_b, ' is out ', &
+               ' of range'
+          CALL message( 'check_parameters', 'PA0151', 1, 2, 0, 6, 0 )
+       ENDIF
+       IF ( .NOT. ANY( dpdxy /= 0.0 ) )  THEN
+          WRITE( message_string, * )  'dp_external is .TRUE. but dpdxy is zero',&
+               ', i.e. the external pressure gradient & will not be applied'
+          CALL message( 'check_parameters', 'PA0152', 0, 1, 0, 6, 0 )
+       ENDIF
+    ENDIF
+    IF ( ANY( dpdxy /= 0.0 ) .AND. .NOT. dp_external )  THEN
+       WRITE( message_string, * )  'dpdxy is nonzero but dp_external is ', &
+            '.FALSE., i.e. the external pressure gradient & will not be applied'
+       CALL message( 'check_parameters', 'PA0153', 0, 1, 0, 6, 0 )
+    ENDIF
+
+!
 !-- Check &userpar parameters
     CALL user_check_parameters