FUNCTION radec2deg, ra_cent, dec_cent, ra, dec
;Coordonnees du centre du champ transformees en radians
x0 = ra_cent*!PI/180.                
y0 = dec_cent*!PI/180.             
;$x0=2.*pi*(13.+ 6./60. + 11.0/3600.)/24.
;$y0=2.*pi*(29.+ 1./60. + 56.0/3600.)/360.;
		  
;Coordonnees des sources transformees en radians
if ((size(ra))[0] ne 0) then ra = reform(ra)
if ((size(dec))[0] ne 0) then dec = reform(dec)
  
X = !PI*ra/180.                 
Y = !PI*dec/180.                
		  
;$X =2*pi*($HH + $HM/60. + $HSEC/3600.)/24.;
;$Y =2*pi*($sDD + $DM/60. + $DS/3600.)/360.;
		  
DA  = X-x0                         
D1  = sin(y0)*sin(Y)               
D2  = cos(y0)*cos(Y)*cos(DA)       
DET = D1+D2                      
y1  = cos(y0)*sin(Y)               
y2  = sin(y0)*cos(Y)*cos(DA)       
	
;coordonnees a tracer, en radians
xplot = atan(cos(Y)*sin(-DA),DET) 
yplot = atan((y1-y2),DET)  
      
size_xplot = size(xplot)
size_yplot = size(yplot)
;print, 'xplot',size_xplot[1]

if (size_xplot[0] eq 0) then size_ra = 1 $ ;scalaire
  else size_ra = size_xplot[1] ;vecteur
if (size_yplot[0] eq 0) then size_dec = 1 $ ;scalaire
  else size_dec = size_yplot[1] ;vecteur

if (size_ra eq 1 and size_dec eq 1) then size_out = 1

if (size_ra ne size_dec) then size_out = max([size_ra,size_dec])
if (size_ra eq size_dec and size_ra ne 1) then size_out = size_ra

resultats      = dblarr(2,size_out)
resultats[0,*] = xplot*180./!PI
resultats[1,*] = yplot*180./!PI
; les coordonnees sont ici renvoyees en degres
return, resultats	

END