Changeset 57c5b2e in flexpart.git
- Timestamp:
- Sep 13, 2015, 7:54:53 PM (9 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- cbad0f1
- Parents:
- f55fdce
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/releaseparticles_mpi.f90
r5f9d14a r57c5b2e 55 55 use com_mod 56 56 use random_mod, only: ran1 57 use mpi_mod, only: mp_partid, maxpart_mpi, mp_partgroup_np, mp_seed 57 use mpi_mod, only: mp_partid, maxpart_mpi, mp_partgroup_np, mp_seed, mpif_mpi_barrier 58 58 59 59 implicit none … … 78 78 logical :: first_call=.true. 79 79 80 ! Use different seed for each process 80 ! Use different seed for each process. 81 !**************************************************************************** 81 82 if (first_call) then 82 83 idummy=idummy+mp_seed … … 152 153 !********************************************************************** 153 154 155 rfraction=rfraction*average_timecorrect 156 157 rfraction=rfraction+xmasssave(i) ! number to be released at this time 158 159 ! number to be released for one process 160 if (mp_partid.lt.mod(int(rfraction),mp_partgroup_np)) then 161 addone=1 162 else 163 addone=0 164 end if 165 166 numrel=int(rfraction/mp_partgroup_np) + addone 167 168 xmasssave(i)=rfraction-int(rfraction) 169 170 else 171 ! All particles are released in this time interval 172 ! ************************************************ 154 173 if (mp_partid.lt.mod(npart(i),mp_partgroup_np)) then 155 174 addone=1 … … 158 177 end if 159 178 160 rfraction=rfraction*average_timecorrect/real(mp_partgroup_np)161 162 rfraction=rfraction+xmasssave(i) ! number to be released at this time163 numrel=int(rfraction)+addone164 xmasssave(i)=rfraction-real(numrel)165 else166 if (mp_partid.lt.mod(npart(i),mp_partgroup_np)) then167 addone=1168 else169 addone=0170 end if171 172 179 numrel=npart(i)/mp_partgroup_np+addone 173 180 endif 174 181 175 182 xaux=xpoint2(i)-xpoint1(i) 176 183 yaux=ypoint2(i)-ypoint1(i) … … 208 215 xmass1(ipart,k)=xmass(i,k)/real(npart(i)) & 209 216 *timecorrect(k)/average_timecorrect 210 ! write (*,*) 'xmass1: ',xmass1(ipart,k),ipart,k217 ! write (*,*) 'xmass1: ',xmass1(ipart,k),ipart,k 211 218 ! Assign certain properties to particle 212 219 !**************************************
Note: See TracChangeset
for help on using the changeset viewer.