Errors of Unclassifiable statement & Invalid character in name in the fortran code compilling -

i have written parts of fortran code attached. when have compiled it, have received errors. code:

!  ! file part of siesta package. ! ! copyright (c) fundacion general universidad autonoma de madrid: ! e.artacho, j.gale, a.garcia, j.junquera, p.ordejon, d.sanchez-portal ! , j.m.soler, 1996- . !  ! use of software constitutes agreement full conditions ! given in siesta license, signed legitimate users. !       subroutine constr( cell, na, isa, amass, xa, stress, fa, ntcon ) c ***************************************************************** c user-written routine implement specific geometric constraints, c orthogonalizing forces , stress undesired changes. c arguments: c real*8  cell(3,3)    : input lattice vectors (bohr) c integer na           : input number of atoms c integer isa(na)      : input species indexes c real*8  amass(na)    : input atomic masses c real*8  xa(3,na)     : input atomic cartesian coordinates (bohr) c real*8  stress( 3,3) : input/output stress tensor (ry/bohr**3) c real*8  fa(3,na)     : input/output atomic forces (ry/bohr) c integer ntcon        : total number of positions constr. imposed c *****************************************************************       implicit         none       integer           na, isa(na), ntcon       double precision  amass(na), cell(3,3), fa(3,na), fz       double precision  stress(3,3), xa(3,na), tmass          fz = fa(3,1)+fa(3,2)+fa(3,3)+fa(3,4)+fa(3,5)+fa(3,6)+fa(3,7)       $ +fa(3,8)+fa(3,9)+fa(3,10)       $ +fa(3,11)+fa(3,12)+fa(3,13)+fa(3,14)+fa(3,15)+fa(3,16)+fa(3,17)       $ +fa(3,18)+fa(3,19)+fa(3,20)       $ +fa(3,21)+fa(3,22)+fa(3,23)+fa(3,24)+fa(3,25)+fa(3,26)+fa(3,27)       $ +fa(3,28)+fa(3,29)+fa(3,30)       $ +fa(3,31)+fa(3,32)+fa(3,33)+fa(3,34)+fa(3,35)+fa(3,36)+fa(3,37)       $ +fa(3,38)+fa(3,39)+fa(3,40)       $ +fa(3,41)+fa(3,42)+fa(3,43)+fa(3,44)+fa(3,45)+fa(3,46)+fa(3,47)       $ +fa(3,48)+fa(3,49)+fa(3,50)       $ +fa(3,51)+fa(3,52)+fa(3,53)+fa(3,54)+fa(3,55)+fa(3,56)+fa(3,57)       $ +fa(3,58)+fa(3,59)+fa(3,60)         tmass=amass(1)+amass(2)+amass(3)+amass(4)+amass(5)+amass(6)       $ +amass(7)+amass(8)+amass(9)+amass(10)       $ +amass(11)+amass(12)+amass(13)+amass(14)+amass(15)+amass(16)       $ +amass(17)+amass(18)+amass(19)+amass(20)       $ +amass(21)+amass(22)+amass(23)+amass(24)+amass(25)+amass(26)       $ +amass(27)+amass(28)+amass(29)+amass(30)       $ +amass(31)+amass(32)+amass(33)+amass(34)+amass(35)+amass(36)       $ +amass(37)+amass(38)+amass(39)+amass(40)       $ +amass(41)+amass(42)+amass(43)+amass(44)+amass(45)+amass(46)       $ +amass(47)+amass(48)+amass(49)+amass(50)       $ +amass(51)+amass(52)+amass(53)+amass(54)+amass(55)+amass(56)       $ +amass(57)+amass(58)+amass(59)+amass(60)         fa(3,1) = fz * amass(1)/tmass        fa(3,2) = fz * amass(2)/tmass        fa(3,3) = fz * amass(3)/tmass        fa(3,4) = fz * amass(4)/tmass        fa(3,5) = fz * amass(5)/tmass        fa(3,6) = fz * amass(6)/tmass        fa(3,7) = fz * amass(7)/tmass        fa(3,8) = fz * amass(8)/tmass        fa(3,9) = fz * amass(9)/tmass        fa(3,10) = fz * amass(10)/tmass        fa(3,11) = fz * amass(11)/tmass        fa(3,12) = fz * amass(12)/tmass        fa(3,13) = fz * amass(13)/tmass        fa(3,14) = fz * amass(14)/tmass        fa(3,15) = fz * amass(15)/tmass        fa(3,16) = fz * amass(16)/tmass        fa(3,17) = fz * amass(17)/tmass        fa(3,18) = fz * amass(18)/tmass        fa(3,19) = fz * amass(19)/tmass        fa(3,20) = fz * amass(20)/tmass        fa(3,21) = fz * amass(21)/tmass        fa(3,22) = fz * amass(22)/tmass        fa(3,23) = fz * amass(23)/tmass        fa(3,24) = fz * amass(24)/tmass        fa(3,25) = fz * amass(25)/tmass        fa(3,26) = fz * amass(26)/tmass        fa(3,27) = fz * amass(27)/tmass        fa(3,28) = fz * amass(28)/tmass        fa(3,29) = fz * amass(29)/tmass        fa(3,30) = fz * amass(30)/tmass        fa(3,31) = fz * amass(31)/tmass        fa(3,32) = fz * amass(32)/tmass        fa(3,33) = fz * amass(33)/tmass        fa(3,34) = fz * amass(34)/tmass        fa(3,35) = fz * amass(35)/tmass        fa(3,36) = fz * amass(36)/tmass        fa(3,37) = fz * amass(37)/tmass        fa(3,38) = fz * amass(38)/tmass        fa(3,39) = fz * amass(39)/tmass        fa(3,40) = fz * amass(40)/tmass        fa(3,41) = fz * amass(41)/tmass        fa(3,42) = fz * amass(42)/tmass        fa(3,43) = fz * amass(43)/tmass        fa(3,44) = fz * amass(44)/tmass        fa(3,45) = fz * amass(45)/tmass        fa(3,46) = fz * amass(46)/tmass        fa(3,47) = fz * amass(47)/tmass        fa(3,48) = fz * amass(48)/tmass        fa(3,49) = fz * amass(49)/tmass        fa(3,50) = fz * amass(50)/tmass        fa(3,51) = fz * amass(51)/tmass        fa(3,52) = fz * amass(52)/tmass        fa(3,53) = fz * amass(53)/tmass        fa(3,54) = fz * amass(54)/tmass        fa(3,55) = fz * amass(55)/tmass        fa(3,56) = fz * amass(56)/tmass        fa(3,57) = fz * amass(57)/tmass        fa(3,58) = fz * amass(58)/tmass        fa(3,59) = fz * amass(59)/tmass        fa(3,60) = fz * amass(60)/tmass         fa(1,1) = 0.0d0        fa(1,2) = 0.0d0        fa(1,3) = 0.0d0        fa(1,4) = 0.0d0        fa(1,5) = 0.0d0        fa(1,6) = 0.0d0        fa(1,7) = 0.0d0        fa(1,8) = 0.0d0        fa(1,9) = 0.0d0        fa(1,10) = 0.0d0        fa(1,11) = 0.0d0        fa(1,12) = 0.0d0        fa(1,13) = 0.0d0        fa(1,14) = 0.0d0        fa(1,15) = 0.0d0        fa(1,16) = 0.0d0        fa(1,17) = 0.0d0        fa(1,18) = 0.0d0        fa(1,19) = 0.0d0        fa(1,20) = 0.0d0        fa(1,21) = 0.0d0        fa(1,22) = 0.0d0        fa(1,23) = 0.0d0        fa(1,24) = 0.0d0        fa(1,25) = 0.0d0        fa(1,26) = 0.0d0        fa(1,27) = 0.0d0        fa(1,28) = 0.0d0        fa(1,29) = 0.0d0        fa(1,30) = 0.0d0        fa(1,31) = 0.0d0        fa(1,32) = 0.0d0        fa(1,33) = 0.0d0        fa(1,34) = 0.0d0        fa(1,35) = 0.0d0        fa(1,36) = 0.0d0        fa(1,37) = 0.0d0        fa(1,38) = 0.0d0        fa(1,39) = 0.0d0        fa(1,40) = 0.0d0        fa(1,41) = 0.0d0        fa(1,42) = 0.0d0        fa(1,43) = 0.0d0        fa(1,44) = 0.0d0        fa(1,45) = 0.0d0        fa(1,46) = 0.0d0        fa(1,47) = 0.0d0        fa(1,48) = 0.0d0        fa(1,49) = 0.0d0        fa(1,50) = 0.0d0        fa(1,51) = 0.0d0        fa(1,52) = 0.0d0        fa(1,53) = 0.0d0        fa(1,54) = 0.0d0        fa(1,55) = 0.0d0        fa(1,56) = 0.0d0        fa(1,57) = 0.0d0        fa(1,58) = 0.0d0        fa(1,59) = 0.0d0        fa(1,60) = 0.0d0         fa(2,1) = 0.0d0        fa(2,2) = 0.0d0        fa(2,3) = 0.0d0        fa(2,4) = 0.0d0        fa(2,5) = 0.0d0        fa(2,6) = 0.0d0        fa(2,7) = 0.0d0        fa(2,8) = 0.0d0        fa(2,9) = 0.0d0        fa(2,10) = 0.0d0        fa(2,11) = 0.0d0        fa(2,12) = 0.0d0        fa(2,13) = 0.0d0        fa(2,14) = 0.0d0        fa(2,15) = 0.0d0        fa(2,16) = 0.0d0        fa(2,17) = 0.0d0        fa(2,18) = 0.0d0        fa(2,19) = 0.0d0        fa(2,20) = 0.0d0        fa(2,21) = 0.0d0        fa(2,22) = 0.0d0        fa(2,23) = 0.0d0        fa(2,24) = 0.0d0        fa(2,25) = 0.0d0        fa(2,26) = 0.0d0        fa(2,27) = 0.0d0        fa(2,28) = 0.0d0        fa(2,29) = 0.0d0        fa(2,30) = 0.0d0        fa(2,31) = 0.0d0        fa(2,32) = 0.0d0        fa(2,33) = 0.0d0        fa(2,34) = 0.0d0        fa(2,35) = 0.0d0        fa(2,36) = 0.0d0        fa(2,37) = 0.0d0        fa(2,38) = 0.0d0        fa(2,39) = 0.0d0        fa(2,40) = 0.0d0        fa(2,41) = 0.0d0        fa(2,42) = 0.0d0        fa(2,43) = 0.0d0        fa(2,44) = 0.0d0        fa(2,45) = 0.0d0        fa(2,46) = 0.0d0        fa(2,47) = 0.0d0        fa(2,48) = 0.0d0        fa(2,49) = 0.0d0        fa(2,50) = 0.0d0        fa(2,51) = 0.0d0        fa(2,52) = 0.0d0        fa(2,53) = 0.0d0        fa(2,54) = 0.0d0        fa(2,55) = 0.0d0        fa(2,56) = 0.0d0        fa(2,57) = 0.0d0        fa(2,58) = 0.0d0        fa(2,59) = 0.0d0        fa(2,60) = 0.0d0        fz = fa(3,151)+fa(3,152)+fa(3,153)+fa(3,154)+fa(3,155)+fa(3,156)       $+fa(3,157)+fa(3,158)+fa(3,159)+fa(3,160)       $+fa(3,161)+fa(3,162)+fa(3,163)+fa(3,164)+fa(3,165)+fa(3,166)       $+fa(3,167)+fa(3,168)+fa(3,169)+fa(3,170)       $+fa(3,171)+fa(3,172)+fa(3,173)+fa(3,174)+fa(3,175)+fa(3,176)       $+fa(3,177)+fa(3,178)+fa(3,179)+fa(3,180)       $+fa(3,181)+fa(3,182)+fa(3,183)+fa(3,184)+fa(3,185)+fa(3,186)       $+fa(3,187)+fa(3,188)+fa(3,189)+fa(3,190)       $+fa(3,191)+fa(3,192)+fa(3,193)+fa(3,194)+fa(3,195)+fa(3,196)       $+fa(3,197)+fa(3,198)+fa(3,199)+fa(3,200)       $+fa(3,201)+fa(3,202)+fa(3,203)+fa(3,204)+fa(3,205)+fa(3,206)       $+fa(3,207)+fa(3,208)+fa(3,209)+fa(3,210)        tmass=amass(151)+amass(152)+amass(153)+amass(154)+amass(155)       $ +amass(156)+amass(157)+amass(158)+amass(159)+amass(160)       $ +amass(161)+amass(162)+amass(163)+amass(164)+amass(165)       $ +amass(166)+amass(167)+amass(168)+amass(169)+amass(170)       $ +amass(171)+amass(172)+amass(173)+amass(174)+amass(175)       $ +amass(176)+amass(177)+amass(178)+amass(179)+amass(180)       $ +amass(181)+amass(182)+amass(183)+amass(184)+amass(185)       $ +amass(186)+amass(187)+amass(188)+amass(189)+amass(190)       $ +amass(191)+amass(192)+amass(193)+amass(194)+amass(195)       $ +amass(196)+amass(197)+amass(198)+amass(199)+amass(200)       $ +amass(201)+amass(202)+amass(203)+amass(204)+amass(205)       $ +amass(206)+amass(207)+amass(208)+amass(209)+amass(210)         fa(3,151) = fz * amass(151)/tmass        fa(3,152) = fz * amass(152)/tmass        fa(3,153) = fz * amass(153)/tmass        fa(3,154) = fz * amass(154)/tmass        fa(3,155) = fz * amass(155)/tmass        fa(3,156) = fz * amass(156)/tmass        fa(3,157) = fz * amass(157)/tmass        fa(3,158) = fz * amass(158)/tmass        fa(3,159) = fz * amass(159)/tmass        fa(3,160) = fz * amass(160)/tmass        fa(3,161) = fz * amass(161)/tmass        fa(3,162) = fz * amass(162)/tmass        fa(3,163) = fz * amass(163)/tmass        fa(3,164) = fz * amass(164)/tmass        fa(3,165) = fz * amass(165)/tmass        fa(3,166) = fz * amass(166)/tmass        fa(3,167) = fz * amass(167)/tmass        fa(3,168) = fz * amass(168)/tmass        fa(3,169) = fz * amass(169)/tmass        fa(3,170) = fz * amass(170)/tmass        fa(3,171) = fz * amass(171)/tmass        fa(3,172) = fz * amass(172)/tmass        fa(3,173) = fz * amass(173)/tmass        fa(3,174) = fz * amass(174)/tmass        fa(3,175) = fz * amass(175)/tmass        fa(3,176) = fz * amass(176)/tmass        fa(3,177) = fz * amass(177)/tmass        fa(3,178) = fz * amass(178)/tmass        fa(3,179) = fz * amass(179)/tmass        fa(3,180) = fz * amass(180)/tmass        fa(3,181) = fz * amass(181)/tmass        fa(3,182) = fz * amass(182)/tmass        fa(3,183) = fz * amass(183)/tmass        fa(3,184) = fz * amass(184)/tmass        fa(3,185) = fz * amass(185)/tmass        fa(3,186) = fz * amass(186)/tmass        fa(3,187) = fz * amass(187)/tmass        fa(3,188) = fz * amass(188)/tmass        fa(3,189) = fz * amass(189)/tmass        fa(3,190) = fz * amass(190)/tmass        fa(3,191) = fz * amass(191)/tmass        fa(3,192) = fz * amass(192)/tmass        fa(3,193) = fz * amass(193)/tmass        fa(3,194) = fz * amass(194)/tmass        fa(3,195) = fz * amass(195)/tmass        fa(3,196) = fz * amass(196)/tmass        fa(3,197) = fz * amass(197)/tmass        fa(3,198) = fz * amass(198)/tmass        fa(3,199) = fz * amass(199)/tmass        fa(3,200) = fz * amass(200)/tmass        fa(3,201) = fz * amass(201)/tmass        fa(3,202) = fz * amass(202)/tmass        fa(3,203) = fz * amass(203)/tmass        fa(3,204) = fz * amass(204)/tmass        fa(3,205) = fz * amass(205)/tmass        fa(3,206) = fz * amass(206)/tmass        fa(3,207) = fz * amass(207)/tmass        fa(3,208) = fz * amass(208)/tmass        fa(3,209) = fz * amass(209)/tmass        fa(3,210) = fz * amass(210)/tmass         fa(1,151) = 0.0d0        fa(1,152) = 0.0d0        fa(1,153) = 0.0d0        fa(1,154) = 0.0d0        fa(1,155) = 0.0d0        fa(1,156) = 0.0d0        fa(1,157) = 0.0d0        fa(1,158) = 0.0d0        fa(1,159) = 0.0d0        fa(1,160) = 0.0d0        fa(1,161) = 0.0d0        fa(1,162) = 0.0d0        fa(1,163) = 0.0d0        fa(1,164) = 0.0d0        fa(1,165) = 0.0d0        fa(1,166) = 0.0d0        fa(1,167) = 0.0d0        fa(1,168) = 0.0d0        fa(1,169) = 0.0d0        fa(1,170) = 0.0d0        fa(1,171) = 0.0d0        fa(1,172) = 0.0d0        fa(1,173) = 0.0d0        fa(1,174) = 0.0d0        fa(1,175) = 0.0d0        fa(1,176) = 0.0d0        fa(1,177) = 0.0d0        fa(1,178) = 0.0d0        fa(1,179) = 0.0d0        fa(1,180) = 0.0d0        fa(1,181) = 0.0d0        fa(1,182) = 0.0d0        fa(1,183) = 0.0d0        fa(1,184) = 0.0d0        fa(1,185) = 0.0d0        fa(1,186) = 0.0d0        fa(1,187) = 0.0d0        fa(1,188) = 0.0d0        fa(1,189) = 0.0d0        fa(1,190) = 0.0d0        fa(1,191) = 0.0d0        fa(1,192) = 0.0d0        fa(1,193) = 0.0d0        fa(1,194) = 0.0d0        fa(1,195) = 0.0d0        fa(1,196) = 0.0d0        fa(1,197) = 0.0d0        fa(1,198) = 0.0d0        fa(1,199) = 0.0d0        fa(1,200) = 0.0d0        fa(1,201) = 0.0d0        fa(1,202) = 0.0d0        fa(1,203) = 0.0d0        fa(1,204) = 0.0d0        fa(1,205) = 0.0d0        fa(1,206) = 0.0d0        fa(1,207) = 0.0d0        fa(1,208) = 0.0d0        fa(1,209) = 0.0d0        fa(1,210) = 0.0d0         fa(2,151) = 0.0d0        fa(2,152) = 0.0d0        fa(2,153) = 0.0d0        fa(2,154) = 0.0d0        fa(2,155) = 0.0d0        fa(2,156) = 0.0d0        fa(2,157) = 0.0d0        fa(2,158) = 0.0d0        fa(2,159) = 0.0d0        fa(2,160) = 0.0d0        fa(2,161) = 0.0d0        fa(2,162) = 0.0d0        fa(2,163) = 0.0d0        fa(2,164) = 0.0d0        fa(2,165) = 0.0d0        fa(2,166) = 0.0d0        fa(2,167) = 0.0d0        fa(2,168) = 0.0d0        fa(2,169) = 0.0d0        fa(2,170) = 0.0d0        fa(2,171) = 0.0d0        fa(2,172) = 0.0d0        fa(2,173) = 0.0d0        fa(2,174) = 0.0d0        fa(2,175) = 0.0d0        fa(2,176) = 0.0d0        fa(2,177) = 0.0d0        fa(2,178) = 0.0d0        fa(2,179) = 0.0d0        fa(2,180) = 0.0d0        fa(2,181) = 0.0d0        fa(2,182) = 0.0d0        fa(2,183) = 0.0d0        fa(2,184) = 0.0d0        fa(2,185) = 0.0d0        fa(2,186) = 0.0d0        fa(2,187) = 0.0d0        fa(2,188) = 0.0d0        fa(2,189) = 0.0d0        fa(2,190) = 0.0d0        fa(2,191) = 0.0d0        fa(2,192) = 0.0d0        fa(2,193) = 0.0d0        fa(2,194) = 0.0d0        fa(2,195) = 0.0d0        fa(2,196) = 0.0d0        fa(2,197) = 0.0d0        fa(2,198) = 0.0d0        fa(2,199) = 0.0d0        fa(2,200) = 0.0d0        fa(2,201) = 0.0d0        fa(2,202) = 0.0d0        fa(2,203) = 0.0d0        fa(2,204) = 0.0d0        fa(2,205) = 0.0d0        fa(2,206) = 0.0d0        fa(2,207) = 0.0d0        fa(2,208) = 0.0d0        fa(2,209) = 0.0d0        fa(2,210) = 0.0d0         ntcon=358        end subroutine 

errors:

c ****************************** *********************************** 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:13:  c user-written routine implement specific geometric constraints, 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:14:  c orthogonalizing forces , stress undesired changes. 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:15:  c arguments: 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:16:  c real*8  cell(3,3)    : input lattice vectors (bohr) 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:17:  c integer na           : input number of atoms 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:18:  c integer isa(na)      : input species indexes 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:19:  c real*8  amass(na)    : input atomic masses 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:20:  c real*8  xa(3,na)     : input atomic cartesian coordinates (bohr) 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:21:  c real*8  stress( 3,3) : input/output stress tensor (ry/bohr**3) 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:22:  c real*8  fa(3,na)     : input/output atomic forces (ry/bohr) 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:23:  c integer ntcon        : total number of positions constr. imposed 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:24:  c ***************************************************************** 1 error: unclassifiable statement @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:27.55:        double precision  amass(na), cell(3,3), fa(3,na),                                                        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:32.7:        $ +fa(3,8)+fa(3,9)+fa(3,10)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:33.7:        $ +fa(3,11)+fa(3,12)+fa(3,13)+fa(3,14)+fa(3,15)+fa(3,16)+fa(3,17)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:34.7:        $ +fa(3,18)+fa(3,19)+fa(3,20)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:35.7:        $ +fa(3,21)+fa(3,22)+fa(3,23)+fa(3,24)+fa(3,25)+fa(3,26)+fa(3,27)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:36.7:        $ +fa(3,28)+fa(3,29)+fa(3,30)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:37.7:        $ +fa(3,31)+fa(3,32)+fa(3,33)+fa(3,34)+fa(3,35)+fa(3,36)+fa(3,37)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:38.7:        $ +fa(3,38)+fa(3,39)+fa(3,40)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:39.7:        $ +fa(3,41)+fa(3,42)+fa(3,43)+fa(3,44)+fa(3,45)+fa(3,46)+fa(3,47)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:40.7:        $ +fa(3,48)+fa(3,49)+fa(3,50)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:41.7:        $ +fa(3,51)+fa(3,52)+fa(3,53)+fa(3,54)+fa(3,55)+fa(3,56)+fa(3,57)        1 error: invalid character in name @ (1) /home/salami/constr_siesta/siesta-3.2/src/constr.f90:42.7:        $ +fa(3,58)+fa(3,59)+fa(3,60)        1 error: invalid character in name @ (1) fatal error: error count reached limit of 25. make: *** [constr.o] error 1 

could me solve problem, please?

your program fixed form fortran source, have in file suffix .f90. suffix used free-form sources. rename file suffix .f. option pass gfortran option -ffixed-form tells consider file in fixed-form regardless of file suffix.

be sure have line continuations $ @ sixth column, otherwise error.