修改了两处不严谨的地方，gfortran环境下会warning，1：get_non_uniformgrid.f90中，IF判断中IF(SourceLength/GridSize.gt.51)改SourceLength/GridSize为整型SourceGridNum；2：替换旧写法PAUSE语句为READ(*,*)

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
+	if(SourceGridNum.gt.51)then
-     do ii=nxs-(SourceLength/GridSize-1)/2,nxs+(SourceLength/GridSize-1)/2,1
+     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
+	if(SourceGridNum.gt.51)then
-		do ii=nxs-(SourceLength/GridSize-1)/2,nxs+1+(SourceLength/GridSize-1)/2,1
+		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
+	IF(SourceGridNum.gt.51)then
-		do ii=nys-(SourceLength/GridSize-1)/2,nys+(SourceLength/GridSize-1)/2,1
+		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
+	IF(SourceGridNum.gt.51)then
-		do ii=nys-(SourceLength/GridSize-1)/2,nys+1+(SourceLength/GridSize-1)/2,1
+		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

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

tem3dfdtd/module/VTK_fortran.f90

@@ -204,7 +204,7 @@ Module VTK_Fortran
                     do i=1,NXB
                         
                         temp = temp + 1
-                        write(unit, fmt='(3f)', iostat=error) Coord_x(i), Coord_y(j), Coord_z(k)
+                        write(unit, fmt='(3F25.16)', iostat=error) Coord_x(i), Coord_y(j), Coord_z(k)
                         
                     end do
                 end do
@@ -250,7 +250,7 @@ Module VTK_Fortran
                 do j=1, NY
                     do i=1,NX
                         temp = temp + 1
-                        write(unit, fmt='(9i)', iostat=error)nums,&
+                        write(unit, fmt='(9I12)', iostat=error)nums,&
                             ((i  ) + (j-1)*NXB + (k-1)*NXB*NYB) -1,  &
                             ((i+1) + (j-1)*NXB + (k-1)*NXB*NYB) -1,  &
                             ((i  ) + (j  )*NXB + (k-1)*NXB*NYB) -1,  &
@@ -266,7 +266,7 @@ Module VTK_Fortran
             write(unit=unit,fmt='(A,I0,A,I0)', iostat=error) "CELL_TYPES ", cells
             
             do i=1, temp
-                write(unit, fmt='(i)', iostat=error)num_cell_types
+                write(unit, fmt='(I12)', iostat=error)num_cell_types
             end do
             write(unit=unit,fmt='(2A,I0)', iostat=error) Value_site_cell," ", cells
 
@@ -311,7 +311,7 @@ Module VTK_Fortran
                 do j=1, NY
                     do i=1,NX
                         temp = temp + 1
-                        write(unit, fmt='(f)', iostat=error)values(i,j,k)
+                        write(unit, fmt='(F25.16)', iostat=error)values(i,j,k)
                     end do
                 end do
             end do