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

Modularization of all bulk cloud physics code components

File:
1 edited

Legend:

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

    r3045 r3274  
    2525! -----------------
    2626! $Id$
     27! Modularization of all bulk cloud physics code components
     28!
     29! 3045 2018-05-28 07:55:41Z Giersch
    2730! Code adjusted according to PALM coding standards
    2831!
     
    131134
    132135    USE arrays_3d,                                                             &
    133         ONLY:  p, pt, q, ql, s, ts, u, us, v, w, zu
    134        
    135     USE cloud_parameters,                                                      &
    136         ONLY:  l_d_cp, pt_d_t
    137        
    138     USE constants,                                                             &
    139         ONLY:  pi
    140        
     136        ONLY:  p, pt, q, ql, s, ts, u, us, v, w, zu, d_exner
     137
     138    USE basic_constants_and_equations_mod,                                     &
     139        ONLY:  pi, lv_d_cp
     140
    141141    USE control_parameters,                                                    &
    142         ONLY:  cloud_droplets, cloud_physics, do2d, do3d, humidity, ibc_uv_b,  &
     142        ONLY:  cloud_droplets, do2d, do3d, humidity, ibc_uv_b,                 &
    143143               message_string, nz_do3d, passive_scalar, simulated_time,        &
    144144               threshold
     
    160160       
    161161    USE kinds
     162
     163    USE bulk_cloud_model_mod,                                                  &
     164        ONLY:  bulk_cloud_model
    162165
    163166    USE pegrid
     
    335338
    336339                CASE ( 'pt', 'pt_xy', 'pt_xz', 'pt_yz' )
    337                    IF ( .NOT. cloud_physics ) THEN
     340                   IF ( .NOT. bulk_cloud_model ) THEN
    338341                      DO  i = nxl_dvrp, nxr_dvrp+1
    339342                         DO  j = nys_dvrp, nyn_dvrp+1
     
    347350                         DO  j = nys_dvrp, nyn_dvrp+1
    348351                            DO  k = nzb, nz_do3d
    349                                local_pf(i,j,k) = pt(k,j,i) + l_d_cp *          &
    350                                                  pt_d_t(k) * ql(k,j,i)
     352                               local_pf(i,j,k) = pt(k,j,i) + lv_d_cp *         &
     353                                                 d_exner(k) * ql(k,j,i)
    351354                            ENDDO
    352355                         ENDDO
     
    372375             
    373376                CASE ( 'ql', 'ql_xy', 'ql_xz', 'ql_yz' )
    374                    IF ( cloud_physics  .OR.  cloud_droplets )  THEN
     377                   IF ( bulk_cloud_model  .OR.  cloud_droplets )  THEN
    375378                      DO  i = nxl_dvrp, nxr_dvrp+1
    376379                         DO  j = nys_dvrp, nyn_dvrp+1
     
    381384                      ENDDO
    382385                   ELSE                     
    383                       message_string = 'if cloud_physics = .FALSE. and ' //    &
     386                      message_string = 'if bulk_cloud_model = .FALSE. and ' // &
    384387                                  'cloud_droplets = .FALSE. '
    385388                                  'output of ' // TRIM( output_variable) //    &
Note: See TracChangeset for help on using the changeset viewer.