← Back to branch summary

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

« back to all changes in this revision

Viewing changes to stdred/echelle/proc/necresp.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:
stdred/echelle/proc/necresp.prg
 
1
! @(#)necresp.prg       19.1 (ES0-DMD) 02/25/03 14:19:42
 
2
! @(#)necresp.prg       19.1  (ESO)  02/25/03  14:19:42
 
3
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
4
!
 
5
!.COPYRIGHT   (C) 1991 European Southern Observatory
 
6
!.IDENT       echresp.prg
 
7
!.AUTHOR      Pascal Ballester,  ESO - Garching
 
8
!.KEYWORDS    Spectroscopy, Echelle, 
 
9
!.PURPOSE     Command RESPONSE/ECHELLE [std] [fluxtab] [response]
 
10
!.VERSION     1.0    Creation    
 
11
!.VERSION     1.1    change `\' to ` ',  KB
 
12
!
 
13
!-------------------------------------------------------
 
14
!
 
15
DEFINE/PARAM  P1  {STD}      I  "Raw standard star:"
 
16
DEFINE/PARAM  P2  {FLUXTAB}  I  "Flux table:"
 
17
DEFINE/PARAM  P3  {RESPONSE} I  "Output response:"
 
18
 
 
19
VERIFY/ECHELLE {P1}
 
20
SET/ECH STD={P1}  FLUXTAB={P2}   RESPONSE={P3}
 
21
VERIFY/SPEC    {FLUXTAB}  MID_STANDARD FLUXTAB
 
22
 
 
23
 
 
24
DEFINE/LOCAL    LOOP/I/1/1    0
 
25
DEFINE/LOCAL    PIXPOS/I/1/1  0   
 
26
DEFINE/LOCAL    NROW/I/1/1    0
 
27
 
 
28
REDUCE/SIMPLY   {STD}  &m  &j  !  Simplified reduction
 
29
FILTER/MEDIAN   &m   &p  {FILTMED(1)},{FILTMED(2)},{FILTMED(3)}
 
30
FILTER/SMOOTH   &p   &n  {FILTSMO(1)},{FILTSMO(2)},{FILTSMO(3)}
 
31
 
 
32
WRITE/OUT       "Compute flux calibration: Read columns :WAVE, :FLUX_W -
 
33
from {FLUXTAB}"
 
34
CONVERT/TABLE   &o = {FLUXTAB} :WAVE :FLUX_W  &n SPECTR
 
35
VERIFY/ECHELLE  &o   REBIN
 
36
VERIFY/ECHELLE  &o AGAINST &n
 
37
COMPUTE/IMAGE   &p = &o/&n
 
38
 
 
39
WRITE/IMAGE     &p  [<,<:@{PIXNUL(1)},>] 0.  ALL  {SESSOUTV}
 
40
 
 
41
NROW = {middummp,NPIX(2)}
 
42
DO LOOP = 1  {NROW}
 
43
   PIXPOS = {middummm.bdf,NPTOT({LOOP})} - PIXNUL(2)
 
44
   WRITE/IMAGE &p [@{PIXPOS},@{LOOP}:>,@{LOOP}]  0.  ALL  {SESSOUTV}
 
45
ENDDO
 
46
 
 
47
COPY/II          middummp  {RESPONSE}
 
48
COPY/DD          {STD}  *,3  {RESPONSE}
 
49
 
 
50
!@ creifnot 2 heat
 
51
!DISP/CHAN       0
 
52
!LOAD/IMAGE      {RESPONSE}  0 8,8
 
53
 
 
54
WRITE/OUT       OUTPUT RESPONSE
 
55
WRITE/OUT       "---------------"
 
56
READ/DESC       {RESPONSE}
 
57
WRITE/OUT       "------------------------------------------------------"
 
58
WRITE/OUT       "Check if RESPONSE looks reasonable..."
 
59
WRITE/OUT       "To correct it use REPEAT/ECHELLE"
 
60
WRITE/OUT       "------------------------------------------------------"
 
61
WRITE/OUT       "Apply RESPONSE to Standard Star"
 
62
 
 
63
VERIFY/ECHELLE  &m AGAINST {RESPONSE}
 
64
COMPUTE/IMAGE   &a = &m*{RESPONSE}
 
65
 
 
66
COPY/DD         {STD}  *,3  &a
 
67
 
 
68
WRITE/OUT       "Use MERGE/ECHELLE &a name 3.0 AVERAGE"
 
69
WRITE/OUT       "to inspect the calibrated standard star"
 
70
 
 
71
RETURN
 
72
 
 
73
 
 
74
 
 
75
 
 
76