Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/check_parameters.f90

    r240 r241  
    28402840       CALL message( 'check_parameters', 'PA0153', 0, 1, 0, 6, 0 )
    28412841    ENDIF
     2842    IF ( conserve_volume_flow )  THEN
     2843       IF ( TRIM( conserve_volume_flow_mode ) == 'default' )  THEN
     2844          IF ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' )  THEN
     2845             conserve_volume_flow_mode = 'inflow_profile'
     2846          ELSE
     2847             conserve_volume_flow_mode = 'initial_profiles'
     2848          ENDIF
     2849       ELSEIF ( TRIM( conserve_volume_flow_mode ) /= 'initial_profiles' .AND.  &
     2850            TRIM( conserve_volume_flow_mode ) /= 'inflow_profile' .AND.  &
     2851            TRIM( conserve_volume_flow_mode ) /= 'bulk_velocity' )  THEN
     2852          WRITE( message_string, * )  'unknown conserve_volume_flow_mode: ', &
     2853               conserve_volume_flow_mode
     2854          CALL message( 'check_parameters', 'PA0154', 1, 2, 0, 6, 0 )
     2855       ENDIF
     2856       IF ( ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' ) .AND. &
     2857            TRIM( conserve_volume_flow_mode ) /= 'inflow_profile' )  THEN
     2858          WRITE( message_string, * )  'noncyclic boundary conditions ', &
     2859               'require & conserve_volume_flow_mode = ''inflow_profile'''
     2860          CALL message( 'check_parameters', 'PA0155', 1, 2, 0, 6, 0 )
     2861       ENDIF
     2862       IF ( bc_lr == 'cyclic'  .AND.  bc_ns == 'cyclic'  .AND.  &
     2863            TRIM( conserve_volume_flow_mode ) == 'inflow_profile' )  THEN
     2864          WRITE( message_string, * )  'cyclic boundary conditions ', &
     2865               'require & conserve_volume_flow_mode = ''initial_profiles''', &
     2866               ' or ''bulk_velocity'''
     2867          CALL message( 'check_parameters', 'PA0156', 1, 2, 0, 6, 0 )
     2868       ENDIF
     2869    ENDIF
     2870    IF ( ( u_bulk /= 0.0 .OR. v_bulk /= 0.0 ) .AND.  &
     2871         ( .NOT. conserve_volume_flow .OR.  &
     2872         TRIM( conserve_volume_flow_mode ) /= 'bulk_velocity' ) )  THEN
     2873       WRITE( message_string, * )  'nonzero bulk velocity requires ', &
     2874            'conserve_volume_flow = .T. and & ', &
     2875            'conserve_volume_flow_mode = ''bulk_velocity'''
     2876       CALL message( 'check_parameters', 'PA0157', 1, 2, 0, 6, 0 )
     2877    ENDIF
    28422878
    28432879!
Note: See TracChangeset for help on using the changeset viewer.