← Back to branch summary

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

« back to all changes in this revision

Viewing changes to contrib/invent/libsrc/rmflst.for

  • 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/invent/libsrc/rmflst.for
 
1
C @(#)rmflst.for        19.1 (ES0-DMD) 02/25/03 13:25:38
 
2
C===========================================================================
 
3
C Copyright (C) 1995 European Southern Observatory (ESO)
 
4
C
 
5
C This program is free software; you can redistribute it and/or 
 
6
C modify it under the terms of the GNU General Public License as 
 
7
C published by the Free Software Foundation; either version 2 of 
 
8
C the License, or (at your option) any later version.
 
9
C
 
10
C This program is distributed in the hope that it will be useful,
 
11
C but WITHOUT ANY WARRANTY; without even the implied warranty of
 
12
C MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
13
C GNU General Public License for more details.
 
14
C
 
15
C You should have received a copy of the GNU General Public 
 
16
C License along with this program; if not, write to the Free 
 
17
C Software Foundation, Inc., 675 Massachusetss Ave, Cambridge, 
 
18
C MA 02139, USA.
 
19
C
 
20
C Corresponding concerning ESO-MIDAS should be addressed as follows:
 
21
C       Internet e-mail: midas@eso.org
 
22
C       Postal address: European Southern Observatory
 
23
C                       Data Management Division 
 
24
C                       Karl-Schwarzschild-Strasse 2
 
25
C                       D 85748 Garching bei Muenchen 
 
26
C                       GERMANY
 
27
C===========================================================================
 
28
C
 
29
C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
30
C
 
31
C
 
32
C
 
33
C----------------------------------------------------------------------
 
34
      SUBROUTINE RMFLST(M, L0, L1, NREG, LSTP, NCAT)
 
35
C
 
36
      IMPLICIT  NONE
 
37
      INCLUDE   'MID_REL_INCL:INVENT.INC/NOLIST'
 
38
C
 
39
      INTEGER   M
 
40
      INTEGER   L0
 
41
      INTEGER   L1
 
42
      INTEGER   NREG
 
43
      INTEGER   LSTP(0:4,0:NREG)
 
44
      INTEGER   NCAT(NIPAR,MAXCNT)
 
45
C
 
46
      INTEGER   LM, LL
 
47
      INTEGER   IREG , LOBJ
 
48
C     INTEGER   NCT1(NIPAR) , NCT2(NIPAR) 
 
49
      INTEGER   NOBJ , NP
 
50
C
 
51
C     REAL   PMT1(NRPAR) , PMT2(NRPAR)
 
52
C     REAL   PRC1(0:MAXSUB) , PRC2(0:MAXSUB)
 
53
C
 
54
      LOGICAL   BUFF1 , BUFF2
 
55
C
 
56
      NP = MAXCNT
 
57
C
 
58
C ******      Find first object in the region LOBJ.
 
59
C
 
60
C      IF ( M .GT. L0 .AND. M .LE. L1 ) THEN
 
61
          BUFF1 = .TRUE.
 
62
          LM = M - L0
 
63
          IREG = NCAT(7,LM)
 
64
C      ELSE
 
65
C          BUFF1 = .FALSE.
 
66
C          READ ( ISF , REC=M ) NCT1 , PMT1 , PRC1
 
67
C          IREG = NCT1(7)
 
68
C      ENDIF
 
69
      LOBJ = LSTP(0,IREG)
 
70
      IF ( LOBJ .EQ. M ) THEN
 
71
C
 
72
C ***     M is the first object on list.
 
73
C
 
74
C          IF ( BUFF1 ) THEN
 
75
              NOBJ = NCAT(8,LM)
 
76
C          ELSE
 
77
C              NOBJ = NCT1(8)
 
78
C          ENDIF
 
79
          LSTP(0,IREG) = NOBJ
 
80
C          IF ( BUFF1 ) THEN
 
81
              NCAT(7,LM) = 0
 
82
              NCAT(8,LM) = 0
 
83
C          ELSE
 
84
C              NCT1(7) = 0
 
85
C              NCT1(8) = 0
 
86
C              WRITE ( ISF , REC=M ) NCT1 , PMT1 , PRC1
 
87
C          ENDIF
 
88
      ELSE IF ( LOBJ .EQ. 0 ) THEN
 
89
C          IF ( BUFF1 ) THEN
 
90
              NCAT(7,LM) = 0
 
91
              NCAT(8,LM) = 0
 
92
C          ELSE
 
93
C              NCT1(7) = 0
 
94
C              NCT1(8) = 0
 
95
C              WRITE ( ISF , REC=M ) NCT1 , PMT1 , PRC1
 
96
C          ENDIF
 
97
      ELSE
 
98
C
 
99
C *** M is in the middle or on the end of list.
 
100
C
 
101
   10     CONTINUE
 
102
C              IF ( LOBJ .GT. L0 .AND. LOBJ .LE. L1 ) THEN
 
103
                  LL = LOBJ - L0
 
104
                  NOBJ = NCAT(8,LL)
 
105
                  BUFF2 = .TRUE.
 
106
C              ELSE
 
107
C                  READ ( ISF , REC=LOBJ ) NCT2 , PMT2 , PRC2
 
108
C                  NOBJ = NCT2(8)
 
109
C                  BUFF2 = .FALSE.
 
110
C              ENDIF
 
111
              IF ( NOBJ .EQ. M ) THEN
 
112
C                  IF ( BUFF1 ) THEN
 
113
                      NOBJ = NCAT(8,LM)
 
114
                      NCAT(7,LM) = 0
 
115
                      NCAT(8,LM) = 0 
 
116
C                  ELSE
 
117
C                      NOBJ = NCT1(8)
 
118
C                      NCT1(7) = 0
 
119
C                      NCT1(8) = 0
 
120
C                      WRITE ( ISF , REC=M ) NCT1 , PMT1 , PRC1
 
121
C                  ENDIF
 
122
C                  IF ( BUFF2 ) THEN
 
123
                      NCAT(8,LL) = NOBJ
 
124
C                  ELSE
 
125
C                      NCT2(8) = NOBJ
 
126
C                      WRITE ( ISF , REC=LOBJ ) NCT2 , PMT2 , PRC2
 
127
C                  ENDIF
 
128
              ELSE IF ( NOBJ .EQ. 0 ) THEN
 
129
C                  IF ( BUFF1 ) THEN
 
130
                      NCAT(7,LM) = 0
 
131
                      NCAT(8,LM) = 0
 
132
C                  ELSE
 
133
C                      NCT1(7) = 0
 
134
C                      NCT1(8) = 0
 
135
C                      WRITE ( ISF , REC=M ) NCT1 , PMT1 , PRC1
 
136
C                  ENDIF
 
137
              ELSE
 
138
                  LOBJ = NOBJ
 
139
                  GOTO 10
 
140
              ENDIF
 
141
          CONTINUE
 
142
      ENDIF
 
143
C
 
144
      IF ( LSTP(3,0) .EQ. M ) THEN
 
145
          LSTP(3,0) = M - 1
 
146
      ENDIF
 
147
C
 
148
      RETURN
 
149
C
 
150
      END
 
151
C