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/data_output_2d.f90

    r4 r72  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Output of precipitation amount/rate and roughness length
    77!
    88! Former revisions:
     
    257257                ENDIF
    258258
     259             CASE ( 'pra*_xy' )        ! 2d-array / integral quantity => no av
     260                CALL exchange_horiz_2d( precipitation_amount )
     261                DO  i = nxl-1, nxr+1
     262                   DO  j = nys-1, nyn+1
     263                      local_pf(i,j,nzb+1) =  precipitation_amount(j,i)
     264                   ENDDO
     265                ENDDO
     266                precipitation_amount = 0.0   ! reset for next integ. interval
     267                resorted = .TRUE.
     268                two_d = .TRUE.
     269                level_z(nzb+1) = zu(nzb+1)
     270
     271             CASE ( 'prr*_xy' )        ! 2d-array
     272                IF ( av == 0 )  THEN
     273                   CALL exchange_horiz_2d( precipitation_rate )
     274                   DO  i = nxl-1, nxr+1
     275                      DO  j = nys-1, nyn+1
     276                         local_pf(i,j,nzb+1) =  precipitation_rate(j,i)
     277                      ENDDO
     278                   ENDDO
     279                ELSE
     280                   CALL exchange_horiz_2d( precipitation_rate_av )
     281                   DO  i = nxl-1, nxr+1
     282                      DO  j = nys-1, nyn+1
     283                         local_pf(i,j,nzb+1) =  precipitation_rate_av(j,i)
     284                      ENDDO
     285                   ENDDO
     286                ENDIF
     287                resorted = .TRUE.
     288                two_d = .TRUE.
     289                level_z(nzb+1) = zu(nzb+1)
     290
    259291             CASE ( 'pt_xy', 'pt_xz', 'pt_yz' )
    260292                IF ( av == 0 )  THEN
     
    419451                ENDIF
    420452                IF ( mode == 'xy' )  level_z = zw
     453
     454             CASE ( 'z0*_xy' )        ! 2d-array
     455                IF ( av == 0 ) THEN
     456                   DO  i = nxl-1, nxr+1
     457                      DO  j = nys-1, nyn+1
     458                         local_pf(i,j,nzb+1) =  z0(j,i)
     459                      ENDDO
     460                   ENDDO
     461                ELSE
     462                   DO  i = nxl-1, nxr+1
     463                      DO  j = nys-1, nyn+1
     464                         local_pf(i,j,nzb+1) =  z0_av(j,i)
     465                      ENDDO
     466                   ENDDO
     467                ENDIF
     468                resorted = .TRUE.
     469                two_d = .TRUE.
     470                level_z(nzb+1) = zu(nzb+1)
    421471
    422472             CASE DEFAULT
Note: See TracChangeset for help on using the changeset viewer.