Ignore:
Timestamp:
Oct 28, 2019 4:03:29 PM (23 months ago)
Author:
schwenkel
Message:

modularize lpm code components of time integration

File:
1 edited

Legend:

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

    r4275 r4276  
    2525! -----------------
    2626! $Id$
     27! Further modularization of lpm code components
     28!
     29! 4275 2019-10-28 15:34:55Z schwenkel
    2730! Move call oft lpm to the end of intermediate timestep loop
    28 ! 
     31!
    2932! 4268 2019-10-17 11:29:38Z schwenkel
    3033! Removing module specific boundary conditions an put them into their modules
    31 ! 
     34!
    3235! 4227 2019-09-10 18:04:34Z gronemeier
    3336! implement new palm_date_time_mod
     
    177180    USE arrays_3d,                                                                                 &
    178181        ONLY:  diss, diss_p, dzu, e, e_p, nc, nc_p, nr, nr_p, prho, pt, pt_p, pt_init, q_init, q,  &
    179                qc, qc_p, ql, ql_c, ql_v, ql_vp, qr, qr_p, q_p, ref_state, rho_ocean, s, s_p, sa_p, &
     182               qc, qc_p, qr, qr_p, q_p, ref_state, rho_ocean, s, s_p, sa_p, &
    180183               tend, u, u_p, v, vpt, v_p, w, w_p
    181184
     
    207210               averaging_interval, averaging_interval_pr, bc_lr_cyc, bc_ns_cyc, bc_pt_t_val,       &
    208211               bc_q_t_val, biometeorology, call_psolver_at_all_substeps,  child_domain,            &
    209                cloud_droplets, constant_flux_layer, constant_heatflux, create_disturbances,        &
     212               constant_flux_layer, constant_heatflux, create_disturbances,        &
    210213               dopr_n, constant_diffusion, coupling_mode, coupling_start_time,                     &
    211214               current_timestep_number, disturbance_created, disturbance_energy_limit, dist_range, &
     
    254257
    255258    USE lagrangian_particle_model_mod,                                                             &
    256         ONLY:  lpm_data_output_ptseries, lpm_interaction_droplets_ptq
     259        ONLY:  lpm_data_output_ptseries
    257260
    258261    USE lsf_nudging_mod,                                                                           &
     
    588591       IF ( nesting_offline )  CALL nesting_offl_input
    589592!
    590 !--    Execute all other module actions routunes
     593!--    Execute all other module actions routines
    591594       CALL module_interface_actions( 'before_timestep' )
    592595
     
    975978!--       ### particle model should be moved before prognostic_equations, in order
    976979!--       to regard droplet interactions directly
    977           IF ( particle_advection  .AND.  time_since_reference_point >= particle_advection_start   &
    978                .AND.  intermediate_timestep_count == intermediate_timestep_count_max )             &
    979           THEN
    980              CALL module_interface_actions( 'after_pressure_solver' )
    981              first_call_lpm = .FALSE.
    982           ENDIF
    983 
    984           IF ( cloud_droplets )  THEN
    985              CALL exchange_horiz( ql, nbgp )
    986              CALL exchange_horiz( ql_c, nbgp )
    987              CALL exchange_horiz( ql_v, nbgp )
    988              CALL exchange_horiz( ql_vp, nbgp )
    989           ENDIF
     980
     981          CALL module_interface_actions( 'after_pressure_solver' )
    990982!
    991983!--       Interaction of droplets with temperature and mixing ratio.
    992984!--       Droplet condensation and evaporation is calculated within
    993985!--       advec_particles.
    994           IF ( cloud_droplets  .AND.  intermediate_timestep_count == intermediate_timestep_count_max ) &
    995           THEN
    996              CALL lpm_interaction_droplets_ptq
    997              CALL exchange_horiz( pt, nbgp )
    998              CALL exchange_horiz( q, nbgp )
    999           ENDIF
    1000986!
    1001987!--       If required, compute liquid water content
Note: See TracChangeset for help on using the changeset viewer.