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

Timestamp:

Nov 2, 2016 1:47:31 PM (8 years ago)

Author:

suehring

Message:

Bugfix, read and write restart data for wall fluxes; additional checks for wall fluxes

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Additional checks for wall_heatflux, wall_humidityflux and wall_scalarflux.
+! Bugfix, check for constant_scalarflux.
 ! 
 ! Former revisions:
@@ -1735 +1736 @@
     ENDIF
 
+    IF ( ANY( wall_heatflux /= 0.0_wp )  .AND.                        &
+         surface_heatflux == 9999999.9_wp )  THEN
+       message_string = 'wall_heatflux additionally requires ' //     &
+                        'setting of surface_heatflux'
+       CALL message( 'check_parameters', 'PA0443', 1, 2, 0, 6, 0 )
+    ENDIF
+
 !
 !   This IF clause needs revision, got too complex!!
@@ -1854 +1862 @@
 !-- Set boundary conditions for total water content
     IF ( humidity )  THEN
+
+       IF ( ANY( wall_humidityflux /= 0.0_wp )  .AND.                        &
+            surface_waterflux == 9999999.9_wp )  THEN
+          message_string = 'wall_humidityflux additionally requires ' //     &
+                           'setting of surface_waterflux'
+          CALL message( 'check_parameters', 'PA0444', 1, 2, 0, 6, 0 )
+       ENDIF
+
        CALL set_bc_scalars( 'q', bc_q_b, bc_q_t, ibc_q_b, ibc_q_t,           &
                             'PA0071', 'PA0072' )
@@ -1883 +1899 @@
     
     IF ( passive_scalar )  THEN
+
+       IF ( ANY( wall_scalarflux /= 0.0_wp )  .AND.                        &
+            surface_scalarflux == 9999999.9_wp )  THEN
+          message_string = 'wall_scalarflux additionally requires ' //     &
+                           'setting of surface_scalarflux'
+          CALL message( 'check_parameters', 'PA0445', 1, 2, 0, 6, 0 )
+       ENDIF
+
+       IF ( surface_scalarflux == 9999999.9_wp )  constant_scalarflux = .FALSE.
+
        CALL set_bc_scalars( 's', bc_s_b, bc_s_t, ibc_s_b, ibc_s_t,             &
                             'PA0071', 'PA0072' )