diff --git a/cat/dat/sim_0051.nl b/cat/dat/sim_0051.nl
index 570a43075997bcc4429a2736c064c9c6274a75f9..77274c7df4af9fc8ccb3d52e7cd07b869a33b7a4 100644
--- a/cat/dat/sim_0051.nl
+++ b/cat/dat/sim_0051.nl
@@ -2,7 +2,7 @@
  ysim    = "fjet01"
  w1      = 8
  w2      = 2
- scl     = 1
+ hscl    = 1
  qmax    = 0.00001
  snsteps = 2000000
  sdt     = 1.0d-3
diff --git a/cat/dat/sim_0052.nl b/cat/dat/sim_0052.nl
index cbdc90412c6791a64dd601be13dda16af7daf748..79abeea14ba09a9dfa58d779a0d08818ea8aeea3 100644
--- a/cat/dat/sim_0052.nl
+++ b/cat/dat/sim_0052.nl
@@ -2,7 +2,7 @@
  ysim    = "fjet01"
  w1      = 8
  w2      = 2
- scl     = 1
+ hscl    = 1
  qmax    = 0.00001
  snsteps = 2000000
  sdt     = 1.0d-3
diff --git a/cat/dat/sim_0053.nl b/cat/dat/sim_0053.nl
new file mode 100644
index 0000000000000000000000000000000000000000..200d943ddadda17f7926c8a40e3f39d2870574d4
--- /dev/null
+++ b/cat/dat/sim_0053.nl
@@ -0,0 +1,12 @@
+ $sim_nl
+ ysim    = "fjet01"
+ w1      = 8
+ w2      = 2
+ hscl    = 1
+ qmax    = 0.00001
+ snsteps = 2000000
+ sdt     = 1.0d-3
+ sngui   = 250
+ snforc  = 1
+ snpert  = 5
+ /
diff --git a/cat/dat/sim_0061.nl b/cat/dat/sim_0061.nl
index 8307319c494fd2ed7eaa888f0ade10016619a71d..096f1e5a5e03cd5a86cb3586b3011080ba443f5a 100644
--- a/cat/dat/sim_0061.nl
+++ b/cat/dat/sim_0061.nl
@@ -2,7 +2,7 @@
  ysim    = "fjet02"
  w1      = 8
  w2      = 2
- scl     = 1
+ hscl    = 1
  qmax    = 0.00001
  snsteps = 2000000
  sdt     = 1.0d-3
diff --git a/cat/dat/sim_0062.nl b/cat/dat/sim_0062.nl
index 0ac6730db1e99307445ecda34b291a4c8f75f825..5ff17f4c159fa3605505e391d0e320393fa397a9 100644
--- a/cat/dat/sim_0062.nl
+++ b/cat/dat/sim_0062.nl
@@ -2,7 +2,7 @@
  ysim    = "fjet02"
  w1      = 8
  w2      = 2
- scl     = 1
+ hscl    = 1
  qmax    = 0.00001
  snsteps = 2000000
  sdt     = 1.0d-3
diff --git a/cat/src/simmod.f90 b/cat/src/simmod.f90
index 9d2bf130244412463c5ee7478a893bbe6e06f2ae..62b2b35be8f75dd21a38ada8189549d3f44bc45e 100644
--- a/cat/src/simmod.f90
+++ b/cat/src/simmod.f90
@@ -36,7 +36,7 @@ character (256) :: ysim = "djet01" ! type of predefined simulation
 !--- parameters of djet01 (initial top hat jet)
 integer :: w1     = 8      ! half width of jet center (in grid points)
 integer :: w2     = 4      ! width of vortex sheet (in grid points)
-integer :: scl    = 1      ! horizontal scale of jet
+integer :: hscl   = 1      ! horizontal scale of jet
 real(8) :: qmax   = 1.0    ! amplitude of vortex sheets
 
 
@@ -83,7 +83,7 @@ implicit none
 
 !--- define sim_namelist
 namelist /sim_nl/ ysim       ,                        &
-                  qmax       ,w1      ,w2     ,scl   ,& 
+                  qmax       ,w1      ,w2     ,hscl  ,& 
                   sdt        ,snsteps ,sngui  ,       &
                   snpert     ,snforc  ,snpost
 
@@ -131,12 +131,12 @@ select case(ysim)
   case("djet01")
      gpvar(:,:) = 0.0
      do jy = 1, ngy
-        if ( jy .ge. ngy/2+1-scl*(w1+w2) .and. & 
-           jy .le. ngy/2-scl*w1 ) then
+        if ( jy .ge. ngy/2+1-hscl*(w1+w2) .and. & 
+           jy .le. ngy/2-hscl*w1 ) then
            gpvar(:,jy) = -qmax
         endif
-        if ( jy .ge. ngy/2+1+scl*w1 .and. & 
-           jy .le. ngy/2+scl*(w1+w2) ) then
+        if ( jy .ge. ngy/2+1+hscl*w1 .and. & 
+           jy .le. ngy/2+hscl*(w1+w2) ) then
            gpvar(:,jy) =  qmax
         endif
      enddo
@@ -150,12 +150,12 @@ select case(ysim)
   case("fjet01")
      gpvar(:,:) = 0.0
      do jy = 1, ngy
-        if ( jy .ge. ngy/2+1-scl*(w1+w2) .and. & 
-           jy .le. ngy/2-scl*w1 ) then
+        if ( jy .ge. ngy/2+1-hscl*(w1+w2) .and. & 
+           jy .le. ngy/2-hscl*w1 ) then
            gpvar(:,jy) = -qmax
         endif
-        if ( jy .ge. ngy/2+1+scl*w1 .and. & 
-           jy .le. ngy/2+scl*(w1+w2) ) then
+        if ( jy .ge. ngy/2+1+hscl*w1 .and. & 
+           jy .le. ngy/2+hscl*(w1+w2) ) then
            gpvar(:,jy) =  qmax
         endif
      enddo
@@ -164,14 +164,14 @@ select case(ysim)
   case("fjet02")
      gpvar(:,:) = 0.0
      do jx = 1, ngx
-        if ( jx .ge. ngx/2+1-scl*(w1+w2) .and. & 
-           jx .le. ngx/2-scl*w1 ) then
-           gpvar(jx,:) =  qmax
-        endif
-        if ( jx .ge. ngx/2+1+scl*w1 .and. & 
-           jx .le. ngx/2+scl*(w1+w2) ) then
+        if ( jx .ge. ngx/2+1-hscl*(w1+w2) .and. & 
+           jx .le. ngx/2-hscl*w1 ) then
            gpvar(jx,:) = -qmax
         endif
+        if ( jx .ge. ngx/2+1+hscl*w1 .and. & 
+           jx .le. ngx/2+hscl*(w1+w2) ) then
+           gpvar(jx,:) =  qmax
+        endif
      enddo
      call sim_wrtgp(gpvar,qfrccde,1)
   case default
@@ -262,6 +262,11 @@ subroutine simstep
 use simmod
 implicit none
 
+
+
+
+
+
 return
 end subroutine simstep