Changeset db712a8 in flexpart.git for src/interpol_rain_nests.f90

Timestamp:

Mar 3, 2016, 12:34:56 PM (8 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:

38b7917

Parents:

b0434e1

Message:

Completed handling of nested wind fields with cloud water (for wet deposition).

File:

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

src/interpol_rain_nests.f90

@@ -88 +88 @@
   !***********************************************************
 
-  if (xt.ge.(real(nxn(ngrid)-1)-0.0001)) &
-       xt=real(nxn(ngrid)-1)-0.0001
-  if (yt.ge.(real(nyn(ngrid)-1)-0.0001)) &
-       yt=real(nyn(ngrid)-1)-0.0001
+  ! if (xt.ge.(real(nxn(ngrid)-1)-0.0001)) &
+  !      xt=real(nxn(ngrid)-1)-0.0001
+  ! if (yt.ge.(real(nyn(ngrid)-1)-0.0001)) &
+  !      yt=real(nyn(ngrid)-1)-0.0001
+
+! ESO make it consistent with interpol_rain
+  if (xt.ge.(real(nxn(ngrid)-1))) xt=real(nxn(ngrid)-1)-0.00001
+  if (yt.ge.(real(nyn(ngrid)-1))) yt=real(nyn(ngrid)-1)-0.00001
 
 
@@ -105 +109 @@
   ix=int(xt)
   jy=int(yt)
+
   ixp=ix+1
   jyp=jy+1