Día Juliano en Fortran

Desarrollado para el curso de Geofísica. Calcula el día juliano para: TANNO, TMES, TDIA, THORA, TMINUTO.

program GREJD
integer TMES, TANNO, TDIA, JDIA
real THORA,TMINUTO
TMINUTO=25
THORA=5                
TDIA=28                
TMES=7            
TANNO=1821
THORA=(THORA+(TMINUTO/60))/24+0.5
if (TMES > 2) then
   TMES=TMES-3                 
else
   TMES=TMES+9
TANNO=TANNO-1
endif
JDIA = (TANNO / 4000) * 1460969;
TANNO = MOD(TANNO,4000);
JDIA=JDIA +(((TANNO / 100)*146097)/4) +((MOD(TANNO,100) * 1461) / 4)+(((153 * TMES) + 2) / 5) +TDIA +1721119;
write (*,*) 'DIA JULIANO:'
PRINT *,'PARTE ENTERA',JDIA
PRINT *,'PARTE DECIMAL',dble(THORA)
end