--- MITgcm_contrib/plumes/plume2dyn.F	2004/05/25 18:08:51	1.2
+++ MITgcm_contrib/plumes/plume2dyn.F	2004/05/25 23:12:50	1.3
@@ -76,12 +76,16 @@
        do j = jdim1,jdim2
        do i = idim1,idim2
        do L = 1,Lmplume
-        sqrtarg = (qplumeav(i,j,L)*qplumeav(i,j,L)) /
+        if(vref(i,j,L,bi,bj).ne.0.) then
+         sqrtarg = (qplumeav(i,j,L)*qplumeav(i,j,L)) /
      .  ( ( (uref(i,j,L,bi,bj)*uref(i,j,L,bi,bj)) /
      .      (vref(i,j,L,bi,bj)*vref(i,j,L,bi,bj)) ) + 1. )
-        qdyn2(i,j,L,bi,bj) = sqrt(sqrtarg)
-        qdyn1(i,j,L,bi,bj) = qdyn2(i,j,L,bi,bj) *
+         qdyn2(i,j,L,bi,bj) = sqrt(sqrtarg)
+         qdyn1(i,j,L,bi,bj) = qdyn2(i,j,L,bi,bj) *
      .                            (uref(i,j,L,bi,bj)/vref(i,j,L,bi,bj))
+        else
+         qdyn1(i,j,L,bi,bj) = qplumeav(i,j,L)
+         qdyn2(i,j,L,bi,bj) = 0.
        enddo
        enddo
        enddo