Changeset 403 for palm/trunk/SOURCE

Timestamp:

Oct 22, 2009 1:57:16 PM (15 years ago)

Author:

franke

Message:

Bugfix in calculation of precipitation_rate

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (1 diff)
• calc_precipitation.f90 (modified) (6 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -16 +16 @@
 Bugfix: calculation of time-averaged surface heatfluxes (sum_up_3d_data)
 
-sum_up_3d_data, time_integration
+Bugfix: calculation of precipitation_rate (calc_precipitation)
+
+sum_up_3d_data, time_integration, calc_precipitation
 
 

palm/trunk/SOURCE/calc_precipitation.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Bugfix in calculation of precipitation_rate(j,i)
 !
 ! Former revisions:
@@ -40 +40 @@
        MODULE PROCEDURE calc_precipitation_ij
     END INTERFACE calc_precipitation
- 
+
  CONTAINS
 
@@ -60 +60 @@
        REAL    ::  dqdt_precip
 
+       precipitation_rate = 0.0
 
-       precipitation_rate = 0.0
- 
        DO  i = nxl, nxr
           DO  j = nys, nyn
@@ -104 +103 @@
        USE control_parameters
        USE indices
-   
+
        IMPLICIT NONE
 
@@ -110 +109 @@
        REAL    ::  dqdt_precip
 
+       precipitation_rate(j,i) = 0.0
 
-       precipitation_rate(j,i) = 0.0
- 
 !
 !--    Ghostpoints are included (although not needed for tend) to avoid a later
@@ -126 +124 @@
 
 !
-!--       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
+!--       Precipitation rate in kg / m**2 / s (= mm/s)
+          precipitation_rate(j,i) = precipitation_rate(j,i) + dqdt_precip * &
+                                                              dzw(k)
 
        ENDDO
 
 !
-!--    Sum up the precipitation amount (unit kg * 0.001 / m**2)
+!--    Sum up the precipitation amount , unit kg / m**2 (= mm)
        IF ( intermediate_timestep_count == intermediate_timestep_count_max     &
             .AND. ( dt_do2d_xy-time_do2d_xy ) < precipitation_amount_interval )&