Ignore:
Timestamp:
Sep 23, 2009 9:40:33 AM (12 years ago)
Author:
raasch
Message:

in-situ AND potential density are calculated and used in the ocean version

File:
1 edited

Legend:

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

    r376 r388  
    44! Actual revisions:
    55! -----------------
     6! Check profiles fpr prho and hyp.
    67! Bugfix: output of averaged 2d/3d quantities requires that an avaraging
    78! interval has been set, respective error message is included
     
    933934
    934935!
    935 !--    If required compute the profile of leaf area density used in the plant canopy model
     936!--    If required compute the profile of leaf area density used in the plant
     937!--    canopy model
    936938       IF ( plant_canopy ) THEN
    937939       
     
    967969
    968970!
    969 !--       In case of no given leaf area density gradients, choose a vanishing gradient
     971!--       In case of no given leaf area density gradients, choose a vanishing
     972!--       gradient
    970973          IF ( lad_vertical_gradient_level(1) == -9999999.9 ) THEN
    971974             lad_vertical_gradient_level(1) = 0.0
     
    20892092
    20902093          CASE ( 'rho' )
    2091              dopr_index(i) = 64
    2092              dopr_unit(i)  = 'kg/m3'
    2093              hom(:,2,64,:) = SPREAD( zu, 2, statistic_regions+1 )
     2094             IF ( .NOT. ocean ) THEN
     2095                message_string = 'data_output_pr = ' // &
     2096                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
     2097                                 'lemented for ocean = .FALSE.'
     2098                CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
     2099             ELSE
     2100                dopr_index(i) = 64
     2101                dopr_unit(i)  = 'kg/m3'
     2102                hom(:,2,64,:) = SPREAD( zu, 2, statistic_regions+1 )
     2103             ENDIF
    20942104
    20952105          CASE ( 'w"sa"' )
     
    21502160                hom(:,2,70,:) = SPREAD( zu, 2, statistic_regions+1 )
    21512161             ENDIF
     2162
     2163          CASE ( 'prho' )
     2164             IF ( .NOT. ocean ) THEN
     2165                message_string = 'data_output_pr = ' // &
     2166                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
     2167                                 'lemented for ocean = .FALSE.'
     2168                CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
     2169             ELSE
     2170                dopr_index(i) = 71
     2171                dopr_unit(i)  = 'kg/m3'
     2172                hom(:,2,71,:) = SPREAD( zu, 2, statistic_regions+1 )
     2173             ENDIF
     2174
     2175          CASE ( 'hyp' )
     2176             dopr_index(i) = 72
     2177             dopr_unit(i)  = 'kPa'
     2178             hom(:,2,72,:) = SPREAD( zu, 2, statistic_regions+1 )
    21522179
    21532180          CASE DEFAULT
     
    29292956!-- Check pressure gradient conditions
    29302957    IF ( dp_external .AND. conserve_volume_flow )  THEN
    2931        WRITE( message_string, * )  'Both dp_external and conserve_volume_flow', &
    2932             ' are .TRUE. but one of them must be .FALSE.'
     2958       WRITE( message_string, * )  'Both dp_external and conserve_volume_flo', &
     2959            'w are .TRUE. but one of them must be .FALSE.'
    29332960       CALL message( 'check_parameters', 'PA0150', 1, 2, 0, 6, 0 )
    29342961    ENDIF
     
    29402967       ENDIF
    29412968       IF ( .NOT. ANY( dpdxy /= 0.0 ) )  THEN
    2942           WRITE( message_string, * )  'dp_external is .TRUE. but dpdxy is zero',&
    2943                ', i.e. the external pressure gradient & will not be applied'
     2969          WRITE( message_string, * )  'dp_external is .TRUE. but dpdxy is ze', &
     2970               'ro, i.e. the external pressure gradient & will not be applied'
    29442971          CALL message( 'check_parameters', 'PA0152', 0, 1, 0, 6, 0 )
    29452972       ENDIF
Note: See TracChangeset for help on using the changeset viewer.