LES ETANGS PALAVASIENS :

MODELISATION D'UN SYSTEME ECOLOGIQUE

Méthode de Runge-Kutta : Programme Fortran


C Last change: R&V 9 Feb 99 1:36 pm

program BIO

IMPLICIT NONE

INTEGER i,j,tf,MAX

REAL*8 dt,tfinal

PARAMETER(dt=100.,MAX=3)

REAL*8 TAB1,TAB2

dimension TAB1(MAX),TAB2(MAX)

C Constitution du fichier resultat

OPEN(1,FILE='c:\Robinson\biotruc.dat')

C Initialisation de la condition initiale

TAB1(1)=64.E-1

TAB1(2)=5.

TAB1(3)=4.

 

C Debut de la boucle temporelle

WRITE(*,*) 'ENTREZ LE TEMPS FINAL'

READ(*,*) tfinal

tf=INT(tfinal/dt)

 

C Affectation de boucle temporelle

do i=1,tf

TAB2(1)=dt*(3.E-6/(7.E-1)**2*(64E-1-TAB1(1))-273.E-7*TAB1(3)

& -52.E-1*7349.E-10*TAB1(2))+TAB1(1)

TAB2(2)=dt*(-7349.E-10*TAB1(2))+TAB1(2)

TAB2(3)=dt*(-273.E-7*TAB1(3))+TAB1(3)

do j=1,3

if (TAB2(j).GT.0) then

TAB1(j)=TAB2(j)

else

TAB1(j)=0.

endif

END do

WRITE(1,*) TAB1(1),TAB1(2),TAB1(3)

WRITE(*,*) TAB1(1),TAB1(2),TAB1(3)

end do

 

C Fermeture du fichier solution

 

CLOSE(1)

End