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

Timestamp:

Jan 15, 2014 1:40:41 PM (10 years ago)

Author:

heinze

Message:

Usage of Dirichlet bottom boundary condition for scalars in conjunction with large scale forcing enabled

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! Use LSF_DATA also in case of Dirichlet bottom boundary condition for scalars
 !
 ! Former revisions:
@@ -1613 +1613 @@
     ENDIF
 
-    IF ( surface_heatflux == 9999999.9  .AND.  .NOT. large_scale_forcing )  THEN
-        constant_heatflux     = .FALSE.
+    IF ( surface_heatflux == 9999999.9  )  THEN
+       constant_heatflux = .FALSE.
+       IF ( large_scale_forcing )  THEN
+          IF ( ibc_pt_b == 0 )  THEN
+             constant_heatflux = .FALSE.
+          ELSEIF ( ibc_pt_b == 1 )  THEN
+             constant_heatflux = .TRUE.
+             IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
+                surface_heatflux = shf_surf(1)
+             ENDIF
+          ENDIF
+       ENDIF
     ELSE
-        constant_heatflux     = .TRUE.
+        constant_heatflux = .TRUE.
         IF ( TRIM( initializing_actions ) /= 'read_restart_data' .AND. &
                large_scale_forcing ) THEN
-           surface_heatflux      = shf_surf(1)
+           surface_heatflux = shf_surf(1)
         ENDIF
     ENDIF
@@ -1742 +1752 @@
        ENDIF
 
-       IF ( surface_waterflux == 9999999.9  .AND.  &
-                             .NOT. large_scale_forcing )  THEN
-           constant_waterflux     = .FALSE.
+       IF ( surface_waterflux == 9999999.9  )  THEN
+          constant_waterflux = .FALSE.
+          IF ( large_scale_forcing )  THEN
+             IF ( ibc_q_b == 0 )  THEN
+                constant_waterflux = .FALSE.
+             ELSEIF ( ibc_q_b == 1 )  THEN
+                constant_waterflux = .TRUE.
+                IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
+                   surface_waterflux = qsws_surf(1)
+                ENDIF
+             ENDIF
+          ENDIF
        ELSE
-           constant_waterflux     = .TRUE.
-           IF ( TRIM( initializing_actions ) /= 'read_restart_data' .AND. &
-               large_scale_forcing ) THEN 
-              surface_waterflux      = qsws_surf(1)
-           ENDIF
+          constant_waterflux = .TRUE.
+          IF ( TRIM( initializing_actions ) /= 'read_restart_data' .AND. &
+                 large_scale_forcing ) THEN
+             surface_waterflux = qsws_surf(1)
+          ENDIF
        ENDIF