Changeset 2703 for palm/trunk/SOURCE


Ignore:
Timestamp:
Dec 15, 2017 8:12:38 PM (7 years ago)
Author:
maronga
Message:

workaround to keep aerodyn. resistances positive in urban surface model

File:
1 edited

Legend:

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

    r2696 r2703  
    2121! Current revisions:
    2222! ------------------
    23 !
     23! Workaround for calculation of r_a
    2424!
    2525! Former revisions:
     
    68906890!--        pt, us, ts are not available for the prognostic time step,
    68916891!--        data from the last time step is used here.
    6892                
    6893            r_a = ( surf_usm_h%pt1(m) - t_surf_h(m) / exn(k) ) /                              &
    6894                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
    6895            r_a_window = ( surf_usm_h%pt1(m) - t_surf_window_h(m) / exn(k) ) /                &
    6896                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
    6897            r_a_green = ( surf_usm_h%pt1(m) - t_surf_green_h(m) / exn(k) ) /                  &
    6898                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
     6892
     6893!--        Workarounf to keep resistances positive
     6894           r_a = ( surf_usm_h%pt1(m) - surf_usm_h%pt_surface(m) ) /                            &
     6895                 ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )   
     6896           r_a_window = r_a
     6897           r_a_green  = r_a
     6898!            r_a = ( surf_usm_h%pt1(m) - t_surf_h(m) / exn(k) ) /                              &
     6899!                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
     6900!            r_a_window = ( surf_usm_h%pt1(m) - t_surf_window_h(m) / exn(k) ) /                &
     6901!                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
     6902!            r_a_green = ( surf_usm_h%pt1(m) - t_surf_green_h(m) / exn(k) ) /                  &
     6903!                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
    68996904               
    69006905!--        make sure that the resistance does not drop to zero
Note: See TracChangeset for help on using the changeset viewer.