; routine to mock up a VOTable containing Chandra and EIS attributes for matched sources

; read in the two input catalogues

; first the EIS data

fmt1='(a1,4x,a20,2x,f5.3,1x,a11,1x,a11,1x,i6,1x,f7.4,1x,f6.4,1x,f7.4,1x,f6.4,1x,f7.4,1x,f6.4,1x,f7.4,1x,f6.4,1x,f7.4,1x,f6.4,1x,f7.4,1x,f6.4,1x,f7.4,1x,f6.4,1x,i1,1x,f10.7,1x,f11.7)'

readfmt,'/home/rgm/idl/my_idl/vo/chandra_eis_new.txt',fmt1,rank,matchstring1,dist,eis_ra,eis_dec,eis_catid,umag,e_umag,bmag,e_bmag,$
vmag,e_vmag,rmag,e_rmag,imag,e_imag,jmag,e_jmag,kmag,e_kmag,ext_flag,eis_raj2000,eis_dej2000

; then the Chandra data

fmt2='(a20,3x,f9.5,1x,f10.6,2x,f5.3,1x,a11,1x,a11,1x,i3,1x,i3,1x,a24,1x,e7.1,1x,e7.1,1x,e7.1,1x,e7.1,1x,f5.2,1x,f4.2,1x,f4.2)'

readfmt,'/home/rgm/idl/my_idl/vo/chandra_chandra.txt',fmt2,matchstring2,chandra_raj2000,chandra_dej2000,xdist,$
chandra_ra,chandra_dec,recno,chandra_catid,chandra_radec_orig,sflux,e_sflux,hflux,e_hflux,hr,e_hr,poserr

; check that matching strings do match
;.r
;for i=0,n_elements(matchstring1)-1 do begin
;	if (matchstring1[i] ne matchstring2[i]) then print,i,' does not match'
;endfor

; looks OK

; remove the VOTable file (if it's already there)
spawn,'rm /home/rgm/idl/my_idl/vo/chandra_eis_vot.xml'

; open the VOTable file and write in all the header info

openw,u,'/home/rgm/idl/my_idl/vo/chandra_eis_vot.xml',/get_lun

printf,u,'<?xml version="1.0" encoding="UTF-8"?>'
printf,u,'<VOTABLE xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://vizier.u-strasbg.fr/xml/VOTable.xsd">'
printf,u,'<DESCRIPTION />'
printf,u,'<DEFINITIONS>'
printf,u,'  <COOSYS ID="J2000" system="eq_FK5" equinox="J2000"/>'
printf,u,'</DEFINITIONS>'
printf,u,'<RESOURCE ID="XXXXXX" name="nnnnnn">'
printf,u,'  <DESCRIPTION>Chandra Deep Field South. 1 Ms catalog (Giacconi+, 2002)</DESCRIPTION>'
printf,u,'  <TABLE ID="J_ApJS_139_369_table2" name="J/ApJS/139/369/table2">'
printf,u,'    <DESCRIPTION>Main Catalog</DESCRIPTION>'
printf,u,'    <FIELD name="ID" ucd="ID_NUMBER" datatype="short">'
printf,u,'      <DESCRIPTION>Unique Detection identification number</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Chandra_RAJ2000" ucd="POS_EQ_RA_MAIN" ref="J2000" datatype="char"  unit="&quot;h:m:s&quot;">'
printf,u,'      <DESCRIPTION>Right ascension (FK5) Equinox=J2000. (computed by VizieR, not part of the original data)</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Chandra_DEJ2000" ucd="POS_EQ_DEC_MAIN" ref="J2000" datatype="char" unit="&quot;d:m:s&quot;">'
printf,u,'      <DESCRIPTION>Declination (FK5) Equinox=J2000. (computed by VizieR, not part of the original data)</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="_RAJ2000_Chandra" ucd="POS_EQ_RA_MAIN" ref="J2000" datatype="float"  unit="deg">'
printf,u,'    <DESCRIPTION>Right ascension (FK5) Equinox=J2000. (computed by VizieR, not part of the original data)</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="_DEJ2000_Chandra" ucd="POS_EQ_DEC_MAIN" ref="J2000" datatype="double"  unit="deg">'
printf,u,'    <DESCRIPTION>Declination (FK5) Equinox=J2000. (computed by VizieR, not part of the original data)</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="SFlux" ucd="PHOT_FLUX_X" datatype="float" unit="mW/m2">'
printf,u,'      <DESCRIPTION>Soft band flux in units of erg/cm^2^/s</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_SFlux" ucd="ERROR" datatype="float"  unit="mW/m2">'
printf,u,'      <DESCRIPTION>? Error in SFlux</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="HFlux" ucd="PHOT_FLUX_X" datatype="float" unit="mW/m2">'
printf,u,'      <DESCRIPTION>Hard band flux in units of erg/cm^2^/s</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_HFlux" ucd="ERROR" datatype="float"  unit="mW/m2">'
printf,u,'      <DESCRIPTION>? Error in HFlux</DESCRIPTION>'
printf,u,'   </FIELD>'
printf,u,'    <FIELD name="HR" ucd="SPECT_HARDNESS-RATIO" datatype="float">'
printf,u,'      <DESCRIPTION>Hardness ratio (5)</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_HR" ucd="ERROR" datatype="float">'
printf,u,'      <DESCRIPTION>Error in HR</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="PosErr" ucd="ERROR" datatype="float" unit="arcsec">'
printf,u,'      <DESCRIPTION>Error on the X-ray source position as calculated in the text</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Match_dist" ucd="POS_ANG_DIST_GENERAL" datatype="float" unit="arcsec">'
printf,u,'    <DESCRIPTION>Distance from between Chandra and EIS position</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="EIS_RAJ2000" ucd="POS_EQ_RA_MAIN" ref="J2000" datatype="char"  unit="&quot;h:m:s&quot;">'
printf,u,'      <DESCRIPTION>Right ascension (FK5) Equinox=J2000. (computed by VizieR, not part of the original data)</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="EIS_DEJ2000" ucd="POS_EQ_DEC_MAIN" ref="J2000" datatype="char" unit="&quot;d:m:s&quot;">'
printf,u,'      <DESCRIPTION>Declination (FK5) Equinox=J2000. (computed by VizieR, not part of the original data)</DESCRIPTION>'
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="_RAJ2000_EIS" ucd="POS_EQ_RA_MAIN" ref="J2000" datatype="float"  unit="deg">'
printf,u,'    <DESCRIPTION>Right ascension (FK5) Equinox=J2000. (computed by VizieR, not part of the original data)</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="_DEJ2000_EIS" ucd="POS_EQ_DEC_MAIN" ref="J2000" datatype="double"  unit="deg">'
printf,u,'    <DESCRIPTION>Declination (FK5) Equinox=J2000. (computed by VizieR, not part of the original data)</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Umag" ucd="PHOT_JHN_U" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>U magnitude in Vega system</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_Umag" ucd="ERROR" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>? error in Umag</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Bmag" ucd="PHOT_JHN_B" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>B magnitude in Vega system</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_Bmag" ucd="ERROR" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>? error in Bmag</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Vmag" ucd="PHOT_JHN_V" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>V magnitude in Vega system</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_Vmag" ucd="ERROR" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>? error in Vmag</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Rmag" ucd="PHOT_JHN_R" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>R magnitude in Vega system</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_Rmag" ucd="ERROR" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>? error in Rmag</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Imag" ucd="PHOT_JHN_I" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>I magnitude in Vega system</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_Imag" ucd="ERROR" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>? error in Imag</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Jmag" ucd="PHOT_JHN_J" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>J magnitude in Vega system</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_Jmag" ucd="ERROR" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>? error in Jmag</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="Kmag" ucd="PHOT_JHN_K" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>K magnitude in Vega system</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,'    <FIELD name="e_Kmag" ucd="ERROR" datatype="float"  unit="mag">'
printf,u,'    <DESCRIPTION>? error in Kmag</DESCRIPTION>' 
printf,u,'    </FIELD>'
printf,u,' '
printf,u,'<DATA><TABLEDATA>'

for i=0,n_elements(matchstring1)-1 do begin

	if (rank[i] eq 1) then begin

		printf,u,'      <TR>'
		str='<TD>'+string(i+1)+'</TD>'
		printf,u,'         ',strcompress(str,/rem)
		printf,u,'         <TD>',chandra_ra[i],'</TD>'
		printf,u,'         <TD>',chandra_dec[i],'</TD>'
		str='<TD>'+string(chandra_raj2000[i])+'</TD>'
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(chandra_dej2000[i])+'</TD>'
		printf,u,'         ',strcompress(str,/rem)	
		str='<TD>'+string(sflux[i])+'</TD>'
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(e_sflux[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(hflux[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(e_hflux[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(hr[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(e_hr[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(poserr[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(dist[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		printf,u,'         <TD>',eis_ra[i],'</TD>'	
		printf,u,'         <TD>',eis_dec[i],'</TD>'	
		str='<TD>'+string(eis_raj2000[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(eis_dej2000[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(umag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)	
		str='<TD>'+string(e_umag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(bmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(e_bmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(vmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(e_vmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(rmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(e_rmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(imag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(e_imag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(jmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(e_jmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(kmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
		str='<TD>'+string(e_kmag[i])+'</TD>'	
		printf,u,'         ',strcompress(str,/rem)
 		printf,u,'      </TR>'

	endif

endfor

printf,u,'</TABLEDATA>'
printf,u,'</DATA>'
printf,u,'  </TABLE>'
printf,u,'</RESOURCE>'
printf,u,'</VOTABLE>'


close,u & free_lun,u

end