weird behaviour
Éric Depagne
edepagne@lcogt.net
Fri Aug 3 18:46:00 GMT 2007
Le Friday 03 August 2007 11:03:14 Martin Tee, vous avez écrit :
> Hi,
Hi Martin,
good catch!
I fixed this, and now it works.
Thanks a lot,
Éric.
>
> Looks like your arguments don't match; 9 actual and 10 dummy.
>
> Regards,
> Martin
>
> Éric Depagne wrote:
> > Le Thursday 02 August 2007 21:00:58 Brooks Moses, vous avez écrit :
> >> Éric Depagne wrote:
> >>> That's the way I call the function:
> >>>
> >>> contrib=LINE(wave,nlo,nup,hlambda,t(k),ne(k),nh1(k),nhe1(k),frac)
> >>
> >> How are these variables declared?
> >
> > Unfortunately, most of them are implicit.
> >
> >>> And that's the first lines of the function:
> >>> FUNCTION
> >>> LINE(WAVE,NBLO,NBUP,WAVEH,T,XNE,HNORM,H1FRC,HE1FRC,DOPPLE)
> >>>
> >>> REAL KAP0,KAP0RED,KAPRED,KAP0BLUE,KAPBLUE
> >>> REAL*8 WAVEH,WAVE,WCON,EMERGE,WMERGE,WSHIFT,WTAIL
> >>> REAL*8 BLUECUT,WLBLUE,REDCUT,WLRED,EMERGEH
> >>> REAL*8 EHYD(100),CONTH(15)
> >>> REAL*4 DOPPLE
> >>> REAL*8 VACAIR
> >>> real pg
> >>> LOGICAL FIRST
> >>> SAVE
> >>
> >> In general, it's a good idea to put IMPLICIT NONE at the beginning of
> >> all declaration blocks, to catch undeclared variables.
> >
> > I agree, but this code is not mine, and it is rougly 20000 lines long,
> > splitted in almost 90 files.
> > And some files are written in pure F77, thus rewriting them so everything
> > is readable is not ( unfrtunately) an option.
> >
> >> In particular, NBLO, NBUP, T, XNE, HNORM, H1FRC, and HE1FRC are not
> >> declared here, and thus will get an implicit type. My guess is that
> >> this does not match the type that they have in the caller, and that the
> >> mismatch is such that the function's location for DOPPLE is off the end
> >> of the call stack that the caller supplies.
> >>
> >> In any case, reducing this to a minimal complete example should enable
> >> you to find the problem. (Steve's suggestions for warning options may
> >> find the problem as well, though I'm not sure of that in this case.)
> >
> > I'll try his suggestion, hoping it will show the things missing.
> >
> > Éric.
> >
> >> - Brooks
--
Un clavier azerty en vaut deux
----------------------------------------------------------
Éric Depagne edepagne@lcogt.net
Las Cumbres Observatory
6740 Cortona Dr
Goleta CA, 93117
----------------------------------------------------------
More information about the Fortran
mailing list