Changeset 72 for palm/trunk/SOURCE/impact_of_latent_heat.f90

Timestamp:

Mar 19, 2007 8:20:46 AM (17 years ago)

Author:

raasch

Message:

preliminary changes for precipitation output

File:

• palm/trunk/SOURCE/impact_of_latent_heat.f90 (modified) (5 diffs)

palm/trunk/SOURCE/impact_of_latent_heat.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! precipitation_rate renamed dqdt_precip
 !
 ! Former revisions:
@@ -52 +52 @@
 
        INTEGER ::  i, j, k
-       REAL    :: precipitation_rate
+       REAL    ::  dqdt_precip
 
  
@@ -60 +60 @@
 
                 IF ( ql(k,j,i) > ql_crit )  THEN
-                   precipitation_rate = prec_time_const * &
-                                        ( ql(k,j,i) - ql_crit )
+                   dqdt_precip = prec_time_const * ( ql(k,j,i) - ql_crit )
                 ELSE
-                   precipitation_rate = 0.0
+                   dqdt_precip = 0.0
                 ENDIF
-                tend(k,j,i) = tend(k,j,i) + precipitation_rate * l_d_cp * &
-                                            pt_d_t(k)
+                tend(k,j,i) = tend(k,j,i) + dqdt_precip * l_d_cp * pt_d_t(k)
 
              ENDDO
@@ -87 +85 @@
        IMPLICIT NONE
 
-       INTEGER :: i, j, k
-       REAL    :: precipitation_rate
+       INTEGER ::  i, j, k
+       REAL    ::  dqdt_precip
 
 
@@ -94 +92 @@
 
           IF ( ql(k,j,i) > ql_crit )  THEN
-             precipitation_rate = prec_time_const * ( ql(k,j,i) - ql_crit )
+             dqdt_precip = prec_time_const * ( ql(k,j,i) - ql_crit )
           ELSE
-             precipitation_rate = 0.0
+             dqdt_precip = 0.0
           ENDIF
-          tend(k,j,i) = tend(k,j,i) + precipitation_rate * l_d_cp * pt_d_t(k)
+          tend(k,j,i) = tend(k,j,i) + dqdt_precip * l_d_cp * pt_d_t(k)
 
        ENDDO