diff --git a/tem3dfdtd/lib/get_non_uniformgrid.f90 b/tem3dfdtd/lib/get_non_uniformgrid.f90
index 69ede10..85b7726 100644
--- a/tem3dfdtd/lib/get_non_uniformgrid.f90
+++ b/tem3dfdtd/lib/get_non_uniformgrid.f90
@@ -30,17 +30,17 @@ SUBROUTINE GET_NON_UNIFORMGRID
   ! ----------------------------end of mesh------------------------------------!
   ! -------------------------------mesh x----------------------------------------!
   IF(Logi_Sourcelenth) THEN    !The number of grids occupied by the source is odd
-	if(SourceLength/GridSize.gt.51)then
-     do ii=nxs-(SourceLength/GridSize-1)/2,nxs+(SourceLength/GridSize-1)/2,1
+	if(SourceGridNum.gt.51)then
+     do ii=nxs-(SourceGridNum-1)/2,nxs+(SourceGridNum-1)/2,1
         Cdelx(ii)=GridSize
      end do
-     do ii=nxs-(SourceLength/GridSize-1)/2-1,1,-1
+     do ii=nxs-(SourceGridNum-1)/2-1,1,-1
         Cdelx(ii)=Cdelx(ii+1)*scale_par
         if(Cdelx(ii).gt.GridSize_MAX)then
            Cdelx(ii)=GridSize_MAX
         end if
      end do
-     do ii=nxs+(SourceLength/GridSize-1)/2+1,nx,1
+     do ii=nxs+(SourceGridNum-1)/2+1,nx,1
         Cdelx(ii)=Cdelx(ii-1)*scale_par
         if(Cdelx(ii).gt.GridSize_MAX)then
            Cdelx(ii)=GridSize_MAX
@@ -64,17 +64,17 @@ SUBROUTINE GET_NON_UNIFORMGRID
      end do
 	ENDIF
   ELSE    !The number of grids occupied by the source is even
-	if(SourceLength/GridSize.gt.51)then
-		do ii=nxs-(SourceLength/GridSize-1)/2,nxs+1+(SourceLength/GridSize-1)/2,1
+	if(SourceGridNum.gt.51)then
+		do ii=nxs-(SourceGridNum-1)/2,nxs+1+(SourceGridNum-1)/2,1
 			Cdelx(ii)=GridSize
 		end do
-		do ii=nxs-(SourceLength/GridSize-1)/2-1,1,-1
+		do ii=nxs-(SourceGridNum-1)/2-1,1,-1
 			Cdelx(ii)=Cdelx(ii+1)*scale_par
 			if(Cdelx(ii).gt.GridSize_MAX)then
 				Cdelx(ii)=GridSize_MAX
 			end if
 		end do
-		do ii=nxs+(SourceLength/GridSize-1)/2+2,nx,1
+		do ii=nxs+(SourceGridNum-1)/2+2,nx,1
 			Cdelx(ii)=Cdelx(ii-1)*scale_par
 			if(Cdelx(ii).gt.GridSize_MAX)then
 				Cdelx(ii)=GridSize_MAX
@@ -101,17 +101,17 @@ SUBROUTINE GET_NON_UNIFORMGRID
   ! -----------------------------end of mesh------------------------------------!
   ! --------------------------------mesh y----------------------------------------!
   IF(Logi_Sourcelenth) THEN    !The number of grids occupied by the source is odd
-	IF(SourceLength/GridSize.gt.51)then
-		do ii=nys-(SourceLength/GridSize-1)/2,nys+(SourceLength/GridSize-1)/2,1
+	IF(SourceGridNum.gt.51)then
+		do ii=nys-(SourceGridNum-1)/2,nys+(SourceGridNum-1)/2,1
 			Cdely(ii)=GridSize
 		enddo
-		do ii=nys-(SourceLength/GridSize-1)/2-1,1,-1
+		do ii=nys-(SourceGridNum-1)/2-1,1,-1
 			Cdely(ii)=Cdely(ii+1)*scale_par
 			if(Cdely(ii).gt.GridSize_MAX)then
 				Cdely(ii)=GridSize_MAX
 			end if
 		end do
-		do ii=nys+(SourceLength/GridSize-1)/2+1,ny,1
+		do ii=nys+(SourceGridNum-1)/2+1,ny,1
 			Cdely(ii)=Cdely(ii-1)*scale_par
 			if(Cdely(ii).gt.GridSize_MAX)then
 				Cdely(ii)=GridSize_MAX
@@ -135,17 +135,17 @@ SUBROUTINE GET_NON_UNIFORMGRID
 		end do
 	end if
   ELSE	!The number of grids occupied by the source is even
-	IF(SourceLength/GridSize.gt.51)then
-		do ii=nys-(SourceLength/GridSize-1)/2,nys+1+(SourceLength/GridSize-1)/2,1
+	IF(SourceGridNum.gt.51)then
+		do ii=nys-(SourceGridNum-1)/2,nys+1+(SourceGridNum-1)/2,1
 			Cdely(ii)=GridSize
 		enddo
-		do ii=nys-(SourceLength/GridSize-1)/2-1,1,-1
+		do ii=nys-(SourceGridNum-1)/2-1,1,-1
 			Cdely(ii)=Cdely(ii+1)*scale_par
 			if(Cdely(ii).gt.GridSize_MAX)then
 				Cdely(ii)=GridSize_MAX
 			end if
 		end do
-		do ii=nys+(SourceLength/GridSize-1)/2+2,ny,1
+		do ii=nys+(SourceGridNum-1)/2+2,ny,1
 			Cdely(ii)=Cdely(ii-1)*scale_par
 			if(Cdely(ii).gt.GridSize_MAX)then
 				Cdely(ii)=GridSize_MAX
diff --git a/tem3dfdtd/lib/time-serious.f90 b/tem3dfdtd/lib/time-serious.f90
index e150da7..5554b8d 100644
--- a/tem3dfdtd/lib/time-serious.f90
+++ b/tem3dfdtd/lib/time-serious.f90
@@ -62,7 +62,7 @@ SUBROUTINE TIME_SERIOUS
      WRITE(10005,*)'The number of iteration steps exceeds the range given in the input.dat, please change it. now the Nstop value is determined by Max_off_time.'
      print*,'The number of iteration steps exceeds the range given in the input.dat, please change it. now the Nstop value is determined by Max_off_time.'
      print*,'I give you a pause here, you should decide to continue or to quit'
-     pause
+     READ(*,*)
      Nstop=Nstop_temp
   ENDIF
   RETURN