[Numpy-discussion] Converting from FORTRAN Equivalence to a Numpy alternative

Rob rob at pythonemproject.com
Sat Feb 16 10:04:03 EST 2002


I hoped there might be a combo FORTRAN/Numpy hacker that might help me
with this.  I have the following statements in a FORTRAN routine that I
am converting to Numpy.  I am wondering how to handle the EQUIVALENCE
statement.  As I understand it, ARL1[1] should be equal to the first
item in the memory block for AR1.  But not knowing how FORTRAN allocates
arrays, its hard to tell what to do in Numpy.

I tried ARL1=ravel(AR1) but that didn't work.  Similarly ARL=AR[:,1,1}
didn't work.  I'm lost.

Thanks in advance for any help.  Rob.


    #COMMON /GGRID/ AR1(11,10,4),AR2(17,5,4),AR3(9,8,4),EPSCF,DXA(3),&
    #            &DYA(3),XSA(3),YSA(3),NXA(3),NYA(3

    #DIMENSION ARL1(1), ARL2(1), ARL3(1)
    #EQUIVALENCE (ARL1,AR1), (ARL2,AR2), (ARL3,AR3), (XS2,XSA(2)),
(YS3,YSA(3))


-- 
The Numeric Python EM Project

www.pythonemproject.com





More information about the NumPy-Discussion mailing list