From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- noao/onedspec/fortran/trans.f | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) create mode 100644 noao/onedspec/fortran/trans.f (limited to 'noao/onedspec/fortran/trans.f') diff --git a/noao/onedspec/fortran/trans.f b/noao/onedspec/fortran/trans.f new file mode 100644 index 00000000..038384ba --- /dev/null +++ b/noao/onedspec/fortran/trans.f @@ -0,0 +1,21 @@ + SUBROUTINE TRANS(Y,A,X) + DIMENSION Y(10), A(10), X(10) + COMMON /NLCPAR/XC(10), N, FIXSEP + LOGICAL FIXSEP +C----- TRANSOFRMATION FOR GAUSSIAN LINES +C +C----- 'N' GAUSSIAN LINES +C + Y(1)=EXP(-0.5*((X(1)-XC(1)-A(2))/A(1))**2) + DO 1000 I=2,N + IF(FIXSEP) THEN + DELTA=A(2) + ELSE + DELTA=A(2*I) + ENDIF + Y(1)=Y(1)+ABS(A(2*I-1)*EXP(-0.5*((X(1)-XC(I)-DELTA)/ + * A(1))**2)) +1000 CONTINUE +C + RETURN + END -- cgit