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

Timestamp:

Feb 19, 2009 4:08:56 PM (13 years ago)

Author:

letzel

Message:

• Option to predefine a target bulk velocity for conserve_volume_flow

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (1 diff)

palm/trunk/SOURCE/check_parameters.f90

@@ -2840 +2840 @@
        CALL message( 'check_parameters', 'PA0153', 0, 1, 0, 6, 0 )
     ENDIF
+    IF ( conserve_volume_flow )  THEN
+       IF ( TRIM( conserve_volume_flow_mode ) == 'default' )  THEN
+          IF ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' )  THEN
+             conserve_volume_flow_mode = 'inflow_profile'
+          ELSE
+             conserve_volume_flow_mode = 'initial_profiles'
+          ENDIF
+       ELSEIF ( TRIM( conserve_volume_flow_mode ) /= 'initial_profiles' .AND.  &
+            TRIM( conserve_volume_flow_mode ) /= 'inflow_profile' .AND.  &
+            TRIM( conserve_volume_flow_mode ) /= 'bulk_velocity' )  THEN
+          WRITE( message_string, * )  'unknown conserve_volume_flow_mode: ', &
+               conserve_volume_flow_mode
+          CALL message( 'check_parameters', 'PA0154', 1, 2, 0, 6, 0 )
+       ENDIF
+       IF ( ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' ) .AND. &
+            TRIM( conserve_volume_flow_mode ) /= 'inflow_profile' )  THEN
+          WRITE( message_string, * )  'noncyclic boundary conditions ', &
+               'require & conserve_volume_flow_mode = ''inflow_profile'''
+          CALL message( 'check_parameters', 'PA0155', 1, 2, 0, 6, 0 )
+       ENDIF
+       IF ( bc_lr == 'cyclic'  .AND.  bc_ns == 'cyclic'  .AND.  &
+            TRIM( conserve_volume_flow_mode ) == 'inflow_profile' )  THEN
+          WRITE( message_string, * )  'cyclic boundary conditions ', &
+               'require & conserve_volume_flow_mode = ''initial_profiles''', &
+               ' or ''bulk_velocity'''
+          CALL message( 'check_parameters', 'PA0156', 1, 2, 0, 6, 0 )
+       ENDIF
+    ENDIF
+    IF ( ( u_bulk /= 0.0 .OR. v_bulk /= 0.0 ) .AND.  &
+         ( .NOT. conserve_volume_flow .OR.  &
+         TRIM( conserve_volume_flow_mode ) /= 'bulk_velocity' ) )  THEN
+       WRITE( message_string, * )  'nonzero bulk velocity requires ', &
+            'conserve_volume_flow = .T. and & ', &
+            'conserve_volume_flow_mode = ''bulk_velocity'''
+       CALL message( 'check_parameters', 'PA0157', 1, 2, 0, 6, 0 )
+    ENDIF
 
 !