Changeset c7e771d in flexpart.git for src/get_wetscav.f90

Timestamp:

Mar 1, 2017, 1:06:56 PM (7 years ago)

Author:

Sabine <sabine.eckhardt@…>

Branches:

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

Children:

5f53d0e

Parents:

b85b020

Message:

interpolate the rain only in space and not in time

File:

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

src/get_wetscav.f90

@@ -123 +123 @@
     endif
 
-
 ! Interpolate large scale precipitation, convective precipitation and
 ! total cloud cover
 ! Note that interpolated time refers to itime-0.5*ltsample [PS]
 !********************************************************************
-    interp_time=nint(itime-0.5*ltsample)
+    interp_time=nint(itime-0.5*ltsample) 
+
+    n=memind(2)
+    if (abs(memtime(1)-interp_time).lt.abs(memtime(2)-interp_time)) &
+         n=memind(1)
 
     if (ngrid.eq.0) then
       call interpol_rain(lsprec,convprec,tcc,nxmax,nymax, &
-           1,nx,ny,memind,real(xtra1(jpart)),real(ytra1(jpart)),1, &
+           1,nx,ny,n,real(xtra1(jpart)),real(ytra1(jpart)),1, &
            memtime(1),memtime(2),interp_time,lsp,convp,cc)
     else
       call interpol_rain_nests(lsprecn,convprecn,tccn, &
-           nxmaxn,nymaxn,1,maxnests,ngrid,nxn,nyn,memind,xtn,ytn,1, &
+           nxmaxn,nymaxn,1,maxnests,ngrid,nxn,nyn,n,xtn,ytn,1, &
            memtime(1),memtime(2),interp_time,lsp,convp,cc)
     endif
@@ -153 +156 @@
 !26  continue
 
-    n=memind(2)
-    if (abs(memtime(1)-interp_time).lt.abs(memtime(2)-interp_time)) &
-         n=memind(1)
 
     if (ngrid.eq.0) then
@@ -227 +227 @@
         act_temp=tt(ix,jy,hz,n)
       endif
-
 
 !***********************
@@ -339 +338 @@
       endif !incloud
 
+
 20  continue