Changeset 07c3e71 in flexpart.git for src/part0.f90

Timestamp:

Jul 24, 2019, 12:16:51 PM (5 years ago)

Author:

Ignacio Pisso <ip@…>

Branches:

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

Children:

25b4532

Parents:

95a45d3

Message:

check dsigma for dry deposition velocity

File:

• src/part0.f90 (modified) (4 diffs)

src/part0.f90

@@ -103 +103 @@
     x01=alog(d01/dquer)/xdummy
     x02=alog(d02/dquer)/xdummy
-
+    !print*,'part0:: d02=' , d02 , 'd01=', d01
 
   ! Area under Gauss-function is calculated and gives mass fraction of interval
@@ -109 +109 @@
 
     fract(i)=0.5*(erf(x01)-erf(x02))
-
+    !print*,'part0:: fract(',i,')', fract(i)
+    !print*,'part0:: fract', fract(i), x01, x02, erf(x01), erf(x02)
 
   ! Geometric mean diameter of interval in [m]
@@ -115 +116 @@
 
     dmean=1.E-6*exp(0.5*alog(d01*d02))
-
+    !print*,'part0:: dmean=', dmean
 
   ! Calculation of time independent parameters of each interval
@@ -132 +133 @@
     vsh(i)=ga*density*dmean*dmean*cun/(18.*myl)
 
+    !print*,'part0:: vsh(',i,')', vsh(i)
+
   end do
 
+  !stop 'part0' 
+
 end subroutine part0