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

Timestamp:

Oct 26, 2015 4:17:44 PM (9 years ago)

Author:

maronga

Message:

various bugfixes and modifications of the atmosphere-land-surface-radiation interaction. Completely re-written routine to calculate surface fluxes (surface_layer_fluxes.f90) that replaces prandtl_fluxes. Minor formatting corrections and renamings

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -14 +14 @@
 ! PALM. If not, see <http://www.gnu.org/licenses/>.
 !
-! Copyright 1997-2014 Leibniz Universitaet Hannover
+! Copyright 1997-2015 Leibniz Universitaet Hannover
 !--------------------------------------------------------------------------------!
 !
 ! Current revisions:
 ! -----------------
-! 
+! Added output of Obukhov length (ol) and radiative heating rates for RRTMG. 
+! Added checks for use of radiation / lsm with topography.
 ! 
 ! Former revisions:
@@ -307 +308 @@
     USE transpose_indices
 
+
     IMPLICIT NONE
 
@@ -650 +652 @@
           WRITE( action, '(A)' )  'cloud_droplets = .TRUE.'
        ENDIF
-       IF ( .NOT. prandtl_layer )  THEN
-          WRITE( action, '(A)' )  'prandtl_layer = .FALSE.'
+       IF ( .NOT. constant_flux_layer )  THEN
+          WRITE( action, '(A)' )  'constant_flux_layer = .FALSE.'
        ENDIF
        IF ( action /= ' ' )  THEN
@@ -1007 +1009 @@
        ENDIF
 
-       IF ( .NOT. prandtl_layer )  THEN
+       IF ( .NOT. constant_flux_layer )  THEN
           message_string = 'lsm requires '//                                   &
-                           'prandtl_layer = .T.'
+                           'constant_flux_layer = .T.'
           CALL message( 'check_parameters', 'PA0400', 1, 2, 0, 6, 0 )
+       ENDIF
+
+       IF ( topography /= 'flat' )  THEN
+          message_string = 'lsm cannot be used ' //  & 
+                           'in combination with  topography /= "flat"'
+          CALL message( 'check_parameters', 'PA0415', 1, 2, 0, 6, 0 )
        ENDIF
 
@@ -1196 +1204 @@
           CALL message( 'check_parameters', 'PA0411', 1, 2, 0, 6, 0 )
        ENDIF
-
-    ENDIF
-
+       IF ( topography /= 'flat' )  THEN
+          message_string = 'radiation scheme cannot be used ' //  & 
+                           'in combination with  topography /= "flat"'
+          CALL message( 'check_parameters', 'PA0414', 1, 2, 0, 6, 0 )
+       ENDIF
+    ENDIF
 
     IF ( .NOT. ( loop_optimization == 'cache'  .OR.                            &
@@ -1729 +1740 @@
 !-- In case of using a prandtl-layer, calculated (or prescribed) surface
 !-- fluxes have to be used in the diffusion-terms
-    IF ( prandtl_layer )  use_surface_fluxes = .TRUE.
+    IF ( constant_flux_layer )  use_surface_fluxes = .TRUE.
 
 !
@@ -1796 +1807 @@
     ELSEIF ( bc_e_b == '(u*)**2+neumann' )  THEN
        ibc_e_b = 2
-       IF ( .NOT. prandtl_layer )  THEN
+       IF ( .NOT. constant_flux_layer )  THEN
           bc_e_b = 'neumann'
           ibc_e_b = 1
@@ -2003 +2014 @@
        IF ( surface_waterflux == 9999999.9_wp  )  THEN
           constant_waterflux = .FALSE.
-          IF ( large_scale_forcing )  THEN
+          IF ( large_scale_forcing .OR. land_surface )  THEN
              IF ( ibc_q_b == 0 )  THEN
                 constant_waterflux = .FALSE.
@@ -2045 +2056 @@
     ELSEIF ( bc_uv_b == 'neumann' )  THEN
        ibc_uv_b = 1
-       IF ( prandtl_layer )  THEN
+       IF ( constant_flux_layer )  THEN
           message_string = 'boundary condition: bc_uv_b = "' // &
-               TRIM( bc_uv_b ) // '" is not allowed with prandtl_layer = .TRUE.'
+               TRIM( bc_uv_b ) // '" is not allowed with constant_flux_layer'  &
+               // ' = .TRUE.'
           CALL message( 'check_parameters', 'PA0075', 1, 2, 0, 6, 0 )
        ENDIF
@@ -2356 +2368 @@
              dopr_unit(i)  = 'm2/s2'
              hom(:,2,12,:) = SPREAD( zw, 2, statistic_regions+1 )
-             IF ( prandtl_layer )  hom(nzb,2,12,:) = zu(1)
+             IF ( constant_flux_layer )  hom(nzb,2,12,:) = zu(1)
 
           CASE ( 'w*u*' )
@@ -2367 +2379 @@
              dopr_unit(i)  = 'm2/s2'
              hom(:,2,14,:) = SPREAD( zw, 2, statistic_regions+1 )
-             IF ( prandtl_layer )  hom(nzb,2,14,:) = zu(1)
+             IF ( constant_flux_layer )  hom(nzb,2,14,:) = zu(1)
 
           CASE ( 'w*v*' )
@@ -2393 +2405 @@
              dopr_unit(i)  = 'm2/s2'
              hom(:,2,19,:) = SPREAD( zw, 2, statistic_regions+1 )
-             IF ( prandtl_layer )  hom(nzb,2,19,:) = zu(1)
+             IF ( constant_flux_layer )  hom(nzb,2,19,:) = zu(1)
 
           CASE ( 'wv' )
@@ -2399 +2411 @@
              dopr_unit(i)  = 'm2/s2'
              hom(:,2,20,:) = SPREAD( zw, 2, statistic_regions+1 )
-             IF ( prandtl_layer )  hom(nzb,2,20,:) = zu(1)
+             IF ( constant_flux_layer )  hom(nzb,2,20,:) = zu(1)
 
           CASE ( 'w*pt*BC' )
@@ -3153 +3165 @@
              ENDIF
 
+          CASE ( 'rad_lw_cs_hr' )
+             IF ( (.NOT. radiation) .OR. radiation_scheme /= 'rrtmg' )  THEN
+                message_string = 'data_output_pr = ' //                        &
+                                 TRIM( data_output_pr(i) ) // ' is not ava' // &
+                                 'lable for radiation = .FALSE. or ' //        &
+                                 'radiation_scheme /= "rrtmg"'
+                CALL message( 'check_parameters', 'PA0413', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 106
+                dopr_unit(i)  = 'K/h'
+                hom(:,2,106,:)  = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'rad_lw_hr' )
+             IF ( (.NOT. radiation) .OR. radiation_scheme /= 'rrtmg' )  THEN
+                message_string = 'data_output_pr = ' //                        &
+                                 TRIM( data_output_pr(i) ) // ' is not ava' // &
+                                 'lable for radiation = .FALSE. or ' //        &
+                                 'radiation_scheme /= "rrtmg"'
+                CALL message( 'check_parameters', 'PA0413', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 107
+                dopr_unit(i)  = 'K/h'
+                hom(:,2,107,:)  = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'rad_sw_cs_hr' )
+             IF ( (.NOT. radiation) .OR. radiation_scheme /= 'rrtmg' )  THEN
+                message_string = 'data_output_pr = ' //                        &
+                                 TRIM( data_output_pr(i) ) // ' is not ava' // &
+                                 'lable for radiation = .FALSE. or ' //        &
+                                 'radiation_scheme /= "rrtmg"'
+                CALL message( 'check_parameters', 'PA0413', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 108
+                dopr_unit(i)  = 'K/h'
+                hom(:,2,108,:)  = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'rad_sw_hr' )
+             IF ( (.NOT. radiation) .OR. radiation_scheme /= 'rrtmg' )  THEN
+                message_string = 'data_output_pr = ' //                        &
+                                 TRIM( data_output_pr(i) ) // ' is not ava' // &
+                                 'lable for radiation = .FALSE. or ' //        &
+                                 'radiation_scheme /= "rrtmg"'
+                CALL message( 'check_parameters', 'PA0413', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 109
+                dopr_unit(i)  = 'K/h'
+                hom(:,2,109,:)  = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
           CASE DEFAULT
 
@@ -3351 +3415 @@
              unit = 'kg/kg'
 
-
-          CASE ( 'rad_lw_in', 'rad_lw_out', 'rad_sw_in', 'rad_sw_out' )
+          CASE ( 'rad_lw_in', 'rad_lw_out', 'rad_lw_cs_hr', 'rad_lw_hr',       &
+                 'rad_sw_in', 'rad_sw_out', 'rad_sw_cs_hr', 'rad_sw_hr' )
              IF ( .NOT. radiation .OR. radiation_scheme /= 'rrtmg' )  THEN
-                message_string = '"output of "' // TRIM( var ) // '" requi' //  &
+                message_string = '"output of "' // TRIM( var ) // '" requi' // &
                                  'res radiation = .TRUE. and ' //              &
                                  'radiation_scheme = "rrtmg"'
@@ -3395 +3459 @@
 
           CASE ( 'c_liq*', 'c_soil*', 'c_veg*', 'ghf_eb*', 'lai*', 'lwp*',     &
-                 'm_liq_eb*', 'pra*', 'prr*', 'qsws*', 'qsws_eb*',             &
+                 'm_liq_eb*', 'ol*', 'pra*', 'prr*', 'qsws*', 'qsws_eb*',      &
                  'qsws_liq_eb*', 'qsws_soil_eb*', 'qsws_veg_eb*', 'rad_net*',  &
                  'rrtm_aldif*', 'rrtm_aldir*', 'rrtm_asdif*', 'rrtm_asdir*',   &
@@ -3514 +3578 @@
              IF ( TRIM( var ) == 'ghf_eb*')  unit = 'W/m2'
              IF ( TRIM( var ) == 'lai*'   )  unit = 'none'
-             IF ( TRIM( var ) == 'lwp*'   )  unit = 'kg/kg*m'
+             IF ( TRIM( var ) == 'lwp*'   )  unit = 'kg/m2'
+             IF ( TRIM( var ) == 'm_liq_eb*'     )  unit = 'm'
+             IF ( TRIM( var ) == 'ol*'   )   unit = 'm'
              IF ( TRIM( var ) == 'pra*'   )  unit = 'mm'
              IF ( TRIM( var ) == 'prr*'   )  unit = 'mm/s'
@@ -3788 +3854 @@
           constant_diffusion = .TRUE.
 
-          IF ( prandtl_layer )  THEN
-             message_string = 'prandtl_layer is not allowed with fixed ' //    &
-                              'value of km'
+          IF ( constant_flux_layer )  THEN
+             message_string = 'constant_flux_layer is not allowed with fixed ' &
+                              // 'value of km'
              CALL message( 'check_parameters', 'PA0123', 1, 2, 0, 6, 0 )
           ENDIF
@@ -3816 +3882 @@
 
 !
-!-- Check value range for rif
-    IF ( rif_min >= rif_max )  THEN
-       WRITE( message_string, * )  'rif_min = ', rif_min, ' must be less ',    &
-                                   'than rif_max = ', rif_max
+!-- Check value range for zeta = z/L
+    IF ( zeta_min >= zeta_max )  THEN
+       WRITE( message_string, * )  'zeta_min = ', zeta_min, ' must be less ',  &
+                                   'than zeta_max = ', zeta_max
        CALL message( 'check_parameters', 'PA0125', 1, 2, 0, 6, 0 )
     ENDIF
@@ -4235 +4301 @@
 
 !
+!-- Check for valid setting of most_method
+    IF ( TRIM( most_method ) /= 'circular'  .AND.                              &
+         TRIM( most_method ) /= 'newton'  .AND.                                &
+         TRIM( most_method ) /= 'lookup' )  THEN
+       message_string = 'most_method = "' // TRIM( most_method ) //      &
+                        '" is unknown'
+       CALL message( 'check_parameters', 'PA0416', 1, 2, 0, 6, 0 )
+    ENDIF
+
+!
 !-- Check &userpar parameters
     CALL user_check_parameters