Ignore:
Timestamp:
Sep 16, 2008 3:29:03 PM (13 years ago)
Author:
raasch
Message:

further adjustments for SGI and other small changes

File:
1 edited

Legend:

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

    r142 r197  
    44! Actual revisions:
    55! -----------------
     6! Spline timeseries splptx etc. removed, timeseries w'u', w'v', w'q' (k=0)
     7! added
    68! Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr)
    79! (like other scalars)
     
    109111       sums_l(nzb+9,pr_palm,0)  = sums_divold_l(sr)  ! old divergence from pres
    110112       sums_l(nzb+10,pr_palm,0) = sums_divnew_l(sr)  ! new divergence from pres
    111 !--    WARNING: next four lines still may have to be adjusted for OpenMP
    112        sums_l(nzb:nzb+2,pr_palm-1,0)    = sums_up_fraction_l(1,1:3,sr)! upstream
    113        sums_l(nzb+3:nzb+5,pr_palm-1,0)  = sums_up_fraction_l(2,1:3,sr)! parts
    114        sums_l(nzb+6:nzb+8,pr_palm-1,0)  = sums_up_fraction_l(3,1:3,sr)! from
    115        sums_l(nzb+9:nzb+11,pr_palm-1,0) = sums_up_fraction_l(4,1:3,sr)! spline
    116113
    117114!
     
    380377             sums_l(nzb+3,pr_palm,tn) = sums_l(nzb+3,pr_palm,tn) + &
    381378                                        ts(j,i)   * rmask(j,i,sr)
     379             IF ( humidity )  THEN
     380                sums_l(nzb+12,pr_palm,tn) = sums_l(nzb+12,pr_palm,tn) + &
     381                                            qs(j,i)   * rmask(j,i,sr)
     382             ENDIF
    382383          ENDDO
    383384       ENDDO
     
    807808!--    above the topography, they are being divided by ngp_2dh(sr)
    808809       sums(nzb:nzb+3,pr_palm)    = sums(nzb:nzb+3,pr_palm)    / &
     810                                    ngp_2dh(sr)
     811       sums(nzb+12,pr_palm)       = sums(nzb+12,pr_palm)       / &    ! qs
    809812                                    ngp_2dh(sr)
    810813!--    eges, e*
     
    882885       hom(:,1,68,sr) = sums(:,68)     ! w*p*
    883886       hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho
     887       hom(:,1,70,sr) = sums(:,70)     ! q*2
    884888
    885889       hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1)
     
    9991003       ts_value(17,sr) = hom(nzb,1,4,sr)            ! pt(0)
    10001004       ts_value(18,sr) = hom(nzb+1,1,4,sr)          ! pt(zp)
    1001        ts_value(19,sr) = hom(nzb+9,1,pr_palm-1,sr)  ! splptx
    1002        ts_value(20,sr) = hom(nzb+10,1,pr_palm-1,sr) ! splpty
    1003        ts_value(21,sr) = hom(nzb+11,1,pr_palm-1,sr) ! splptz
     1005       ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr)    ! u'w'    at k=0
     1006       ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr)    ! v'w'    at k=0
     1007       ts_value(21,sr) = hom(nzb+12,1,pr_palm,sr)   ! w'q'    at k=0
     1008
    10041009       IF ( ts_value(5,sr) /= 0.0 )  THEN
    10051010          ts_value(22,sr) = ts_value(4,sr)**2 / &
Note: See TracChangeset for help on using the changeset viewer.