← Back to branch summary

~ubuntu-branches/debian/jessie/eso-midas/jessie

« back to all changes in this revision

Viewing changes to contrib/tsa/proc/tsaort.prg

  • Committer: Package Import Robot
  • Author(s): Ole Streicher
  • Date: 2014-04-22 14:44:58 UTC
  • Revision ID: package-import@ubuntu.com-20140422144458-okiwi1assxkkiz39
Tags: upstream-13.09pl1.2+dfsg
Import upstream version 13.09pl1.2+dfsg

Show diffs side-by-side

added added

removed removed

Lines of Context:
contrib/tsa/proc/tsaort.prg
 
1
! @(#)tsaort.prg        19.1 (ESO-IPG) 02/25/03 13:33:13
 
2
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
 
3
!.COPYRIGHT (c)  1992     European Southern Observatory and Warsaw Observatory
 
4
!.IDENT     tsaort.prg
 
5
!.AUTHOR    Alex Schwarzenberg-Czerny, ESO and Warsaw Observatory 
 
6
!.KEYWORD   TSA, MIDAS, Time Series Analysis
 
7
!.LANGUAGE  MIDAS
 
8
!.PURPOSE   Compute multiharmonic Fourier spectrum by orthogonal projections
 
9
!.VERSION   0.0               June 1992
 
10
!.RETURNS   None  
 
11
!.ENVIRON   TSA context
 
12
!-----------------------------------------------------------------------------
 
13
!
 
14
!   Get parameters
 
15
!
 
16
CROSSREF INTAB OUTIMA STARTTSA STEPTSA NSTEPS ORDERTSA  
 
17
! 
 
18
DEFINE/PARAM P1 {INROOT1}  TABLE  "Enter input table: "
 
19
DEFINE/PARAM P2 {OUTROOT}  IMAGE  "Enter output image: "
 
20
DEFINE/PARAM P3 {STARTTSA}    NUMBER "Enter frequency start: "
 
21
DEFINE/PARAM P4 {STEPTSA}     NUMBER "Enter frequency step: "
 
22
DEFINE/PARAM P5 {NSTEPS}   NUMBER "Enter number of frequencies: "
 
23
DEFINE/PARAM P6 {ORDERTSA}    NUMBER "Enter order (2*no of harmonics): "
 
24
!
 
25
WRITE/KEYW INROOT1/C/1/60 {P1}
 
26
WRITE/KEYW OUTROOT/C/1/60 {P2}
 
27
WRITE/KEYW STARTTSA/D/1/1    {P3}
 
28
WRITE/KEYW STEPTSA/D/1/1     {P4}
 
29
WRITE/KEYW NSTEPS/I/1/1   {P5}
 
30
WRITE/KEYW ORDERTSA/I/1/1    {P6}
 
31
IF ORDERTSA .LT. 1 THEN
 
32
  WRITE/KEYW ORDERTSA/I/1/1 2
 
33
  WRITE/OUT *** Wrong ORDERTSA, set to 2
 
34
ELSEIF ORDERTSA .GT. 100 THEN
 
35
  WRITE/KEYW ORDERTSA/I/1/1 100
 
36
  WRITE/OUT *** Wrong ORDERTSA, set to 100
 
37
ENDIF
 
38
!
 
39
SELECT/TABLE {INROOT1}     SELECT.AND.:TIME.NE.NULL.AND.:VALUE.NE.NULL
 
40
OUTPUTI(2) = ORDERTSA
 
41
IF OUTPUTI(1) .LE. OUTPUTI(2) THEN
 
42
  WRITE/OUT *** Too few selected/correct data or too many bins
 
43
  RETURN/EXIT
 
44
ENDIF
 
45
COPY/TABLE  {INROOT1}      ZZMID1WORK
 
46
!
 
47
!   Call routine
 
48
!
 
49
WRITE/KEYW IN_A/C/1/60    ZZMID1WORK
 
50
WRITE/KEYW OUT_A/C/1/60   {OUTROOT}
 
51
WRITE/KEYW   HISTORY    "ORT/TSA"
 
52
RUN CON_EXE:TSAORT
 
53
DELETE/TABLE  ZZMID1WORK NO