diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
| commit | fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch) | |
| tree | bdda434976bc09c864f2e4fa6f16ba1952b1e555 /noao/onedspec/fortran/trans.f | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'noao/onedspec/fortran/trans.f')
| -rw-r--r-- | noao/onedspec/fortran/trans.f | 21 |
1 files changed, 21 insertions, 0 deletions
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 |