Ignore:
Timestamp:
Mar 22, 2007 9:54:05 AM (14 years ago)
Author:
raasch
Message:

preliminary update for changes concerning non-cyclic boundary conditions

File:
1 edited

Legend:

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

    r51 r75  
    55! -----------------
    66! Collection of time series quantities moved from routine flow_statistics to
    7 ! here, routine user_statistics is called for each statistic region
     7! here, routine user_statistics is called for each statistic region,
     8! moisture renamed humidity
    89!
    910! Former revisions:
     
    115116!--    total water content, specific humidity and liquid water potential
    116117!--    temperature
    117        IF ( moisture )  THEN
     118       IF ( humidity )  THEN
    118119          !$OMP DO
    119120          DO  i = nxl, nxr
     
    164165             sums_l(:,2,0) = sums_l(:,2,0) + sums_l(:,2,i)
    165166             sums_l(:,4,0) = sums_l(:,4,0) + sums_l(:,4,i)
    166              IF ( moisture )  THEN
     167             IF ( humidity )  THEN
    167168                sums_l(:,41,0) = sums_l(:,41,0) + sums_l(:,41,i)
    168169                sums_l(:,44,0) = sums_l(:,44,0) + sums_l(:,44,i)
     
    187188       CALL MPI_ALLREDUCE( sums_l(nzb,4,0), sums(nzb,4), nzt+2-nzb, MPI_REAL, &
    188189                           MPI_SUM, comm2d, ierr )
    189        IF ( moisture ) THEN
     190       IF ( humidity ) THEN
    190191          CALL MPI_ALLREDUCE( sums_l(nzb,44,0), sums(nzb,44), nzt+2-nzb, &
    191192                              MPI_REAL, MPI_SUM, comm2d, ierr )
     
    208209       sums(:,2) = sums_l(:,2,0)
    209210       sums(:,4) = sums_l(:,4,0)
    210        IF ( moisture ) THEN
     211       IF ( humidity ) THEN
    211212          sums(:,44) = sums_l(:,44,0)
    212213          sums(:,41) = sums_l(:,41,0)
     
    231232!
    232233!--    Humidity and cloud parameters
    233        IF ( moisture ) THEN
     234       IF ( humidity ) THEN
    234235          sums(:,44) = sums(:,44) / ngp_2dh_outer(:,sr)
    235236          sums(:,41) = sums(:,41) / ngp_2dh_outer(:,sr)
     
    365366!
    366367!--             Buoyancy flux, water flux (humidity flux) w"q"
    367                 IF ( moisture ) THEN
     368                IF ( humidity ) THEN
    368369                   sums_l(k,45,tn) = sums_l(k,45,tn)                           &
    369370                                         - 0.5 * ( kh(k,j,i) + kh(k+1,j,i) )   &
     
    407408                sums_l(nzb,61,tn) = sums_l(nzb,61,tn) + &
    408409                                    0.0 * rmask(j,i,sr)           ! v"pt"
    409                 IF ( moisture )  THEN
     410                IF ( humidity )  THEN
    410411                   sums_l(nzb,48,tn) = sums_l(nzb,48,tn) + &
    411412                                       qsws(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
     
    437438                sums_l(nzt,61,tn) = sums_l(nzt,61,tn) + &
    438439                                    0.0 * rmask(j,i,sr)           ! v"pt"
    439                 IF ( moisture )  THEN
     440                IF ( humidity )  THEN
    440441                   sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
    441442                                       qswst(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
     
    499500!--             Buoyancy flux, water flux, humidity flux and liquid water
    500501!--             content
    501                 IF ( moisture )  THEN
     502                IF ( humidity )  THEN
    502503                   pts = 0.5 * ( vpt(k,j,i)   - hom(k,1,44,sr) + &
    503504                                 vpt(k+1,j,i) - hom(k+1,1,44,sr) )
Note: See TracChangeset for help on using the changeset viewer.