Changeset 1788 for palm/trunk/SOURCE/radiation_model.f90

Timestamp:

Mar 10, 2016 11:01:04 AM (8 years ago)

Author:

maronga

Message:

added support for water and paved surfaced in land surface model / minor changes

File:

• palm/trunk/SOURCE/radiation_model.f90 (modified) (9 diffs)

palm/trunk/SOURCE/radiation_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added new albedo class for pavements / roads.
 !
 ! Former revisions:
@@ -112 +112 @@
 
 #if defined ( __rrtmg )
+
+!     USE netcdf_interface,                                                      &
+!         ONLY:  nc_stat, netcdf_handle_error
+
     USE parrrsw,                                                               &
         ONLY:  naerec, nbndsw
@@ -139 +143 @@
 !
 !-- Predefined Land surface classes (albedo_type) after Briegleb (1992)
-    CHARACTER(37), DIMENSION(0:16), PARAMETER :: albedo_type_name = (/      &
+    CHARACTER(37), DIMENSION(0:17), PARAMETER :: albedo_type_name = (/      &
                                    'user defined                         ', & !  0 
                                    'ocean                                ', & !  1
@@ -156 +160 @@
                                    'land ice                             ', & ! 14 
                                    'sea ice                              ', & ! 15 
-                                   'snow                                 '  & ! 16
+                                   'snow                                 ', & ! 16
+                                   'pavement/roads                       '  & ! 17
                                                          /)
 
@@ -163 +168 @@
                     day_init     = 172,  & !< day of the year at model start (21/06)
                     dots_rad     = 0       !< starting index for timeseries output 
-
-
-
-
-
 
     LOGICAL ::  unscheduled_radiation_calls = .TRUE., & !< flag parameter indicating whether additional calls of the radiation code are allowed
@@ -214 +214 @@
 !-- Land surface albedos for solar zenith angle of 60° after Briegleb (1992)     
 !-- (shortwave, longwave, broadband):   sw,      lw,      bb,
-    REAL(wp), DIMENSION(0:2,1:16), PARAMETER :: albedo_pars = RESHAPE( (/& 
+    REAL(wp), DIMENSION(0:2,1:17), PARAMETER :: albedo_pars = RESHAPE( (/& 
                                    0.06_wp, 0.06_wp, 0.06_wp,            & !  1
                                    0.09_wp, 0.28_wp, 0.19_wp,            & !  2
@@ -230 +230 @@
                                    0.90_wp, 0.65_wp, 0.77_wp,            & ! 14
                                    0.90_wp, 0.65_wp, 0.77_wp,            & ! 15
-                                   0.95_wp, 0.70_wp, 0.82_wp             & ! 16
-                                 /), (/ 3, 16 /) )
+                                   0.95_wp, 0.70_wp, 0.82_wp,            & ! 16
+                                   0.08_wp, 0.08_wp, 0.08_wp             & ! 17
+                                 /), (/ 3, 17 /) )
 
     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: &
@@ -275 +276 @@
                     rrtm_iaer = 0,     & !< aerosol option flag (0: no aerosol layers, for lw only: 6 (requires setting of rrtm_sw_ecaer), 10: one or more aerosol layers (not implemented)
                     rrtm_idrv = 1        !< longwave upward flux calculation option (0,1)
+
+    INTEGER(iwp) :: nc_stat !< local variable for storin the result of netCDF calls for error message handling
 
     LOGICAL :: snd_exists = .FALSE.      !< flag parameter to check whether a user-defined input files exists
@@ -1030 +1033 @@
                  rrtm_aldir(0,:,:) = aldif
                  rrtm_asdir(0,:,:) = asdif
+
+!
+!--        Asphalt
+           ELSEIF ( albedo_type == 17 )  THEN
+                 rrtm_aldir(0,:,:) = aldif
+                 rrtm_asdir(0,:,:) = asdif
 !
 !--        Land surfaces