Il Linguaggio Fortran 90/95

92 Istruzioni di controllo
READ(*,*) radicando
IF (radicando>=0.0) THEN
WRITE(*,*) " La sua radice quadrata vale: "
WRITE(*,*) SQRT(radicando)
ELSE
WRITE(*,*) " Errore: inserito valore negativo! "
END IF
STOP
END PROGRAM radq
Il programma seguente, invece, valuta il giorno dell’anno in cui cade la Pasqua sfruttando
l’algoritmo proposto da Oudin e Tondering. Al fine di comprendere le istruzioni contenute
nel programma giova ricordare che la Pasqua cade sempre la domenica seguente il plenilunio
successivo all’equinozio di primavera (21 marzo). Si ricorda, inoltre, che si chiama età della
luna a un dato momento, il numero di giorni trascorsi dall’ultimo novilunio fino a quel momento
(per tutti i giorni dell’anno si può dunque esprimere l’età della luna, in generale mediante un
numero intero di giorni). Si chiama, invece, epatta relativa a un determinato anno l’età della
luna al primo gennaio di quell’anno. Convenedo di chiamare 0 il giorno in cui la Luna è nuova,
l’epatta può variare tra 0 e 29. Infine, sovrapponendo al ciclo solare annuale il ciclo metonico
di 19 anni ne consegue che ogni 19 anni, prescindendo da errori di computo, i fenomeni celesti
solari e lunari devono ripetersi identicamente.
PROGRAM Pasqua
! Scopo: calcolare la data in cui cade la Pasqua in un dato anno
IMPLICIT NONE
INTEGER :: year, metcyc, century, error1, error2, day
INTEGER :: epact, luna, temp
!
WRITE(*,*) " Inserire l’anno (successivo al 1582): "
READ(*,*) year
! Calcolo del ciclo metonico
metcyc = MOD(year,19)+1
century = (year/100)+1
error1 = (3*century/4)-12
error2 = ((8*century+5)/25)-5
day = (5*Year/4)-error1-10
temp = 11*metcyc+20+error2-error1
! Calcolo dell’epatta
epact = MOD(temp,30)
IF(epact11).OR.epact==24) THEN
epact = epact+1