Ignore:
Timestamp:
Apr 12, 2018 2:47:44 PM (3 years ago)
Author:
suehring
Message:

Minor revision of static input file checks, bugfix in initialization of surface-fractions in LSM; minor bugfix in initialization of albedo at window-surfaces; for clearer access of albedo and emissivity introduce index for vegetation/wall, pavement/green-wall and water/window surfaces

File:
1 edited

Legend:

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

    r2817 r2963  
    2525! -----------------
    2626! $Id$
     27! Introduce index for vegetation/wall, pavement/green-wall and water/window
     28! surfaces, for clearer access of surface fraction, albedo, emissivity, etc. .
     29!
     30! 2817 2018-02-19 16:32:21Z knoop
    2731! Preliminary gust module interface implemented
    2832!
     
    263267
    264268    USE surface_mod,                                                           &
    265         ONLY:  surf_def_h, surf_lsm_h, surf_usm_h
     269        ONLY:  ind_pav_green, ind_veg_wall, ind_wat_win, surf_def_h,           &
     270               surf_lsm_h, surf_usm_h
    266271
    267272    USE turbulence_closure_mod,                                                &
     
    568573                      i                   = surf_usm_h%i(m)           
    569574                      j                   = surf_usm_h%j(m)
    570                       local_pf(i,j,nzb+1) = surf_usm_h%frac(0,m)     *          &
     575                      local_pf(i,j,nzb+1) = surf_usm_h%frac(ind_veg_wall,m)  *  &
    571576                                            surf_usm_h%wghf_eb(m)        +      &
    572                                             surf_usm_h%frac(1,m)     *          &
     577                                            surf_usm_h%frac(ind_pav_green,m) *  &
    573578                                            surf_usm_h%wghf_eb_green(m)  +      &
    574                                             surf_usm_h%frac(2,m)     *          &
     579                                            surf_usm_h%frac(ind_wat_win,m)   *  &
    575580                                            surf_usm_h%wghf_eb_window(m)
    576581                   ENDDO
     
    915920                      j   = surf_usm_h%j(m)
    916921                      local_pf(i,j,nzb+1) =                                          &
    917                                  ( surf_usm_h%frac(0,m) * surf_usm_h%r_a(m)       +  &
    918                                    surf_usm_h%frac(1,m) * surf_usm_h%r_a_green(m) +  &
    919                                    surf_usm_h%frac(2,m) * surf_usm_h%r_a_window(m) )
     922                                 ( surf_usm_h%frac(ind_veg_wall,m)  *                &
     923                                   surf_usm_h%r_a(m)       +                         &
     924                                   surf_usm_h%frac(ind_pav_green,m) *                &
     925                                   surf_usm_h%r_a_green(m) +                         &
     926                                   surf_usm_h%frac(ind_wat_win,m)   *                &
     927                                   surf_usm_h%r_a_window(m) )
    920928                   ENDDO
    921929                ELSE
Note: See TracChangeset for help on using the changeset viewer.