Changeset 4858 for palm/trunk/SOURCE/radiation_model_mod.f90

Timestamp:

Jan 29, 2021 3:27:59 PM (4 years ago)

Author:

raasch

Message:

rrtmg input file names changed, default value for unscheduled calls of the radiation model set to FALSE

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! default value for unscheduled calls of the radiation model set to FALSE,
+! rrtmg input file names changed
+!
+! 4843 2021-01-15 15:22:11Z raasch
 ! local namelist parameter added to switch off the module although the respective module namelist
 ! appears in the namelist file
@@ -664 +668 @@
     INTEGER(iwp) ::  day_of_year                    !< day of the current year
 
-    LOGICAL ::  unscheduled_radiation_calls = .TRUE., &  !< flag parameter indicating whether additional calls
-                                                         !< of the radiation code are allowed
-                constant_albedo = .FALSE.,            &  !< flag parameter indicating whether the albedo may
-                                                         !< change depending on zenith
-                force_radiation_call = .FALSE.,       &  !< flag parameter for unscheduled radiation calls
-                lw_radiation = .TRUE.,                &  !< flag parameter indicating whether longwave radiation shall be calculated
-                radiation = .FALSE.,                  &  !< flag parameter indicating whether the radiation model is used
-                sun_up    = .TRUE.,                   &  !< flag parameter indicating whether the sun is up or down
-                sw_radiation = .TRUE.,                &  !< flag parameter indicating whether shortwave
-                                                         !< radiation shall be calculated
-                sun_direction = .FALSE.,              &  !< flag parameter indicating whether solar direction shall be calculated
-                average_radiation = .FALSE.,          &  !< flag to set the calculation of radiation averaging for the domain
-                radiation_interactions = .FALSE.,     &  !< flag to activiate RTM (TRUE only if vertical
-                                                         !< urban/land surface and trees exist)
-                surface_reflections = .TRUE.,         &  !< flag to switch the calculation of radiation
-                                                         !<interaction between surfaces.
-                                                         !< When it switched off, only the effect of buildings and trees shadow
-                                                         !< will be considered. However fewer SVFs are expected.
-                radiation_interactions_on = .TRUE.       !< namelist flag to force RTM activiation regardless
-                                                         !<to vertical urban/land surface and trees
+    LOGICAL ::  unscheduled_radiation_calls = .FALSE., &  !< flag parameter indicating whether additional calls
+                                                          !< of the radiation code are allowed
+                constant_albedo = .FALSE.,             &  !< flag parameter indicating whether the albedo may
+                                                          !< change depending on zenith
+                force_radiation_call = .FALSE.,        &  !< flag parameter for unscheduled radiation calls
+                lw_radiation = .TRUE.,                 &  !< flag parameter indicating whether longwave radiation shall be calculated
+                radiation = .FALSE.,                   &  !< flag parameter indicating whether the radiation model is used
+                sun_up    = .TRUE.,                    &  !< flag parameter indicating whether the sun is up or down
+                sw_radiation = .TRUE.,                 &  !< flag parameter indicating whether shortwave
+                                                          !< radiation shall be calculated
+                sun_direction = .FALSE.,               &  !< flag parameter indicating whether solar direction shall be calculated
+                average_radiation = .FALSE.,           &  !< flag to set the calculation of radiation averaging for the domain
+                radiation_interactions = .FALSE.,      &  !< flag to activiate RTM (TRUE only if vertical
+                                                          !< urban/land surface and trees exist)
+                surface_reflections = .TRUE.,          &  !< flag to switch the calculation of radiation
+                                                          !<interaction between surfaces.
+                                                          !< When it switched off, only the effect of buildings and trees shadow
+                                                          !< will be considered. However fewer SVFs are expected.
+                radiation_interactions_on = .TRUE.        !< namelist flag to force RTM activiation regardless
+                                                          !<to vertical urban/land surface and trees
 
     REAL(wp), PARAMETER ::  emissivity_atm_clsky = 0.8_wp  !< emissivity of the clear-sky atmosphere
@@ -1922 +1926 @@
  SUBROUTINE radiation_init
 
+#if defined( __rrtmg )
+    USE control_parameters,                                                                        &
+        ONLY:  run_identifier
+#endif
+
     IMPLICIT NONE
 
@@ -3104 +3113 @@
 !
 !--    Initialize RRTMG, before check if files are existent
-       INQUIRE( FILE = 'rrtmg_lw.nc', EXIST = lw_exists )
+       INQUIRE( FILE = 'RRTMG_LW', EXIST = lw_exists )
        IF ( .NOT. lw_exists )  THEN
-          message_string = 'Input file rrtmg_lw.nc &for rrtmg missing. ' //                        &
-                         '&Please provide <jobname>_lsw file in the INPUT directory.'
+          message_string = 'Input file RRTMG_LW for rrtmg model missing.& Please provide ' //      &
+                           TRIM( run_identifier ) // '_lsw file in the INPUT directory.'
           CALL message( 'radiation_init', 'PA0583', 1, 2, 0, 6, 0 )
        ENDIF
-       INQUIRE( FILE = 'rrtmg_sw.nc', EXIST = sw_exists )
+       INQUIRE( FILE = 'RRTMG_SW', EXIST = sw_exists )
        IF ( .NOT. sw_exists )  THEN
-          message_string = 'Input file rrtmg_sw.nc &for rrtmg missing. ' //                        &
-                         '&Please provide <jobname>_rsw file in the INPUT directory.'
+          message_string = 'Input file RRTMG_SW for rrtmg model missing.& Please provide ' //      &
+                           TRIM( run_identifier ) // '_rsw file in the INPUT directory.'
           CALL message( 'radiation_init', 'PA0584', 1, 2, 0, 6, 0 )
        ENDIF
@@ -3124 +3133 @@
        INQUIRE( FILE = 'RAD_SND_DATA', EXIST = snd_exists )
        IF ( .NOT. snd_exists )  THEN
-          rrtm_input_file = 'rrtmg_lw.nc'
+          rrtm_input_file = 'RRTMG_LW'
        ENDIF