Changeset 5f9d14a in flexpart.git for src/conccalc_mpi.f90

Timestamp:

Apr 8, 2015, 2:23:27 PM (9 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:

1585284

Parents:

cd85138

Message:

Updated wet depo scheme

File:

• src/conccalc_mpi.f90 (modified) (3 diffs, 1 prop)

src/conccalc_mpi.f90

@@ -62 +62 @@
   real :: xl,yl,wx,wy,w
   real,parameter :: factor=.596831, hxmax=6.0, hymax=4.0, hzmax=150.
-
+  integer :: mind2
+  ! mind2        eso: pointer to 2nd windfield in memory
 
   ! For forward simulations, make a loop over the number of species;
@@ -70 +71 @@
 
   if (mp_measure_time) call mpif_mtime('conccalc',0)
+
+  mind2=memind(2)
 
   do i=1,numpart
@@ -131 +134 @@
   !*****************************************************************************
       do ind=indz,indzp
-        rhoprof(ind-indz+1)=p1*rho(ix ,jy ,ind,2) &
-             +p2*rho(ixp,jy ,ind,2) &
-             +p3*rho(ix ,jyp,ind,2) &
-             +p4*rho(ixp,jyp,ind,2)
+        ! rhoprof(ind-indz+1)=p1*rho(ix ,jy ,ind,2) &
+        !      +p2*rho(ixp,jy ,ind,2) &
+        !      +p3*rho(ix ,jyp,ind,2) &
+        !      +p4*rho(ixp,jyp,ind,2)
+        rhoprof(ind-indz+1)=p1*rho(ix ,jy ,ind,mind2) &
+             +p2*rho(ixp,jy ,ind,mind2) &
+             +p3*rho(ix ,jyp,ind,mind2) &
+             +p4*rho(ixp,jyp,ind,mind2)
+
       end do
       rhoi=(dz1*rhoprof(2)+dz2*rhoprof(1))*dz