Changeset 6a678e3 in flexpart.git for src/timemanager_mpi.f90

Timestamp:

Jan 25, 2016, 2:17:34 PM (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:

26f6039

Parents:

fddc6ec

Message:

Added option to use double precision for calculating the deposition fields

File:

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

src/timemanager_mpi.f90

@@ -114 +114 @@
   real :: decfact
 
-  real :: drydeposit(maxspec),gridtotalunc,wetgridtotalunc
-  real :: drygridtotalunc,xold,yold,zold,xmassfract
+  real(sp) :: gridtotalunc
+  real(dep_prec) :: drygridtotalunc,wetgridtotalunc,drydeposit(maxspec)
+  real :: xold,yold,zold,xmassfract
   real, parameter :: e_inv = 1.0/exp(1.0)
 !double precision xm(maxspec,maxpointspec_act),
@@ -261 +262 @@
 ! TODO: z0(7) changes with time, so should be dimension (numclass,2) to
 ! allow transfer of the future value in the background
-        call MPI_Bcast(z0,numclass,mp_pp,id_read,MPI_COMM_WORLD,mp_ierr)
+        call MPI_Bcast(z0,numclass,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)
         call mpif_gf_request
       end if
@@ -756 +757 @@
                    
                    !CGZ-lifetime: Check mass fraction left/save lifetime
-                   ! if(lroot.and.real(npart(npoint(j)))*xmass1(j,ks)/xmass(npoint(j),ks).lt.inv_e.and.checklifetime(j,ks).eq.0.)then
+                   ! if(lroot.and.real(npart(npoint(j)))*xmass1(j,ks)/xmass(npoint(j),ks).lt.e_inv.and.checklifetime(j,ks).eq.0.)then
                        !Mass below 1% of initial >register lifetime
                    !     checklifetime(j,ks)=abs(itra1(j)-itramem(j))