need help very urgent

This is my final year project...I have to run the program using silverfrost ftn95 in order to complete my project.But as I compiled it shows to much of errors.I have no idea about it.So, could you guys please help me...Your assistance is highly appreciated.Thank you.( I need to submit the result by this Thursday)

so this is the code.

c####################################################################### c c SPH: Smoothed Particle Hydrodynamics c c####################################################################### c program SPH c implicit real*8(a-h,o-z) c include 'parameter.f' c parameter( & TOTAL_TIME=5.0D-6, ! ƒVƒ_{ƒ…ƒŒ[ƒVƒ‡ƒ“‚̍ő厞ŠÔ & MAXLOOP=100000) ! ƒVƒ}ƒ…ƒŒ[ƒVƒ‡ƒ“‚̍őåŒJ‚è•Ô‚µ” c common/model/ DISA,XMAX,XMIN,XSIZE,YMAX,YMIN,YSIZE, & NUMX1,NUMY1,NUMX2,NUMY2,NUMP,NP1,NP2 c common/particle/ X(NDIM,MAXP), X0(NDIM,MAXP) c common/property/ RHO1,RHO2,CS1,CS2,G1,G2, & RHO(MAXP),PMASS(MAXP),CS(MAXP) c common/composite/ CMAT(3,3,2),MATER(MAXP) c common/boundary/ & V_BC_GIVEN(NDIM,MAXBC),V_LOAD_GIVEN(NDIM,MAXBC), & I_BC_GIVEN(NDIM,MAXBC),I_LOAD_GIVEN(NDIM,MAXBC), & N_BC_GIVEN(NDIM),N_LOAD_GIVEN(NDIM),SPEED c common/calc/ V(NDIM,MAXP),VISC(NEIGHBOR,MAXP),P(MAXP),E(MAXP) c common/param/ H0,H,DT,PDT c common/neighbor_table/ & NEIGH(NEIGHBOR,MAXP),NNEIGH(MAXP), & NEIGH0(NEIGHBOR,MAXP),NNEIGH0(MAXP),JCONNECT0(NEIGHBOR,MAXP), & RIJ(NDIM,NEIGHBOR,MAXP),RIJ1(NEIGHBOR,MAXP),RIJ2(NEIGHBOR,MAXP), & RIJ0(NEIGHBOR,MAXP) c common/kernel/ GRADW(NDIM,NEIGHBOR,MAXP) c common/stress/ & SIGM(3,MAXP),HENSAS(3,MAXP),HENSASOLD(3,MAXP),STRNDOT(3,MAXP) c common/vardot/ DRHODT(MAXP),DVDT(NDIM,MAXP),DEDT(MAXP) c common/rkgill/ DX(NDIM,MAXP),DV(NDIM,MAXP),DR(MAXP),DE(MAXP) c common/yield/ SGMY(MAXP),EPSP(MAXP),T(MAXP),JYIELD(MAXP) c common/damage/ DMG(MAXP) c c==== output files open(2,file='rst.txt') open(3,file='out_damage.plt') open(4,file='rst_damage.csv') open(7,file='out_stress.plt') open(8,file='out_pstress.plt') open(100,file='check_coordinates.csv') c c==== input call inp(angle_inc) c call set_param c call set_constitutive_matrix c c==== initialize do ip=1,NUMP X(1:NDIM,ip) = X0(1:NDIM,ip) ! À•W E(ip) = 0.0D0 ! ƒGƒlƒ‹ƒM[ do k=1,3 SIGM(k,ip) = 0.0D0 ! ‰ž—Í HENSAS(k,ip) = 0.0D0 ! •Î·‰ž—Í HENSASOLD(k,ip) = 0.0D0 STRNDOT(k,ip) = 0.0D0 ! ‚Ђ¸‚Ý‘¬“x enddo c DX(1:NDIM,ip) = 0.0D0 ! ƒ‹ƒ“ƒQEƒNƒbƒ^–@ DV(1:NDIM,ip) = 0.0D0 ! ‚Å‚ÌŽžŠÔÏ•ª‚ÅŽg‚¤ DR(ip) = 0.0D0 DE(ip) = 0.0D0 c SGMY(ip) = YIELD0 ! ~•š‰ž—Í EPSP(ip) = 0.0D0 ! ‘Š“–‘Y«‚Ђ¸‚Ý T(ip) = 300.0D0 ! ‰·“x c DMG(ip) = 0.0D0 ! ‘¹ƒpƒ‰ƒ[ƒ^ enddo c c---- define boundary conditions call definebc c V(1:NDIM,1:NUMP) = 0.0D0 ! ‘¬“x‚̏‰Šú‰» c do ip=NP1+1,NUMP ! ”òãđ̂̏‰‘¬ V(1,ip) = SPEED * dcos(angle_inc) V(2,ip) = SPEED * dsin(angle_inc) enddo

Tidied up a bit and unicode characters removed.

c####################################################################### 
c 
c SPH: Smoothed Particle Hydrodynamics 
c 
c####################################################################### 
c 
      program SPH 
c 
      implicit real*8(a-h,o-z) 
c 
c      include 'parameter.f' 
c??????????????????????????????????????????????????????????????????????????????????????????????????
c parameter statement to satisfy missing values, assumed to be in 'parameter.f' commented out above
      parameter( maxp=1000,ndim=100,maxbc=200,neighbor=10) 
c 
      parameter( TOTAL_TIME=5.0D-6,MAXLOOP=100000) 
c 
      common/model/ DISA,XMAX,XMIN,XSIZE,YMAX,YMIN,YSIZE, 
     & NUMX1,NUMY1,NUMX2,NUMY2,NUMP,NP1,NP2 
c 
      common/particle/ X(NDIM,MAXP), X0(NDIM,MAXP) 
c 
      common/property/ RHO1,RHO2,CS1,CS2,G1,G2, 
     & RHO(MAXP),PMASS(MAXP),CS(MAXP) 
c 
      common/composite/ CMAT(3,3,2),MATER(MAXP) 
c 
      common/boundary/ 
     & V_BC_GIVEN(NDIM,MAXBC),V_LOAD_GIVEN(NDIM,MAXBC), 
     & I_BC_GIVEN(NDIM,MAXBC),I_LOAD_GIVEN(NDIM,MAXBC), 
     & N_BC_GIVEN(NDIM),N_LOAD_GIVEN(NDIM),SPEED 
c 
      common/calc/ V(NDIM,MAXP),VISC(NEIGHBOR,MAXP),P(MAXP),E(MAXP) 
c 
      common/param/ H0,H,DT,PDT 
c 
      common/neighbor_table/ 
     & NEIGH(NEIGHBOR,MAXP),NNEIGH(MAXP), 
     & NEIGH0(NEIGHBOR,MAXP),NNEIGH0(MAXP),JCONNECT0(NEIGHBOR,MAXP), 
     & RIJ(NDIM,NEIGHBOR,MAXP),RIJ1(NEIGHBOR,MAXP),RIJ2(NEIGHBOR,MAXP), 
     & RIJ0(NEIGHBOR,MAXP) 
c 
      common/kernel/ GRADW(NDIM,NEIGHBOR,MAXP) 
c 
      common/stress/ 
     & SIGM(3,MAXP),HENSAS(3,MAXP),HENSASOLD(3,MAXP),STRNDOT(3,MAXP) 
c 
      common/vardot/ DRHODT(MAXP),DVDT(NDIM,MAXP),DEDT(MAXP) 
c 
      common/rkgill/ DX(NDIM,MAXP),DV(NDIM,MAXP),DR(MAXP),DE(MAXP) 
c 
      common/yield/ SGMY(MAXP),EPSP(MAXP),T(MAXP),JYIELD(MAXP) 
c 
      common/damage/ DMG(MAXP) 
c 
c==== output files 
c??????????????????????????????????????????????????????????????????????????????????????????????????
c better to use higher numbers for files to avoid conflict with predefined unit numbers
c use 12, 13, 14, 17 & 18
      open(12,file='rst.txt') 
      open(13,file='out_damage.plt') 
      open(14,file='rst_damage.csv') 
      open(17,file='out_stress.plt') 
      open(18,file='out_pstress.plt') 
      open(100,file='check_coordinates.csv') 
c??????????????????????????????????????????????????????????????????????????????????????????????????
c  assign a value to 'YIELDO'
      YIELD0 = 415d0 !MPa
c 
c==== input 
      call inp(angle_inc) 
c 
      call set_param 
c 
      call set_constitutive_matrix 
c 
c==== initialize 
      do ip=1,NUMP 
        X(1:NDIM,ip) = X0(1:NDIM,ip) 
        E(ip) = 0.0D0
        do k=1,3 
          SIGM(k,ip) = 0.0D0 
          HENSAS(k,ip) = 0.0D0  
          HENSASOLD(k,ip) = 0.0D0 
          STRNDOT(k,ip) = 0.0D0 
        enddo 
c 
        DX(1:NDIM,ip) = 0.0D0 
        DV(1:NDIM,ip) = 0.0D0
        DR(ip) = 0.0D0 
        DE(ip) = 0.0D0 
c 
        SGMY(ip) = YIELD0
        EPSP(ip) = 0.0D0 
        T(ip) = 300.0D0 
c 
        DMG(ip) = 0.0D0 
      enddo 

Continued in next post ...

c 
c---- define boundary conditions 
      call definebc 
c 
      V(1:NDIM,1:NUMP) = 0.0D0  
c 
      do ip=NP1+1,NUMP 
        V(1,ip) = SPEED * dcos(angle_inc) 
        V(2,ip) = SPEED * dsin(angle_inc) 
      enddo

      end