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

Timestamp:

Dec 16, 2019 8:17:03 AM (4 years ago)

Author:

Giersch

Message:

Topography closed channel flow with symmetric boundaries implemented, ID tag in radiation module corrected

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -25 +25 @@
 ! -----------------
 ! 4172 2019-08-20 11:55:33Z oliver.maas
+! Checks for closed channel flow implemented
+! 
+! 11:55:33Z oliver.maas
 ! Move 2-m potential temperature output to diagnostic_output_quantities
 ! 
@@ -498 +501 @@
           WRITE( action, '(A)' )  'cloud_droplets = .TRUE.'
        ENDIF
-       IF ( .NOT. constant_flux_layer )  THEN
+       IF ( .NOT. constant_flux_layer .AND. topography /= 'closed_channel' )   &
+       THEN
           WRITE( action, '(A)' )  'constant_flux_layer = .FALSE.'
        ENDIF
        IF ( action /= ' ' )  THEN
-          message_string = 'a non-flat topography does not allow ' //          &
+          message_string = 'The specified topography does not allow ' //       &
                            TRIM( action )
           CALL message( 'check_parameters', 'PA0014', 1, 2, 0, 6, 0 )
        ENDIF
-
+!
+!--    Check illegal/untested parameter combinations for closed channel
+       If ( topography == 'closed_channel' ) THEN 
+          symmetry_flag = 1
+          message_string = 'Bottom and top boundary are treated equal'
+          CALL message( 'check_parameters', 'PA0676', 0, 0, 0, 6, 0 )
+       
+          IF ( dz(1) /= dz(COUNT( dz /= -1.0_wp )) .OR.                        &
+               dz_stretch_level /= -9999999.9_wp) THEN
+             WRITE( message_string, * )  'dz should be equal close to the ' // &
+                                         'boundaries due to symmetrical problem'
+             CALL message( 'check_parameters', 'PA0313', 1, 2, 0, 6, 0 )
+          ENDIF
+       
+          IF ( constant_flux_layer ) THEN
+             WRITE( message_string, * )  'A constant flux layer is not '//     &
+                                         'allowed if a closed channel '//      &
+                                         'shall be used'
+             CALL message( 'check_parameters', 'PA0688', 1, 2, 0, 6, 0 )
+          ENDIF
+       
+          IF ( ocean_mode ) THEN
+             WRITE( message_string, * )  'The ocean mode is not allowed if '// &
+                                         'a closed channel shall be used'
+             CALL message( 'check_parameters', 'PA0689', 1, 2, 0, 6, 0 )
+          ENDIF
+       
+          IF ( momentum_advec /= 'ws-scheme' .OR.                              &
+               scalar_advec /= 'ws-scheme' ) THEN
+             WRITE( message_string, * )  'A closed channel require the '//     &
+                                         'upwind scheme of Wicker and ' //     &
+                                         'Skamarock as the advection scheme'
+             CALL message( 'check_parameters', 'PA0690', 1, 2, 0, 6, 0 )
+          ENDIF
+       ENDIF
     ENDIF