-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathdmcond.f
88 lines (88 loc) · 2.2 KB
/
dmcond.f
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
SUBROUTINE DM_COND ( iflno, isrch, irow, icol, sflag, iret )
C************************************************************************
C* DM_COND *
C* *
C* This subroutine checks whether the search conditions for a single *
C* search are met. *
C* *
C* DM_COND ( IFLNO, ISRCH, IROW, ICOL, SFLAG, IRET ) *
C* *
C* Input parameters: *
C* IFLNO INTEGER File number *
C* ISRCH INTEGER Search number *
C* IROW INTEGER Row to check *
C* ICOL INTEGER Column to check *
C* *
C* Output parameters: *
C* SFLAG LOGICAL Result of search *
C* IRET INTEGER Return code *
C* 0 = normal return *
C** *
C* Log: *
C* M. desJardins/GSFC 4/87 *
C************************************************************************
INCLUDE 'GEMPRM.PRM'
INCLUDE 'dmcmn.cmn'
C
LOGICAL sflag
C-----------------------------------------------------------------------
iret = 0
C
C* Compute which header corresponds to this column.
C
jcol = krow ( iflno ) + icol
C
C* Check that the row and column headers are defined.
C
sflag = .true.
IF ( kheadr ( 0, irow, iflno ) .eq. IMISSD ) THEN
sflag = .false.
RETURN
END IF
C
IF ( ( icol .ne. 0 ) .and. ( kheadr ( 0, jcol, iflno ) .eq.
+ IMISSD ) ) THEN
sflag = .false.
RETURN
END IF
C
C* Check that this is a valid search.
C
IF ( isrch .eq. 0 ) THEN
IF ( .not. srcflg ( iflno ) ) RETURN
ELSE IF ( isrch .gt. nsrch ( iflno ) ) THEN
RETURN
END IF
C
C* Check row conditions.
C
DO i = 1, ksnrow ( isrch, iflno )
loc = kslrow ( i, isrch, iflno )
IF (( kheadr (loc,irow,iflno) .lt. ksrlov (i,isrch,iflno))
+ .or.
+ ( kheadr (loc,irow,iflno) .gt. ksrhiv (i,isrch,iflno)))
+ THEN
sflag = .false.
RETURN
END IF
END DO
C
C* Only check column conditions if column is not 0.
C
IF ( icol .eq. 0 ) RETURN
C
C* Check column conditions.
C
DO i = 1, ksncol ( isrch, iflno )
loc = kslcol ( i, isrch, iflno )
IF (( kheadr (loc,jcol,iflno) .lt. ksclov (i,isrch,iflno))
+ .or.
+ ( kheadr (loc,jcol,iflno) .gt. kschiv (i,isrch,iflno)))
+ THEN
sflag = .false.
RETURN
END IF
END DO
C*
RETURN
END