source: palm/trunk/SOURCE/pmc_particle_interface.f90 @ 3129

Last change on this file since 3129 was 3123, checked in by suehring, 6 years ago

Correct working precision for REAL and INTEGER numbers

  • Property svn:keywords set to Id
File size: 44.1 KB
Line 
1MODULE pmc_particle_interface
2
3!------------------------------------------------------------------------------!
4! This file is part of PALM.
5!
6! PALM is free software: you can redistribute it and/or modify it under the
7! terms of the GNU General Public License as published by the Free Software
8! Foundation, either version 3 of the License, or (at your option) any later
9! version.
10!
11! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
12! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
13! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
14!
15! You should have received a copy of the GNU General Public License along with
16! PALM. If not, see <http://www.gnu.org/licenses/>.
17!
18! Copyright 1997-2018 Leibniz Universitaet Hannover
19!------------------------------------------------------------------------------!
20!
21! Current revisions:
22! ------------------
23!
24!
25! Former revisions:
26! -----------------!
27! $Id: pmc_particle_interface.f90 3123 2018-07-12 16:21:53Z gronemeier $
28! Correct working precision for REAL numbers
29!
30! 3049 2018-05-29 13:52:36Z Giersch
31! Error messages revised
32!
33! 2967 2018-04-13 11:22:08Z Giersch
34! bugfix: missing parallel cpp-directives added
35!
36! 2884 2018-03-14 08:33:20Z scharf
37! B: corrected KIND of variable "parsize" for some MPI calls
38!
39! 2841 2018-02-27 15:02:57Z knoop
40! Bugfix: wrong placement of include 'mpif.h' corrected
41!
42! 2809 2018-02-15 09:55:58Z schwenkel
43! Bugfix for gfortran: Replace the function C_SIZEOF with STORAGE_SIZE.
44
45!
46! Initial Version (by Klaus Ketelsen)
47!
48!
49! Description:
50! ------------
51! Introduce particle transfer in nested models.
52! Interface to palm lpm model to handel particle transfer between parent and
53! child model.
54!------------------------------------------------------------------------------!
55
56   USE, INTRINSIC ::  ISO_C_BINDING
57
58#if defined( __parallel )  &&  !defined( __mpifh )
59   USE MPI
60#endif
61
62   USE kinds
63
64   USE pegrid,                                                                 &
65       ONLY: myidx,myidy
66
67   USE indices,                                                                &
68       ONLY: nx, ny, nxl, nxr, nys, nyn, nzb, nzt, nxlg, nxrg, nysg, nyng, nbgp
69
70   USE grid_variables,                                                         &
71       ONLY:  dx, dy
72
73   USE arrays_3d,                                                              &
74        ONLY: zw
75
76   USE control_parameters,                                                     &
77       ONLY:  message_string
78
79   USE particle_attributes,                                                    &
80       ONLY:  prt_count, particles, grid_particles,                            &
81              particle_type, number_of_particles, zero_particle,               &
82              ibc_par_t, ibc_par_lr, ibc_par_ns, alloc_factor
83
84   USE lpm_pack_and_sort_mod
85
86   USE lpm_exchange_horiz_mod,                                                 &
87       ONLY: realloc_particles_array
88
89#if defined( __parallel )
90   USE pmc_general,                                                            &
91       ONLY: pedef
92
93   USE pmc_parent,                                                             &
94       ONLY: children, pmc_s_fillbuffer, pmc_s_getdata_from_buffer,            &
95             pmc_s_get_child_npes
96
97   USE pmc_child,                                                              &
98       ONLY:  me, pmc_c_getbuffer, pmc_c_putbuffer
99
100   USE pmc_interface,                                                          &
101       ONLY: cpl_id, get_number_of_childs, nr_part,  part_adr, nested_run,     &
102             get_childid, get_child_edges, nr_partc,  part_adrc,               &
103             coarse_bound, coord_x, coord_y, cg,   get_child_gridspacing,      &
104             lower_left_coord_x, lower_left_coord_y
105
106    USE pmc_handle_communicator,                                               &
107        ONLY:  pmc_parent_for_child
108
109   USE pmc_mpi_wrapper,                                                        &
110       ONLY: pmc_send_to_parent, pmc_recv_from_child
111
112#endif
113
114   IMPLICIT NONE
115
116#if defined( __parallel )  &&  defined( __mpifh )
117   INCLUDE "mpif.h"
118#endif
119
120   PRIVATE
121   SAVE
122
123   TYPE  coarse_particle_def
124      INTEGER(iwp) ::  nr_particle
125     
126      TYPE(particle_type),ALLOCATABLE,DIMENSION(:) ::  parent_particles
127   END TYPE  coarse_particle_def
128
129   INTEGER(iwp),PARAMETER ::  Min_particles_per_column = 100      !<
130   INTEGER(iwp),PARAMETER ::  max_nr_particle_in_rma_win = 100000 !<
131
132   INTEGER(iwp) ::  nr_fine_in_coarse  !< Number of fine grid cells in coarse grid (one direction)
133   INTEGER(iwp) ::  particle_win_child !<
134
135   INTEGER(iwp),ALLOCATABLE,DIMENSION(:) ::  particle_win_parent !<
136   
137   TYPE(C_PTR), ALLOCATABLE,DIMENSION(:) ::  buf_ptr !<
138   
139   TYPE(particle_type), DIMENSION(:),POINTER ::  particle_in_win !<
140
141   TYPE(coarse_particle_def),ALLOCATABLE,DIMENSION(:,:) ::  coarse_particles !<
142
143   
144
145   INTERFACE pmcp_g_init
146      MODULE PROCEDURE pmcp_g_init
147   END  INTERFACE pmcp_g_init
148
149   INTERFACE pmcp_g_alloc_win
150      MODULE PROCEDURE pmcp_g_alloc_win
151   END  INTERFACE pmcp_g_alloc_win
152
153   INTERFACE pmcp_c_get_particle_from_parent
154      MODULE PROCEDURE pmcp_c_get_particle_from_parent
155   END  INTERFACE pmcp_c_get_particle_from_parent
156
157   INTERFACE pmcp_c_send_particle_to_parent
158      MODULE PROCEDURE pmcp_c_send_particle_to_parent
159   END  INTERFACE pmcp_c_send_particle_to_parent
160
161   INTERFACE pmcp_p_fill_particle_win
162      MODULE PROCEDURE pmcp_p_fill_particle_win
163   END  INTERFACE pmcp_p_fill_particle_win
164
165   INTERFACE pmcp_p_empty_particle_win
166      MODULE PROCEDURE pmcp_p_empty_particle_win
167   END  INTERFACE pmcp_p_empty_particle_win
168
169   INTERFACE pmcp_g_print_number_of_particles
170      MODULE PROCEDURE pmcp_g_print_number_of_particles
171   END  INTERFACE pmcp_g_print_number_of_particles
172
173   INTERFACE pmcp_p_delete_particles_in_fine_grid_area
174      MODULE PROCEDURE pmcp_p_delete_particles_in_fine_grid_area
175   END  INTERFACE pmcp_p_delete_particles_in_fine_grid_area
176
177   PUBLIC pmcp_g_init, pmcp_g_alloc_win, pmcp_c_get_particle_from_parent
178   PUBLIC pmcp_c_send_particle_to_parent, pmcp_p_fill_particle_win, pmcp_g_print_number_of_particles
179   PUBLIC pmcp_p_empty_particle_win, pmcp_p_delete_particles_in_fine_grid_area
180
181 CONTAINS
182
183!------------------------------------------------------------------------------!
184! Description:
185! ------------
186!> general routine:
187!> Initializing actions of the particle interface
188!> check particle boundary conditions for the child models
189!------------------------------------------------------------------------------!
190 SUBROUTINE pmcp_g_init
191 
192    IMPLICIT NONE
193 
194    INTEGER(iwp) ::  nr_childs !< Number of child models of the current model
195
196#if defined( __parallel )
197
198    nr_childs = get_number_of_childs()
199!
200!-- Check if the current model has child models
201    IF ( nr_childs > 0 )   THEN
202       ALLOCATE( nr_part(nysg:nyng, nxlg:nxrg) )
203       ALLOCATE( part_adr(nysg:nyng, nxlg:nxrg) )
204       nr_part = 0
205       part_adr   = 0
206    ENDIF
207
208!
209!-- Set the boundary conditions to nested for all non root (i.e child) models
210    IF ( cpl_id > 1 )  THEN
211   
212       IF ( ibc_par_t /= 3 )   THEN
213          ibc_par_t  = 3
214          message_string = 'In Child model:  ibc_par_t is automatically set to nested '
215          CALL message( 'pmcp_g_init ', 'PA0477', 0, 1, 0, 6, 0 )
216       ENDIF
217   
218       IF ( ibc_par_lr /= 3 )   THEN
219          ibc_par_lr = 3
220          message_string = 'In Child model:  ibc_par_lr is automatically set to nested '
221          CALL message( 'pmcp_g_init ', 'PA0478', 0, 1, 0, 6, 0 )
222       ENDIF
223       
224       IF ( ibc_par_ns /= 3 )   THEN
225          ibc_par_ns = 3
226          message_string = 'In Child model:  ibc_par_ns is automatically set to nested '
227          CALL message( 'pmcp_g_init ', 'PA0479', 0, 1, 0, 6, 0 )
228       ENDIF
229       
230    ENDIF
231
232#endif
233 END SUBROUTINE pmcp_g_init
234!------------------------------------------------------------------------------!
235! Description:
236! ------------
237!> general routine:
238!> allocate the MPI windows
239!------------------------------------------------------------------------------!
240 SUBROUTINE pmcp_g_alloc_win
241 
242    IMPLICIT NONE
243   
244    INTEGER(iwp) ::  m         !< loop index
245    INTEGER(iwp) ::  ierr      !< error code
246    INTEGER(iwp) ::  icl       !< left boundary in coarse(parent) index space
247    INTEGER(iwp) ::  icr       !< right boundary in coarse(parent) index space
248    INTEGER(iwp) ::  jcs       !< south boundary in coarse(parent) index space
249    INTEGER(iwp) ::  jcn       !< north boundary in coarse(parent) index space
250    INTEGER(iwp) ::  child_id  !< Id of a child model
251    INTEGER(iwp) ::  nr_childs !< Number of child models of the current model
252   
253    INTEGER ::  parsize !<
254    TYPE(C_PTR), SAVE ::  ptr !<
255   
256    TYPE(particle_type),DIMENSION(:),POINTER ::  win_buffer !<
257   
258    INTEGER(iwp),DIMENSION(1) ::  buf_shape !<
259
260#if defined( __parallel )
261    INTEGER(KIND=MPI_ADDRESS_KIND) ::  parsize_mpi_address_kind !<
262    INTEGER(KIND=MPI_ADDRESS_KIND) ::  winsize !<
263
264!
265!-- If the model has a parent model prepare the structures for transfer
266    IF ( cpl_id > 1 )  THEN
267
268       parsize_mpi_address_kind = STORAGE_SIZE(zero_particle)/8
269
270       CALL MPI_ALLOC_MEM( parsize_mpi_address_kind , MPI_INFO_NULL, ptr, ierr )
271       parsize = parsize_mpi_address_kind
272       buf_shape(1) = 1
273       CALL C_F_POINTER( ptr, win_buffer, buf_shape )
274       CALL MPI_WIN_CREATE( win_buffer, parsize_mpi_address_kind, parsize,     &
275                            MPI_INFO_NULL, me%intra_comm, particle_win_child,  &
276                            ierr )
277
278!
279!--    Child domain boundaries in the parent index space
280       icl = coarse_bound(1)
281       icr = coarse_bound(2)
282       jcs = coarse_bound(3)
283       jcn = coarse_bound(4)
284
285       ALLOCATE( coarse_particles(jcs:jcn, icl:icr) )
286
287       coarse_particles(:,:)%nr_particle = 0
288    ENDIF
289
290!
291!-- If the model has child models prepare the structures for transfer
292    nr_childs = get_number_of_childs()
293    IF ( nr_childs > 0 )   THEN
294       ALLOCATE( particle_win_parent(nr_childs) )
295       ALLOCATE( buf_ptr(nr_childs) )
296       DO  m = 1, nr_childs
297          child_id = get_childid(m)
298          parsize_mpi_address_kind = STORAGE_SIZE(zero_particle)/8
299          parsize = parsize_mpi_address_kind
300
301          winsize = max_nr_particle_in_rma_win * parsize_mpi_address_kind
302          CALL MPI_ALLOC_MEM( winsize , MPI_INFO_NULL, buf_ptr(m), ierr )
303          buf_shape(1) = max_nr_particle_in_rma_win
304          CALL C_F_POINTER( buf_ptr(m), win_buffer, buf_shape )
305          CALL MPI_WIN_CREATE( win_buffer, winsize, parsize, MPI_INFO_NULL,    &
306                               children(child_id)%intra_comm,                  &
307                               particle_win_parent(m), ierr )
308          ENDDO
309    ENDIF
310
311#endif
312 END SUBROUTINE pmcp_g_alloc_win
313
314
315!------------------------------------------------------------------------------!
316! Description:
317! ------------
318!> child routine:
319!> Read/get particles out of the parent MPI window
320!------------------------------------------------------------------------------!
321 SUBROUTINE pmcp_c_get_particle_from_parent
322 
323    IMPLICIT NONE
324
325    INTEGER(iwp) ::  i    !< x grid index
326    INTEGER(iwp) ::  icl  !< left boundary in coarse(parent) index space
327    INTEGER(iwp) ::  ierr !< error code
328    INTEGER(iwp) ::  ij   !< combined xy index for the buffer array
329    INTEGER(iwp) ::  ip   !< loop index (child PEs)
330    INTEGER(iwp) ::  j    !< y grid index
331    INTEGER(iwp) ::  jcs  !< south boundary in coarse(parent) index space
332    INTEGER(iwp) ::  nr   !< number of particles to receive from a parent box
333   
334    INTEGER ::  parsize !<
335
336#if defined( __parallel )
337    TYPE(pedef), POINTER ::  ape !< TO_DO Klaus: give a description and better name of the variable
338
339    INTEGER(KIND=MPI_ADDRESS_KIND) ::  parsize_mpi_address_kind !<
340    INTEGER(KIND=MPI_ADDRESS_KIND) ::  target_disp !<
341
342    IF ( cpl_id > 1 )  THEN
343   
344       CALL pmc_c_getbuffer( particle_transfer = .TRUE. ) !Get number of particle/column and offset in RMA window xx
345
346!
347!--    Wait for buffer to fill.
348!
349!--    The parent side (in pmc_s_fillbuffer) is filling the buffer in the MPI RMA window
350!--    When the filling is complete, a MPI_BARRIER is called.
351!--    The child is not allowd to access the parent-buffer before it is completely filled
352!--    Synchronization is done implicitely in pmc_c_getbuffer and pmc_s_fillbuffer on the parent side
353
354       icl = coarse_bound(1)
355       jcs = coarse_bound(3)
356
357       DO  ip = 1, me%inter_npes
358
359          ape => me%pes(ip)
360
361          DO  ij = 1, ape%nrele
362              j = ape%locind(ij)%j + jcs - 1
363              i = ape%locind(ij)%i + icl - 1
364              nr = nr_partc(j,i)
365              IF ( nr > 0 )  THEN
366
367                 CALL check_and_alloc_coarse_particle (i, j, nr)
368                 parsize_mpi_address_kind = STORAGE_SIZE(zero_particle)/8
369                 parsize = parsize_mpi_address_kind
370                 target_disp = part_adrc(j,i) - 1
371                 CALL MPI_WIN_LOCK( MPI_LOCK_SHARED , ip - 1, 0,               &
372                                    particle_win_child, ierr )
373                 CALL MPI_GET( coarse_particles(j,i)%parent_particles,         &
374                               nr * parsize, MPI_BYTE, ip - 1, target_disp,    &
375                               nr * parsize, MPI_BYTE, particle_win_child, ierr )
376                 CALL MPI_WIN_UNLOCK( ip - 1, particle_win_child, ierr )
377             ENDIF
378             coarse_particles(j,i)%nr_particle = nr
379          ENDDO
380       ENDDO
381
382       CALL c_copy_particle_to_child_grid
383    ENDIF
384
385#endif
386 END SUBROUTINE pmcp_c_get_particle_from_parent
387
388
389!------------------------------------------------------------------------------!
390! Description:
391! ------------
392!> child routine:
393!> Write/put particles into the parent MPI window
394!------------------------------------------------------------------------------!
395  SUBROUTINE pmcp_c_send_particle_to_parent
396 
397    IMPLICIT NONE
398   
399    INTEGER(iwp) ::  disp_offset            !<
400    INTEGER(iwp) ::  i                      !< x loop index
401    INTEGER(iwp) ::  icl                    !< left boundary in coarse(parent) index space
402    INTEGER(iwp) ::  icr                    !< right boundary in coarse(parent) index space
403    INTEGER(iwp) ::  ierr                   !< error code
404    INTEGER(iwp) ::  ij                     !< combined xy index for the buffer array
405    INTEGER(iwp) ::  ip                     !< loop index (child PEs)
406    INTEGER(iwp) ::  j                      !< y loop index
407    INTEGER(iwp) ::  jcs                    !< south boundary in coarse(parent) index space
408    INTEGER(iwp) ::  jcn                    !< north boundary in coarse(parent) index space
409    INTEGER(iwp) ::  max_nr_particle_per_pe !< maximum number of particles per PE (depending on grid apect ratio)
410    INTEGER(iwp) ::  n                      !< shorter variable name for nr_fine_in_coarse
411    INTEGER(iwp) ::  nr                     !< shorter variabel name for nr_partc
412    INTEGER(iwp) ::  pe_offset              !< offset index of the current PE
413   
414    INTEGER ::  parsize !<
415   
416    REAL(wp) ::  eps=0.00001 !< used in calculations to avoid rounding errors
417    REAL(wp) ::  xx          !< number of fine grid cells inside a coarse grid cell in x-direction
418    REAL(wp) ::  yy          !< number of fine grid cells inside a coarse grid cell in y-direction
419
420 !   TYPE(particle_type) ::  dummy_part !< dummy particle (needed for size calculations)
421
422#if defined( __parallel )
423    TYPE(pedef), POINTER ::  ape !< TO_DO Klaus: give a description and better name of the variable
424
425    INTEGER(KIND=MPI_ADDRESS_KIND) ::  parsize_mpi_address_kind !<
426    INTEGER(KIND=MPI_ADDRESS_KIND) ::  target_disp !<
427   
428   
429    IF ( cpl_id > 1 )  THEN
430       CALL c_copy_particle_to_coarse_grid
431
432!
433!--    Child domain boundaries in the parent index space
434
435       icl = coarse_bound(1)
436       icr = coarse_bound(2)
437       jcs = coarse_bound(3)
438       jcn = coarse_bound(4)
439
440       nr_partc = 0
441       
442       DO i = icl, icr
443          DO j = jcs, jcn
444             nr_partc(j,i) = coarse_particles(j,i)%nr_particle
445          ENDDO
446       ENDDO
447       part_adrc = 0
448
449!
450!--    compute number of fine grid cells in coarse grid (one direction)
451       xx = (cg%dx + eps) / dx ! +eps to avoid rounding error
452       yy = (cg%dy + eps) / dy
453       nr_fine_in_coarse = MAX( INT(xx), INT(yy) )
454
455       IF ( MOD( coord_x(0), cg%dx ) /= 0.0 .OR. MOD( coord_y(0), cg%dy ) /= 0.0 ) THEN
456          nr_fine_in_coarse = nr_fine_in_coarse + 1
457       ENDIF
458
459!
460!--    Assign a number to my child PE to select different areas in the RMA window on server side
461!--    With this number a square of child PEs is defined which share the same coarse grid cells
462
463       n           = nr_fine_in_coarse ! local variable n to make folloing statements shorter
464       pe_offset   = MOD( myidx, n ) * n + MOD( myidy, n)
465       max_nr_particle_per_pe = max_nr_particle_in_rma_win / ( n * n )
466       disp_offset            = pe_offset * max_nr_particle_per_pe
467       parsize_mpi_address_kind = STORAGE_SIZE(zero_particle)/8
468       parsize = parsize_mpi_address_kind
469       DO  ip = 1, me%inter_npes
470
471          ape => me%pes(ip)
472
473          target_disp = disp_offset
474          DO ij = 1, ape%nrele
475             j  = ape%locind(ij)%j + jcs - 1
476             i  = ape%locind(ij)%i + icl - 1
477             nr = nr_partc(j,i)
478             IF( nr > 0 ) THEN
479                IF ( target_disp + nr - disp_offset >= max_nr_particle_per_pe ) THEN
480                   WRITE(9,*) 'RMA window too small on child ',                &
481                              target_disp + nr - disp_offset,                  &
482                              max_nr_particle_per_pe, max_nr_particle_in_rma_win
483                   message_string = 'RMA window too small on child'
484                   CALL message( 'pmci_create_child_arrays', 'PA0480', 3, 2, 0, 6, 0 )
485                ENDIF
486                CALL MPI_WIN_LOCK( MPI_LOCK_SHARED , ip - 1, 0,                &
487                                   particle_win_child, ierr )
488                CALL MPI_PUT( coarse_particles(j,i)%parent_particles,          &
489                              nr * parsize, MPI_BYTE, ip - 1, target_disp,     &
490                              nr * parsize, MPI_BYTE,  particle_win_child, ierr )
491                CALL MPI_WIN_UNLOCK( ip - 1, particle_win_child, ierr )
492                part_adrc(j,i) = target_disp + 1
493                target_disp    = target_disp + nr
494             ENDIF
495          ENDDO
496       ENDDO
497
498       CALL pmc_c_putbuffer ( particle_transfer = .TRUE. )   !Send new number of particle/column and offset to parent
499
500    ENDIF
501
502#endif
503 END SUBROUTINE pmcp_c_send_particle_to_parent
504
505
506!------------------------------------------------------------------------------!
507! Description:
508! ------------
509!> parent routine:
510!> write particles into the MPI window
511!------------------------------------------------------------------------------!
512 SUBROUTINE pmcp_p_fill_particle_win
513 
514    IMPLICIT NONE
515
516    LOGICAL      ::  active_particle !< Particles located in the fine/child grid area are marked as active (to be transferred)
517    LOGICAL,SAVE ::  lfirst = .TRUE. !<
518   
519    INTEGER(iwp) ::  child_id            !< id of the child model
520    INTEGER(iwp) ::  i                   !< x grid box index
521    INTEGER(iwp) ::  ierr                !< error code
522    INTEGER(iwp) ::  ij                  !< combined xy index for the buffer array
523    INTEGER(iwp) ::  ip                  !< loop index (child PEs)
524    INTEGER(iwp) ::  j                   !< y grid box index
525    INTEGER(iwp) ::  k                   !< z grid box index
526    INTEGER(iwp) ::  m                   !< loop index (number of childs)
527    INTEGER(iwp) ::  n                   !< loop index (number of particles)
528    INTEGER(iwp) ::  nr_part_col         !< Number of particles to transfer per column
529    INTEGER(iwp) ::  number_of_particles !<
530    INTEGER(iwp) ::  pindex              !<
531    INTEGER(iwp) ::  tot_particle_count  !< Total number of particles per child
532   
533    REAL(wp) ::  dx_child   !< child grid spacing
534    REAL(wp) ::  dy_child   !< child grid spacing
535    REAL(wp) ::  dz_child   !< child grid spacing
536    REAL(wp) ::  ny_coord   !< north coordinate of child grid
537    REAL(wp) ::  ny_coord_b !< north coordinate of child grid boundary
538    REAL(wp) ::  lx_coord   !< left coordinate of child grid
539    REAL(wp) ::  lx_coord_b !< left coordinate of child grid boundary
540    REAL(wp) ::  rx_coord   !< right coordinate of child grid
541    REAL(wp) ::  rx_coord_b !< right coordinate of child grid boundary
542    REAL(wp) ::  sy_coord   !< south coordinate of child grid
543    REAL(wp) ::  sy_coord_b !< south coordinate of child grid boundary
544    REAL(wp) ::  uz_coord   !< upper coordinate of child grid
545    REAL(wp) ::  uz_coord_b !< upper coordinate of child grid boundary
546    REAL(wp) ::  x          !< particle position
547    REAL(wp) ::  xo         !< origin of particle
548    REAL(wp) ::  y          !< particle position
549    REAL(wp) ::  yo         !< origin of particle
550    REAL(wp) ::  z          !< particle position
551   
552    INTEGER(iwp),DIMENSION(1) ::  buf_shape !<
553   
554#if defined( __parallel )
555    TYPE(pedef), POINTER ::  ape !< TO_DO Klaus: give a description and better name of the variable
556
557    DO  m = 1, get_number_of_childs()
558
559       child_id = get_childid(m)
560
561       CALL get_child_edges( m, lx_coord, lx_coord_b, rx_coord, rx_coord_b,    &
562                                sy_coord, sy_coord_b, ny_coord, ny_coord_b,    &
563                                uz_coord, uz_coord_b) 
564
565       CALL get_child_gridspacing( m, dx_child, dy_child, dz_child )
566
567       IF ( lfirst )   THEN
568          WRITE(9,'(a,5f10.2)') 'edges          ',lx_coord,rx_coord,sy_coord,ny_coord,uz_coord
569          WRITE(9,'(a,5f10.2)') 'edges boundary ',lx_coord_b,rx_coord_b,sy_coord_b,ny_coord_b,uz_coord_b
570          WRITE(9,'(a,5f10.2)') 'child spacing  ',dx_child, dy_child, dz_child,lower_left_coord_x,lower_left_coord_y
571       ENDIF
572!
573!--    reset values for every child
574       tot_particle_count = 0 
575       nr_part  = 0           
576       part_adr = 0
577       pindex   = 1
578
579       buf_shape(1) = max_nr_particle_in_rma_win
580       CALL C_F_POINTER( buf_ptr(m), particle_in_win , buf_shape )
581
582       DO  ip = 1, children(child_id)%inter_npes
583
584          ape => children(child_id)%pes(ip)
585
586          nr_part_col   = 0 
587         
588          DO  ij = 1, ape%nrele
589             
590!
591!--          Inside the PMC adressing of 3d arrays starts with 1
592             i = ape%locind(ij)%i + nxl - nbgp - 1
593             j = ape%locind(ij)%j + nys - nbgp - 1
594             nr_part_col = 0 ! Number of particles to transfer per column
595             part_adr(j,i) = pindex
596             
597             DO  k = nzb + 1, nzt
598                number_of_particles = prt_count(k,j,i)
599               
600                IF ( number_of_particles <= 0 )  CYCLE
601               
602                particles => grid_particles(k,j,i)%particles(1:number_of_particles)
603
604                ! Select particles within boundary area
605
606                DO n = 1, number_of_particles
607                   x = particles(n)%x
608                   y = particles(n)%y
609                   z = particles(n)%z
610!
611!--                check if the particle is located in the fine grid area
612                   active_particle = ((x > lx_coord .AND. x < rx_coord) .AND.  & 
613                                      (y > sy_coord .AND. y < ny_coord) .AND.  &
614                                      (z > 0.000000001 .AND. z < uz_coord))
615                   IF ( active_particle .AND. particles(n)%particle_mask )  THEN
616                     
617                      particle_in_win(pindex) = particles(n)
618!
619!--                   Change particle positions and origin relative to global origin
620                      particle_in_win(pindex)%x = particle_in_win(pindex)%x + lower_left_coord_x
621                      particle_in_win(pindex)%y = particle_in_win(pindex)%y + lower_left_coord_y
622                      particle_in_win(pindex)%origin_x = particle_in_win(pindex)%origin_x + lower_left_coord_x
623                      particle_in_win(pindex)%origin_y = particle_in_win(pindex)%origin_y + lower_left_coord_y
624
625                      tot_particle_count = tot_particle_count + 1
626                      nr_part_col        = nr_part_col + 1
627                      pindex             = pindex + 1
628                      IF ( pindex > max_nr_particle_in_rma_win ) THEN
629                         WRITE(9,*) 'RMA window too small on parent ',pindex, max_nr_particle_in_rma_win
630                         message_string = 'RMA window too small on parent'
631                         CALL message( 'pmci_create_child_arrays', 'PA0481', 3, 2, 0, 6, 0 )   ! PA number has to be adjusted
632                     ENDIF
633                   END IF
634                ENDDO
635             ENDDO
636             nr_part(j,i) = nr_part_col
637          ENDDO
638       ENDDO
639
640       CALL pmc_s_fillbuffer( child_id, particle_transfer = .TRUE. )
641    ENDDO
642
643    lfirst = .FALSE.
644
645#endif
646 END SUBROUTINE pmcp_p_fill_particle_win
647
648 
649!------------------------------------------------------------------------------!
650! Description:
651! ------------
652!> parent routine:
653!> delete particles from the MPI window
654!------------------------------------------------------------------------------!
655 SUBROUTINE pmcp_p_empty_particle_win
656
657    IMPLICIT NONE
658
659    INTEGER(iwp) ::  child_id           !< model id of the child
660    INTEGER(iwp) ::  ip                 !< loop index (child PEs)
661    INTEGER(iwp) ::  m                  !< loop index (number of childs)
662
663    INTEGER(iwp),DIMENSION(1) ::  buf_shape !<
664
665#if defined( __parallel )
666    DO  m = 1, get_number_of_childs()
667
668       child_id = get_childid(m)
669
670       buf_shape(1) = max_nr_particle_in_rma_win
671       CALL C_F_POINTER( buf_ptr(m), particle_in_win , buf_shape )
672
673!
674!--    In some cells of the coarse grid, there are contributions from more than one child process
675!--    Therfore p_copy_particle_to_org_grid is done for one child process per call
676
677       DO ip = 1, pmc_s_get_child_npes( child_id )
678          nr_part  = 0
679          part_adr = 0
680
681          CALL pmc_s_getdata_from_buffer( child_id, particle_transfer = .TRUE., child_process_nr = ip )
682
683          CALL p_copy_particle_to_org_grid( m, child_id )
684       ENDDO
685
686    ENDDO
687
688#endif
689 END SUBROUTINE pmcp_p_empty_particle_win
690
691
692!------------------------------------------------------------------------------!
693! Description:
694! ------------
695!> parent routine:
696!> After the transfer mark all parent particles that are still inside on of the
697!> child areas for deletion.
698!------------------------------------------------------------------------------!
699 SUBROUTINE pmcp_p_delete_particles_in_fine_grid_area
700
701    IMPLICIT NONE
702
703    LOGICAL ::  to_delete !< particles outside of model domain are marked as to_delete
704   
705    INTEGER(iwp) ::  i !< loop index (x grid)
706    INTEGER(iwp) ::  j !< loop index (y grid)
707    INTEGER(iwp) ::  k !< loop index (z grid)
708    INTEGER(iwp) ::  m !< loop index (number of particles)
709    INTEGER(iwp) ::  n !< loop index (number of childs)
710   
711    REAL(wp) ::  dx_child   !< child grid spacing
712    REAL(wp) ::  dy_child   !< child grid spacing
713    REAL(wp) ::  dz_child   !< child grid spacing
714    REAL(wp) ::  ny_coord   !< north coordinate of child grid
715    REAL(wp) ::  ny_coord_b !< north coordinate of child grid boundary
716    REAL(wp) ::  lx_coord   !< left coordinate of child grid
717    REAL(wp) ::  lx_coord_b !< left coordinate of child grid boundary
718    REAL(wp) ::  rx_coord   !< right coordinate of child grid
719    REAL(wp) ::  rx_coord_b !< right coordinate of child grid boundary
720    REAL(wp) ::  sy_coord   !< south coordinate of child grid
721    REAL(wp) ::  sy_coord_b !< south coordinate of child grid boundary
722    REAL(wp) ::  uz_coord   !< upper coordinate of child grid
723    REAL(wp) ::  uz_coord_b !< upper coordinate of child grid boundary
724    REAL(wp) ::  x          !< particle position
725    REAL(wp) ::  y          !< particle position
726    REAL(wp) ::  z          !< particle position
727   
728#if defined( __parallel )
729    DO  m = 1, get_number_of_childs()
730       CALL get_child_edges( m, lx_coord, lx_coord_b, rx_coord, rx_coord_b,    &
731                                sy_coord, sy_coord_b, ny_coord, ny_coord_b,    &
732                                uz_coord, uz_coord_b )
733
734
735       CALL get_child_gridspacing( m, dx_child, dy_child, dz_child ) 
736
737       DO i = nxl, nxr
738          DO j = nys, nyn
739             DO k = nzb, nzt
740                number_of_particles = prt_count(k,j,i)
741
742                IF ( number_of_particles == 0 ) CYCLE
743
744                particles => grid_particles(k,j,i)%particles(1:number_of_particles)
745
746                DO n = 1, number_of_particles
747                   x = particles(n)%x
748                   y = particles(n)%y
749                   z = particles(n)%z
750
751                   to_delete = ((x > lx_coord .AND. x < rx_coord) .AND.        &
752                                (y > sy_coord .AND. y < ny_coord) .AND.        &
753                                (z > 0.000000001 .AND. z < uz_coord))
754
755                   IF ( to_delete ) THEN
756                      particles(n)%particle_mask = .FALSE.
757                   ENDIF
758                ENDDO
759             ENDDO
760          ENDDO
761       ENDDO
762
763    ENDDO
764
765#endif
766 END SUBROUTINE pmcp_p_delete_particles_in_fine_grid_area
767
768
769!------------------------------------------------------------------------------!
770! Description:
771! ------------
772!> general routine:
773!> print the total number of of the current model and its child models into a file
774!------------------------------------------------------------------------------!
775 SUBROUTINE pmcp_g_print_number_of_particles (my_time,local_nr_particles)
776 
777    USE pegrid,                                                                &
778        ONLY: myid
779 
780    IMPLICIT NONE
781
782    INTEGER(iwp),INTENT(IN) ::  local_nr_particles !<
783   
784    REAL(wp),INTENT(IN) ::  my_time !<
785     
786    LOGICAL, SAVE ::  is_file_open=.FALSE. !<
787   
788    INTEGER(iwp) ::  child_id           !<
789    INTEGER(iwp) ::  child_nr_particles !< total number of particles in all child models
790    INTEGER(iwp) ::  ierr               !< error code
791    INTEGER(iwp) ::  m                  !< loop index (number of childs)
792   
793    CHARACTER(LEN=16) ::  fname !< filename
794   
795    INTEGER(iwp),DIMENSION(2) ::  ivalr !< integer value to be received
796    INTEGER(iwp),DIMENSION(2) ::  ivals !< integer value to be send
797   
798#if defined( __parallel )
799    child_nr_particles = 0
800    IF ( myid == 0 )  THEN
801       IF ( cpl_id > 1 )  THEN
802          ivals(1) = local_nr_particles
803          CALL pmc_send_to_parent( ivals, 1, 0, 400, ierr )
804       ENDIF
805       DO  m = 1, SIZE( pmc_parent_for_child ) - 1
806
807          child_id = pmc_parent_for_child(m)
808          CALL pmc_recv_from_child( child_id, ivalr, 1, 0, 400, ierr )
809          child_nr_particles = child_nr_particles + ivalr(1)
810       ENDDO
811
812       IF ( SIZE( pmc_parent_for_child ) > 1 ) THEN
813          IF ( .NOT. is_file_open )  THEN !kk muss noch auf file_open umgestellt werden
814             WRITE(fname,'(a,i2.2)') 'nr_particles_',cpl_id
815             OPEN (333,file = fname)
816             is_file_open = .true.
817          ENDIF
818          WRITE(333,'(f12.2,3i10)') my_time,local_nr_particles + child_nr_particles,local_nr_particles,child_nr_particles
819       ENDIF
820    ENDIF
821
822#endif
823 END SUBROUTINE pmcp_g_print_number_of_particles
824
825
826!------------------------------------------------------------------------------!
827!------------------------------------------------------------------------------!
828! Private subroutines
829!------------------------------------------------------------------------------!
830!------------------------------------------------------------------------------!
831
832!------------------------------------------------------------------------------!
833! Description:
834! ------------
835!> child routine
836!> update the size of the local buffer (coarse_particles)
837!------------------------------------------------------------------------------!
838 SUBROUTINE check_and_alloc_coarse_particle (ic,jc,nr,with_copy)
839   
840    IMPLICIT NONE
841   
842    INTEGER(iwp),INTENT(IN) ::  ic !< coarse x grid index
843    INTEGER(iwp),INTENT(IN) ::  jc !< coarse y grid index
844    INTEGER(iwp),INTENT(IN) ::  nr !< number of particles to be transferred/stored in local buffer
845   
846    LOGICAL,INTENT(IN),OPTIONAL ::  with_copy !< copy particles in buffer? or reallocate empty buffer
847
848    LOGICAL :: with_copy_lo !< local variable of with copy
849   
850    INTEGER(iwp) ::  new_size !< new size of the local buffer
851    INTEGER(iwp) ::  old_size !< old size of the local buffer
852   
853    TYPE(particle_type), DIMENSION(:), ALLOCATABLE ::  tmp_particles_d !<
854
855#if defined( __parallel )
856    with_copy_lo = .FALSE.
857    IF ( PRESENT( with_copy ) ) with_copy_lo = with_copy
858
859    IF ( .NOT. ALLOCATED( coarse_particles(jc,ic)%parent_particles ) ) THEN
860       new_size = MAX( nr, Min_particles_per_column )
861       ALLOCATE( coarse_particles(jc,ic)%parent_particles(new_size) )
862    ELSEIF ( nr > SIZE( coarse_particles(jc,ic)%parent_particles ) ) THEN
863
864       old_size = SIZE( coarse_particles(jc,ic)%parent_particles )
865       new_size = old_size * ( 1.0_wp + alloc_factor / 100.0_wp )
866       new_size = MAX( nr, new_size, old_size + Min_particles_per_column )
867
868!
869!--    Copy existing particles to new particle buffer
870       IF ( with_copy_lo ) THEN
871          ALLOCATE( tmp_particles_d(old_size) )
872          tmp_particles_d(1:old_size) = coarse_particles(jc,ic)%parent_particles
873
874          DEALLOCATE( coarse_particles(jc,ic)%parent_particles )
875          ALLOCATE( coarse_particles(jc,ic)%parent_particles(new_size) )
876
877          coarse_particles(jc,ic)%parent_particles(1:old_size)          = tmp_particles_d(1:old_size)
878          coarse_particles(jc,ic)%parent_particles(old_size+1:new_size) = zero_particle
879
880          DEALLOCATE( tmp_particles_d )
881!
882!--    allocate or reallocate an empty buffer
883       ELSE
884          DEALLOCATE( coarse_particles(jc,ic)%parent_particles )
885          ALLOCATE( coarse_particles(jc,ic)%parent_particles(new_size) )
886       ENDIF
887    ENDIF
888
889#endif
890 END SUBROUTINE check_and_alloc_coarse_particle
891
892
893!------------------------------------------------------------------------------!
894! Description:
895! ------------
896!> child routine:
897!> copy/sort particles out of the local buffer into the respective grid boxes
898!------------------------------------------------------------------------------!
899 SUBROUTINE c_copy_particle_to_child_grid
900 
901    IMPLICIT NONE
902 
903    INTEGER(iwp) ::  ic  !< coarse x grid index
904    INTEGER(iwp) ::  icl !< left boundary in coarse(parent) index space
905    INTEGER(iwp) ::  icr !< right boundary in coarse(parent) index space
906    INTEGER(iwp) ::  ip  !< x grid index
907    INTEGER(iwp) ::  jc  !< coarse y grid index
908    INTEGER(iwp) ::  jcn !< north boundary in coarse(parent) index space
909    INTEGER(iwp) ::  jcs !< south boundary in coarse(parent) index space
910    INTEGER(iwp) ::  jp  !< y grid index
911    INTEGER(iwp) ::  kp  !< z grid index
912    INTEGER(iwp) ::  n   !< loop index (number of particles)
913    INTEGER(iwp) ::  nr  !< short variable for name number or particles
914   
915    REAL(wp) ::  xc  !< child x coordinate
916    REAL(wp) ::  xoc !< child x origin
917    REAL(wp) ::  yc  !< child y coordinate
918    REAL(wp) ::  yoc !< child y origin
919    REAL(wp) ::  zc  !< child z coordinate
920
921#if defined( __parallel )
922!
923!-- Child domain boundaries in the parent index space
924    icl = coarse_bound(1)
925    icr = coarse_bound(2)
926    jcs = coarse_bound(3)
927    jcn = coarse_bound(4)
928
929    DO ic = icl, icr
930       DO jc = jcs, jcn
931          nr = coarse_particles(jc,ic)%nr_particle
932
933          IF ( nr > 0 ) THEN
934             DO n = 1, nr
935                xc =  coarse_particles(jc,ic)%parent_particles(n)%x-lower_left_coord_x
936                yc =  coarse_particles(jc,ic)%parent_particles(n)%y-lower_left_coord_y
937                zc =  coarse_particles(jc,ic)%parent_particles(n)%z
938                xoc = coarse_particles(jc,ic)%parent_particles(n)%origin_x-lower_left_coord_x
939                yoc = coarse_particles(jc,ic)%parent_particles(n)%origin_y-lower_left_coord_y
940                ip = xc / dx
941                jp = yc / dy
942                kp = nzt
943                DO WHILE ( zw(kp-1) > zc .AND. kp > nzb + 1 )         ! kk search loop has to be optimzed !!!
944                   kp = kp - 1
945                ENDDO 
946
947                IF ( ip >= nxl .AND. ip <= nxr .AND. jp >= nys .AND. jp <= nyn ) THEN
948                   prt_count(kp,jp,ip) = prt_count(kp,jp,ip) + 1
949                   IF ( prt_count(kp,jp,ip) > SIZE( grid_particles(kp,jp,ip)%particles ) ) THEN
950                      CALL realloc_particles_array( ip, jp, kp )
951                   ENDIF
952                   coarse_particles(jc,ic)%parent_particles(n)%x = xc                   ! Adjust coordinates to child grid
953                   coarse_particles(jc,ic)%parent_particles(n)%y = yc
954                   coarse_particles(jc,ic)%parent_particles(n)%origin_x = xoc           ! Adjust origins to child grid
955                   coarse_particles(jc,ic)%parent_particles(n)%origin_y = yoc
956                   grid_particles(kp,jp,ip)%particles(prt_count(kp,jp,ip)) = coarse_particles(jc,ic)%parent_particles(n)
957                ENDIF
958             ENDDO
959          ENDIF
960       ENDDO
961    ENDDO
962
963#endif
964 END SUBROUTINE c_copy_particle_to_child_grid
965
966
967!------------------------------------------------------------------------------!
968! Description:
969! ------------
970!> child routine:
971!> copy particles which left the model area into the local buffer
972!------------------------------------------------------------------------------!
973 SUBROUTINE c_copy_particle_to_coarse_grid
974 
975    IMPLICIT NONE
976   
977    LOGICAL ::  boundary_particle !<
978   
979    INTEGER(iwp) ::  i    !< loop index (x grid)
980    INTEGER(iwp) ::  ic   !< loop index (coarse x grid)
981    INTEGER(iwp) ::  icl  !< left boundary in coarse(parent) index space
982    INTEGER(iwp) ::  icr  !< left boundary in coarse(parent) index space
983    INTEGER(iwp) ::  ierr !< error code
984    INTEGER(iwp) ::  j    !< loop index (y grid)
985    INTEGER(iwp) ::  jc   !< loop index (coarse y grid)
986    INTEGER(iwp) ::  jcs  !< south boundary in coarse(parent) index space
987    INTEGER(iwp) ::  jcn  !< north boundary in coarse(parent) index space
988    INTEGER(iwp) ::  k    !< loop index (z grid)
989    INTEGER(iwp) ::  n    !< loop index (number of particles)
990   
991    REAL(wp) ::  x       !< x coordinate
992    REAL(wp) ::  xo      !< x origin
993    REAL(wp) ::  x_left  !< absolute left boundary
994    REAL(wp) ::  x_right !< absolute right boundary
995    REAL(wp) ::  y       !< left boundary in coarse(parent) index space
996    REAL(wp) ::  yo      !< y origin
997    REAL(wp) ::  y_north !< absolute north boundary
998    REAL(wp) ::  y_south !< absoulte south boundary
999    REAL(wp) ::  z       !< z coordinate
1000
1001#if defined( __parallel )
1002!
1003!-- Child domain boundaries in the parent index space
1004
1005    icl = coarse_bound(1)
1006    icr = coarse_bound(2)
1007    jcs = coarse_bound(3)
1008    jcn = coarse_bound(4)
1009
1010!
1011!-- absolute coordinates
1012    x_left  = coord_x(0)
1013    x_right = coord_x(nx) + dx
1014    y_south = coord_y(0)
1015    y_north = coord_y(ny) + dy
1016
1017!   Clear Particle Buffer
1018
1019    DO ic = icl, icr
1020       DO jc = jcs, jcn
1021          coarse_particles(jc,ic)%nr_particle = 0
1022       ENDDO
1023    ENDDO
1024
1025!
1026!-- Determine particles which leave the inner area in east or west dirextion
1027!-- Compute only first (nxl) and last (nxr) loop iterration.
1028    DO  i = nxl, nxr
1029       DO  j = nys, nyn
1030          DO  k = nzb + 1, nzt
1031             number_of_particles = prt_count(k,j,i)
1032             IF ( number_of_particles <= 0 )  CYCLE
1033             particles => grid_particles(k,j,i)%particles(1:number_of_particles)
1034             DO  n = 1, number_of_particles
1035                IF ( particles(n)%particle_mask )  THEN
1036                   x =  particles(n)%x+ lower_left_coord_x 
1037                   y =  particles(n)%y+ lower_left_coord_y 
1038                   xo = particles(n)%origin_x + lower_left_coord_x
1039                   yo = particles(n)%origin_y + lower_left_coord_y
1040                   z =  particles(n)%z
1041                   
1042                   boundary_particle = .FALSE.
1043                   boundary_particle = boundary_particle .OR. (x < x_left  .OR. x > x_right)
1044                   boundary_particle = boundary_particle .OR. (y < y_south .OR. y > y_north)
1045                   boundary_particle = boundary_particle .OR. (z > zw(nzt))
1046
1047                   IF ( boundary_particle ) THEN                     
1048                      ic = x / cg%dx                     !TODO anpassen auf Mehrfachnesting
1049                      jc = y / cg%dy
1050
1051                      IF ( ic >= icl .AND. ic <= icr .AND. jc >= jcs .AND. jc <= jcn ) THEN
1052                         coarse_particles(jc,ic)%nr_particle = coarse_particles(jc,ic)%nr_particle + 1
1053                         CALL check_and_alloc_coarse_particle( ic, jc, coarse_particles(jc,ic)%nr_particle, with_copy=.TRUE. )
1054
1055                         coarse_particles(jc,ic)%parent_particles(coarse_particles(jc,ic)%nr_particle)   = particles(n)
1056                         coarse_particles(jc,ic)%parent_particles(coarse_particles(jc,ic)%nr_particle)%x = x   !adapt to global coordinates
1057                         coarse_particles(jc,ic)%parent_particles(coarse_particles(jc,ic)%nr_particle)%y = y
1058                         coarse_particles(jc,ic)%parent_particles(coarse_particles(jc,ic)%nr_particle)%origin_x = xo
1059                         coarse_particles(jc,ic)%parent_particles(coarse_particles(jc,ic)%nr_particle)%origin_y = yo
1060!--                      Mark particle as deleted after copying it to the transfer buffer
1061                         grid_particles(k,j,i)%particles(n)%particle_mask = .FALSE.
1062                      ELSE
1063                         WRITE(9,'(a,10i6)') 'This should not happen ',i,j,k,ic,jc,icl,icr,jcs,jcn
1064                         CALL MPI_Abort( MPI_COMM_WORLD, 9999, ierr )
1065                      ENDIF
1066                   ENDIF
1067                ENDIF
1068             ENDDO
1069          ENDDO
1070       ENDDO
1071    ENDDO
1072
1073!
1074!- Pack particles (eliminate those marked for deletion),
1075!- determine new number of particles
1076   CALL lpm_sort_in_subboxes
1077
1078#endif
1079 END SUBROUTINE c_copy_particle_to_coarse_grid
1080
1081
1082!------------------------------------------------------------------------------!
1083! Description:
1084! ------------
1085!> parent routine:
1086!> copy/sort particles from the MPI window into the respective grid boxes
1087!------------------------------------------------------------------------------!
1088 SUBROUTINE p_copy_particle_to_org_grid( m, child_id )
1089
1090    IMPLICIT NONE
1091
1092    INTEGER(iwp),INTENT(IN) ::  child_id !<
1093    INTEGER(iwp),INTENT(IN) ::  m        !<
1094
1095    INTEGER(iwp) ::  i      !< loop index (x grid)
1096    INTEGER(iwp) ::  j      !< loop index (y grid)
1097    INTEGER(iwp) ::  k      !< loop index (z grid)
1098    INTEGER(iwp) ::  n      !< loop index (nr part)
1099    INTEGER(iwp) ::  nr     !< short variable name for nr_part
1100    INTEGER(iwp) ::  pindex !< short variable name part_adr
1101   
1102    REAL(wp) ::  x  !< x coordinate
1103    REAL(wp) ::  xo !< x origin
1104    REAL(wp) ::  y  !< y coordinate
1105    REAL(wp) ::  yo !< y origin
1106    REAL(wp) ::  z  !< z coordinate
1107   
1108    INTEGER(iwp),DIMENSION(1) ::  buf_shape !<
1109
1110#if defined( __parallel )
1111    buf_shape(1) = max_nr_particle_in_rma_win
1112    CALL C_F_POINTER( buf_ptr(m), particle_in_win , buf_shape )
1113
1114    DO i = nxl, nxr
1115       DO j = nys, nyn
1116          nr = nr_part(j,i)
1117          IF ( nr > 0 ) THEN
1118             pindex = part_adr(j,i)
1119             DO n = 1, nr
1120                x = particle_in_win(pindex)%x-lower_left_coord_x
1121                y = particle_in_win(pindex)%y-lower_left_coord_y
1122                z = particle_in_win(pindex)%z
1123                xo = particle_in_win(pindex)%origin_x-lower_left_coord_x
1124                yo = particle_in_win(pindex)%origin_y-lower_left_coord_y
1125                k = nzt + 1
1126                DO WHILE ( zw(k-1) > z .AND. k > nzb + 1 )                               ! kk search loop has to be optimzed !!!
1127                   k = k - 1
1128                END DO
1129
1130                prt_count(k,j,i) = prt_count(k,j,i) + 1
1131                IF ( prt_count(k,j,i) > SIZE( grid_particles(k,j,i)%particles ) ) THEN
1132                   CALL realloc_particles_array( i, j, k )
1133                ENDIF
1134                grid_particles(k,j,i)%particles(prt_count(k,j,i)) = particle_in_win(pindex)
1135               
1136!
1137!--             Update particle positions and origins relative to parent domain
1138                grid_particles(k,j,i)%particles(prt_count(k,j,i))%x = x
1139                grid_particles(k,j,i)%particles(prt_count(k,j,i))%y = y
1140                grid_particles(k,j,i)%particles(prt_count(k,j,i))%origin_x = xo
1141                grid_particles(k,j,i)%particles(prt_count(k,j,i))%origin_y = yo
1142                pindex = pindex + 1
1143             ENDDO
1144          ENDIF
1145       ENDDO
1146    ENDDO
1147
1148#endif
1149 END SUBROUTINE p_copy_particle_to_org_grid
1150
1151
1152END MODULE pmc_particle_interface
Note: See TracBrowser for help on using the repository browser.