Changeset 2801 for palm/trunk/SOURCE/lpm.f90

Timestamp:

Feb 14, 2018 4:01:55 PM (6 years ago)

Author:

thiele

Message:

Introduce particle transfer in nested models

File:

• palm/trunk/SOURCE/lpm.f90 (modified) (8 diffs)

palm/trunk/SOURCE/lpm.f90

@@ +1 @@
+MODULE lpm_mod
+
 !> @file lpm.f90
 !------------------------------------------------------------------------------!
@@ -25 +27 @@
 ! -----------------
 ! $Id$
+! Changed lpm from subroutine to module.
+! Introduce particle transfer in nested models.
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -139 +145 @@
 !> Particle advection
 !------------------------------------------------------------------------------!
- SUBROUTINE lpm
  
 
@@ -163 +168 @@
         ONLY: lpm_create_particle, PHASE_RELEASE
 
-    USE lpm_pack_and_sort_mod,                                                   &
-        ONLY:  lpm_sort_in_subboxes, lpm_sort_timeloop_done
+    USE lpm_pack_and_sort_mod
 
     USE particle_attributes,                                                   &
@@ -179 +183 @@
     USE pegrid
 
+ USE pmc_particle_interface,                                                &
+     ONLY: pmcp_c_get_particle_from_parent, pmcp_p_fill_particle_win,       &
+           pmcp_c_send_particle_to_parent, pmcp_p_empty_particle_win,       &
+           pmcp_p_delete_particles_in_fine_grid_area
+
+    USE pmc_interface,                                                      &
+       ONLY: nested_run
+
+ IMPLICIT NONE
+ PRIVATE
+ SAVE
+
+ INTERFACE lpm
+    MODULE PROCEDURE lpm
+ END INTERFACE lpm
+
+ PUBLIC lpm
+
+CONTAINS
+ SUBROUTINE lpm
     IMPLICIT NONE
 
@@ -424 +448 @@
 !--    Horizontal boundary conditions including exchange between subdmains
        CALL lpm_exchange_horiz
-!
-!--    Pack particles (eliminate those marked for deletion),
-!--    determine new number of particles
-       CALL lpm_sort_in_subboxes
-!
-!--    Initialize variables for the next (sub-) timestep, i.e., for marking 
-!--    those particles to be deleted after the timestep
-       deleted_particles = 0
+
+       IF ( .NOT. dt_3d_reached .OR. .NOT. nested_run )   THEN    ! IF .FALSE., lpm_sort_in_subboxes is done inside pcmp
+!
+!--       Pack particles (eliminate those marked for deletion),
+!--       determine new number of particles
+          CALL lpm_sort_in_subboxes
+!
+!--       Initialize variables for the next (sub-) timestep, i.e., for marking
+!--       those particles to be deleted after the timestep
+          deleted_particles = 0
+       ENDIF
 
        IF ( dt_3d_reached )  EXIT
@@ -437 +464 @@
        first_loop_stride = .FALSE.
     ENDDO   ! timestep loop
+!    
+!-- in case of nested runs do the transfer of particles after every full model time step
+    IF ( nested_run )   THEN
+       CALL particles_from_parent_to_child
+       CALL particles_from_child_to_parent
+       CALL pmcp_p_delete_particles_in_fine_grid_area
+
+       CALL lpm_sort_in_subboxes
+
+       deleted_particles = 0
+    ENDIF
 
 !
@@ -469 +507 @@
 
  END SUBROUTINE lpm
+
+ SUBROUTINE particles_from_parent_to_child
+    IMPLICIT NONE
+
+    CALL pmcp_c_get_particle_from_parent                         ! Child actions
+    CALL pmcp_p_fill_particle_win                                ! Parent actions
+
+    RETURN
+ END SUBROUTINE particles_from_parent_to_child
+
+ SUBROUTINE particles_from_child_to_parent
+    IMPLICIT NONE
+
+    CALL pmcp_c_send_particle_to_parent                         ! Child actions
+    CALL pmcp_p_empty_particle_win                              ! Parent actions
+
+    RETURN
+ END SUBROUTINE particles_from_child_to_parent
+
+
+END MODULE lpm_mod