Changeset 3633 for palm/trunk/SOURCE

Timestamp:

Dec 17, 2018 4:17:57 PM (6 years ago)

Author:

schwenkel

Message:

your commit message

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Include check for rrtmg files
+! 
+! 3630 2018-12-17 11:04:17Z knoop
 ! - fix initialization of date and time after calling zenith
 ! - fix a bug in radiation_solar_pos
@@ -798 +801 @@
 
     LOGICAL :: snd_exists = .FALSE.      !< flag parameter to check whether a user-defined input files exists
+    LOGICAL :: sw_exists = .FALSE.       !< flag parameter to check whether that required rrtmg sw file exists
+    LOGICAL :: lw_exists = .FALSE.       !< flag parameter to check whether that required rrtmg lw file exists
+
 
     REAL(wp), PARAMETER :: mol_mass_air_d_wv = 1.607793_wp !< molecular weight dry air / water vapor
@@ -2595 +2601 @@
 
 !
-!--       Initialize RRTMG
+!--       Initialize RRTMG, before check if files are existent
+          INQUIRE( FILE='rrtmg_lw.nc' // TRIM( coupling_char ), EXIST=lw_exists )
+          IF ( .NOT. lw_exists )  THEN
+             message_string = 'Input file rrtmg_lw.nc' //                &
+                            TRIM( coupling_char ) // ' for rrtmg missing. ' // &
+                            '&Please provide <jobname>_lsw file in the INPUT directory.'
+             CALL message( 'radiation_init', 'PA0583', 1, 2, 0, 6, 0 )
+          ENDIF         
+          INQUIRE( FILE='rrtmg_sw.nc' // TRIM( coupling_char ), EXIST=sw_exists )
+          IF ( .NOT. sw_exists )  THEN
+             message_string = 'Input file rrtmg_sw.nc' //                &
+                            TRIM( coupling_char ) // ' for rrtmg missing. ' // &
+                            '&Please provide <jobname>_rsw file in the INPUT directory.'
+             CALL message( 'radiation_init', 'PA0584', 1, 2, 0, 6, 0 )
+          ENDIF          
+          
           IF ( lw_radiation )  CALL rrtmg_lw_ini ( c_p )
           IF ( sw_radiation )  CALL rrtmg_sw_ini ( c_p )
-
+          
 !
 !--       Set input files for RRTMG