Changeset 4198 for palm/trunk/SOURCE

Timestamp:

Aug 29, 2019 3:17:48 PM (5 years ago)

Author:

gronemeier

Message:

Add check into radiation_check_parameters if rotation angle is set to 0. Using the radiation model for a rotated model domain is not allowed due to missing implementation within the radiation model.

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Prohibit execution of radiation model if rotation_angle is not zero
+! 
+! 4197 2019-08-29 14:33:32Z suehring
 ! Revise steering of surface albedo initialization when albedo_pars is provided
 ! 
@@ -178 +181 @@
 !>               to INT broke raytracing      -- P. Krc
 !> @todo Optimize radiation_tendency routines
+!> @todo Consider rotated model domains (rotation_angle/=0.0)
 !>
 !> @note Many variables have a leading dummy dimension (0:0) in order to
@@ -1352 +1356 @@
 
        USE control_parameters,                                                 &
-           ONLY: land_surface, message_string, urban_surface
+           ONLY: land_surface, message_string, rotation_angle, urban_surface
 
        USE netcdf_data_input_mod,                                              &
@@ -1446 +1450 @@
           message_string = 'dt_radiation must be > 0.0'  
           CALL message( 'check_parameters', 'PA0591', 1, 2, 0, 6, 0 ) 
+       ENDIF
+!
+!--    Check rotation angle
+       !> @todo Remove this limitation
+       IF ( rotation_angle /= 0.0 )  THEN
+          message_string = 'rotation of the model domain is not considered in the radiation ' //   &
+                           'model.&Using rotation_angle /= 0.0 is not allowed in combination ' //  &
+                           'with the radiation model at the moment!'
+          CALL message( 'check_parameters', 'PA0675', 1, 2, 0, 6, 0 ) 
        ENDIF
  
@@ -5318 +5331 @@
  SUBROUTINE radiation_interaction
 
+    USE control_parameters,                                                    &
+        ONLY:  rotation_angle
+
      IMPLICIT NONE
 
@@ -5333 +5349 @@
      REAL(wp), DIMENSION(0:nsurf_type) ::  costheta           !< direct irradiance factor of solar angle
      REAL(wp), DIMENSION(nz_urban_b:nz_urban_t) ::  pchf_prep          !< precalculated factor for canopy temperature tendency
-     REAL(wp), PARAMETER               :: alpha = 0._wp      !< grid rotation (TODO: synchronize with rotation_angle
-                                                             !< from netcdf_data_input_mod)
      REAL(wp)                          ::  pc_box_area, pc_abs_frac, pc_abs_eff
      REAL(wp)                          ::  asrc               !< area of source face
@@ -5373 +5387 @@
      vnorm(3,:) = idir(:)
      mrot(1, :) = (/ 1._wp,  0._wp,      0._wp      /)
-     mrot(2, :) = (/ 0._wp,  COS(alpha), SIN(alpha) /)
-     mrot(3, :) = (/ 0._wp, -SIN(alpha), COS(alpha) /)
+     mrot(2, :) = (/ 0._wp,  COS(rotation_angle), SIN(rotation_angle) /)
+     mrot(3, :) = (/ 0._wp, -SIN(rotation_angle), COS(rotation_angle) /)
      sunorig = (/ cos_zenith, sun_dir_lat, sun_dir_lon /)
      sunorig = MATMUL(mrot, sunorig)