pass itime to init_domainfill to enable start which is not time 0

src/init_domainfill.f90

@@ -19,7 +19,7 @@
 ! along with FLEXPART.  If not, see <http://www.gnu.org/licenses/>.   *
 !**********************************************************************
 
-subroutine init_domainfill
+subroutine init_domainfill(itime)
 !
 !*****************************************************************************
 !                                                                            *
@@ -52,6 +52,8 @@ subroutine init_domainfill
 
   implicit none
 
+  integer :: itime
+
   integer :: j,ix,jy,kz,ncolumn,numparttot
   real :: gridarea(0:nymax-1),pp(nzmax),ylat,ylatp,ylatm,hzone
   real :: cosfactm,cosfactp,deltacol,dz1,dz2,dz,pnew,fractus
@@ -272,8 +274,12 @@ subroutine init_domainfill
                 numparticlecount=numparticlecount+1
                 npoint(numpart+jj)=numparticlecount
                 idt(numpart+jj)=mintime
-                itra1(numpart+jj)=0
-                itramem(numpart+jj)=0
+                if (WATERCYCLE) then
+                   itra1(numpart+jj)=itime
+                else
+                   itra1(numpart+jj)=0
+                endif
+                itramem(numpart+jj)=itime !sec, changed from 0 to itime for watercycle, is 0 when called in standard FP
                 itrasplit(numpart+jj)=itra1(numpart+jj)+ldirect* &
                      itsplit
                 xmass1(numpart+jj,1)=colmass(ix,jy)/real(ncolumn)
@@ -309,7 +315,7 @@ subroutine init_domainfill
 
 
   xmassperparticle=colmasstotal/real(numparttot)
-
+  write(*,*) 'Initdomainfill: ',numparttot,xmassperparticle
 
 ! Make sure that all particles are within domain
 !***********************************************