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

Timestamp:

Sep 27, 2017 10:36:13 AM (7 years ago)

Author:

maronga

Message:

lsm now allows for moist soil and roots below paved surfaces

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Updates pavement types and albedo parameters
+! 
+! 2328 2017-08-03 12:34:22Z maronga
 ! Emissivity can now be set individually for each pixel.
 ! Albedo type can be inferred from land surface model.
@@ -216 +219 @@
 !
 !-- Predefined Land surface classes (albedo_type) after Briegleb (1992)
-    CHARACTER(37), DIMENSION(0:18), PARAMETER :: albedo_type_name = (/      &
+    CHARACTER(37), DIMENSION(0:33), PARAMETER :: albedo_type_name = (/      &
                                    'user defined                         ', & !  0 
                                    'ocean                                ', & !  1
@@ -234 +237 @@
                                    'sea ice                              ', & ! 15 
                                    'snow                                 ', & ! 16
-                                   'pavement/roads                       ', & ! 17
-                                   'bare soil                            '  & ! 18
+                                   'bare soil                            ', & ! 17
+                                   'asphalt/concrete mix                 ', & ! 18
+                                   'asphalt (asphalt concrete)           ', & ! 19
+                                   'concrete (Portland concrete)         ', & ! 20
+                                   'sett                                 ', & ! 21
+                                   'paving stones                        ', & ! 22
+                                   'cobblestone                          ', & ! 23
+                                   'metal                                ', & ! 24
+                                   'wood                                 ', & ! 25
+                                   'gravel                               ', & ! 26
+                                   'fine gravel                          ', & ! 27
+                                   'pebblestone                          ', & ! 28
+                                   'woodchips                            ', & ! 29
+                                   'tartan (sports)                      ', & ! 30
+                                   'artifical turf (sports)              ', & ! 31
+                                   'clay (sports)                        ', & ! 32
+                                   'building (dummy)                     '  & ! 33
                                                          /)
+
 
     INTEGER(iwp) :: albedo_type  = 9999999, & !< Albedo surface type
@@ -292 +311 @@
 !-- Land surface albedos for solar zenith angle of 60° after Briegleb (1992)     
 !-- (shortwave, longwave, broadband):   sw,      lw,      bb,
-    REAL(wp), DIMENSION(0:2,1:18), PARAMETER :: albedo_pars = RESHAPE( (/& 
+    REAL(wp), DIMENSION(0:2,1:33), PARAMETER :: albedo_pars = RESHAPE( (/& 
                                    0.06_wp, 0.06_wp, 0.06_wp,            & !  1
                                    0.09_wp, 0.28_wp, 0.19_wp,            & !  2
@@ -310 +329 @@
                                    0.95_wp, 0.70_wp, 0.82_wp,            & ! 16
                                    0.08_wp, 0.08_wp, 0.08_wp,            & ! 17
-                                   0.17_wp, 0.17_wp, 0.17_wp             & ! 18
-                                 /), (/ 3, 18 /) )
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 18
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 19
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 20
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 21
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 22
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 23
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 24
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 25
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 26
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 27
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 28
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 29
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 30
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 31
+                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 32
+                                   0.17_wp, 0.17_wp, 0.17_wp             & ! 33
+                                 /), (/ 3, 33 /) )
 
     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: &
@@ -1739 +1773 @@
 
 !
-!--        Asphalt
+!--        Bare soil
            ELSEIF ( albedo_type == 17 )  THEN
                  rrtm_aldir(0,:,:) = aldif
@@ -1745 +1779 @@
 
 !
-!--        Bare soil
-           ELSEIF ( albedo_type == 18 )  THEN
+!--        For impermeable surfaces, use values from the lookup table
+           ELSEIF ( albedo_type > 17 )  THEN
                  rrtm_aldir(0,:,:) = aldif
                  rrtm_asdir(0,:,:) = asdif
-
 !
 !--        Land surfaces
            ELSE
-              SELECT CASE ( albedo_type )
+               SELECT CASE ( albedo_type )
 
 !