Home

Context Navigation

source: palm/trunk/SOURCE/pmc_parent_mod.f90 @ 4558

Last change on this file since 4558 was 4360, checked in by suehring, 5 years ago
Bugfix in output of time-averaged plant-canopy quanities; Output of plant-canopy data only where tall canopy is defined; land-surface model: fix wrong location strings; tests: update urban test case; all source code files: copyright update
Property svn:keywords set to `Id`
File size: 34.4 KB

Rev	Line
[4246]	1	MODULE pmc_parent
	2
	3	!------------------------------------------------------------------------------!
	4	! This file is part of the PALM model system.
	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	!
[4360]	18	! Copyright 1997-2020 Leibniz Universitaet Hannover
[4246]	19	!------------------------------------------------------------------------------!
	20	!
	21	! Current revisions:
	22	! ------------------
[4213]	23	!
	24	!
[4246]	25	! Former revisions:
	26	! -----------------
	27	! $Id: pmc_parent_mod.f90 4360 2020-01-07 11:25:50Z moh.hefny $
	28	!
	29	!
[4245]	30	! 4213 2019-09-02 14:25:56Z suehring
[4246]	31	! Allocate array for index_list_2d also for zero-size arrays, in order to
	32	! avoid errors when array bound checks are enabled
	33	!
	34	! 4212 2019-09-02 14:23:05Z suehring
	35	! Corrected "Former revisions" section
	36	!
	37	! 3962 2019-05-08 19:40:33Z suehring
	38	! Bugfixes in initial settings of child and parent communication patterns.
	39	!
	40	! 3655 2019-01-07 16:51:22Z knoop
	41	! explicit kind settings
	42	!
	43	! 1762 2016-02-25 12:31:13Z hellstea
	44	! Initial revision by K. Ketelsen
	45	!
	46	! Description:
	47	! ------------
	48	!> Parent part of Palm Model Coupler
	49	!------------------------------------------------------------------------------!
	50
	51	#if defined( __parallel )
	52	USE, INTRINSIC :: ISO_C_BINDING
	53
	54	#if !defined( __mpifh )
	55	USE MPI
	56	#endif
	57	USE kinds
	58	USE pmc_general, &
	59	ONLY: arraydef, childdef, da_namedef, da_namelen, pedef, &
	60	pmc_g_setname, pmc_max_array, pmc_max_models, pmc_sort
	61
	62	USE pmc_handle_communicator, &
	63	ONLY: m_model_comm,m_model_rank,m_model_npes, m_to_child_comm, &
	64	m_world_rank, pmc_parent_for_child
	65
	66	USE pmc_mpi_wrapper, &
	67	ONLY: pmc_alloc_mem, pmc_bcast, pmc_time
	68
	69	IMPLICIT NONE
	70
	71	#if defined( __mpifh )
	72	INCLUDE "mpif.h"
	73	#endif
	74
	75	PRIVATE
	76	SAVE
	77
	78	TYPE childindexdef
	79	INTEGER :: nrpoints !<
	80	INTEGER, DIMENSION(:,:), ALLOCATABLE :: index_list_2d !<
	81	END TYPE childindexdef
	82
	83	TYPE(childdef), DIMENSION(pmc_max_models),PUBLIC :: children !<
	84	TYPE(childindexdef), DIMENSION(pmc_max_models) :: indchildren !<
	85
	86	INTEGER :: next_array_in_list = 0 !<
	87
	88
	89	PUBLIC pmc_parent_for_child
	90
	91
	92	INTERFACE pmc_parentinit
	93	MODULE PROCEDURE pmc_parentinit
	94	END INTERFACE pmc_parentinit
	95
	96	INTERFACE pmc_s_set_2d_index_list
	97	MODULE PROCEDURE pmc_s_set_2d_index_list
	98	END INTERFACE pmc_s_set_2d_index_list
	99
	100	INTERFACE pmc_s_clear_next_array_list
	101	MODULE PROCEDURE pmc_s_clear_next_array_list
	102	END INTERFACE pmc_s_clear_next_array_list
	103
	104	INTERFACE pmc_s_getnextarray
	105	MODULE PROCEDURE pmc_s_getnextarray
	106	END INTERFACE pmc_s_getnextarray
	107
	108	INTERFACE pmc_s_set_dataarray
	109	MODULE PROCEDURE pmc_s_set_dataarray_2d
	110	MODULE PROCEDURE pmc_s_set_dataarray_3d
	111	MODULE PROCEDURE pmc_s_set_dataarray_ip2d
	112	END INTERFACE pmc_s_set_dataarray
	113
	114	INTERFACE pmc_s_setind_and_allocmem
	115	MODULE PROCEDURE pmc_s_setind_and_allocmem
	116	END INTERFACE pmc_s_setind_and_allocmem
	117
	118	INTERFACE pmc_s_fillbuffer
	119	MODULE PROCEDURE pmc_s_fillbuffer
	120	END INTERFACE pmc_s_fillbuffer
	121
	122	INTERFACE pmc_s_getdata_from_buffer
	123	MODULE PROCEDURE pmc_s_getdata_from_buffer
	124	END INTERFACE pmc_s_getdata_from_buffer
	125
	126	INTERFACE pmc_s_set_active_data_array
	127	MODULE PROCEDURE pmc_s_set_active_data_array
	128	END INTERFACE pmc_s_set_active_data_array
	129
	130	INTERFACE pmc_s_get_child_npes
	131	MODULE PROCEDURE pmc_s_get_child_npes
	132	END INTERFACE pmc_s_get_child_npes
	133
	134	PUBLIC pmc_parentinit, pmc_s_clear_next_array_list, pmc_s_fillbuffer, &
	135	pmc_s_getdata_from_buffer, pmc_s_getnextarray, &
	136	pmc_s_setind_and_allocmem, pmc_s_set_active_data_array, &
	137	pmc_s_set_dataarray, pmc_s_set_2d_index_list, &
	138	pmc_s_get_child_npes
	139
	140	CONTAINS
	141
	142
	143	SUBROUTINE pmc_parentinit
	144
	145	IMPLICIT NONE
	146
	147	INTEGER(iwp) :: childid !<
	148	INTEGER(iwp) :: i !<
	149	INTEGER(iwp) :: j !<
	150	INTEGER(iwp) :: istat !<
	151
	152	DO i = 1, SIZE( pmc_parent_for_child )-1
	153
	154	childid = pmc_parent_for_child( i )
	155
	156	children(childid)%model_comm = m_model_comm
	157	children(childid)%inter_comm = m_to_child_comm(childid)
	158
	159	!
	160	!-- Get rank and size
	161	CALL MPI_COMM_RANK( children(childid)%model_comm, &
	162	children(childid)%model_rank, istat )
	163	CALL MPI_COMM_SIZE( children(childid)%model_comm, &
	164	children(childid)%model_npes, istat )
	165	CALL MPI_COMM_REMOTE_SIZE( children(childid)%inter_comm, &
	166	children(childid)%inter_npes, istat )
	167	!
	168	!-- Intra communicator is used for MPI_GET
	169	CALL MPI_INTERCOMM_MERGE( children(childid)%inter_comm, .FALSE., &
	170	children(childid)%intra_comm, istat )
	171	CALL MPI_COMM_RANK( children(childid)%intra_comm, &
	172	children(childid)%intra_rank, istat )
	173
	174	ALLOCATE( children(childid)%pes(children(childid)%inter_npes))
	175	!
	176	!-- Allocate array of TYPE arraydef for all child PEs to store information
	177	!-- of the transfer array
	178	DO j = 1, children(childid)%inter_npes
	179	ALLOCATE( children(childid)%pes(j)%array_list(pmc_max_array) )
	180	ENDDO
	181
	182	CALL get_da_names_from_child (childid)
	183
	184	ENDDO
	185
	186	END SUBROUTINE pmc_parentinit
	187
	188
	189
	190	SUBROUTINE pmc_s_set_2d_index_list( childid, index_list )
	191
	192	IMPLICIT NONE
	193
	194	INTEGER(iwp), INTENT(IN) :: childid !<
	195	INTEGER(iwp), DIMENSION(:,:), INTENT(INOUT) :: index_list !<
	196
	197	INTEGER(iwp) :: ian !<
	198	INTEGER(iwp) :: ie !<
	199	INTEGER(iwp) :: ip !<
	200	INTEGER(iwp) :: is !<
	201	INTEGER(iwp) :: istat !<
	202
	203
	204	IF ( m_model_rank == 0 ) THEN
	205	!
	206	!-- Sort to ascending parent process order
	207	CALL pmc_sort( index_list, 6 )
	208	is = 1
	209	DO ip = 0, m_model_npes-1
	210	!
	211	!-- Split into parent processes
	212	ie = is - 1
	213	!
	214	!-- There may be no entry for this process
	215	IF ( is <= SIZE( index_list,2 ) .AND. ie >= 0 ) THEN
	216	DO WHILE ( index_list(6,ie+1 ) == ip )
	217	ie = ie + 1
	218	IF ( ie == SIZE( index_list,2 ) ) EXIT
	219	ENDDO
	220	ian = ie - is + 1
	221	ELSE
	222	is = -1
	223	ie = -2
	224	ian = 0
	225	ENDIF
	226	!
	227	!-- Send data to other parent processes
	228	IF ( ip == 0 ) THEN
	229	indchildren(childid)%nrpoints = ian
	230	!
	231	!-- Allocate array for index_list_2d. Note, the array will also be
	232	!-- allocated in case ian=0, in order to avoid errors when
	233	!-- array bounds are checked.
	234	ALLOCATE( indchildren(childid)%index_list_2d(6,1:ian) )
	235	IF ( ian > 0) THEN
	236	indchildren(childid)%index_list_2d(:,1:ian) = &
	237	index_list(:,is:ie)
	238	ENDIF
	239	ELSE
	240	CALL MPI_SEND( ian, 1, MPI_INTEGER, ip, 1000, m_model_comm, &
	241	istat )
	242	IF ( ian > 0) THEN
	243	CALL MPI_SEND( index_list(1,is), 6*ian, MPI_INTEGER, ip, &
	244	1001, m_model_comm, istat )
	245	ENDIF
	246	ENDIF
	247	is = ie + 1
	248	ENDDO
	249	ELSE
	250	CALL MPI_RECV( indchildren(childid)%nrpoints, 1, MPI_INTEGER, 0, 1000, &
	251	m_model_comm, MPI_STATUS_IGNORE, istat )
	252	ian = indchildren(childid)%nrpoints
	253	!
	254	!-- Allocate array for index_list_2d. Note, the array will also be
	255	!-- allocated in case ian=0, in order to avoid errors when
	256	!-- array bounds are checked.
	257	ALLOCATE( indchildren(childid)%index_list_2d(6,1:ian) )
	258	IF ( ian > 0 ) THEN
	259	CALL MPI_RECV( indchildren(childid)%index_list_2d, 6*ian, &
	260	MPI_INTEGER, 0, 1001, m_model_comm, &
	261	MPI_STATUS_IGNORE, istat)
	262	ENDIF
	263	ENDIF
	264	CALL set_pe_index_list( children(childid), &
	265	indchildren(childid)%index_list_2d, &
	266	indchildren(childid)%nrpoints )
	267
	268	END SUBROUTINE pmc_s_set_2d_index_list
	269
	270
	271
	272	SUBROUTINE pmc_s_clear_next_array_list
	273
	274	IMPLICIT NONE
	275
	276	next_array_in_list = 0
	277
	278	END SUBROUTINE pmc_s_clear_next_array_list
	279
	280
	281
	282	LOGICAL FUNCTION pmc_s_getnextarray( childid, myname )
	283
	284	!
	285	!-- Althoug there are no linked lists any more in PMC, this call still looks like working with a list
	286
	287	CHARACTER(LEN=*), INTENT(OUT) :: myname !<
	288	INTEGER(iwp), INTENT(IN) :: childid !<
	289
	290	TYPE(arraydef), POINTER :: ar
	291	TYPE(pedef), POINTER :: ape
	292
	293	next_array_in_list = next_array_in_list + 1
	294	!
	295	!-- Array names are the same on all children processes, so take first
	296	!-- process to get the name
	297	ape => children(childid)%pes(1)
	298
	299	IF ( next_array_in_list > ape%nr_arrays ) THEN
	300	!
	301	!-- All arrays are done
	302	pmc_s_getnextarray = .FALSE.
	303	RETURN
	304	ENDIF
	305
	306	ar => ape%array_list(next_array_in_list)
	307	myname = ar%name
	308	!
	309	!-- Return true if there is still an array in the list
	310
	311	pmc_s_getnextarray = .TRUE.
	312
	313	END FUNCTION pmc_s_getnextarray
	314
	315
	316
	317	SUBROUTINE pmc_s_set_dataarray_2d( childid, array, array_2 )
	318
	319	IMPLICIT NONE
	320
	321	INTEGER(iwp), INTENT(IN) :: childid !<
	322
	323	REAL(wp), INTENT(IN), DIMENSION(:,:), POINTER :: array !<
	324	REAL(wp), INTENT(IN), DIMENSION(:,:), POINTER, OPTIONAL :: array_2 !<
	325
	326	INTEGER(iwp) :: nrdims !<
	327	INTEGER(iwp), DIMENSION(4) :: dims !<
	328	TYPE(C_PTR) :: array_adr !<
	329	TYPE(C_PTR) :: second_adr !<
	330
	331
	332	dims = 1
	333	nrdims = 2
	334	dims(1) = SIZE( array,1 )
	335	dims(2) = SIZE( array,2 )
	336	array_adr = C_LOC( array )
	337
	338	IF ( PRESENT( array_2 ) ) THEN
	339	second_adr = C_LOC(array_2)
	340	CALL pmc_s_setarray( childid, nrdims, dims, array_adr, &
	341	second_adr = second_adr)
	342	ELSE
	343	CALL pmc_s_setarray( childid, nrdims, dims, array_adr )
	344	ENDIF
	345
	346	END SUBROUTINE pmc_s_set_dataarray_2d
	347
	348	SUBROUTINE pmc_s_set_dataarray_ip2d( childid, array )
	349
	350	IMPLICIT NONE
	351
	352	INTEGER(iwp),INTENT(IN) :: childid !<
	353
	354	INTEGER(idp), INTENT(IN), DIMENSION(:,:), POINTER :: array !<
	355
	356	INTEGER(iwp) :: nrdims !<
	357	INTEGER(iwp), DIMENSION(4) :: dims !<
	358	TYPE(C_PTR) :: array_adr !<
	359
	360
	361	dims = 1
	362	nrdims = 2
	363	dims(1) = SIZE( array,1 )
	364	dims(2) = SIZE( array,2 )
	365	array_adr = C_LOC( array )
	366
	367	CALL pmc_s_setarray( childid, nrdims, dims, array_adr , dimkey=22)
	368
	369	END SUBROUTINE pmc_s_set_dataarray_ip2d
	370
	371
	372	SUBROUTINE pmc_s_set_dataarray_3d( childid, array, nz_cl, nz, array_2 )
	373
	374	IMPLICIT NONE
	375
	376	INTEGER(iwp), INTENT(IN) :: childid !<
	377	INTEGER(iwp), INTENT(IN) :: nz !<
	378	INTEGER(iwp), INTENT(IN) :: nz_cl !<
	379
	380	REAL(wp), INTENT(IN), DIMENSION(:,:,:), POINTER :: array !<
	381	REAL(wp), INTENT(IN), DIMENSION(:,:,:), POINTER, OPTIONAL :: array_2 !<
	382
	383	INTEGER(iwp) :: nrdims !<
	384	INTEGER(iwp), DIMENSION(4) :: dims !<
	385	TYPE(C_PTR) :: array_adr !<
	386	TYPE(C_PTR) :: second_adr !<
	387
	388	nrdims = 3
	389	dims(1) = SIZE( array,1 )
	390	dims(2) = SIZE( array,2 )
	391	dims(3) = SIZE( array,3 )
	392	dims(4) = nz_cl+dims(1)-nz ! works for first dimension 1:nz and 0:nz+1
	393
	394	array_adr = C_LOC(array)
	395	!
	396	!-- In PALM's pointer version, two indices have to be stored internally.
	397	!-- The active address of the data array is set in swap_timelevel.
	398	IF ( PRESENT( array_2 ) ) THEN
	399	second_adr = C_LOC( array_2 )
	400	CALL pmc_s_setarray( childid, nrdims, dims, array_adr, &
	401	second_adr = second_adr)
	402	ELSE
	403	CALL pmc_s_setarray( childid, nrdims, dims, array_adr )
	404	ENDIF
	405
	406	END SUBROUTINE pmc_s_set_dataarray_3d
	407
	408
	409
	410	SUBROUTINE pmc_s_setind_and_allocmem( childid )
	411
	412	USE control_parameters, &
	413	ONLY: message_string
	414
	415	IMPLICIT NONE
	416
	417	!
	418	!-- Naming convention for appendices: _pc -> parent to child transfer
	419	!-- _cp -> child to parent transfer
	420	!-- send -> parent to child transfer
	421	!-- recv -> child to parent transfer
	422	INTEGER(iwp), INTENT(IN) :: childid !<
	423
	424	INTEGER(iwp) :: arlen !<
	425	INTEGER(iwp) :: i !<
	426	INTEGER(iwp) :: ierr !<
	427	INTEGER(iwp) :: j !<
	428	INTEGER(iwp) :: lo_nr_arrays !< store number of arrays in local variiab le
	429	INTEGER(iwp) :: myindex !<
	430	INTEGER(iwp) :: total_npes !< Total Number of PEs Parent and Child
	431
	432	INTEGER(idp) :: bufsize !< size of MPI data window
	433	INTEGER(KIND=MPI_ADDRESS_KIND) :: winsize !<
	434
	435	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: myindex_s
	436	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: myindex_r
	437
	438	TYPE(C_PTR) :: base_ptr !<
	439	TYPE(pedef), POINTER :: ape !<
	440	TYPE(arraydef), POINTER :: ar !<
	441
	442	REAL(wp),DIMENSION(:), POINTER, SAVE :: base_array_pc !< base array for parent to child transfer
	443	REAL(wp),DIMENSION(:), POINTER, SAVE :: base_array_cp !< base array for child to parent transfer
	444
	445	call MPI_COMM_SIZE (children(childid)%intra_comm, total_npes, ierr)
	446	!
	447	!-- Parent to child direction
	448	myindex = 1
	449	bufsize = 8
	450	!
	451	!-- All Child processes get the same number of arrays
	452	!-- Therfore the number of arrays form the first Child process can be used for Dimension.
	453	lo_nr_arrays = children(childid)%pes(1)%nr_arrays
	454
	455	ALLOCATE( myindex_s(lo_nr_arrays,0:total_npes-1) )
	456	ALLOCATE( myindex_r(lo_nr_arrays,0:total_npes-1) )
	457
	458	myindex_s = 0
	459
	460	!
	461	!-- First stride: compute size and set index
	462	DO i = 1, children(childid)%inter_npes
	463
	464	ape => children(childid)%pes(i)
	465
	466	DO j = 1, ape%nr_arrays
	467
	468	ar => ape%array_list(j)
	469	IF ( ar%nrdims == 2 ) THEN
	470	arlen = ape%nrele
	471	ELSEIF ( ar%nrdims == 3 ) THEN
	472	arlen = ape%nrele * ar%a_dim(4)
	473	ELSE
	474	arlen = -1
	475	ENDIF
	476	ar%sendindex = myindex
	477	!
	478	! Using intra communicator for MPU_Alltoall, the numbers of the child processes are after the paremt ones
	479
	480	myindex_s(j,i-1+children(childid)%model_npes) = myindex
	481
	482	myindex = myindex + arlen
	483	bufsize = bufsize + arlen
	484	ar%sendsize = arlen
	485	ENDDO
	486
	487	ENDDO
	488	!
	489	!-- Using MPI_Alltoall to send indices from Parent to Child
	490	!-- The data comming back from the child processes are ignored.
	491
	492	CALL MPI_ALLTOALL( myindex_s, lo_nr_arrays, MPI_INTEGER, &
	493	myindex_r, lo_nr_arrays, MPI_INTEGER, &
	494	children(childid)%intra_comm, ierr )
	495
	496	!
	497	!-- Using MPI_Alltoall to receive indices from Child
	498	myindex_s = 0
	499	myindex_r = 0
	500
	501	CALL MPI_ALLTOALL( myindex_s, lo_nr_arrays, MPI_INTEGER, &
	502	myindex_r, lo_nr_arrays, MPI_INTEGER, &
	503	children(childid)%intra_comm, ierr )
	504	!
	505	!-- Create RMA (One Sided Communication) window for data buffer parent to
	506	!-- child transfer.
	507	!-- The buffer of MPI_GET (counterpart of transfer) can be PE-local, i.e.
	508	!-- it can but must not be part of the MPI RMA window. Only one RMA window is
	509	!-- required to prepare the data for
	510	!-- parent -> child transfer on the parent side
	511	!-- and for
	512	!-- child -> parent transfer on the child side
	513	CALL pmc_alloc_mem( base_array_pc, bufsize )
	514	children(childid)%totalbuffersize = bufsize * wp
	515
	516	winsize = bufsize * wp
	517	CALL MPI_WIN_CREATE( base_array_pc, winsize, wp, MPI_INFO_NULL, &
	518	children(childid)%intra_comm, &
	519	children(childid)%win_parent_child, ierr )
	520	!
	521	!-- Open window to set data
	522	CALL MPI_WIN_FENCE( 0, children(childid)%win_parent_child, ierr )
	523	!
	524	!-- Second stride: set buffer pointer
	525	DO i = 1, children(childid)%inter_npes
	526
	527	ape => children(childid)%pes(i)
	528
	529	DO j = 1, ape%nr_arrays
	530
	531	ar => ape%array_list(j)
	532	ar%sendbuf = C_LOC( base_array_pc(ar%sendindex) )
	533
	534	IF ( ar%sendindex + ar%sendsize > bufsize ) THEN
	535	WRITE( message_string, '(a,i4,4i7,1x,a)' ) &
	536	'parent buffer too small ',i, &
	537	ar%sendindex,ar%sendsize,ar%sendindex+ar%sendsize, &
	538	bufsize,trim(ar%name)
	539	CALL message( 'pmc_s_setind_and_allocmem', 'PA0429', 3, 2, 0, 6, 0 )
	540	ENDIF
	541	ENDDO
	542	ENDDO
	543	!
	544	!-- Child to parent direction
	545	bufsize = 8
	546	!
	547	!-- First stride: compute size and set index
	548	DO i = 1, children(childid)%inter_npes
	549	ape => children(childid)%pes(i)
	550	DO j = 1, ape%nr_arrays
	551	ar => ape%array_list(j)
	552	!
	553	!-- Receive index from child
	554	IF ( ar%nrdims == 3 ) THEN
	555	bufsize = MAX( bufsize, &
	556	INT( ape%nrele * ar%a_dim(4), MPI_ADDRESS_KIND ) )
	557	ELSE
	558	bufsize = MAX( bufsize, INT( ape%nrele, MPI_ADDRESS_KIND ) )
	559	ENDIF
	560	ar%recvindex = myindex_r(j,i-1+children(childid)%model_npes)
	561	ENDDO
	562	ENDDO
	563
	564	DEALLOCATE( myindex_s )
	565	DEALLOCATE( myindex_r )
	566
	567	!
	568	!-- Create RMA (one sided communication, RMA = Remote Memory Access) data buffer.
	569	!-- The buffer for MPI_GET can be PE local, i.e. it can but must not be part of
	570	!-- the MPI RMA window
	571	CALL pmc_alloc_mem( base_array_cp, bufsize, base_ptr )
	572	children(childid)%totalbuffersize = bufsize * wp
	573
	574	CALL MPI_BARRIER( children(childid)%intra_comm, ierr )
	575	!
	576	!-- Second stride: set buffer pointer
	577	DO i = 1, children(childid)%inter_npes
	578	ape => children(childid)%pes(i)
	579	DO j = 1, ape%nr_arrays
	580	ar => ape%array_list(j)
	581	ar%recvbuf = base_ptr
	582	ENDDO
	583	ENDDO
	584
	585	END SUBROUTINE pmc_s_setind_and_allocmem
	586
	587
	588
	589	SUBROUTINE pmc_s_fillbuffer( childid, waittime, particle_transfer )
	590
	591	IMPLICIT NONE
	592
	593	INTEGER(iwp), INTENT(IN) :: childid !<
	594
	595	REAL(wp), INTENT(OUT), OPTIONAL :: waittime !<
	596	LOGICAL, INTENT(IN), OPTIONAL :: particle_transfer !<
	597
	598
	599	INTEGER(iwp) :: ierr !<
	600	INTEGER(iwp) :: ij !<
	601	INTEGER(iwp) :: ip !<
	602	INTEGER(iwp) :: j !<
	603	INTEGER(iwp) :: myindex !<
	604
	605	LOGICAL :: lo_ptrans
	606
	607	INTEGER(iwp), DIMENSION(1) :: buf_shape
	608
	609	REAL(wp) :: t1 !<
	610	REAL(wp) :: t2 !<
	611	REAL(wp), POINTER, DIMENSION(:) :: buf !<
	612	REAL(wp), POINTER, DIMENSION(:,:) :: data_2d !<
	613	REAL(wp), POINTER, DIMENSION(:,:,:) :: data_3d !<
	614	INTEGER(idp), POINTER, DIMENSION(:) :: ibuf !<
	615	INTEGER(idp), POINTER, DIMENSION(:,:) :: idata_2d !<
	616
	617	TYPE(pedef), POINTER :: ape !<
	618	TYPE(arraydef), POINTER :: ar !<
	619
	620	!
	621	!-- Synchronization of the model is done in pmci_synchronize.
	622	!-- Therefor the RMA window can be filled without
	623	!-- sychronization at this point and a barrier is not necessary.
	624	!-- Please note that waittime has to be set in pmc_s_fillbuffer AND
	625	!-- pmc_c_getbuffer
	626	IF ( PRESENT( waittime) ) THEN
	627	t1 = pmc_time()
	628	CALL MPI_BARRIER( children(childid)%intra_comm, ierr )
	629	t2 = pmc_time()
	630	waittime = t2- t1
	631	ENDIF
	632
	633	lo_ptrans = .FALSE.
	634	IF ( PRESENT( particle_transfer)) lo_ptrans = particle_transfer
	635
	636	DO ip = 1, children(childid)%inter_npes
	637	ape => children(childid)%pes(ip)
	638	DO j = 1, ape%nr_arrays
	639	ar => ape%array_list(j)
	640	myindex = 1
	641
	642	IF ( ar%dimkey == 2 .AND. .NOT.lo_ptrans ) THEN ! PALM 2D REAL*8 Array
	643
	644	buf_shape(1) = ape%nrele
	645	CALL C_F_POINTER( ar%sendbuf, buf, buf_shape )
	646	CALL C_F_POINTER( ar%data, data_2d, ar%a_dim(1:2) )
	647	DO ij = 1, ape%nrele
	648	buf(myindex) = data_2d(ape%locind(ij)%j,ape%locind(ij)%i)
	649	myindex = myindex + 1
	650	ENDDO
	651
	652	ELSEIF ( ar%dimkey == 3 .AND. .NOT.lo_ptrans ) THEN ! PALM 3D REAL*8 Array
	653
	654	buf_shape(1) = ape%nrele*ar%a_dim(4)
	655	CALL C_F_POINTER( ar%sendbuf, buf, buf_shape )
	656	CALL C_F_POINTER( ar%data, data_3d, ar%a_dim(1:3) )
	657	DO ij = 1, ape%nrele
	658	buf(myindex:myindex+ar%a_dim(4)-1) = &
	659	data_3d(1:ar%a_dim(4),ape%locind(ij)%j,ape%locind(ij)%i)
	660	myindex = myindex + ar%a_dim(4)
	661	ENDDO
	662	ELSEIF ( ar%dimkey == 22 .AND. lo_ptrans ) THEN ! 2D INTEGER*8 Array for particle Transfer
	663
	664	buf_shape(1) = ape%nrele
	665	CALL C_F_POINTER( ar%sendbuf, ibuf, buf_shape )
	666	CALL C_F_POINTER( ar%data, idata_2d, ar%a_dim(1:2) )
	667	DO ij = 1, ape%nrele
	668	ibuf(myindex) = idata_2d(ape%locind(ij)%j,ape%locind(ij)%i)
	669	myindex = myindex + 1
	670	ENDDO
	671	ENDIF
	672	ENDDO
	673	ENDDO
	674	!
	675	!-- Buffer is filled
	676	CALL MPI_BARRIER( children(childid)%intra_comm, ierr )
	677
	678	END SUBROUTINE pmc_s_fillbuffer
	679
	680
	681
	682	SUBROUTINE pmc_s_getdata_from_buffer( childid, waittime , particle_transfer, child_process_nr )
	683
	684	IMPLICIT NONE
	685
	686	INTEGER(iwp), INTENT(IN) :: childid !<
	687	REAL(wp), INTENT(OUT), OPTIONAL :: waittime !<
	688	LOGICAL, INTENT(IN), OPTIONAL :: particle_transfer !<
	689	INTEGER(iwp), INTENT(IN), OPTIONAL :: child_process_nr !<
	690
	691	INTEGER(iwp) :: ierr !<
	692	INTEGER(iwp) :: ij !<
	693	INTEGER(iwp) :: ip !<
	694	INTEGER(iwp) :: ip_start !<
	695	INTEGER(iwp) :: ip_end !<
	696	INTEGER(iwp) :: j !<
	697	INTEGER(iwp) :: myindex !<
	698	INTEGER(iwp) :: nr !<
	699	INTEGER(iwp) :: target_pe !<
	700	INTEGER(KIND=MPI_ADDRESS_KIND) :: target_disp !<
	701
	702	LOGICAL :: lo_ptrans
	703
	704	INTEGER(iwp), DIMENSION(1) :: buf_shape !<
	705
	706	REAL(wp) :: t1 !<
	707	REAL(wp) :: t2 !<
	708	REAL(wp), POINTER, DIMENSION(:) :: buf !<
	709	REAL(wp), POINTER, DIMENSION(:,:) :: data_2d !<
	710	REAL(wp), POINTER, DIMENSION(:,:,:) :: data_3d !<
	711	INTEGER(idp), POINTER, DIMENSION(:) :: ibuf !<
	712	INTEGER(idp), POINTER, DIMENSION(:,:) :: idata_2d !<
	713
	714	TYPE(pedef), POINTER :: ape !<
	715	TYPE(arraydef), POINTER :: ar !<
	716
	717
	718	t1 = pmc_time()
	719
	720	IF(PRESENT(child_process_nr)) then
	721	ip_start = child_process_nr
	722	ip_end = child_process_nr
	723	ELSE
	724	ip_start = 1
	725	ip_end = children(childid)%inter_npes
	726	END IF
	727
	728	lo_ptrans = .FALSE.
	729	IF ( PRESENT( particle_transfer)) lo_ptrans = particle_transfer
	730
	731	IF(ip_start == 1) THEN
	732	!
	733	!-- Wait for child to fill buffer
	734	CALL MPI_BARRIER( children(childid)%intra_comm, ierr )
	735	t2 = pmc_time() - t1
	736	IF ( PRESENT( waittime ) ) waittime = t2
	737
	738	CALL MPI_BARRIER( children(childid)%intra_comm, ierr )
	739	ENDIF
	740
	741	DO ip = ip_start,ip_end
	742	ape => children(childid)%pes(ip)
	743	DO j = 1, ape%nr_arrays
	744	ar => ape%array_list(j)
	745
	746	IF ( ar%recvindex < 0 ) CYCLE
	747
	748	IF ( ar%dimkey == 2 .AND. .NOT.lo_ptrans ) THEN
	749	nr = ape%nrele
	750	ELSEIF ( ar%dimkey == 3 .AND. .NOT.lo_ptrans ) THEN
	751	nr = ape%nrele * ar%a_dim(4)
	752	ELSE IF ( ar%dimkey == 22 .AND. lo_ptrans) THEN
	753	nr = ape%nrele
	754	ELSE
	755	CYCLE !particle array are not transfered here
	756	ENDIF
	757	buf_shape(1) = nr
	758	IF(lo_ptrans) THEN
	759	CALL C_F_POINTER( ar%recvbuf, ibuf, buf_shape )
	760	ELSE
	761	CALL C_F_POINTER( ar%recvbuf, buf, buf_shape )
	762	ENDIF
	763
	764	!
	765	!-- MPI passive target RMA
	766	IF ( nr > 0 ) THEN
	767	target_disp = ar%recvindex - 1
	768	!
	769	!-- Child processes are located behind parent process
	770	target_pe = ip - 1 + m_model_npes
	771	CALL MPI_WIN_LOCK( MPI_LOCK_SHARED, target_pe, 0, &
	772	children(childid)%win_parent_child, ierr )
	773	IF(lo_ptrans) THEN
	774	CALL MPI_GET( ibuf, nr8, MPI_BYTE, target_pe, target_disp, nr8, & !There is no MPI_INTEGER8 datatype
	775	MPI_BYTE, children(childid)%win_parent_child, ierr )
	776	ELSE
	777	CALL MPI_GET( buf, nr, MPI_REAL, target_pe, target_disp, nr, &
	778	MPI_REAL, children(childid)%win_parent_child, ierr )
	779	ENDIF
	780	CALL MPI_WIN_UNLOCK( target_pe, &
	781	children(childid)%win_parent_child, ierr )
	782	ENDIF
	783	myindex = 1
	784	IF ( ar%dimkey == 2 .AND. .NOT.lo_ptrans ) THEN
	785
	786	CALL C_F_POINTER( ar%data, data_2d, ar%a_dim(1:2) )
	787	DO ij = 1, ape%nrele
	788	data_2d(ape%locind(ij)%j,ape%locind(ij)%i) = buf(myindex)
	789	myindex = myindex + 1
	790	ENDDO
	791
	792	ELSEIF ( ar%dimkey == 3 .AND. .NOT.lo_ptrans ) THEN
	793
	794	CALL C_F_POINTER( ar%data, data_3d, ar%a_dim(1:3))
	795	DO ij = 1, ape%nrele
	796	data_3d(1:ar%a_dim(4),ape%locind(ij)%j,ape%locind(ij)%i) = &
	797	buf(myindex:myindex+ar%a_dim(4)-1)
	798	myindex = myindex + ar%a_dim(4)
	799	ENDDO
	800
	801	ELSE IF ( ar%dimkey == 22 .AND. lo_ptrans) THEN
	802
	803	CALL C_F_POINTER( ar%data, idata_2d, ar%a_dim(1:2) )
	804	DO ij = 1, ape%nrele
	805	idata_2d(ape%locind(ij)%j,ape%locind(ij)%i) = ibuf(myindex)
	806	myindex = myindex + 1
	807	ENDDO
	808
	809	ENDIF
	810	ENDDO
	811	ENDDO
	812
	813	END SUBROUTINE pmc_s_getdata_from_buffer
	814
	815
	816
	817	SUBROUTINE get_da_names_from_child( childid )
	818
	819	!
	820	!-- Get data array description and name from child
	821	IMPLICIT NONE
	822
	823	INTEGER(iwp), INTENT(IN) :: childid !<
	824
	825	TYPE(da_namedef) :: myname !<
	826
	827	DO
	828	CALL pmc_bcast( myname%couple_index, 0, comm=m_to_child_comm(childid) )
	829
	830	IF ( myname%couple_index == -1 ) EXIT
	831
	832	CALL pmc_bcast( myname%parentdesc, 0, comm=m_to_child_comm(childid) )
	833	CALL pmc_bcast( myname%nameonparent, 0, comm=m_to_child_comm(childid) )
	834	CALL pmc_bcast( myname%childdesc, 0, comm=m_to_child_comm(childid) )
	835	CALL pmc_bcast( myname%nameonchild, 0, comm=m_to_child_comm(childid) )
	836
	837	CALL pmc_g_setname( children(childid), myname%couple_index, &
	838	myname%nameonparent )
	839	ENDDO
	840
	841	END SUBROUTINE get_da_names_from_child
	842
	843
	844
	845	SUBROUTINE pmc_s_setarray(childid, nrdims, dims, array_adr, second_adr, dimkey )
	846
	847	!
	848	!-- Set array for child inter process 0
	849	IMPLICIT NONE
	850
	851	INTEGER(iwp), INTENT(IN) :: childid !<
	852	INTEGER(iwp), INTENT(IN) :: nrdims !<
	853	INTEGER(iwp), INTENT(IN), DIMENSION(:) :: dims !<
	854
	855	TYPE(C_PTR), INTENT(IN) :: array_adr !<
	856	TYPE(C_PTR), INTENT(IN), OPTIONAL :: second_adr !<
	857	INTEGER(iwp), INTENT(IN), OPTIONAL :: dimkey !<
	858
	859	INTEGER(iwp) :: i !< local counter
	860
	861	TYPE(pedef), POINTER :: ape !<
	862	TYPE(arraydef), POINTER :: ar !<
	863
	864
	865	DO i = 1, children(childid)%inter_npes
	866	ape => children(childid)%pes(i)
	867	ar => ape%array_list(next_array_in_list)
	868	ar%nrdims = nrdims
	869	ar%dimkey = nrdims
	870	IF(PRESENT(dimkey)) ar%dimkey = dimkey
	871	ar%a_dim = dims
	872	ar%data = array_adr
	873	IF ( PRESENT( second_adr ) ) THEN
	874	ar%po_data(1) = array_adr
	875	ar%po_data(2) = second_adr
	876	ELSE
	877	ar%po_data(1) = C_NULL_PTR
	878	ar%po_data(2) = C_NULL_PTR
	879	ENDIF
	880	ENDDO
	881
	882	END SUBROUTINE pmc_s_setarray
	883
	884
	885
	886	SUBROUTINE pmc_s_set_active_data_array( childid, iactive )
	887
	888	IMPLICIT NONE
	889
	890	INTEGER(iwp), INTENT(IN) :: childid !<
	891	INTEGER(iwp), INTENT(IN) :: iactive !<
	892
	893	INTEGER(iwp) :: ip !<
	894	INTEGER(iwp) :: j !<
	895
	896	TYPE(pedef), POINTER :: ape !<
	897	TYPE(arraydef), POINTER :: ar !<
	898
	899	DO ip = 1, children(childid)%inter_npes
	900	ape => children(childid)%pes(ip)
	901	DO j = 1, ape%nr_arrays
	902	ar => ape%array_list(j)
	903	if(mod(ar%dimkey,10) == 2) CYCLE !Not for 2D array
	904	IF ( iactive == 1 .OR. iactive == 2 ) THEN
	905	ar%data = ar%po_data(iactive)
	906	ENDIF
	907	ENDDO
	908	ENDDO
	909
	910	END SUBROUTINE pmc_s_set_active_data_array
	911
	912	INTEGER FUNCTION pmc_s_get_child_npes (child_id)
	913	IMPLICIT NONE
	914
	915	INTEGER(iwp),INTENT(IN) :: child_id
	916
	917	pmc_s_get_child_npes = children(child_id)%inter_npes
	918
	919	RETURN
	920	END FUNCTION pmc_s_get_child_npes
	921
	922
	923	SUBROUTINE set_pe_index_list( mychild, index_list, nrp )
	924
	925	IMPLICIT NONE
	926
	927	INTEGER(iwp), INTENT(IN), DIMENSION(:,:) :: index_list !<
	928	INTEGER(iwp), INTENT(IN) :: nrp !<
	929
	930	TYPE(childdef), INTENT(INOUT) :: mychild !<
	931
	932	INTEGER(iwp) :: i !<
	933	INTEGER(iwp) :: ierr !<
	934	INTEGER(iwp) :: ind !<
	935	INTEGER(iwp) :: indwin !<
	936	INTEGER(iwp) :: indwin2 !<
	937	INTEGER(iwp) :: i2 !<
	938	INTEGER(iwp) :: j !<
	939	INTEGER(iwp) :: rempe !<
	940	INTEGER(KIND=MPI_ADDRESS_KIND) :: winsize !<
	941
	942	INTEGER(iwp), DIMENSION(mychild%inter_npes) :: remind !<
	943
	944	INTEGER(iwp), DIMENSION(:), POINTER :: remindw !<
	945	INTEGER(iwp), DIMENSION(:), POINTER :: rldef !<
	946
	947	TYPE(pedef), POINTER :: ape !<
	948
	949	!
	950	!-- First, count entries for every remote child process
	951	DO i = 1, mychild%inter_npes
	952	ape => mychild%pes(i)
	953	ape%nrele = 0
	954	ENDDO
	955	!
	956	!-- Loop over number of coarse grid cells
	957	DO j = 1, nrp
	958	rempe = index_list(5,j) + 1 ! process number on remote process
	959	ape => mychild%pes(rempe)
	960	ape%nrele = ape%nrele + 1 ! Increment number of elements for this child process
	961	ENDDO
	962
	963	DO i = 1, mychild%inter_npes
	964	ape => mychild%pes(i)
	965	ALLOCATE( ape%locind(ape%nrele) )
	966	ENDDO
	967
	968	remind = 0
	969	!
	970	!-- Second, create lists
	971	!-- Loop over number of coarse grid cells
	972	DO j = 1, nrp
	973	rempe = index_list(5,j) + 1
	974	ape => mychild%pes(rempe)
	975	remind(rempe) = remind(rempe)+1
	976	ind = remind(rempe)
	977	ape%locind(ind)%i = index_list(1,j)
	978	ape%locind(ind)%j = index_list(2,j)
	979	ENDDO
	980	!
	981	!-- Prepare number of elements for children processes
	982	CALL pmc_alloc_mem( rldef, mychild%inter_npes*2 )
	983	!
	984	!-- Number of child processes * size of INTEGER (i just arbitrary INTEGER)
	985	winsize = mychild%inter_npesSTORAGE_SIZE(i)/82
	986
	987	CALL MPI_WIN_CREATE( rldef, winsize, iwp, MPI_INFO_NULL, &
	988	mychild%intra_comm, indwin, ierr )
	989	!
	990	!-- Open window to set data
	991	CALL MPI_WIN_FENCE( 0, indwin, ierr )
	992
	993	rldef(1) = 0 ! index on remote process 0
	994	rldef(2) = remind(1) ! number of elements on remote process 0
	995	!
	996	!-- Reserve buffer for index array
	997	DO i = 2, mychild%inter_npes
	998	i2 = (i-1) * 2 + 1
	999	rldef(i2) = rldef(i2-2) + rldef(i2-1) * 2 ! index on remote process
	1000	rldef(i2+1) = remind(i) ! number of elements on remote process
	1001	ENDDO
	1002	!
	1003	!-- Close window to allow child to access data
	1004	CALL MPI_WIN_FENCE( 0, indwin, ierr )
	1005	!
	1006	!-- Child has retrieved data
	1007	CALL MPI_WIN_FENCE( 0, indwin, ierr )
	1008
	1009	i2 = 2 * mychild%inter_npes - 1
	1010	winsize = ( rldef(i2) + rldef(i2+1) ) * 2
	1011	!
	1012	!-- Make sure, MPI_ALLOC_MEM works
	1013	winsize = MAX( winsize, INT( 1, MPI_ADDRESS_KIND ) )
	1014
	1015	CALL pmc_alloc_mem( remindw, INT( winsize ) )
	1016
	1017	CALL MPI_BARRIER( m_model_comm, ierr )
	1018	CALL MPI_WIN_CREATE( remindw, winsize*STORAGE_SIZE(i)/8, iwp, MPI_INFO_NULL, &
	1019	mychild%intra_comm, indwin2, ierr )
	1020	!
	1021	!-- Open window to set data
	1022	CALL MPI_WIN_FENCE( 0, indwin2, ierr )
	1023	!
	1024	!-- Create the 2D index list
	1025	DO j = 1, nrp
	1026	rempe = index_list(5,j) + 1 ! process number on remote process
	1027	ape => mychild%pes(rempe)
	1028	i2 = rempe * 2 - 1
	1029	ind = rldef(i2) + 1
	1030	remindw(ind) = index_list(3,j)
	1031	remindw(ind+1) = index_list(4,j)
	1032	rldef(i2) = rldef(i2)+2
	1033	ENDDO
	1034	!
	1035	!-- All data are set
	1036	CALL MPI_WIN_FENCE( 0, indwin2, ierr )
	1037	!
	1038	!-- Don't know why, but this barrier is necessary before windows can be freed
	1039	!-- TODO: find out why this is required
	1040	CALL MPI_BARRIER( mychild%intra_comm, ierr )
	1041
	1042	CALL MPI_WIN_FREE( indwin, ierr )
	1043	CALL MPI_WIN_FREE( indwin2, ierr )
	1044
	1045	!
	1046	!-- TODO: check if the following idea needs to be done
	1047	!-- Sollte funktionieren, Problem mit MPI implementation
	1048	!-- https://www.lrz.de/services/software/parallel/mpi/onesided
	1049	!-- CALL MPI_Free_mem (remindw, ierr)
	1050
	1051	END SUBROUTINE set_pe_index_list
	1052
	1053	#endif
	1054	END MODULE pmc_parent

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |