Ignore:
Timestamp:
Oct 16, 2017 12:41:56 PM (6 years ago)
Author:
schwenkel
Message:

extended by cloud_droplets option

File:
1 edited

Legend:

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

    r2321 r2547  
    2525! -----------------
    2626! $Id$
     27! extended by cloud_droplets option
     28!
     29! 2321 2017-07-24 15:57:07Z schwenkel
    2730! Bugfix: Correct index in lookup table for Obukhov length
    2831!
     
    201204
    202205    USE control_parameters,                                                    &
    203         ONLY:  cloud_physics, constant_heatflux, constant_scalarflux,          &
    204                constant_waterflux, coupling_mode, g, humidity, ibc_e_b,        &
    205                ibc_pt_b, initializing_actions, kappa,                          &
     206        ONLY:  cloud_droplets, cloud_physics, constant_heatflux,               &
     207               constant_scalarflux, constant_waterflux, coupling_mode, g,      &
     208               humidity, ibc_e_b, ibc_pt_b, initializing_actions, kappa,       &
    206209               intermediate_timestep_count, intermediate_timestep_count_max,   &
    207210               land_surface, large_scale_forcing, lsf_surf,                    &
     
    291294!--    temperature and specific humidity at first grid level from the fields pt
    292295!--    and q
    293        IF ( cloud_physics )  THEN
     296       IF ( cloud_physics  .OR.  cloud_droplets )  THEN
    294297!
    295298!--       First call for horizontal default-type surfaces (l=0 - upward facing,
     
    13961399
    13971400             ENDDO
    1398           ELSEIF ( cloud_physics )  THEN
     1401          ELSEIF ( cloud_physics  .OR.  cloud_droplets )  THEN
    13991402             !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo )
    14001403             DO  m = 1, surf%ns
     
    15291532          k   = surf%k(m)
    15301533
    1531           surf%pt1(m) = pt(k,j,i) + l_d_cp * pt_d_t(k) * ql(k,j,i)
    1532           surf%qv1(m) = q(k,j,i) - ql(k,j,i)
     1534          IF ( cloud_physics ) THEN
     1535             surf%pt1(m) = pt(k,j,i) + l_d_cp * pt_d_t(k) * ql(k,j,i)
     1536             surf%qv1(m) = q(k,j,i) - ql(k,j,i)
     1537          ELSEIF( cloud_droplets ) THEN
     1538             surf%pt1(m) = pt(k,j,i) + l_d_cp * pt_d_t(k) * ql(k,j,i)
     1539             surf%qv1(m) = q(k,j,i)
     1540          ENDIF
    15331541
    15341542       ENDDO
     
    15841592          ENDIF
    15851593
    1586           IF ( cloud_physics )  THEN
     1594          IF ( cloud_physics  .OR.  cloud_droplets )  THEN
    15871595             !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo )
    15881596             DO  m = 1, surf%ns   
     
    16851693             ENDIF
    16861694
    1687              IF ( cloud_physics )  THEN
     1695             IF ( cloud_physics  .OR.  cloud_droplets )  THEN
    16881696               !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo )
    16891697                DO  m = 1, surf%ns   
Note: See TracChangeset for help on using the changeset viewer.