Ignore:
Timestamp:
Sep 24, 2018 3:42:55 PM (3 years ago)
Author:
knoop
Message:

Modularization of all bulk cloud physics code components

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/interaction_droplets_ptq.f90

    r3040 r3274  
    2525! -----------------
    2626! $Id$
     27! Modularization of all bulk cloud physics code components
     28!
     29! 3040 2018-05-25 10:22:08Z schwenkel
    2730! Changed the name specific humidity to mixing ratio
    2831!
     
    7477!
    7578! 799 2011-12-21 17:48:03Z franke
    76 ! Bugfix: pt_d_t(k) was missing in calculation of pt_p
     79! Bugfix: d_exner(k) was missing in calculation of pt_p
    7780!
    7881! RCS Log replace by Id keyword, revision history cleaned up
     
    8992 MODULE interaction_droplets_ptq_mod
    9093 
     94
     95    USE arrays_3d,                                                             &
     96        ONLY:  pt_p, ql_c, q_p, d_exner
     97
     98    USE basic_constants_and_equations_mod,                                     &
     99        ONLY:  lv_d_cp
     100
     101    USE indices,                                                               &
     102        ONLY:  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
     103
     104    USE kinds
     105
     106    USE pegrid
    91107
    92108    PRIVATE
     
    111127    SUBROUTINE i_droplets_ptq
    112128
    113        USE arrays_3d,                                                          &
    114            ONLY:  pt_p, ql_c, q_p
    115            
    116        USE cloud_parameters,                                                   &
    117            ONLY:  l_d_cp, pt_d_t
    118            
    119        USE indices,                                                            &
    120            ONLY:  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
    121            
    122        USE kinds
    123 
    124        USE pegrid
    125 
    126129       IMPLICIT NONE
    127130
     
    140143
    141144                q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i) * flag
    142                 pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k)  &
     145                pt_p(k,j,i) = pt_p(k,j,i) + lv_d_cp * ql_c(k,j,i) * d_exner(k) &
    143146                                                        * flag
    144147             ENDDO
     
    155158!------------------------------------------------------------------------------!
    156159    SUBROUTINE i_droplets_ptq_ij( i, j )
    157 
    158        USE arrays_3d,                                                          &
    159            ONLY:  pt_p, ql_c, q_p
    160 
    161        USE cloud_parameters,                                                   &
    162            ONLY:  l_d_cp, pt_d_t
    163 
    164        USE indices,                                                            &
    165            ONLY:  nzb, nzt, wall_flags_0
    166 
    167        USE kinds,                                                              &
    168            ONLY:  iwp, wp
    169 
    170        USE pegrid
    171160
    172161       IMPLICIT NONE
     
    185174
    186175          q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i) * flag
    187           pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k) * flag
     176          pt_p(k,j,i) = pt_p(k,j,i) + lv_d_cp * ql_c(k,j,i) * d_exner(k) * flag
    188177       ENDDO
    189178
Note: See TracChangeset for help on using the changeset viewer.