OB.DAAC Logo
NASA Logo
Ocean Color Science Software

ocssw V2022
dirini.f
Go to the documentation of this file.
1  SUBROUTINE dirini
2  i (numstr, numclm, ifxclm,
3  o mrkclm, idfclm, nrfclm, mapclm, skyclm, mrkstr)
4 C-----------------------------------------------------------------------
5 C MODULE NAME: STDIRINI
6 C
7 C
8 C PURPOSE: TO INITIALIZE DIRECT MATCH ALGORITHM BY INITIALIZING
9 C SEVERAL FIELDS IN THE CLUMP AND OBSERVATION DATA.
10 C
11 C
12 C ARGUMENT LIST:
13 C ARGUMENT I/O TYPE DIMENS DESCRIPTION
14 C -------- --- ---- ------ -----------
15 C NUMSTR I I*4 NUMBER OF OBSERVATIONS
16 C NUMCLM I I*4 NUMBER OF CLUMPS
17 C IFXCLM I I*4 MAX # OF REF STAR MATCHES
18 C (FIRST DIM. FOR MAPCLM & SKYCLM ARRAYS)
19 C MRKCLM I O I*4 * CLUMP STATUS FLAGS
20 C IDFCLM I O I*4 * CLUMP IDENTIFICATION FLAGS
21 C NRFCLM I O I*4 * NUMBER OF REFERENCE STAR MATCHES PER CLUMP
22 C MAPCLM I O I*4 10,* LIST OF REFERENCE STAR MATCH CATALOG #'S
23 C SKYCLM I O R*4 10,3,* LIST OF CORRECTED REFERENCE STAR POSITIONS
24 C MRKSTR I O I*4 * OBSERVATION STATUS FLAGS
25 C
26 C
27 C COMMON BLOCK VARIABLES USED:
28 C COMMON VAR I/O VAR I/O VAR I/O VAR I/O
29 C ------ --- --- --- --- --- --- --- ---
30 C CMDEBG LEVDBG I LUDBUG I
31 C
32 C ++INCLUDE STCMDEBG
33  INTEGER*4 LEVDBG(8),LUDBUG
34  COMMON /cmdebg/levdbg,ludbug
35 C
36 C EXTERNAL FILES REFERENCED:
37 C FILENAME OPERATION FORTRAN UNIT ID
38 C -------- --------- ---------------
39 C NONE
40 C
41 C EXTERNAL REFERENCES:
42 C --------------------------------------------------------------------
43 C NONE
44 C
45 C SUBROUTINE CALLED FROM:
46 C --------------------------------------------------------------------
47 C STDIRECT - DIRECT MATCH ALGORITHM
48 C
49 C
50 C CONSTRAINTS, RESTRICTIONS, MESSAGES, NOTES:
51 C --------------------------------------------------------------------
52 C NONE
53 C
54 C REQUIREMENTS REFERENCES:
55 C --------------------------------------------------------------------
56 C NONE
57 C
58 C DEVELOPMENT HISTORY:
59 C DATE AUTHOR DESCRIPTION
60 C -------- ------ -----------
61 C 8/ 8/88 R.J. BURLEY DESIGN
62 C 5/11/89 R.J. BURLEY CODED
63 C-----------------------------------------------------------------------
64 C METHOD:
65 C DO FOR ALL OBSERVATIONS
66 C IF (OBS HAS BEEN DROPPED BECAUSE OF NON-IDENTIFICATION) THEN
67 C SET STATUS FLAG FOR THIS OBSERVATION TO 0
68 C ENDIF
69 C ENDDO FOR
70 C DO FOR ALL CLUMPS
71 C IF ((MRKCLM(CLUMP#) = 0).OR.(MRKCLM(CLUMP#) INDICATES DROPPED
72 C DUE TO NON-IDENTIFICATION)) THEN
73 C SET CLUMP STATUS FLAG TO 0
74 C SET CLUMP IDENTIFICATION FLAG TO 0
75 C SET # OF REFERENCE STAR MATCHES TO 0
76 C SET ALL ELEMENTS OF LIST OF MATCHED SKYMAP #'S TO 0
77 C SET ALL ELEMENTS OF LIST OF CORRECTED REF STAR POSITIONS TO 0
78 C ENDIF
79 C ENDDO FOR
80 C RETURN
81 C-----------------------------------------------------------------------
82 C
83 C * DEFINE PARAMETER VARIABLES
84  real*4 skyclm(10,3,*)
85 C
86  INTEGER*4 NUMSTR , NUMCLM , IFXCLM
87  INTEGER*4 MRKCLM(*), IDFCLM(*), NRFCLM(*), MAPCLM(10,*), MRKSTR(*)
88 C
89 C * DECLARE LOCAL VARIABLES
90  INTEGER*4 IOBS, ICLM, IMATCH
91 C
92 C INITIALIZE ROUTINE
93  IF (levdbg(7) .NE. 0) WRITE (ludbug,1000)
94 C
95 C
96 C UNDROP UNIDENTIFIED OBSERVATIONS
97  DO 100 iobs = 1,numstr
98  IF ((mrkstr(iobs) .EQ. 16).OR.(mrkstr(iobs) .EQ. 17)) THEN
99  mrkstr(iobs) = 0
100  ENDIF
101 100 CONTINUE
102 C
103 C UNDROP UNIDENTIFIED CLUMPS AND
104 C RESET IDENTIFICATION ARRAYS
105  DO 300 iclm = 1,numclm
106  IF ((mrkclm(iclm) .EQ. 0).OR.(mrkclm(iclm) .EQ. 14).OR.
107  1 (mrkclm(iclm) .EQ. 15)) THEN
108  mrkclm(iclm) = 0
109  idfclm(iclm) = 0
110  nrfclm(iclm) = 0
111  DO 200 imatch = 1,ifxclm
112  mapclm(imatch,iclm) = 0
113  skyclm(imatch,1,iclm) = 0.0
114  skyclm(imatch,2,iclm) = 0.0
115  skyclm(imatch,3,iclm) = 0.0
116 200 CONTINUE
117  ENDIF
118 300 CONTINUE
119 C
120 C NORMAL TERMINATION
121  IF (levdbg(7) .NE. 0) WRITE (ludbug,2000)
122  RETURN
123 C
124 C FORMAT SECTION
125 1000 FORMAT(' *** ENTER DIRINI ***')
126 2000 FORMAT(' *** EXIT DIRINI ***')
127  END
#define real
Definition: DbAlgOcean.cpp:26
subroutine dirini(NUMSTR, NUMCLM, IFXCLM, MRKCLM, IDFCLM, NRFCLM, MAPCLM, SKYCLM, MRKSTR)
Definition: dirini.f:4