character*132 line character*40 flname,outname,imagefile,sptname,base,SETDEF character*3 number,oldnum,zero data zero/'000'/ real*4 oscstrt WRITE(6,*)'PROGRAM TO GENERATE REFIX COMMAND FILES' WRITE(6,*)' 930127' WRITE(6,*) write(6,*)' Enter basename of images in series (bcxx_):' read(5,11) base write(6,*)'base=',base write(6,*)' Enter numbers of first and last images in series:' read(5,*) n1,n2 11 format(a40) write(6,*)'enter oscillation CENTER PHI for first image:' read(5,*) oscstrt WRITE(6,*)'ENTER PHI INCREMENT:' READ(5,*)DPHI open(unit=3,file='stills.COM',status='new',carriagecontrol='list') SETDEF='$set def a:[berry.mosf]' write(3,*) SETDEF WRITE(3,*)'$assign c:[berry.dat] D' open(unit=4,file='PART4.COM',status='new',carriagecontrol='list') write(4,*) SETDEF write(4,*)'$run part4' write(4,*)'0.3 !partiality' open(unit=1,file='rrefix.COM',status='new',carriagecontrol='list') write(1,*) SETDEF do 150 i=n1,n2 encode(3,81,number) i 81 format(i3) j=0 101 j=j+1 if(number(j:j).eq.' ') goto 101 number=zero(:j-1)//number(j:) j=index(base,' ')-1 imagefile='D:'//base(:j)//number//'.image' C write(6,*)j,base(:j),number,imagefile SPtname=BASE(:4)//number(2:)//'.SPT' write(3,31) sptname 31 format('$ assign/user ',a30,' SPOTS') write(3,*)'$ run b:[vladimir.mosco]IMSTILLS' write(3,51) imagefile 51 format(' FILE ',A50) write(3,61) oscstrt+(i-n1)*DPHI 61 format(' PHI',f7.2) write(3,*)'PIXEL 0.150' write(3,*)'RADIUS 11.0 60.0' write(3,*)'SELECT MINPIXEL 6' write(3,*)'BINS OSC SPOTS 30' write(3,*)'THRESHOLD 10 ' write(3,*)'SITE DLAB 90.00 90.00' write(3,*)'GO' c write entries to part4.com WRITE(4,*) SPTNAME IF (I.GT.n1+1) THEN SPtname=BASE(:4)//OLDNUM(2:)//'PF.SPT' WRITE(4,*) SPTNAME c write entries to rrefix.com c if ((i.ne.n1).and.(i.ne.n2)) then write(1,*)'$ ASSIGN/USER '//sptname(:6)//'pf.MAT MATRIX' write(1,*)'$ ASSIGN/USER '//sptname(:6)//'pf.lp REFIX.LP' write(1,*)'$ ASSIGN/USER '//sptname(:6)//'pf.spt SPOTS' write(1,*)'$ @refix' c endif ENDIF OLDNUM=NUMBER 150 continue write(3,151) 151 format('$@part4') write(3,*)'$@rrefix' write(3,*)'$!!!!!del bc*.spt;*' write(3,*)'$del bc*.lp;*' write(3,*)'$@rdmat' write(3,*)'$!!!!!del bc*.mat;*' close (unit=3) WRITE(4,*) close (unit=4) close (unit=1) open(unit=4,file='rdmat.COM',status='new',carriagecontrol='list') write(4,*) SETDEF write(4,*)'$run rdmat' write(4,111) BASE(:4) 111 format(a4) WRITE(4,*) 'pf.mat' write(4,*) n1+1,n2-1 close (unit=4) end