Due to the lapse in federal government funding, NASA is not updating this website. We sincerely regret this inconvenience.
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