Changeset 6a678e3 in flexpart.git for src/timemanager.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.f90 (modified) (2 diffs)

src/timemanager.f90

@@ -111 +111 @@
   ! real :: us(maxpart),vs(maxpart),ws(maxpart)
   ! integer(kind=2) :: cbt(maxpart)
-  real :: drydeposit(maxspec),gridtotalunc,wetgridtotalunc
-  real :: drygridtotalunc,xold,yold,zold,xmassfract
+  real(sp) :: gridtotalunc
+  real(dep_prec) :: drydeposit(maxspec),wetgridtotalunc,drygridtotalunc
+  real :: xold,yold,zold,xmassfract
   real, parameter :: e_inv = 1.0/exp(1.0)
   !double precision xm(maxspec,maxpointspec_act),
@@ -595 +596 @@
 !ZHG 2015
                   !CGZ-lifetime: Check mass fraction left/save lifetime
-                   ! if(real(npart(npoint(j)))*xmass1(j,ks)/xmass(npoint(j),ks).lt.inv_e.and.checklifetime(j,ks).eq.0.)then
+                   ! if(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))