Changeset 8a65cb0 in flexpart.git for src/interpol_rain.f90

Timestamp:

Mar 2, 2015, 3:11:55 PM (9 years ago)

Author:

Espen Sollum ATMOS <espen@…>

Branches:

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

Children:

1d207bb

Parents:

60403cd

Message:

Added code, makefile for dev branch

File:

• src/interpol_rain.f90 (modified) (2 diffs)

src/interpol_rain.f90

@@ -66 +66 @@
   !                                                                           *
   !****************************************************************************
+  use par_mod, only: numwfmem
 
   implicit none
 
-  integer :: nx,ny,nxmax,nymax,nzmax,memind(2),m,ix,jy,ixp,jyp
+  integer :: nx,ny,nxmax,nymax,nzmax,memind(numwfmem),m,ix,jy,ixp,jyp
   integer :: itime,itime1,itime2,level,indexh
-  real :: yy1(0:nxmax-1,0:nymax-1,nzmax,2)
-  real :: yy2(0:nxmax-1,0:nymax-1,nzmax,2)
-  real :: yy3(0:nxmax-1,0:nymax-1,nzmax,2)
+  real :: yy1(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
+  real :: yy2(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
+  real :: yy3(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
   real :: ddx,ddy,rddx,rddy,dt1,dt2,dt,y1(2),y2(2),y3(2)
   real :: xt,yt,yint1,yint2,yint3,p1,p2,p3,p4
@@ -115 +116 @@
     indexh=memind(m)
 
+
     y1(m)=p1*yy1(ix ,jy ,level,indexh) &
          + p2*yy1(ixp,jy ,level,indexh) &