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

Timestamp:

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

Author:

raasch

Message:

preliminary changes for precipitation output

File:

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

palm/trunk/SOURCE/calc_precipitation.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Precipitation rate and amount are calculated/stored,
+! + module control_parameters
 !
 ! Former revisions:
@@ -48 +49 @@
        USE cloud_parameters
        USE constants
+       USE control_parameters
        USE indices
 
@@ -53 +55 @@
 
        INTEGER ::  i, j, k
-       REAL    :: precipitation_rate
+       REAL    ::  dqdt_precip
 
+
+       precipitation_rate = 0.0
  
        DO  i = nxl, nxr
@@ -61 +65 @@
 
                 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
+                tend(k,j,i) = tend(k,j,i) - dqdt_precip
+!
+!--             Precipitation rate in (kg * 0.001) / m**2 / s (because 1kg
+!--             gives 1 mm)
+                precipitation_rate(j,i) = precipitation_rate(j,i) + &
+                                          dqdt_precip * dzw(k) * 0.001
 
              ENDDO
+!
+!--          Sum up the precipitation amount (unit kg * 0.001 / m**2)
+             IF ( intermediate_timestep_count ==         &
+                  intermediate_timestep_count_max  .AND. &
+                  ( dt_do2d_xy-time_do2d_xy ) < precipitation_amount_interval )&
+             THEN
+                precipitation_amount(j,i) = precipitation_amount(j,i) + &
+                                            precipitation_rate(j,i) * dt_3d
+             ENDIF
           ENDDO
        ENDDO
@@ -83 +100 @@
        USE cloud_parameters
        USE constants
+       USE control_parameters
        USE indices
    
        IMPLICIT NONE
 
-       INTEGER :: i, j, k
-       REAL    :: precipitation_rate
+       INTEGER ::  i, j, k
+       REAL    ::  dqdt_precip
 
 
+       precipitation_rate(j,i) = 0.0
+ 
+!
+!--    Ghostpoints are included (although not needed for tend) to avoid a later
+!--    exchange of these data for the precipitation amount/rate arrays
        DO  k = nzb_2d(j,i)+1, nzt
 
           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
+          tend(k,j,i) = tend(k,j,i) - dqdt_precip
+
+!
+!--       Precipitation rate in (kg * 0.001) / m**2 / s (because 1kg gives 1 mm)
+!          precipitation_rate(j,i) = precipitation_rate(j,i) + dqdt_precip * &
+!                                                              dzw(k) * 0.001
+          precipitation_rate(j,i) = 1.0
 
        ENDDO
+
+!
+!--    Sum up the precipitation amount (unit kg * 0.001 / m**2)
+       IF ( intermediate_timestep_count == intermediate_timestep_count_max     &
+            .AND. ( dt_do2d_xy-time_do2d_xy ) < precipitation_amount_interval )&
+       THEN
+          precipitation_amount(j,i) = precipitation_amount(j,i) + &
+                                      precipitation_rate(j,i) * dt_3d
+       ENDIF
 
     END SUBROUTINE calc_precipitation_ij