← Back to branch summary

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

« back to all changes in this revision

Viewing changes to contrib/surfphot/proc/bgfit.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/surfphot/proc/bgfit.prg
 
1
! @(#)bgfit.prg 19.1 (ES0-DMD) 02/25/03 13:31:07
 
2
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
3
!.IDENT:    bgfit.prg
 
4
!.PURPOSE:  implement FIT/BACKGROUND
 
5
!.AUTHORS:  E. Held
 
6
!.USE:      FIT/BACK result_frame = input_spec deg,max_iter clip_1,clip_next
 
7
!                    up_lim surface
 
8
!           FIT/BACK input_spec deg,max_iter clip_1,clip_next up_lim surface
 
9
!           FIT/BACK result_frame = in_frame COEF surface
 
10
!           FIT/BACK in_frame COEF surface
 
11
!
 
12
!           with input_spec either in_frame,in_mask (a) or in_frame (b)
 
13
!.VERSION:  840518 rev. 1.2
 
14
!.VERSION:  840614 KB ???
 
15
!.VERSION:  901108 KB ???
 
16
!----------------------------------------------------------------------------
 
17
DEFINE/LOCAL L/I/1/5 0,0,0,0,0
 
18
!
 
19
IF P2(1:1) .EQ. "=" THEN                !handle option result = ...
 
20
   DEFINE/PARAM P1 ? IMA "Enter result frame:"
 
21
   DEFINE/PARAM P4 3,10 N
 
22
   DEFINE/PARAM P5 1.,2.5 N
 
23
   DEFINE/PARAM P6 0.05 N
 
24
   DEFINE/PARAM P7 + IMA
 
25
!
 
26
   WRITE/KEYW OUT_A 'P1'
 
27
   L = M$INDEX(P3,",")
 
28
   IF L .GT. 0 THEN
 
29
      L(2) = L(1) - 1
 
30
      WRITE/KEYW IN_A 'P3(1:{L(2)})'
 
31
      L(2) = L(1) + 1
 
32
      WRITE/KEYW IN_B 'P3({L(2)}:>)'
 
33
   ELSE
 
34
      WRITE/KEYW IN_A 'P3'
 
35
      WRITE/KEYW IN_B +
 
36
   ENDIF
 
37
!
 
38
   IF P4(1:4) .EQ. "COEF" THEN
 
39
      WRITE/KEYW INPUTC BG_COEF         !coefficients are already calculated
 
40
      WRITE/KEYW OUT_B 'P5'
 
41
   ELSE
 
42
      WRITE/KEYW INPUTC ?
 
43
      WRITE/KEYW INPUTI/I/1/2 'P4'
 
44
      WRITE/KEYW INPUTR/R/1/2 'P5'
 
45
      WRITE/KEYW INPUTR/R/3/1 'P6'
 
46
      WRITE/KEYW OUT_B 'P7'
 
47
      WRITE/KEYW BG_COEF/D/1/28 0.0 ALL
 
48
   ENDIF
 
49
!
 
50
ELSE                                    !handle option input_spec ...
 
51
   WRITE/KEYW OUT_A +
 
52
   L = M$INDEX(P1,",")
 
53
   IF L .GT. 0 THEN
 
54
      L(2) = L(1) - 1
 
55
      WRITE/KEYW IN_A 'P1(1:{L(2)})'
 
56
      L(2) = L(1) + 1
 
57
      WRITE/KEYW IN_B 'P1({L(2)}:>)'
 
58
   ELSE
 
59
      WRITE/KEYW IN_A 'P1'
 
60
      WRITE/KEYW IN_B +
 
61
   ENDIF
 
62
   IF P2(1:4) .EQ. "COEF" THEN
 
63
      DEFINE/PARAM P3 + IMA
 
64
      WRITE/KEYW INPUTC BG_COEF         !coefficients are already calculated
 
65
      WRITE/KEYW OUT_B 'P3'
 
66
   ELSE
 
67
      DEFINE/PARAM P2 3,10
 
68
      DEFINE/PARAM P3 1.,2.5
 
69
      DEFINE/PARAM P4 .05
 
70
      DEFINE/PARAM P5 + IMA
 
71
!
 
72
      WRITE/KEYW INPUTC ?
 
73
      WRITE/KEYW INPUTI/I/1/2 'P2'
 
74
      WRITE/KEYW INPUTR/R/1/2 'P3'
 
75
      WRITE/KEYW INPUTR/R/3/1 'P4'
 
76
      WRITE/KEYW OUT_B 'P5'
 
77
      WRITE/KEYW BG_COEF/D/1/28 0.0 ALL
 
78
   ENDIF
 
79
!
 
80
ENDIF
 
81
!
 
82
RUN CON_EXE:BACKGRND                     !now execute the program