Changeset 1970 for palm/trunk/SOURCE

Timestamp:

Jul 18, 2016 2:27:12 PM (8 years ago)

Author:

suehring

Message:

Bugfix, check for ibc_q_b and constant_waterflux in case of land-surface scheme

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix, check for ibc_q_b and constant_waterflux in case of land-surface scheme
 ! 
 ! Former revisions:
@@ -3921 +3921 @@
 !--    A given surface value implies Dirichlet boundary condition for
 !--    the respective quantity. In this case specification of a constant flux is
-!--    forbidden.
-       IF ( ibc_b == 0  .AND.  constant_flux )  THEN
-          message_string = 'boundary condition: bc_' // TRIM( sq ) // '_b ' // &
-                           '= "' // TRIM( bc_b ) // '" is not allowed wi' // &
-                           'th prescribed surface flux'
-          CALL message( 'check_parameters', err_nr_3, 1, 2, 0, 6, 0 )
+!--    forbidden. However, an exception is made for large-scale forcing as well
+!--    as land-surface model.
+       IF ( .NOT. land_surface  .AND.  .NOT. large_scale_forcing )  THEN
+          IF ( ibc_b == 0  .AND.  constant_flux )  THEN
+             message_string = 'boundary condition: bc_' // TRIM( sq ) // '_b ' //  &
+                              '= "' // TRIM( bc_b ) // '" is not allowed with ' // &
+                              'prescribed surface flux'
+             CALL message( 'check_parameters', err_nr_3, 1, 2, 0, 6, 0 )
+          ENDIF
        ENDIF
        IF ( constant_waterflux  .AND.  surface_initial_change /= 0.0_wp )  THEN