Index: convect.f
===================================================================
--- convect.f	(revision 954)
+++ convect.f	(working copy)
@@ -191,6 +191,7 @@
       end
 c-----------------------------------------------------------------------
       subroutine int_vel(c_t,t0,c,n,nc,ct,nid)
+cDEC$ ATTRIBUTES FORCEINLINE :: int_vel
 
 c     Interpolate convecting velocity field c(t_1,...,t_nconv) to
 c     time t0 and return result in c_t.
@@ -204,15 +205,17 @@
       parameter (lwtmax=10)
       real wt(0:lwtmax)
 c
-      if (nc.gt.lwtmax) then
-         write(6,*) nid,'ERROR int_vel: lwtmax too small',lwtmax,m0
-         call exitt
-      endif
+c      if (nc.gt.lwtmax) then
+c         write(6,*) nid,'ERROR int_vel: lwtmax too small',lwtmax,m0
+c         call exitt
+c      endif
 c
       no = nc-1
       call fd_weights_full(t0,ct(0),no,0,wt)  ! interpolation weights
 
-      call rzero(c_t,n)
+c      call rzero(c_t,n)
+      c_t(:)=0.0
+
       do j=1,n
       do i=0,no
          c_t(j) = c_t(j) + wt(i)*c(j,i)