Changeset 5f9d14a in flexpart.git for src/initial_cond_calc.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/initial_cond_calc.f90 (modified) (2 diffs)

src/initial_cond_calc.f90

@@ -44 +44 @@
   real :: ddx,ddy
   real :: rhoprof(2),rhoi,xl,yl,wx,wy,w
+  integer :: mind2
+  ! mind2        eso: pointer to 2nd windfield in memory
 
 
@@ -88 +90 @@
   ! Take density from 2nd wind field in memory (accurate enough, no time interpolation needed)
   !*****************************************************************************
+    mind2=memind(2)
+
     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,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