Ignore:
Timestamp:
Mar 19, 2007 8:20:46 AM (15 years ago)
Author:
raasch
Message:

preliminary changes for precipitation output

File:
1 edited

Legend:

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

    r4 r72  
    44! Actual revisions:
    55! -----------------
    6 !
     6! +sum-up of precipitation rate and roughness length (prr*, z0*)
    77!
    88! Former revisions:
     
    8080                pr_av = 0.0
    8181
     82             CASE ( 'prr*' )
     83                IF ( .NOT. ALLOCATED( precipitation_rate_av ) )  THEN
     84                   ALLOCATE( precipitation_rate_av(nys-1:nyn+1,nxl-1:nxr+1) )
     85                ENDIF
     86                precipitation_rate_av = 0.0
     87
    8288             CASE ( 'pt' )
    8389                IF ( .NOT. ALLOCATED( pt_av ) )  THEN
     
    163169                ENDIF
    164170                w_av = 0.0
     171
     172             CASE ( 'z0*' )
     173                IF ( .NOT. ALLOCATED( z0_av ) )  THEN
     174                   ALLOCATE( z0_av(nys-1:nyn+1,nxl-1:nxr+1) )
     175                ENDIF
     176                z0_av = 0.0
    165177
    166178             CASE DEFAULT
     
    239251             ENDDO
    240252
     253          CASE ( 'pr*' )
     254             DO  i = nxl-1, nxr+1
     255                DO  j = nys-1, nyn+1
     256                   precipitation_rate_av(j,i) = precipitation_rate_av(j,i) + &
     257                                                precipitation_rate(j,i)
     258                ENDDO
     259             ENDDO
     260
    241261          CASE ( 'pt' )
    242262             IF ( .NOT. cloud_physics ) THEN
     
    372392             ENDDO
    373393
     394          CASE ( 'z0*' )
     395             DO  i = nxl-1, nxr+1
     396                DO  j = nys-1, nyn+1
     397                   z0_av(j,i) = z0_av(j,i) + z0(j,i)
     398                ENDDO
     399             ENDDO
     400
    374401          CASE DEFAULT
    375402!
Note: See TracChangeset for help on using the changeset viewer.