Changeset 3256 for palm/trunk


Ignore:
Timestamp:
Sep 17, 2018 12:20:07 PM (6 years ago)
Author:
suehring
Message:

Enable initialization of z0q for vegetation, pavement and water surfaces via namelist input

File:
1 edited

Legend:

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

    r3248 r3256  
    2525! -----------------
    2626! $Id$
     27! Enable initialization of z0q for vegetation, pavement and water surfaces via
     28! namelist input.
     29!
     30! 3248 2018-09-14 09:42:06Z sward
    2731! Minor formating changes
    2832!
     
    31403144          IF ( z0h_vegetation == 9999999.9_wp )  THEN
    31413145             z0h_vegetation = vegetation_pars(ind_v_z0qh,vegetation_type)
    3142           ENDIF   
     3146          ENDIF
     3147         
     3148          IF ( z0q_vegetation == 9999999.9_wp )  THEN
     3149             z0q_vegetation = vegetation_pars(ind_v_z0qh,vegetation_type)
     3150          ENDIF
    31433151         
    31443152          IF ( lambda_surface_stable == 9999999.9_wp )  THEN
     
    31773185             surf_lsm_h%z0(m)               = z0_vegetation
    31783186             surf_lsm_h%z0h(m)              = z0h_vegetation
    3179              surf_lsm_h%z0q(m)              = z0h_vegetation
     3187             surf_lsm_h%z0q(m)              = z0q_vegetation
    31803188             surf_lsm_h%lambda_surface_s(m) = lambda_surface_stable
    31813189             surf_lsm_h%lambda_surface_u(m) = lambda_surface_unstable
     
    32063214                surf_lsm_v(l)%z0(m)               = z0_vegetation
    32073215                surf_lsm_v(l)%z0h(m)              = z0h_vegetation
    3208                 surf_lsm_v(l)%z0q(m)              = z0h_vegetation
     3216                surf_lsm_v(l)%z0q(m)              = z0q_vegetation
    32093217                surf_lsm_v(l)%lambda_surface_s(m) = lambda_surface_stable
    32103218                surf_lsm_v(l)%lambda_surface_u(m) = lambda_surface_unstable
     
    34613469             z0h_water = water_pars(ind_w_z0h,water_type)       
    34623470          ENDIF 
     3471         
     3472          IF ( z0q_water == 9999999.9_wp )  THEN
     3473             z0q_water = water_pars(ind_w_z0h,water_type)       
     3474          ENDIF
    34633475
    34643476          IF ( albedo_type == 9999999  .AND.  albedo == 9999999.9_wp )  THEN
     
    34793491             surf_lsm_h%z0(m)               = z0_water
    34803492             surf_lsm_h%z0h(m)              = z0h_water
    3481              surf_lsm_h%z0q(m)              = z0h_water
     3493             surf_lsm_h%z0q(m)              = z0q_water
    34823494             surf_lsm_h%lambda_surface_s(m) = 1.0E10_wp
    34833495             surf_lsm_h%lambda_surface_u(m) = 1.0E10_wp               
     
    35003512                surf_lsm_v(l)%z0(m)               = z0_water
    35013513                surf_lsm_v(l)%z0h(m)              = z0h_water
    3502                 surf_lsm_v(l)%z0q(m)              = z0h_water
     3514                surf_lsm_v(l)%z0q(m)              = z0q_water
    35033515                surf_lsm_v(l)%lambda_surface_s(m) = 1.0E10_wp
    35043516                surf_lsm_v(l)%lambda_surface_u(m) = 1.0E10_wp               
     
    37123724             z0h_pavement = pavement_pars(ind_p_z0h,pavement_type)
    37133725          ENDIF
     3726         
     3727          IF ( z0q_pavement == 9999999.9_wp )  THEN
     3728             z0q_pavement = pavement_pars(ind_p_z0h,pavement_type)
     3729          ENDIF
    37143730
    37153731          IF ( pavement_heat_conduct == 9999999.9_wp )  THEN
     
    37553771             surf_lsm_h%z0(m)                  = z0_pavement
    37563772             surf_lsm_h%z0h(m)                 = z0h_pavement
    3757              surf_lsm_h%z0q(m)                 = z0h_pavement
     3773             surf_lsm_h%z0q(m)                 = z0q_pavement
    37583774             surf_lsm_h%lambda_surface_s(m)    = pavement_heat_conduct         &
    37593775                                                  * ddz_soil(nzb_soil)         &
     
    37933809                surf_lsm_v(l)%z0(m)                  = z0_pavement
    37943810                surf_lsm_v(l)%z0h(m)                 = z0h_pavement
    3795                 surf_lsm_v(l)%z0q(m)                 = z0h_pavement
     3811                surf_lsm_v(l)%z0q(m)                 = z0q_pavement
    37963812                surf_lsm_v(l)%lambda_surface_s(m)    = pavement_heat_conduct   &
    37973813                                                  * ddz_soil(nzb_soil)         &
Note: See TracChangeset for help on using the changeset viewer.