Another problem is that a namelist that contains all the parameters is not global available to the subroutines.
I define and read the namelist in the main program as follows:
REAL :: alpha, vmaxNP, vmaxPZ, Ks, mP, mZ, Phalf, Pzero
INTEGER :: dt, tend
NAMELIST /NPZNAME/ alpha, vmaxNP, vmaxPZ...
Dear Forum,
I tried to use a Runge-Kutta method to solve a system of three ODEs (dNdt, dPdt, dZdt).
The code is splitted in the main program, the ODE subroutine and a solver subroutine.
I am very new to Fortran (90), therefore my problem could be very basic and resulting from misunderstanding...
This site uses cookies to help personalise content, tailor your experience and to keep you logged in if you register.
By continuing to use this site, you are consenting to our use of cookies.