c+
c name:		variables.for
c purpose:	This program lists stars which have been flagged with a "V*" either in "otype" 
c		or at the beginning of the indentifier
c inputs:	The four starsout files generated by stars2smei.for
c outputs:	Printed list
c calls:	
c mods:		Jan 2003, Andrew Buffington (UCSD; (858)-534-6630; abuffington@ucsd.edu)
c-
	program variables
c
c	
c
	integer*4 i,j,m,irahour,iramin,idecdeg,idecmin,idecsec,iarg,jarg,iv
	real*4 amagb,amagv,rasec,amagsmei
	character*13 starsout
	character*1 identifier(12),v(2)
	character*4 otype,vstar,eclipse(4)
	character*1 obafg,subtype,lum(6)
	data vstar /'  V*'/
	data eclipse /' EB*',' Al*',' bL*',' WU*'/
	v(1)='V'
	v(2)='*'
	iv=0
c
	do m=1,4
	   if(m.eq.1)starsout = 'n_poleout.dat'
	   if(m.eq.2)starsout = 'fieldout+.dat' 
	   if(m.eq.3)starsout = 'fieldout-.dat' 
	   if(m.eq.4)starsout = 's_poleout.dat'
	   open(11,file=starsout)
	   do i=1,3000
	      read(11,11)obafg,subtype,lum,otype,irahour,iramin,rasec,idecdeg,idecmin,idecsec,amagb,amagv,amagsmei,
     &	      iarg,jarg,identifier			
11	      format(8a1,2x,a4,2x,i2,1x,i2,1x,f5.2,2x,i3,1x,i2,1x,i2,8x,f5.2,6x,f5.2,6x,f5.2,2i4,2x,12a1)
	      if(irahour.eq.0.and.iramin.eq.0.and.rasec.eq.0.0)go to 99	!detects EOF
	      if(amagsmei.gt.4.0)go to 98
	      if(identifier(1).eq.v(1).and.identifier(2).eq.v(2))then
	        do j=1,4
	          if(otype.eq.eclipse(j))go to 98
	        enddo
	        iv=iv+1
	      	print 12,iv,identifier,otype,obafg,subtype,lum,irahour,iramin,rasec,idecdeg,idecmin,idecsec,amagb,amagv,amagsmei
     		go to 98
     	      endif
     	      if(otype.eq.vstar)then
     	        iv=iv+1
	      	print 12,iv,identifier,otype,obafg,subtype,lum,irahour,iramin,rasec,idecdeg,idecmin,idecsec,amagb,amagv,amagsmei
12		format(i4,2x,12a1,a4,2x,8a1,2x,i2,1x,i2,1x,f5.2,2x,i3,1x,i2,1x,i2,2x,3f5.2)
     	      endif
98	   continue     	      
     	   enddo
     	   close(11)
99	continue
	enddo
	print *,iv,'non-eclipsing variables brighter than 4th found'
	end