Changeset 3127


Ignore:
Timestamp:
Jul 15, 2018 8:01:25 AM (3 years ago)
Author:
maronga
Message:

bugfix for high vegetation surface temperatures

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r3091 r3127  
    2525! -----------------
    2626! $Id$
     27! Some adjustments for pavement parameters
     28! Limit magnus formula to avoid negative q_s (leads to model crash)
     29!
     30! 3091 2018-06-28 16:20:35Z suehring
    2731! Add check for local roughness length not exceeding surface-layer height and
    2832! limit roughness length where necessary.
     
    865869!--   0.0-0.01, 0.01-0.03, 0.03-0.07, 0.07-0.15, 0.15-0.30, 0.30-0.50,    0.50-1.25,    1.25-3.00
    866870    REAL(wp), DIMENSION(0:7,1:15), PARAMETER :: pavement_subsurface_pars_1 = RESHAPE( (/ &
    867        1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & !  1
    868        1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & !  2
    869        1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & !  3
     871       0.75_wp,   0.75_wp,   0.75_wp,   0.75_wp,   0.75_wp,   0.75_wp, 9999999.9_wp, 9999999.9_wp, & !  1
     872       0.75_wp,   0.75_wp,   0.75_wp,   0.75_wp,   0.75_wp,   0.75_wp, 9999999.9_wp, 9999999.9_wp, & !  2
     873       0.89_wp,   0.89_wp,   0.89_wp,   0.89_wp,   0.89_wp,   0.89_wp, 9999999.9_wp, 9999999.9_wp, & !  3
    870874       1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & !  4
    871875       1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & !  5
     
    873877       1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & !  7
    874878       1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & !  8
    875        1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & !  9
     879       0.70_wp,   0.70_wp,   0.70_wp,   0.70_wp,   0.70_wp,   0.70_wp, 9999999.9_wp, 9999999.9_wp, & !  9
    876880       1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & ! 10
    877881       1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp,   1.00_wp, 9999999.9_wp, 9999999.9_wp, & ! 11
     
    888892       1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 9999999.9_wp, 9999999.9_wp, & !  1
    889893       1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 9999999.9_wp, 9999999.9_wp, & !  2
    890        1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 9999999.9_wp, 9999999.9_wp, & !  3
     894       1.76E6_wp, 1.76E6_wp, 1.76E6_wp, 1.76E6_wp, 1.76E6_wp, 1.76E6_wp, 9999999.9_wp, 9999999.9_wp, & !  3
    891895       1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 9999999.9_wp, 9999999.9_wp, & !  4
    892896       1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 1.94E6_wp, 9999999.9_wp, 9999999.9_wp, & !  5
     
    19231927!
    19241928!--    Calculate water vapour pressure at saturation and convert to hPa
    1925        e_s = 0.01_wp * magnus( surf_t_surface%var_1d(m) )
     1929!--    The magnus formula is limited to temperatures up to 333.15 K to
     1930!--    avoid negative values of q_s
     1931       e_s = 0.01_wp * magnus( MIN(surf_t_surface%var_1d(m), 333.15_wp) )
    19261932
    19271933!
     
    22792285!
    22802286!--       Calculate water vapour pressure at saturation and convert to hPa
    2281           e_s = 0.01_wp * magnus( surf_t_surface_p%var_1d(m) )                   
     2287          e_s = 0.01_wp * magnus( MIN(surf_t_surface%var_1d(m), 333.15_wp) )                   
    22822288
    22832289!
  • palm/trunk/SOURCE/radiation_model_mod.f90

    r3123 r3127  
    2828! -----------------
    2929! $Id$
     30! A few pavement parameters updated.
     31!
     32! 3123 2018-07-12 16:21:53Z suehring
    3033! Correct working precision for INTEGER number
    3134!
     
    566569                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 18
    567570                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 19
    568                                    0.17_wp, 0.17_wp, 0.17_wp,            & ! 20
     571                                   0.30_wp, 0.30_wp, 0.30_wp,            & ! 20
    569572                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 21
    570573                                   0.17_wp, 0.17_wp, 0.17_wp,            & ! 22
Note: See TracChangeset for help on using the changeset viewer.