Ignore:
Timestamp:
Oct 27, 2011 10:56:21 AM (10 years ago)
Author:
heinze
Message:

Output of liquid water potential temperature in case of cloud_physics=.T. enabled

File:
1 edited

Legend:

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

    r668 r771  
    44! Current revisions:
    55! -----------------
     6! lpt_av
    67!
    78! Former revisions:
     
    7172                e_av = 0.0
    7273
     74             CASE ( 'lpt' )
     75                IF ( .NOT. ALLOCATED( lpt_av ) )  THEN
     76                   ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     77                ENDIF
     78                lpt_av = 0.0
     79
    7380             CASE ( 'lwp*' )
    7481                IF ( .NOT. ALLOCATED( lwp_av ) )  THEN
     
    239246                   DO  k = nzb, nzt+1
    240247                      e_av(k,j,i) = e_av(k,j,i) + e(k,j,i)
     248                   ENDDO
     249                ENDDO
     250             ENDDO
     251
     252          CASE ( 'lpt' )
     253             DO  i = nxlg, nxrg
     254                DO  j = nysg, nyng
     255                   DO  k = nzb, nzt+1
     256                      lpt_av(k,j,i) = lpt_av(k,j,i) + pt(k,j,i)
    241257                   ENDDO
    242258                ENDDO
Note: See TracChangeset for help on using the changeset viewer.