Ignore:
Timestamp:
Nov 8, 2018 2:28:46 PM (6 years ago)
Author:
schwenkel
Message:

Minor bugfixes for bcm cache version and initializing prr

File:
1 edited

Legend:

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

    r3452 r3507  
    2525! -----------------
    2626! $Id$
     27! Minor bugfixes for bcm cache version and initializing prr
     28!
     29! 3452 2018-10-30 13:13:34Z schwenkel
    2730! Bugfix for profiles output
    2831!
     
    958961!--       Liquid water content and precipitation amount
    959962!--       are zero at beginning of the simulation
    960           IF ( bulk_cloud_model )  THEN
    961              ql = 0.0_wp
     963          ql = 0.0_wp
    962964! TODO ???
    963              qc = 0.0_wp
    964              precipitation_amount = 0.0_wp
    965           ENDIF
     965          qc = 0.0_wp
     966          precipitation_amount = 0.0_wp
     967          prr = 0.0_wp
    966968
    967969!
     
    34693471!--       Predetermine flag to mask topography
    34703472          flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    3471           nc_auto = MERGE ( nc(k,j,i), nc_const, microphysics_morrison )
     3473          IF ( microphysics_morrison ) THEN
     3474             nc_auto = nc(k,j,i)
     3475          ELSE
     3476             nc_auto = nc_const
     3477          ENDIF
    34723478
    34733479          IF ( qc(k,j,i) > eps_sb  .AND.  nc_auto > eps_mr )  THEN
     
    36193625!--       Predetermine flag to mask topography
    36203626          flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    3621           nc_accr = MERGE ( nc(k,j,i), nc_const, microphysics_morrison )
     3627          IF ( microphysics_morrison ) THEN
     3628             nc_accr = nc(k,j,i)
     3629          ELSE
     3630             nc_accr = nc_const
     3631          ENDIF
    36223632
    36233633          IF ( ( qc(k,j,i) > eps_sb )  .AND.  ( qr(k,j,i) > eps_sb )  .AND.  &
     
    38593869!--       Predetermine flag to mask topography
    38603870          flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    3861           nc_sedi = MERGE( nc(k,j,i), nc_const, microphysics_morrison )
     3871          IF ( microphysics_morrison ) THEN
     3872             nc_sedi = nc(k,j,i)
     3873          ELSE
     3874             nc_sedi = nc_const
     3875          ENDIF
    38623876!
    38633877!--       Sedimentation fluxes for number concentration are only calculated
Note: See TracChangeset for help on using the changeset viewer.