subroutine writed(modvis2,starno) C C C. A. Hummel 27 July 1993 C include 'constnts.inc' include 'data.inc' C real*8 modvis2(mxrec,mxfilt) integer*4 firstid,lastid write(outc,'(a,$)') & ' Write observed or model visibilities ? (1/0) ' read(inc,*)icom open(7,file='binafit.dat') firstid=1 lastid=1 idate=1 zero=0.1 c c Go through the whole data file and ... c DO ID=1,NDATA c c ... work on single night, when date or baseline changes c if(date(id+1).ne.date(id).or.baseline(id+1).ne.baseline(id) & .or.id+1.gt.ndata)then write(7,1)date(id),baseline(id), & lambda0(id,1),lambda0(id,3),lambda0(id,2) 1 format(3x,i6,'.',3x,i2,2x,f5.1,17x,f5.1,17x,f5.1) do n=1,lastid-firstid+1 if(icom.eq.1)then write(7,2)hours(n+firstid-1),starno, & obsvis2(n+firstid-1,1), & obsvis2(n+firstid-1,1),viserr(n+firstid-1,1), & obsvis2(n+firstid-1,3), & obsvis2(n+firstid-1,3),viserr(n+firstid-1,3), & obsvis2(n+firstid-1,2), & obsvis2(n+firstid-1,2),viserr(n+firstid-1,2) else write(7,2)hours(n+firstid-1),starno, & modvis2(n+firstid-1,1), & modvis2(n+firstid-1,1),zero, & modvis2(n+firstid-1,3), & modvis2(n+firstid-1,3),zero, & modvis2(n+firstid-1,2), & modvis2(n+firstid-1,2),zero 2 format(1x,f9.6,1x,i4,3(3(1x,f6.3),1x)) endif enddo c c Now prepare for the next date idate=idate+1 firstid=lastid+1 lastid=firstid-1 c c continue here, if date did not change endif lastid=lastid+1 ENDDO close(7) return end