From fbb2555eae396344c9eeb9b4408e868f14d6c5e9 Mon Sep 17 00:00:00 2001
From: Bill Ladwig <ladwig@ucar.edu>
Date: Thu, 16 Nov 2017 16:22:25 -0700
Subject: [PATCH] Performance improvements

---
 fortran/wrf_vinterp.f90 | 11 +----------
 1 file changed, 1 insertion(+), 10 deletions(-)

diff --git a/fortran/wrf_vinterp.f90 b/fortran/wrf_vinterp.f90
index 7f7c2f1..bbcf168 100644
--- a/fortran/wrf_vinterp.f90
+++ b/fortran/wrf_vinterp.f90
@@ -184,14 +184,6 @@ SUBROUTINE wrf_vintrp(datain, dataout, pres, tk, qvp, ght, terrain,&
         cvcord = 't'
     END IF
 
-    !$OMP PARALLEL DO COLLAPSE(2)
-    DO j = 1,ns
-        DO i = 1,ew
-            tempout(i,j) = rmsg
-        END DO
-    END DO
-    !$OMP END PARALLEL DO
-
     DO nreqlvs = 1,numlevels
         IF (cvcord .EQ. 'z') THEN
             ! Convert rlevel to meters from km
@@ -211,6 +203,7 @@ SUBROUTINE wrf_vintrp(datain, dataout, pres, tk, qvp, ght, terrain,&
         !$OMP vt, tlev, gammam, e, tlcl) REDUCTION (+:log_errcnt, interp_errcnt)
         DO j=1,ns
             DO i=1,ew
+                tempout(i,j) = rmsg
                 ! Get the interpolated value that is within the model domain
                 ifound = 0
                 DO k = 1,nz-1
@@ -411,8 +404,6 @@ SUBROUTINE wrf_vintrp(datain, dataout, pres, tk, qvp, ght, terrain,&
         END DO
         !$OMP END PARALLEL DO
 
-
-
         IF (log_errcnt > 0) THEN
             errstat = ALGERR
             WRITE(errmsg, *) "Pres=0.  Unable to take log of 0."