I've written a fortran 90/95 program to do a parametric study over a
range of 3 variables. This results in a structured program with 4
embedded loops as I am doing a Newton-Raphson iterative scheme for one
of the variables. An outline of my program looks something like the
following:

program ...
! Parameters, variables, and arrays defined here

open(unit=35,file='...',status='...',position='...')

do i = 1, npts + 1
do j = 1, npts + 1
do k= 1, npts + 1
do n = 1, nsteps

! Define initial arrays
theta(i,j,k,n) = ...
T(i,j,k,n)= ...
tau(i,j,k,n) = ...
ftau(i,j,k,n) = ...

I can't post to the j3 mailing list, but there are a couple of articles
I'd like to respond to there. I hesitate to just post them here (though
their context is clear, some members may regard it as unfair to
publicize their internal machinations). I have no need to post
there the vast majority of the time. So I don't want to rock the
boat by requesting the privilege.

--
J. Giles

"I conclude that there are two ways of constructing a software
design: One way is to make it so simple that there are obviously
no deficiencies and the other way is to make it so complicated
that there are no obvious deficiencies."

Hi,

what is the difference between the MOD and the MODULO functions in
Fortran 90? According to the book Fortran 90-95 Explained,
mod(a,p) returns a-int(a/p)*p and
modulo(a,p) returns a-floor(a/p)*p if a and p are real,
or a-floor(a//p)*p if a and p are integers,
where // is "ordinary mathematical division".

I don't see what the difference is. In my tests, both functions always
return the same value.

Thanks,
Michael

CAR CHEAP INSURANCE
http://carcheapinsurance.googlepages.com