Ignore:
Timestamp:
Feb 18, 2009 5:50:38 PM (13 years ago)
Author:
letzel
Message:
  • External pressure gradient (check_parameters, init_3d_model, header, modules, parin, prognostic_equations)
  • New topography case 'single_street_canyon'
File:
1 edited

Legend:

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

    r232 r240  
    55! -----------------
    66! Bugfix: pressure included for profile output
     7! Check pressure gradient conditions
    78!
    89! Former revisions:
     
    28162817
    28172818!
     2819!-- Check pressure gradient conditions
     2820    IF ( dp_external .AND. conserve_volume_flow )  THEN
     2821       WRITE( message_string, * )  'Both dp_external and conserve_volume_flow', &
     2822            ' are .TRUE. but one of them must be .FALSE.'
     2823       CALL message( 'check_parameters', 'PA0150', 1, 2, 0, 6, 0 )
     2824    ENDIF
     2825    IF ( dp_external )  THEN
     2826       IF ( dp_level_b < zu(nzb) .OR. dp_level_b > zu(nzt) )  THEN
     2827          WRITE( message_string, * )  'dp_level_b = ', dp_level_b, ' is out ', &
     2828               ' of range'
     2829          CALL message( 'check_parameters', 'PA0151', 1, 2, 0, 6, 0 )
     2830       ENDIF
     2831       IF ( .NOT. ANY( dpdxy /= 0.0 ) )  THEN
     2832          WRITE( message_string, * )  'dp_external is .TRUE. but dpdxy is zero',&
     2833               ', i.e. the external pressure gradient & will not be applied'
     2834          CALL message( 'check_parameters', 'PA0152', 0, 1, 0, 6, 0 )
     2835       ENDIF
     2836    ENDIF
     2837    IF ( ANY( dpdxy /= 0.0 ) .AND. .NOT. dp_external )  THEN
     2838       WRITE( message_string, * )  'dpdxy is nonzero but dp_external is ', &
     2839            '.FALSE., i.e. the external pressure gradient & will not be applied'
     2840       CALL message( 'check_parameters', 'PA0153', 0, 1, 0, 6, 0 )
     2841    ENDIF
     2842
     2843!
    28182844!-- Check &userpar parameters
    28192845    CALL user_check_parameters
Note: See TracChangeset for help on using the changeset viewer.