corrected reading of beam description file (bad handling of blank lines and of case nalpha=2)

src/gray-externals.f

@@ -898,7 +898,7 @@ c@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
       implicit real*8(a-h,o-z)
 c     
       character*255 filenmeqq,filenmprf,filenmbm
-      character*20 beamname
+      character*20 beamname,beamnameskip
 c     
       integer beamid, nisteer, fdeg, jumprow, nbeam, nalpha, nbeta,
      . iopt, incheck, nxcoord, nycoord, nxest, nyest, lwrk, kwrk,
@@ -989,7 +989,7 @@ c     unused beams' data
         jumprow=0
 c  c====================================================================================
         do i=1,beamid-1
-          read(nfbeam,*) beamname, iox, fghz, nalpha, nbeta
+          read(nfbeam,*) beamnameskip, iox, fghz, nalpha, nbeta
           jumprow = jumprow+nalpha*nbeta
         end do
 c  c====================================================================================
@@ -1000,8 +1000,12 @@ c     beam of interest
 c       
 c  c====================================================================================
 c     unused beams' data grids
-        do i=1,(jumprow + (nbeam - beamid))
+        do i=1,(nbeam - beamid)
-          read(nfbeam,*)
+          read(nfbeam,*) beamnameskip
+        end do
+        do i=1,jumprow
+          read(nfbeam,*) alphast,betast,x00,y00,z00,
+     .     waist1,waist2,rci1,rci2,phi1,phi2
         end do
 c  c====================================================================================
 c       
@@ -1037,7 +1041,7 @@ c     fdeg = 0 alpha, beta free variables
 c            1 alpha free variable
 c            2 beta free variable
 c            3 no free variables
-        fdeg = 2*1/nalpha + 1/nbeta
+        fdeg = 2*(1/nalpha) + 1/nbeta
       end if
 c#######################################################################################
 c