Changeset d005a67 in flexpart.git for src/timemanager_mpi.f90
- Timestamp:
- May 21, 2019, 1:09:18 PM (5 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug
- Children:
- f963113, 0a98afe
- Parents:
- 5d74ed9 (diff), 0c8c7f2 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent. - File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/timemanager_mpi.f90
r20963b1 r0c8c7f2 113 113 integer :: j,ks,kp,l,n,itime=0,nstop,nstop1,memstat=0 114 114 ! integer :: ksp 115 integer :: ip 115 integer :: ip,irec 116 116 integer :: loutnext,loutstart,loutend 117 117 integer :: ix,jy,ldeltat,itage,nage,idummy … … 129 129 ! Measure time spent in timemanager 130 130 if (mp_measure_time) call mpif_mtime('timemanager',0) 131 131 132 132 133 ! First output for time 0 … … 532 533 end if 533 534 534 else ! :TODO: check for zeroing in the netcdf module535 else 535 536 call concoutput_surf_nest(itime,outnum) 536 537 end if … … 593 594 46 format(' Simulated ',f7.1,' hours (',i13,' s), ',i13, ' particles') 594 595 if (ipout.ge.1) then 596 if (mp_measure_time) call mpif_mtime('iotime',0) 597 irec=0 595 598 do ip=0, mp_partgroup_np-1 596 if (ip.eq.mp_partid) call partoutput(itime) ! dump particle positions 599 if (ip.eq.mp_partid) then 600 if (mod(itime,ipoutfac*loutstep).eq.0) call partoutput(itime) ! dump particle positions 601 if (ipout.eq.3) call partoutput_average(itime,irec) ! dump particle positions 602 endif 603 if (ipout.eq.3) irec=irec+npart_per_process(ip) 597 604 call mpif_mpi_barrier 598 605 end do 606 if (mp_measure_time) call mpif_mtime('iotime',1) 599 607 end if 600 608 … … 757 765 if (mp_measure_time) call mpif_mtime('advance',1) 758 766 767 ! Calculate average position for particle dump output 768 !**************************************************** 769 770 if (ipout.eq.3) call partpos_average(itime,j) 771 759 772 760 773 ! Calculate the gross fluxes across layer interfaces … … 895 908 do ip=0, mp_partgroup_np-1 896 909 if (ip.eq.mp_partid) then 897 !if (mp_dbg_mode) write(*,*) 'call partoutput(itime), proc, mp_partid',ip,mp_partid910 if (mp_dbg_mode) write(*,*) 'call partoutput(itime), proc, mp_partid',ip,mp_partid 898 911 call partoutput(itime) ! dump particle positions 899 912 end if
Note: See TracChangeset
for help on using the changeset viewer.