Changeset f251e57 in flexpart.git for src/verttransform_ecmwf.f90

Timestamp:

Jun 25, 2018, 7:44:14 PM (6 years ago)

Author:

pesei <petra seibert at univie ac at>

Branches:

univie

Children:

c0884a8

Parents:

8b3d324

Message:

fix ticket:140 (ref position for vert. grid) in verttransform_gfs, minor changes in both verttransforms

File:

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

src/verttransform_ecmwf.f90

@@ -157 +157 @@
 
     if (ioldref .eq. 0) then ! old reference grid point
-       do jy=0,nymin1
-         do ix=0,nxmin1
-           if (ps(ix,jy,1,n).gt.1000.e2) then
-             ixref=ix
-             jyref=jy
-             goto 3
-           endif
-         end do
-       end do
+      do jy=0,nymin1
+        do ix=0,nxmin1
+          if (ps(ix,jy,1,n).gt.1000.e2) then
+            ixref=ix
+            jyref=jy
+            goto 3
+          endif
+        end do
+      end do
 3     continue
 
@@ -175 +175 @@
 
       psmean = sum( ps(:,:,1,n) ) / (nx*ny)
-      print*,'height: fg psmean',psmean
+!      print*,'height: fg psmean',psmean
       psstd = sqrt(sum( (ps(:,:,1,n) - psmean)**2 ) / (nx*ny))
 
@@ -201 +201 @@
 
       pintref = akz(kz)+bkz(kz)*ps(ixref,jyref,1,n)
-      tv = tth(ixref,jyref,kz,n)*(1.+0.608*qvh(ixref,jyref,kz,n))
+      tv(ixref,jyref) = tth(ixref,jyref,kz,n)*(1.+0.608*qvh(ixref,jyref,kz,n))
 
       if (abs(tv(ixref,jyref)-tvoldref) .gt. 0.2) then
@@ -214 +214 @@
       tvoldref = tv(ixref,jyref)
       poldref = pintref
-      print*,'height=',kz,height(kz),tvoldref,poldref
+!      print*,'height=',kz,height(kz),tvoldref,poldref
 
     end do