← Back to branch summary

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

« back to all changes in this revision

Viewing changes to stdred/feros/locproc/temperature.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/feros/locproc/temperature.prg
 
1
! @(#)temperature.prg   19.1 (ESO-IPG) 02/25/03 14:21:47
 
2
! @@ temperature.prg      
 
3
! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
4
!
 
5
! MIDAS procedure  temperature.prg  
 
6
! A. Kaufer     981124  
 
7
!
 
8
! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
9
!
 
10
define/par p1 ? C "Enter name for file list   : "
 
11
define/par p2 ? N "Enter filenumber start,end : "
 
12
!
 
13
define/local range/i/1/2 {p2}
 
14
define/local infile/c/1/40 " " all
 
15
define/local i/i/1/1 0 
 
16
define/local ind/i/1/1 0 
 
17
define/local num/i/1/1 0 
 
18
define/local ftype/i/1/1 0 
 
19
 
 
20
crea/tab {p1} 3 1000
 
21
crea/col {p1} :num   R*4 F5.2
 
22
crea/col {p1} :JD24  R*8 F13.5
 
23
crea/col {p1} :tspec R*4 F5.2
 
24
crea/col {p1} :troom R*4 F5.2
 
25
crea/col {p1} :rhum  R*4 F7.2
 
26
 
 
27
do num = {range(1)} {range(2)} 1
 
28
!
 
29
 i = i+1
 
30
 ftype = 0
 
31
!
 
32
 ind = M$EXIST("fero{num}.mt")
 
33
 if  ind .eq. 1 then
 
34
  infile = "fero{num}.mt"
 
35
  crea/ima &header 1,1 1,1 NODATA
 
36
  cop/dd {infile} *,3 &header
 
37
  comp/st &header
 
38
  infile = "middummheader.bdf"
 
39
  ftype = 3
 
40
 endif
 
41
! ind = M$EXIST("f{num}ext1.bdf")
 
42
! if  ind .eq. 1 then
 
43
!  infile = "f{num}ext1.bdf"
 
44
!  ftype = 2
 
45
! endif
 
46
 ind = M$EXIST("f{num}1.bdf")
 
47
 if ind .eq. 1 then 
 
48
  infile = "f{num}1.bdf"
 
49
  ftype = 1
 
50
 endif
 
51
!
 
52
 if ftype .eq. 0 goto next
 
53
!
 
54
 {p1},:num,{i}   = {num}
 
55
 copy/dk {infile} O_TIME/D/4/1 jd24/d/1/1
 
56
 {p1},:jd24,{i}  = m$value(jd24)
 
57
 {p1},:tspec,{i} = {{infile},TSPEC}
 
58
 {p1},:troom,{i} = {{infile},TROOM}
 
59
 {p1},:rhum,{i}  = {{infile},RHROOM}
 
60
!
 
61
 read/tab {p1} @{i}
 
62
 next:
 
63
enddo
 
64
!
 
65
set/gra stype=2 xa ya=15,16.5 color=1
 
66
sele/tab {p1} :tspec.gt.0
 
67
plo/tab {p1}  :jd24  :tspec 
 
68
set/gra color=4
 
69
over/tab {p1} :jd24  :troom 3
 
70
set/gra ya color=1
 
71
 
 
72