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

Timestamp:

Oct 24, 2018 6:39:32 PM (6 years ago)

Author:

gronemeier

Message:

new surface-data output; renamed output variables (pt to theta, rho_air to rho, rho_ocean to rho_sea_water)

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Renamed output variables
+! Add checks for surface data output
+! 
+! 3419 2018-10-24 17:27:31Z gronemeier
 ! Remove offline nesting in if clause for pressure top boundary condition
 ! 
@@ -777 +781 @@
 
     USE statistics
+
+    USE surface_output_mod,                                                    &
+        ONLY:  surface_output_check_parameters
 
     USE synthetic_turbulence_generator_mod,                                    &
@@ -1494 +1501 @@
     IF ( urban_surface )        CALL usm_check_parameters
     IF ( wind_turbine )         CALL wtm_check_parameters
+!
+!-- Checks for surface data output
+    IF ( surface_data_output )  CALL surface_output_check_parameters
 
 !
@@ -2363 +2373 @@
              hom(:,2,3,:)  = SPREAD( zw, 2, statistic_regions+1 )
 
-          CASE ( 'pt', '#pt' )
+          CASE ( 'theta', '#theta' )
              IF ( .NOT. bulk_cloud_model ) THEN
                 dopr_index(i) = 4
@@ -2452 +2462 @@
              hom(:,2,15,:) = SPREAD( zw, 2, statistic_regions+1 )
 
-          CASE ( 'w"pt"' )
+          CASE ( 'w"theta"' )
              dopr_index(i) = 16
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,16,:) = SPREAD( zw, 2, statistic_regions+1 )
 
-          CASE ( 'w*pt*' )
+          CASE ( 'w*theta*' )
              dopr_index(i) = 17
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,17,:) = SPREAD( zw, 2, statistic_regions+1 )
 
-          CASE ( 'wpt' )
+          CASE ( 'wtheta' )
              dopr_index(i) = 18
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
@@ -2479 +2489 @@
              IF ( constant_flux_layer )  hom(nzb,2,20,:) = zu(1)
 
-          CASE ( 'w*pt*BC' )
+          CASE ( 'w*theta*BC' )
              dopr_index(i) = 21
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,21,:) = SPREAD( zw, 2, statistic_regions+1 )
 
-          CASE ( 'wptBC' )
+          CASE ( 'wthetaBC' )
              dopr_index(i) = 22
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
@@ -2504 +2514 @@
              hom(:,2,32,:) = SPREAD( zw, 2, statistic_regions+1 )
 
-          CASE ( 'pt*2' )
+          CASE ( 'theta*2' )
              dopr_index(i) = 33
              dopr_unit(i)  = 'K2'
@@ -2603 +2613 @@
              ENDIF
 
-          CASE ( 'lpt', '#lpt' )
+          CASE ( 'thetal', '#thetal' )
              IF ( .NOT. bulk_cloud_model ) THEN
                 message_string = 'data_output_pr = ' //                        &
@@ -2621 +2631 @@
              ENDIF
 
-          CASE ( 'vpt', '#vpt' )
+          CASE ( 'thetav', '#thetav' )
              dopr_index(i) = 44
              dopr_unit(i)  = 'K'
@@ -2632 +2642 @@
              ENDIF
 
-          CASE ( 'w"vpt"' )
+          CASE ( 'w"thetav"' )
              dopr_index(i) = 45
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,45,:) = SPREAD( zw, 2, statistic_regions+1 )
 
-          CASE ( 'w*vpt*' )
+          CASE ( 'w*thetav*' )
              dopr_index(i) = 46
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,46,:) = SPREAD( zw, 2, statistic_regions+1 )
 
-          CASE ( 'wvpt' )
+          CASE ( 'wthetav' )
              dopr_index(i) = 47
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
@@ -2798 +2808 @@
              hom(:,2,57,:) = SPREAD( zu, 2, statistic_regions+1 )
 
-          CASE ( 'u"pt"' )
+          CASE ( 'u"theta"' )
              dopr_index(i) = 58
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,58,:) = SPREAD( zu, 2, statistic_regions+1 )
 
-          CASE ( 'u*pt*' )
+          CASE ( 'u*theta*' )
              dopr_index(i) = 59
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,59,:) = SPREAD( zu, 2, statistic_regions+1 )
 
-          CASE ( 'upt_t' )
+          CASE ( 'utheta_t' )
              dopr_index(i) = 60
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,60,:) = SPREAD( zu, 2, statistic_regions+1 )
 
-          CASE ( 'v"pt"' )
+          CASE ( 'v"theta"' )
              dopr_index(i) = 61
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,61,:) = SPREAD( zu, 2, statistic_regions+1 )
 
-          CASE ( 'v*pt*' )
+          CASE ( 'v*theta*' )
              dopr_index(i) = 62
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
              hom(:,2,62,:) = SPREAD( zu, 2, statistic_regions+1 )
 
-          CASE ( 'vpt_t' )
+          CASE ( 'thetav_t' )
              dopr_index(i) = 63
              dopr_unit(i)  = TRIM ( heatflux_output_unit )
@@ -2855 +2865 @@
              hom(:,2,72,:) = SPREAD( zu, 2, statistic_regions+1 )
 
-          CASE ( 'rho_air' )
+          CASE ( 'rho' )
              dopr_index(i)  = 119
              dopr_unit(i)   = 'kg/m3'
              hom(:,2,119,:) = SPREAD( zu, 2, statistic_regions+1 )
 
-          CASE ( 'rho_air_zw' )
+          CASE ( 'rho_zw' )
              dopr_index(i)  = 120
              dopr_unit(i)   = 'kg/m3'
@@ -3023 +3033 @@
              unit = 'm2/s2'
 
-          CASE ( 'lpt' )
+          CASE ( 'thetal' )
              IF (  .NOT.  bulk_cloud_model )  THEN
                 message_string = 'output of "' // TRIM( var ) // '" requi' //  &
@@ -3041 +3051 @@
              IF ( TRIM( var ) == 'pr' )  unit = 'm'
 
-          CASE ( 'q', 'vpt' )
+          CASE ( 'q', 'thetav' )
              IF (  .NOT.  humidity )  THEN
                 message_string = 'output of "' // TRIM( var ) // '" requi' //  &
@@ -3048 +3058 @@
              ENDIF
              IF ( TRIM( var ) == 'q'   )  unit = 'kg/kg'
-             IF ( TRIM( var ) == 'vpt' )  unit = 'K'
+             IF ( TRIM( var ) == 'thetav' )  unit = 'K'
 
           CASE ( 'ql' )
@@ -3084 +3094 @@
              unit = 'kg/m3'
 
-          CASE ( 'p', 'pt', 'u', 'v', 'w' )
+          CASE ( 'p', 'theta', 'u', 'v', 'w' )
              IF ( TRIM( var ) == 'p'  )  unit = 'Pa'
-             IF ( TRIM( var ) == 'pt' )  unit = 'K'
+             IF ( TRIM( var ) == 'theta' )  unit = 'K'
              IF ( TRIM( var ) == 'u'  )  unit = 'm/s'
              IF ( TRIM( var ) == 'v'  )  unit = 'm/s'
@@ -3093 +3103 @@
 
           CASE ( 'ghf*', 'lwp*', 'ol*', 'qsws*', 'r_a*',                       &
-                 'shf*', 'ssws*', 't*', 'tsurf*', 'u*', 'z0*', 'z0h*', 'z0q*' )
+                 'shf*', 'ssws*', 't*', 'tsurf*', 'us*', 'z0*', 'z0h*', 'z0q*' )
              IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
                 message_string = 'illegal value for data_output: "' //         &
@@ -3141 +3151 @@
              IF ( TRIM( var ) == 't*'     )  unit = 'K'
              IF ( TRIM( var ) == 'tsurf*' )  unit = 'K'  
-             IF ( TRIM( var ) == 'u*'     )  unit = 'm/s'
+             IF ( TRIM( var ) == 'us*'    )  unit = 'm/s'
              IF ( TRIM( var ) == 'z0*'    )  unit = 'm'
              IF ( TRIM( var ) == 'z0h*'   )  unit = 'm'