Changeset 72d3a5a in flexpart.git

Timestamp:

May 31, 2017, 10:11:57 AM (7 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:

e52967c

Parents:

6d420c3

Message:

Bugfix in verttransform (call to function ew)

Location:

src

Files:

• verttransform.f90 (modified) (3 diffs)
• verttransform_nests.f90 (modified) (2 diffs)

src/verttransform.f90

@@ -50 +50 @@
 ! added the variable cloud for use with scavenging - descr. in com_mod
 !*****************************************************************************
+! Date: 2017-05-30 modification of a bug in ew. Don Morton (CTBTO project)   *
+!*****************************************************************************
 !                                                                            *
 ! Variables:                                                                 *
@@ -155 +157 @@
 
 
-    tvold(ixm,jym)=tt2(ixm,jym,1,n)*(1.+0.378*ew*(td2(ixm,jym,1,n))/ &
+    tvold(ixm,jym)=tt2(ixm,jym,1,n)*(1.+0.378*ew(td2(ixm,jym,1,n))/ &
          ps(ixm,jym,1,n))
     pold(ixm,jym)=ps(ixm,jym,1,n)
@@ -215 +217 @@
   do jy=0,nymin1
     do ix=0,nxmin1
-      tvold(ix,jy)=tt2(ix,jy,1,n)*(1.+0.378*ew*(td2(ix,jy,1,n))/ &
+      tvold(ix,jy)=tt2(ix,jy,1,n)*(1.+0.378*ew(td2(ix,jy,1,n))/ &
            ps(ix,jy,1,n))
     enddo

src/verttransform_nests.f90

@@ -56 +56 @@
 !  the actual field dimensions
 !*****************************************************************************
+! Date: 2017-05-30 modification of a bug in ew. Don Morton (CTBTO project)   *
+!*****************************************************************************
 !                                                                            *
 ! Variables:                                                                 *
@@ -107 +109 @@
     do jy=0,nyn(l)-1
       do ix=0,nxn(l)-1
-        tvold(ix,jy)=tt2n(ix,jy,1,n,l)*(1.+0.378*ew*(td2n(ix,jy,1,n,l))/ &
+        tvold(ix,jy)=tt2n(ix,jy,1,n,l)*(1.+0.378*ew(td2n(ix,jy,1,n,l))/ &
              psn(ix,jy,1,n,l))
       end do