diff --git a/components/homme/src/share/cube_mod.F90 b/components/homme/src/share/cube_mod.F90
index 864c24c895d5..e0092c8945ae 100644
--- a/components/homme/src/share/cube_mod.F90
+++ b/components/homme/src/share/cube_mod.F90
@@ -2301,13 +2301,6 @@ subroutine CubeSetupEdgeIndex(Edge)
     integer,allocatable                :: backwardV(:), backwardP(:)
     integer                            :: i,ii
 
-    ii=Edge%tail_face
-
-    !map to correct location - for now all on same nbr side have same wgt, so take the first one
-    ii = Edge%tail%nbrs_ptr(ii)
-  
-    np0 = Edge%tail%nbrs_wgt(ii)
-
     sFace = Edge%tail_face
     dFace = Edge%head_face
     ! Do not reverse the indices
diff --git a/components/homme/src/share/edge_mod.F90 b/components/homme/src/share/edge_mod.F90
index 9ca007502dce..22f97450ac22 100644
--- a/components/homme/src/share/edge_mod.F90
+++ b/components/homme/src/share/edge_mod.F90
@@ -604,7 +604,6 @@ subroutine edgeVpack(edge,v,vlyr,kptr,ielem)
     integer :: i,k,ir,ll,iptr
 
     integer :: is,ie,in,iw
-    integer :: nce
 
     !call t_adj_detailf(+2)
     !call t_startf('edgeVpack')
@@ -688,12 +687,15 @@ subroutine edgeVpack(edge,v,vlyr,kptr,ielem)
     endif
 
     !set the max_corner_elem
-    nce = max_corner_elem
 ! SWEST
     do ll=swest,swest+max_corner_elem-1
         if (edge%putmap(ll,ielem) /= -1) then
             do k=1,vlyr
-                edge%buf(nce*(kptr+k-1)+edge%putmap(ll,ielem)+1)=v(1  ,1 ,k)
+                iptr = (kptr+k-1)+edge%putmap(ll,ielem)+1
+                if (iptr > size(edge%buf)) then
+                    print *,'ERROR 1: ',size(edge%buf),iptr,edge%putmap(ll,ielem)
+                endif
+                edge%buf(iptr)=v(1  ,1 ,k)
             end do
         end if
     end do
@@ -702,7 +704,11 @@ subroutine edgeVpack(edge,v,vlyr,kptr,ielem)
     do ll=swest+max_corner_elem,swest+2*max_corner_elem-1
         if (edge%putmap(ll,ielem) /= -1) then
             do k=1,vlyr
-                edge%buf(nce*(kptr+k-1)+edge%putmap(ll,ielem)+1)=v(np ,1 ,k)
+                iptr = (kptr+k-1)+edge%putmap(ll,ielem)+1
+                if (iptr > size(edge%buf)) then
+                    print *,'ERROR 2: ',size(edge%buf),iptr
+                endif
+                edge%buf(iptr)=v(np ,1 ,k)
 !                edge%buf(kptr+k,edge%putmap(ll,ielem)+1)=v(np ,1 ,k)
             end do
         end if
@@ -712,7 +718,11 @@ subroutine edgeVpack(edge,v,vlyr,kptr,ielem)
     do ll=swest+3*max_corner_elem,swest+4*max_corner_elem-1
         if (edge%putmap(ll,ielem) /= -1) then
             do k=1,vlyr
-                edge%buf(nce*(kptr+k-1)+edge%putmap(ll,ielem)+1)=v(np ,np,k)
+                iptr = (kptr+k-1)+edge%putmap(ll,ielem)+1
+                if (iptr > size(edge%buf)) then
+                    print *,'ERROR 3: ',size(edge%buf),iptr
+                endif
+                edge%buf(iptr)=v(np ,np,k)
 !                edge%buf(kptr+k,edge%putmap(ll,ielem)+1)=v(np ,np,k)
             end do
         end if
@@ -722,7 +732,11 @@ subroutine edgeVpack(edge,v,vlyr,kptr,ielem)
     do ll=swest+2*max_corner_elem,swest+3*max_corner_elem-1
         if (edge%putmap(ll,ielem) /= -1) then
             do k=1,vlyr
-                edge%buf(nce*(kptr+k-1)+edge%putmap(ll,ielem)+1)=v(1  ,np,k)
+                iptr = (kptr+k-1)+edge%putmap(ll,ielem)+1
+                if (iptr > size(edge%buf)) then
+                    print *,'ERROR 4: ',size(edge%buf),iptr
+                endif
+                edge%buf(iptr)=v(1  ,np,k)
 !                edge%buf(kptr+k,edge%putmap(ll,ielem)+1)=v(1  ,np,k)
             end do
         end if
@@ -748,7 +762,6 @@ subroutine edgeSpack(edge,v,vlyr,kptr,ielem)
     integer :: i,k,ir,ll,iptr
 
     integer :: is,ie,in,iw
-    integer :: nce
     real (kind=real_kind) :: tmp
 
 !pw call t_adj_detailf(+2)
@@ -769,13 +782,11 @@ subroutine edgeSpack(edge,v,vlyr,kptr,ielem)
        edge%buf(iptr+in+1)   = v(k) ! North
        edge%buf(iptr+iw+1)   = v(k) ! West
     enddo
-    !set the max_corner_elem
-    nce = max_corner_elem
 ! SWEST
     do ll=swest,swest+max_corner_elem-1
         if (edge%putmap(ll,ielem) /= -1) then
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 edge%buf(iptr+edge%putmap(ll,ielem)+1)=v(k)
             end do
         end if
@@ -785,7 +796,7 @@ subroutine edgeSpack(edge,v,vlyr,kptr,ielem)
     do ll=swest+max_corner_elem,swest+2*max_corner_elem-1
         if (edge%putmap(ll,ielem) /= -1) then
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 edge%buf(iptr+edge%putmap(ll,ielem)+1)=v(k)
             end do
         end if
@@ -795,7 +806,7 @@ subroutine edgeSpack(edge,v,vlyr,kptr,ielem)
     do ll=swest+3*max_corner_elem,swest+4*max_corner_elem-1
         if (edge%putmap(ll,ielem) /= -1) then
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 edge%buf(iptr+edge%putmap(ll,ielem)+1)=v(k)
             end do
         end if
@@ -805,7 +816,7 @@ subroutine edgeSpack(edge,v,vlyr,kptr,ielem)
     do ll=swest+2*max_corner_elem,swest+3*max_corner_elem-1
         if (edge%putmap(ll,ielem) /= -1) then
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 edge%buf(iptr+edge%putmap(ll,ielem)+1)=v(k)
             end do
         end if
@@ -972,7 +983,7 @@ subroutine edgeVunpack(edge,v,vlyr,kptr,ielem)
     !type (EdgeDescriptor_t)            :: desc
 
     ! Local
-    integer :: i,k,ll,iptr,nce
+    integer :: i,k,ll,iptr
     integer :: is,ie,in,iw
     integer :: ks,ke,kblock
     logical :: done
@@ -1014,11 +1025,10 @@ subroutine edgeVunpack(edge,v,vlyr,kptr,ielem)
     enddo
 
 ! SWEST
-    nce = max_corner_elem
     do ll=swest,swest+max_corner_elem-1
         if(edge%getmap(ll,ielem) /= -1) then 
             do k=1,vlyr
-                v(1  ,1 ,k)=v(1 ,1 ,k)+edge%receive(nce*(kptr+k-1)+edge%getmap(ll,ielem)+1)
+                v(1  ,1 ,k)=v(1 ,1 ,k)+edge%receive((kptr+k-1)+edge%getmap(ll,ielem)+1)
             enddo
         endif
     end do
@@ -1027,7 +1037,7 @@ subroutine edgeVunpack(edge,v,vlyr,kptr,ielem)
     do ll=swest+max_corner_elem,swest+2*max_corner_elem-1
         if(edge%getmap(ll,ielem) /= -1) then 
             do k=1,vlyr
-                v(np ,1 ,k)=v(np,1 ,k)+edge%receive(nce*(kptr+k-1)+edge%getmap(ll,ielem)+1)
+                v(np ,1 ,k)=v(np,1 ,k)+edge%receive((kptr+k-1)+edge%getmap(ll,ielem)+1)
             enddo
         endif
     end do
@@ -1036,7 +1046,7 @@ subroutine edgeVunpack(edge,v,vlyr,kptr,ielem)
     do ll=swest+3*max_corner_elem,swest+4*max_corner_elem-1
         if(edge%getmap(ll,ielem) /= -1) then 
             do k=1,vlyr
-                v(np ,np,k)=v(np,np,k)+edge%receive(nce*(kptr+k-1)+edge%getmap(ll,ielem)+1)
+                v(np ,np,k)=v(np,np,k)+edge%receive((kptr+k-1)+edge%getmap(ll,ielem)+1)
             enddo
         endif
     end do
@@ -1045,7 +1055,7 @@ subroutine edgeVunpack(edge,v,vlyr,kptr,ielem)
     do ll=swest+2*max_corner_elem,swest+3*max_corner_elem-1
         if(edge%getmap(ll,ielem) /= -1) then 
             do k=1,vlyr
-                v(1  ,np,k)=v(1 ,np,k)+edge%receive(nce*(kptr+k-1)+edge%getmap(ll,ielem)+1)
+                v(1  ,np,k)=v(1 ,np,k)+edge%receive((kptr+k-1)+edge%getmap(ll,ielem)+1)
             enddo
         endif
     end do
@@ -1338,6 +1348,7 @@ subroutine edgeDGVunpack(edge,v,vlyr,kptr,ielem)
     end do
 
     nce = max_corner_elem
+!   this is probably broken.  nce should be 1?  MT 2016/2/9
     i = swest
     if(edge%getmap(i,ielem) /= -1) then
       do k=1,vlyr
@@ -1387,7 +1398,7 @@ subroutine edgeVunpackMAX(edge,v,vlyr,kptr,ielem)
 
     ! Local
 
-    integer :: i,k,l,iptr,nce
+    integer :: i,k,l,iptr
     integer :: is,ie,in,iw
 
     threadsafe=.false.
@@ -1406,12 +1417,11 @@ subroutine edgeVunpackMAX(edge,v,vlyr,kptr,ielem)
        end do
     end do
 
-    nce = max_corner_elem
 ! SWEST
     do l=swest,swest+max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                v(1  ,1 ,k)=MAX(v(1 ,1 ,k),edge%receive(nce*(kptr+k-1)+edge%getmap(l,ielem)+1))
+                v(1  ,1 ,k)=MAX(v(1 ,1 ,k),edge%receive((kptr+k-1)+edge%getmap(l,ielem)+1))
             enddo
         endif
     end do
@@ -1420,7 +1430,7 @@ subroutine edgeVunpackMAX(edge,v,vlyr,kptr,ielem)
     do l=swest+max_corner_elem,swest+2*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                v(np ,1 ,k)=MAX(v(np,1 ,k),edge%receive(nce*(kptr+k-1)+edge%getmap(l,ielem)+1))
+                v(np ,1 ,k)=MAX(v(np,1 ,k),edge%receive((kptr+k-1)+edge%getmap(l,ielem)+1))
             enddo
         endif
     end do
@@ -1429,7 +1439,7 @@ subroutine edgeVunpackMAX(edge,v,vlyr,kptr,ielem)
     do l=swest+3*max_corner_elem,swest+4*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then
             do k=1,vlyr
-                v(np ,np,k)=MAX(v(np,np,k),edge%receive(nce*(kptr+k-1)+edge%getmap(l,ielem)+1))
+                v(np ,np,k)=MAX(v(np,np,k),edge%receive((kptr+k-1)+edge%getmap(l,ielem)+1))
             enddo
         endif
     end do
@@ -1438,7 +1448,7 @@ subroutine edgeVunpackMAX(edge,v,vlyr,kptr,ielem)
     do l=swest+2*max_corner_elem,swest+3*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                v(1  ,np,k)=MAX(v(1 ,np,k),edge%receive(nce*(kptr+k-1)+edge%getmap(l,ielem)+1))
+                v(1  ,np,k)=MAX(v(1 ,np,k),edge%receive((kptr+k-1)+edge%getmap(l,ielem)+1))
             enddo
         endif
     end do
@@ -1458,7 +1468,7 @@ subroutine edgeSunpackMAX(edge,v,vlyr,kptr,ielem)
 
     ! Local
 
-    integer :: i,k,l,iptr,nce
+    integer :: i,k,l,iptr
     integer :: is,ie,in,iw
 
 !pw call t_startf('edgeSunpack')
@@ -1473,12 +1483,11 @@ subroutine edgeSunpackMAX(edge,v,vlyr,kptr,ielem)
        v(k) = MAX(v(k),edge%receive(iptr+is+1),edge%receive(iptr+ie+1),edge%receive(iptr+in+1),edge%receive(iptr+iw+1))
     end do
 
-    nce = max_corner_elem
 ! SWEST
     do l=swest,swest+max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 v(k)=MAX(v(k),edge%receive(iptr+edge%getmap(l,ielem)+1))
             enddo
         endif
@@ -1488,7 +1497,7 @@ subroutine edgeSunpackMAX(edge,v,vlyr,kptr,ielem)
     do l=swest+max_corner_elem,swest+2*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 v(k)=MAX(v(k),edge%receive(iptr+edge%getmap(l,ielem)+1))
             enddo
         endif
@@ -1498,7 +1507,7 @@ subroutine edgeSunpackMAX(edge,v,vlyr,kptr,ielem)
     do l=swest+3*max_corner_elem,swest+4*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 v(k)=MAX(v(k),edge%receive(iptr+edge%getmap(l,ielem)+1))
             enddo
         endif
@@ -1508,7 +1517,7 @@ subroutine edgeSunpackMAX(edge,v,vlyr,kptr,ielem)
     do l=swest+2*max_corner_elem,swest+3*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 v(k)=MAX(v(k),edge%receive(iptr+edge%getmap(l,ielem)+1))
             enddo
         endif
@@ -1530,7 +1539,7 @@ subroutine edgeSunpackMIN(edge,v,vlyr,kptr,ielem)
 
     ! Local
 
-    integer :: i,k,l,iptr,nce
+    integer :: i,k,l,iptr
     integer :: is,ie,in,iw
 
 !pw call t_startf('edgeSunpack')
@@ -1545,12 +1554,11 @@ subroutine edgeSunpackMIN(edge,v,vlyr,kptr,ielem)
        v(k) = MIN(v(k),edge%receive(iptr+is+1),edge%receive(iptr+ie+1),edge%receive(iptr+in+1),edge%receive(iptr+iw+1))
     end do
 
-    nce = max_corner_elem
 ! SWEST
     do l=swest,swest+max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 v(k)=MiN(v(k),edge%receive(iptr+edge%getmap(l,ielem)+1))
             enddo
         endif
@@ -1560,7 +1568,7 @@ subroutine edgeSunpackMIN(edge,v,vlyr,kptr,ielem)
     do l=swest+max_corner_elem,swest+2*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 v(k)=MIN(v(k),edge%receive(iptr+edge%getmap(l,ielem)+1))
             enddo
         endif
@@ -1570,7 +1578,7 @@ subroutine edgeSunpackMIN(edge,v,vlyr,kptr,ielem)
     do l=swest+3*max_corner_elem,swest+4*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 v(k)=MIN(v(k),edge%receive(iptr+edge%getmap(l,ielem)+1))
             enddo
         endif
@@ -1580,7 +1588,7 @@ subroutine edgeSunpackMIN(edge,v,vlyr,kptr,ielem)
     do l=swest+2*max_corner_elem,swest+3*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                iptr = nce*(kptr+k-1)
+                iptr = (kptr+k-1)
                 v(k)=MIN(v(k),edge%receive(iptr+edge%getmap(l,ielem)+1))
             enddo
         endif
@@ -1602,7 +1610,7 @@ subroutine edgeVunpackMIN(edge,v,vlyr,kptr,ielem)
 
     ! Local
 
-    integer :: i,k,l,iptr,nce
+    integer :: i,k,l,iptr
     integer :: is,ie,in,iw
 
     threadsafe=.false.
@@ -1622,11 +1630,10 @@ subroutine edgeVunpackMIN(edge,v,vlyr,kptr,ielem)
     end do
 
 ! SWEST
-    nce = max_corner_elem
     do l=swest,swest+max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                v(1  ,1 ,k)=MIN(v(1 ,1 ,k),edge%receive(nce*(kptr+k-1)+edge%getmap(l,ielem)+1))
+                v(1  ,1 ,k)=MIN(v(1 ,1 ,k),edge%receive((kptr+k-1)+edge%getmap(l,ielem)+1))
             enddo
         endif
     end do
@@ -1635,7 +1642,7 @@ subroutine edgeVunpackMIN(edge,v,vlyr,kptr,ielem)
     do l=swest+max_corner_elem,swest+2*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                v(np ,1 ,k)=MIN(v(np,1 ,k),edge%receive(nce*(kptr+k-1)+edge%getmap(l,ielem)+1))
+                v(np ,1 ,k)=MIN(v(np,1 ,k),edge%receive((kptr+k-1)+edge%getmap(l,ielem)+1))
             enddo
         endif
     end do
@@ -1644,7 +1651,7 @@ subroutine edgeVunpackMIN(edge,v,vlyr,kptr,ielem)
     do l=swest+3*max_corner_elem,swest+4*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                v(np ,np,k)=MIN(v(np,np,k),edge%receive(nce*(kptr+k-1)+edge%getmap(l,ielem)+1))
+                v(np ,np,k)=MIN(v(np,np,k),edge%receive((kptr+k-1)+edge%getmap(l,ielem)+1))
             enddo
         endif
     end do
@@ -1653,7 +1660,7 @@ subroutine edgeVunpackMIN(edge,v,vlyr,kptr,ielem)
     do l=swest+2*max_corner_elem,swest+3*max_corner_elem-1
         if(edge%getmap(l,ielem) /= -1) then 
             do k=1,vlyr
-                v(1  ,np,k)=MIN(v(1 ,np,k),edge%receive(nce*(kptr+k-1)+edge%getmap(l,ielem)+1))
+                v(1  ,np,k)=MIN(v(1 ,np,k),edge%receive((kptr+k-1)+edge%getmap(l,ielem)+1))
             enddo
         endif
     end do
diff --git a/components/homme/src/share/unit_tests_mod.F90 b/components/homme/src/share/unit_tests_mod.F90
index bb2c7eb7f239..dc4fe759a00c 100644
--- a/components/homme/src/share/unit_tests_mod.F90
+++ b/components/homme/src/share/unit_tests_mod.F90
@@ -431,7 +431,8 @@ subroutine test_subcell_div_fluxes(elem,deriv,nets,nete)
 
       do i=1,intervals
       do j=1,intervals
-        t = ABS(test(i,j)-values(i,j))/MAX(ABS(test(i,j)),ABS(values(i,j)))
+!        t = ABS(test(i,j)-values(i,j))/MAX(ABS(test(i,j)),ABS(values(i,j)))
+        t = ABS(test(i,j)-values(i,j))/MAX(ABS(MAXVAL(test(:,:))),ABS(MAXVAL(values(:,:))))
         if (.0000001<t) then
           print *,__FILE__,__LINE__,ie,i,j,test(i,j),values(i,j),t
           success = .false.
diff --git a/components/homme/test/jw_baroclinic/baro.job b/components/homme/test/jw_baroclinic/baro.job
index e0953defc583..fe9089fdf53f 100755
--- a/components/homme/test/jw_baroclinic/baro.job
+++ b/components/homme/test/jw_baroclinic/baro.job
@@ -25,7 +25,7 @@ set input = $HOMME/test/jw_baroclinic    # input files for test case
 set vdir = $HOMME/test/vcoord            # vertical coordinate files
 set bld = $wdir/bld                      # cmake/build directory
 
-set MACH = $HOMME/cmake/machineFiles/redsky.cmake
+set MACH = $HOMME/cmake/machineFiles/skybride.cmake
 #set MACH = $HOMME/cmake/machineFiles/titan.cmake
 #set MACH = $HOMME/cmake/machineFiles/titan-cuda.cmake    # CMAKE file for Titan GPU support
 #set MACH = $HOMME/cmake/machineFiles/darwin.cmake
@@ -125,7 +125,8 @@ set hypervis_scaling = 0
 
 if ( $ne == 0 ) then
     # CAM defaults
-    set mesh_file="$HOME/work/notes/squadgen/multiscale/wpacific_v2.g"
+    set mesh_file="/ascldap/users/mataylo/codes/mapping/grids/wpacific_v2.g"
+    #set mesh_file="/ascldap/users/mataylo/codes/mapping/grids/mountain_10_x8.g"
     set tstep=75
     set hypervis_scaling = 3.2  
     set nu = 8.0e-8
diff --git a/components/homme/test/sw_conservative/swtc6-tensor-hv.sh b/components/homme/test/sw_conservative/swtc6-tensor-hv.sh
index 6cf190496a6b..3c440c1f0043 100755
--- a/components/homme/test/sw_conservative/swtc6-tensor-hv.sh
+++ b/components/homme/test/sw_conservative/swtc6-tensor-hv.sh
@@ -1,8 +1,9 @@
 #!/bin/bash -f
 set -e
-HOMME=~/codes/acme-dev/models/atm/homme
+HOMME=~/codes/acme-dev/components/homme
+MACH=$HOMME/cmake/machineFiles/skybridge.cmake
+
 wdir=~/scratch1/sweqx
-MACH=$HOMME/cmake/machineFiles/redsky.cmake
 NCPU=0
 
 
@@ -15,6 +16,8 @@ hvscaling=3.2
 test_case=swtc6
 name=${test_case}
 mesh=~/codes/mapping/grids/mountain_10_x8.g
+mesh=~/codes/homme-ncarsvn/test/mesh_refine/grids/uniform_15.g
+
 
 #
 # run/build directoryies
@@ -31,7 +34,7 @@ cd $builddir
 if ! ([ -e  CMakeCache.txt ]) then
   rm -rf CMakeFiles CMakeCache.txt
   #cmake -C $MACH -DSWEQX_USE_PIO=1 -DSWEQX_PLEV=1  -DSWEQX_NP=4 $HOMME
-  cmake -C $MACH -DSWEQX_PLEV=1  -DSWEQX_NP=4 $HOMME
+  cmake -C $MACH -DSWEQX_PLEV=4  -DSWEQX_NP=4 $HOMME
   exit
 fi
 
@@ -46,6 +49,8 @@ rsync $mesh .
 let sfreq=24*3600
 sfreq=`echo "$sfreq / $tstep" | bc`
 
+
+
 	sed s/tstep.\*/"tstep = $tstep"/  $input/swtc6-tensor-hv.nl |\
 	sed s/hypervis_scaling.\*/"hypervis_scaling = $hvscaling"/  |\
 	sed s/nu=.\*/"nu= $nu"/  |\