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/compute_vpt.f90

    r2718 r3274  
    2525! -----------------
    2626! $Id$
     27! Modularization of all bulk cloud physics code components
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    6770
    6871    USE arrays_3d,                                                             &
    69         ONLY:  pt, q, ql, vpt
    70        
     72        ONLY:  pt, q, ql, vpt, d_exner
     73
     74    USE basic_constants_and_equations_mod,                                     &
     75        ONLY:  lv_d_cp
     76
     77    USE control_parameters,                                                    &
     78        ONLY:  cloud_droplets
     79
    7180    USE indices,                                                               &
    7281        ONLY:  nzb, nzt
    73        
    74     USE cloud_parameters,                                                      &
    75         ONLY:  l_d_cp, pt_d_t
    76        
    77     USE control_parameters,                                                    &
    78         ONLY:  cloud_droplets, cloud_physics
    79        
     82
    8083    USE kinds
     84
     85    USE bulk_cloud_model_mod,                                                  &
     86        ONLY:  bulk_cloud_model
    8187
    8288    IMPLICIT NONE
     
    8490    INTEGER(iwp) :: k   !<
    8591
    86     IF ( .NOT. cloud_physics  .AND.  .NOT. cloud_droplets )  THEN
     92    IF ( .NOT. bulk_cloud_model  .AND.  .NOT. cloud_droplets )  THEN
    8793       vpt = pt * ( 1.0_wp + 0.61_wp * q )
    88     ELSE IF (cloud_physics)  THEN
     94    ELSE IF (bulk_cloud_model)  THEN
    8995       DO  k = nzb, nzt+1
    90           vpt(k,:,:) = ( pt(k,:,:) + pt_d_t(k) * l_d_cp * ql(k,:,:) ) *        &
     96          vpt(k,:,:) = ( pt(k,:,:) + d_exner(k) * lv_d_cp * ql(k,:,:) ) *      &
    9197                       ( 1.0_wp + 0.61_wp * q(k,:,:) - 1.61_wp * ql(k,:,:) )
    9298       ENDDO
Note: See TracChangeset for help on using the changeset viewer.