Ignore:
Timestamp:
Jan 11, 2018 12:01:27 PM (4 years ago)
Author:
suehring
Message:

Output of resistance also urban-type surfaces

File:
1 edited

Legend:

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

    r2718 r2735  
    2525! -----------------
    2626! $Id$
     27! output of r_a moved from land-surface to consider also urban-type surfaces
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    844847                IF ( mode == 'xy' )  level_z = zu
    845848
    846 
     849             CASE ( 'r_a*_xy' )        ! 2d-array
     850                IF ( av == 0 )  THEN
     851                   DO  m = 1, surf_lsm_h%ns
     852                      i                   = surf_lsm_h%i(m)           
     853                      j                   = surf_lsm_h%j(m)
     854                      local_pf(i,j,nzb+1) = surf_lsm_h%r_a(m)
     855                   ENDDO
     856
     857                   DO  m = 1, surf_usm_h%ns
     858                      i   = surf_usm_h%i(m)           
     859                      j   = surf_usm_h%j(m)
     860                      local_pf(i,j,nzb+1) =                                          &
     861                                 ( surf_usm_h%frac(0,m) * surf_usm_h%r_a(m)       +  &
     862                                   surf_usm_h%frac(1,m) * surf_usm_h%r_a_green(m) +  &
     863                                   surf_usm_h%frac(2,m) * surf_usm_h%r_a_window(m) )
     864                   ENDDO
     865
     866                ELSE
     867                   DO  i = nxl, nxr
     868                      DO  j = nys, nyn
     869                         local_pf(i,j,nzb+1) = r_a_av(j,i)
     870                      ENDDO
     871                   ENDDO
     872                ENDIF
     873                resorted       = .TRUE.
     874                two_d          = .TRUE.
     875                level_z(nzb+1) = zu(nzb+1)
    847876
    848877             CASE ( 'rho_ocean_xy', 'rho_ocean_xz', 'rho_ocean_yz' )
Note: See TracChangeset for help on using the changeset viewer.