Changeset 32b49c3 in flexpart.git for src/timemanager_mpi.f90

Timestamp:

Apr 13, 2016, 11:57:08 AM (8 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

341f4b7

Parents:

e31b3b5

Message:

Parallel version can now save/restart simulations with IPOUT/IPIN

File:

• src/timemanager_mpi.f90 (modified) (5 diffs)

src/timemanager_mpi.f90

@@ -104 +104 @@
 
   logical :: reqv_state=.false. ! .true. if waiting for a MPI_Irecv to complete
-  integer :: j,ks,kp,l,n,itime=0,nstop,nstop1,memstat=0,mind
+  integer :: j,ks,kp,l,n,itime=0,nstop,nstop1,memstat=0 !,mind
 ! integer :: ksp
   integer :: ip
@@ -115 +115 @@
 
   real(sp) :: gridtotalunc
-  real(dep_prec) :: drygridtotalunc,wetgridtotalunc,drydeposit(maxspec)
+  real(dep_prec) :: drygridtotalunc=0_dep_prec,wetgridtotalunc=0_dep_prec,&
+       & drydeposit(maxspec)=0_dep_prec
   real :: xold,yold,zold,xmassfract
   real, parameter :: e_inv = 1.0/exp(1.0)
@@ -137 +138 @@
 
 !  itime=0
-  if (lroot) then
-  !  write(*,45) itime,numpart*mp_partgroup_np,gridtotalunc,wetgridtotalunc,drygridtotalunc
-    write(*,46) float(itime)/3600,itime,numpart*mp_partgroup_np
+  if (lroot.or.mp_dev_mode) then
+    write(*,45) itime,numpart*mp_partgroup_np,gridtotalunc,wetgridtotalunc,drygridtotalunc
+  !  write(*,46) float(itime)/3600,itime,numpart*mp_partgroup_np
     
     if (verbosity.gt.0) then
@@ -275 +276 @@
 ! For validation and tests: call the function below to set all fields to simple
 ! homogeneous values
-    if (mp_dev_mode) call set_fields_synthetic
+!    if (mp_dev_mode) call set_fields_synthetic
 
 !*******************************************************************************
@@ -554 +555 @@
         
         !CGZ-lifetime: output species lifetime
-        if (lroot) then
+        if (lroot.or.mp_dev_mode) then
         !   write(*,*) 'Overview species lifetime in days', &
         !        real((species_lifetime(:,1)/species_lifetime(:,2))/real(3600.0*24.0))