Changeset 403 for palm/trunk


Ignore:
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:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r402 r403  
    1616Bugfix: calculation of time-averaged surface heatfluxes (sum_up_3d_data)
    1717
    18 sum_up_3d_data, time_integration
     18Bugfix: calculation of precipitation_rate (calc_precipitation)
     19
     20sum_up_3d_data, time_integration, calc_precipitation
    1921
    2022
  • palm/trunk/SOURCE/calc_precipitation.f90

    r77 r403  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix in calculation of precipitation_rate(j,i)
    77!
    88! Former revisions:
     
    4040       MODULE PROCEDURE calc_precipitation_ij
    4141    END INTERFACE calc_precipitation
    42  
     42
    4343 CONTAINS
    4444
     
    6060       REAL    ::  dqdt_precip
    6161
     62       precipitation_rate = 0.0
    6263
    63        precipitation_rate = 0.0
    64  
    6564       DO  i = nxl, nxr
    6665          DO  j = nys, nyn
     
    104103       USE control_parameters
    105104       USE indices
    106    
     105
    107106       IMPLICIT NONE
    108107
     
    110109       REAL    ::  dqdt_precip
    111110
     111       precipitation_rate(j,i) = 0.0
    112112
    113        precipitation_rate(j,i) = 0.0
    114  
    115113!
    116114!--    Ghostpoints are included (although not needed for tend) to avoid a later
     
    126124
    127125!
    128 !--       Precipitation rate in (kg * 0.001) / m**2 / s (because 1kg gives 1 mm)
    129 !          precipitation_rate(j,i) = precipitation_rate(j,i) + dqdt_precip * &
    130 !                                                              dzw(k) * 0.001
    131           precipitation_rate(j,i) = 1.0
     126!--       Precipitation rate in kg / m**2 / s (= mm/s)
     127          precipitation_rate(j,i) = precipitation_rate(j,i) + dqdt_precip * &
     128                                                              dzw(k)
    132129
    133130       ENDDO
    134131
    135132!
    136 !--    Sum up the precipitation amount (unit kg * 0.001 / m**2)
     133!--    Sum up the precipitation amount , unit kg / m**2 (= mm)
    137134       IF ( intermediate_timestep_count == intermediate_timestep_count_max     &
    138135            .AND. ( dt_do2d_xy-time_do2d_xy ) < precipitation_amount_interval )&
Note: See TracChangeset for help on using the changeset viewer.