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

    r39 r72  
    44! Actual revisions:
    55! -----------------
    6 !
     6! precipitation_rate renamed dqdt_precip
    77!
    88! Former revisions:
     
    5252
    5353       INTEGER ::  i, j, k
    54        REAL    :: precipitation_rate
     54       REAL    ::  dqdt_precip
    5555
    5656 
     
    6060
    6161                IF ( ql(k,j,i) > ql_crit )  THEN
    62                    precipitation_rate = prec_time_const * &
    63                                         ( ql(k,j,i) - ql_crit )
     62                   dqdt_precip = prec_time_const * ( ql(k,j,i) - ql_crit )
    6463                ELSE
    65                    precipitation_rate = 0.0
     64                   dqdt_precip = 0.0
    6665                ENDIF
    67                 tend(k,j,i) = tend(k,j,i) + precipitation_rate * l_d_cp * &
    68                                             pt_d_t(k)
     66                tend(k,j,i) = tend(k,j,i) + dqdt_precip * l_d_cp * pt_d_t(k)
    6967
    7068             ENDDO
     
    8785       IMPLICIT NONE
    8886
    89        INTEGER :: i, j, k
    90        REAL    :: precipitation_rate
     87       INTEGER ::  i, j, k
     88       REAL    ::  dqdt_precip
    9189
    9290
     
    9492
    9593          IF ( ql(k,j,i) > ql_crit )  THEN
    96              precipitation_rate = prec_time_const * ( ql(k,j,i) - ql_crit )
     94             dqdt_precip = prec_time_const * ( ql(k,j,i) - ql_crit )
    9795          ELSE
    98              precipitation_rate = 0.0
     96             dqdt_precip = 0.0
    9997          ENDIF
    100           tend(k,j,i) = tend(k,j,i) + precipitation_rate * l_d_cp * pt_d_t(k)
     98          tend(k,j,i) = tend(k,j,i) + dqdt_precip * l_d_cp * pt_d_t(k)
    10199
    102100       ENDDO
Note: See TracChangeset for help on using the changeset viewer.