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

Timestamp:

Sep 23, 2009 9:40:33 AM (15 years ago)

Author:

raasch

Message:

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

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Check profiles fpr prho and hyp.
 ! Bugfix: output of averaged 2d/3d quantities requires that an avaraging
 ! interval has been set, respective error message is included
@@ -933 +934 @@
 
 !
-!--    If required compute the profile of leaf area density used in the plant canopy model 
+!--    If required compute the profile of leaf area density used in the plant
+!--    canopy model 
        IF ( plant_canopy ) THEN
        
@@ -967 +969 @@
 
 !
-!--       In case of no given leaf area density gradients, choose a vanishing gradient
+!--       In case of no given leaf area density gradients, choose a vanishing
+!--       gradient
           IF ( lad_vertical_gradient_level(1) == -9999999.9 ) THEN
              lad_vertical_gradient_level(1) = 0.0
@@ -2089 +2092 @@
 
           CASE ( 'rho' )
-             dopr_index(i) = 64
-             dopr_unit(i)  = 'kg/m3'
-             hom(:,2,64,:) = SPREAD( zu, 2, statistic_regions+1 )
+             IF ( .NOT. ocean ) THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for ocean = .FALSE.'
+                CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 64
+                dopr_unit(i)  = 'kg/m3'
+                hom(:,2,64,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
 
           CASE ( 'w"sa"' )
@@ -2150 +2160 @@
                 hom(:,2,70,:) = SPREAD( zu, 2, statistic_regions+1 )
              ENDIF
+
+          CASE ( 'prho' )
+             IF ( .NOT. ocean ) THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for ocean = .FALSE.'
+                CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 71
+                dopr_unit(i)  = 'kg/m3'
+                hom(:,2,71,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'hyp' )
+             dopr_index(i) = 72
+             dopr_unit(i)  = 'kPa'
+             hom(:,2,72,:) = SPREAD( zu, 2, statistic_regions+1 )
 
           CASE DEFAULT
@@ -2929 +2956 @@
 !-- Check pressure gradient conditions
     IF ( dp_external .AND. conserve_volume_flow )  THEN
-       WRITE( message_string, * )  'Both dp_external and conserve_volume_flow', &
-            ' are .TRUE. but one of them must be .FALSE.'
+       WRITE( message_string, * )  'Both dp_external and conserve_volume_flo', &
+            'w are .TRUE. but one of them must be .FALSE.'
        CALL message( 'check_parameters', 'PA0150', 1, 2, 0, 6, 0 )
     ENDIF
@@ -2940 +2967 @@
        ENDIF
        IF ( .NOT. ANY( dpdxy /= 0.0 ) )  THEN
-          WRITE( message_string, * )  'dp_external is .TRUE. but dpdxy is zero',&
-               ', i.e. the external pressure gradient & will not be applied'
+          WRITE( message_string, * )  'dp_external is .TRUE. but dpdxy is ze', &
+               'ro, i.e. the external pressure gradient & will not be applied'
           CALL message( 'check_parameters', 'PA0152', 0, 1, 0, 6, 0 )
        ENDIF