Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/radiation_model_mod.f90

    r2328 r2504  
    2525! -----------------
    2626! $Id$
     27! Updates pavement types and albedo parameters
     28!
     29! 2328 2017-08-03 12:34:22Z maronga
    2730! Emissivity can now be set individually for each pixel.
    2831! Albedo type can be inferred from land surface model.
     
    216219!
    217220!-- Predefined Land surface classes (albedo_type) after Briegleb (1992)
    218     CHARACTER(37), DIMENSION(0:18), PARAMETER :: albedo_type_name = (/      &
     221    CHARACTER(37), DIMENSION(0:33), PARAMETER :: albedo_type_name = (/      &
    219222                                   'user defined                         ', & !  0
    220223                                   'ocean                                ', & !  1
     
    234237                                   'sea ice                              ', & ! 15
    235238                                   'snow                                 ', & ! 16
    236                                    'pavement/roads                       ', & ! 17
    237                                    'bare soil                            '  & ! 18
     239                                   'bare soil                            ', & ! 17
     240                                   'asphalt/concrete mix                 ', & ! 18
     241                                   'asphalt (asphalt concrete)           ', & ! 19
     242                                   'concrete (Portland concrete)         ', & ! 20
     243                                   'sett                                 ', & ! 21
     244                                   'paving stones                        ', & ! 22
     245                                   'cobblestone                          ', & ! 23
     246                                   'metal                                ', & ! 24
     247                                   'wood                                 ', & ! 25
     248                                   'gravel                               ', & ! 26
     249                                   'fine gravel                          ', & ! 27
     250                                   'pebblestone                          ', & ! 28
     251                                   'woodchips                            ', & ! 29
     252                                   'tartan (sports)                      ', & ! 30
     253                                   'artifical turf (sports)              ', & ! 31
     254                                   'clay (sports)                        ', & ! 32
     255                                   'building (dummy)                     '  & ! 33
    238256                                                         /)
     257
    239258
    240259    INTEGER(iwp) :: albedo_type  = 9999999, & !< Albedo surface type
     
    292311!-- Land surface albedos for solar zenith angle of 60° after Briegleb (1992)     
    293312!-- (shortwave, longwave, broadband):   sw,      lw,      bb,
    294     REAL(wp), DIMENSION(0:2,1:18), PARAMETER :: albedo_pars = RESHAPE( (/&
     313    REAL(wp), DIMENSION(0:2,1:33), PARAMETER :: albedo_pars = RESHAPE( (/&
    295314                                   0.06_wp, 0.06_wp, 0.06_wp,            & !  1
    296315                                   0.09_wp, 0.28_wp, 0.19_wp,            & !  2
     
    310329                                   0.95_wp, 0.70_wp, 0.82_wp,            & ! 16
    311330                                   0.08_wp, 0.08_wp, 0.08_wp,            & ! 17
    312                                    0.17_wp, 0.17_wp, 0.17_wp             & ! 18
    313                                  /), (/ 3, 18 /) )
     331                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 18
     332                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 19
     333                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 20
     334                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 21
     335                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 22
     336                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 23
     337                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 24
     338                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 25
     339                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 26
     340                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 27
     341                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 28
     342                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 29
     343                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 30
     344                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 31
     345                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 32
     346                                   0.17_wp, 0.17_wp, 0.17_wp             & ! 33
     347                                 /), (/ 3, 33 /) )
    314348
    315349    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: &
     
    17391773
    17401774!
    1741 !--        Asphalt
     1775!--        Bare soil
    17421776           ELSEIF ( albedo_type == 17 )  THEN
    17431777                 rrtm_aldir(0,:,:) = aldif
     
    17451779
    17461780!
    1747 !--        Bare soil
    1748            ELSEIF ( albedo_type == 18 )  THEN
     1781!--        For impermeable surfaces, use values from the lookup table
     1782           ELSEIF ( albedo_type > 17 )  THEN
    17491783                 rrtm_aldir(0,:,:) = aldif
    17501784                 rrtm_asdir(0,:,:) = asdif
    1751 
    17521785!
    17531786!--        Land surfaces
    17541787           ELSE
    1755               SELECT CASE ( albedo_type )
     1788               SELECT CASE ( albedo_type )
    17561789
    17571790!
Note: See TracChangeset for help on using the changeset viewer.