Ignore:
Timestamp:
Feb 8, 2018 1:24:35 PM (4 years ago)
Author:
suehring
Message:

Output of ground-heat flux at natural- and urban-type surfaces in one output variable; enable restart data of _av variables that belong to both land- and urban-surface model

File:
1 edited

Legend:

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

    r2743 r2797  
    2525! -----------------
    2626! $Id$
     27! Enable output of ground-heat flux also at urban surfaces.
     28!
     29! 2743 2018-01-12 16:03:39Z suehring
    2730! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2.
    2831!
     
    545548                ENDIF
    546549                IF ( mode == 'xy' )  level_z = zu
     550
     551             CASE ( 'ghf*_xy' )        ! 2d-array
     552                IF ( av == 0 )  THEN
     553                   DO  m = 1, surf_lsm_h%ns
     554                      i                   = surf_lsm_h%i(m)           
     555                      j                   = surf_lsm_h%j(m)
     556                      local_pf(i,j,nzb+1) = surf_lsm_h%ghf(m)
     557                   ENDDO
     558                   DO  m = 1, surf_usm_h%ns
     559                      i                   = surf_usm_h%i(m)           
     560                      j                   = surf_usm_h%j(m)
     561                      local_pf(i,j,nzb+1) = surf_usm_h%frac(0,m)     *          &
     562                                            surf_usm_h%wghf_eb(m)        +      &
     563                                            surf_usm_h%frac(1,m)     *          &
     564                                            surf_usm_h%wghf_eb_green(m)  +      &
     565                                            surf_usm_h%frac(2,m)     *          &
     566                                            surf_usm_h%wghf_eb_window(m)
     567                   ENDDO
     568                ELSE
     569                   DO  i = nxl, nxr
     570                      DO  j = nys, nyn
     571                         local_pf(i,j,nzb+1) = ghf_av(j,i)
     572                      ENDDO
     573                   ENDDO
     574                ENDIF
     575
     576                resorted = .TRUE.
     577                two_d = .TRUE.
     578                level_z(nzb+1) = zu(nzb+1)
    547579
    548580             CASE ( 'ol*_xy' )        ! 2d-array
Note: See TracChangeset for help on using the changeset viewer.