--- MITgcm_contrib/plumes/plume2dyn.F	2004/05/26 03:10:27	1.4
+++ MITgcm_contrib/plumes/plume2dyn.F	2004/06/01 22:31:34	1.6
@@ -46,7 +46,7 @@
       _RL qdyn2(idim1:idim2,jdim1:jdim2,Lmplume,Nsx,Nsy)
 
       integer i,j,L,iplume
-      _RL qplumeav(i1,j2,Lmplume)
+      _RL qplumeav(i2,j2,Lmplume)
       _RL sqrtarg
 
 C First step - compute the average of qplume over Nxplume
@@ -88,6 +88,7 @@
         else
          qdyn1(i,j,L,bi,bj) = qplumeav(i,j,L)
          qdyn2(i,j,L,bi,bj) = 0.
+        endif
        enddo
        enddo
        enddo