C23456===============================================
      subroutine yyddmm_jul(iy,im,id,jd)

C Esta rutina lee el dia mes anho y devuelve el dia 
C juliano. El algoritmo de conversion lo baje de la 
C pagina de Peter Meyer (http://www.hermetic.ch/cal_stud/jdn.htm)
C
C KOBI MOSQUERA
C kobi@chavin.igp.gob.pe

      integer iy,im,id,jd

      jd=( 1461*(iy+4800+ (im-14) / 12 ) ) /4 +
     & ( 367* (im-2-12*((im-14)/12)) )/12 -
     & ( 3*((iy + 4900 + (im-14)/12)/100))/ 4 +
     & id -32075

      end