Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 1 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFAGB Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1 C 2 C 3 C 4 SUBROUTINE WCFAGB 5 # (PKQ,PKLOG,WRCPP,SYSPP,NPK,IRC) 6 C 7 C + + + PURPOSE + + + 8 C WRC FLOOD FREQ ANALYSIS BULL 17 GUIDELINES EXECUTIVE RTNE. 9 C 10 C FLOOD FREQUENCY CALCULATIONS AS DESCRIBED IN U.S. WATER RESOURCES 11 C COUNCIL HYDROLOGY COMMITTEE BULLETIN 17 (1976), GUIDELINES FOR DETER- 12 C MINING FLOOD FLOW FREQUENCY, AND REVISED AS BULL 17-A (6/77). 13 C AND BULL. 17-B, 1981. 14 C 15 C ALTHOUGH THESE CALCULATIONS ARE BELIEVED TO CONFORM TO THE WRC GUIDE- 16 C LINES, THE USER IS RESPONSIBLE FOR ALL INTERPRETATIONS AND APPLI- 17 C CATIONS OF THE RESULTS. 18 C 19 C SEE FORMAT NO 1 FOR CURRENT VERSION NUMBER AND DATE. 20 C WCFAGB VERSION 2.0 5/1/77 BY W.KIRBY, U.S.GEOLOGICAL SURVEY, 21 C NATIONAL CTR STOP 430, RESTON, VA 22092. (703) 860-6947. 22 C VERSION 2.0 INCORPORATES REVISIONS TO WRC BULL 17 THROUGH 4/4/77. 23 C AND PUBLISHED IN BULL 17-A (6/77). 24 C REV 5/79 VER 2.2 - MINOR CHANGES IN -DLO,-FCA,--CSA,-CM0 FOR HAND- 25 C LING NONSTANDARD CONDITIONS. NO CHANGE IN STANDARD 26 C WRC CALCULATIONS. 27 C REV 7/20/79 VER 2.2 - INCREASED PRECISION IN COMPUTING SUMS AND 28 C SUMH IN -ASP, -DHH, -DLO. 29 C 30 C REV 10/18/79 VER 2.3 - FIX ROUND-OFF IN --CWS. USE STUTX FOR EXPEC- 31 C TED PROB IN --FCX. 32 C 33 C REV 1/81 VER 2.4 BULL 17-B. -- NEW WEIGHTED SKEW USING STD ERRORS 34 C OF GENERALIZED AND STATION SKEWS. (WCFCWS, FCA) 35 C -- THREE-POINT FIT PEARSON TYPE III M, S, G AFTER COND PROB ADJ 36 C -- COMPUTE WEIGHTED SKEW AFTER COND PROB ADJ, NOT BEFORE. 37 C -- AUTOMATIC HIGH OUTLIER TEST AND REVISED LOW OUTLIER TEST, 38 C USING GRUBBS-BECK NORMAL OUTLIER TEST. (WCFDHH, DLO, OUTKGB.) 39 C 40 C REV 4/28/81 WK - WCFDHH,DLO - USE NUMBER OF PEAKS ABOVE FLOOD BASE 41 C IN OUTLIER TESTS. (NOT NSYS/HISTPD) 42 C WCFDHH - PRINT HI-OUT TEST CRIT IF SUPERSEDED BY MIN HIST PK. 43 C 44 C REV 11/5/81 WK - WCFFCA CHECK *ABS-VALUE* OF WRCSKW IN TABLE RANGE 45 C 46 C REV 12/83 - K.FLYNN - REVISIONS FOR PRIME - WCFCM0 47 C 48 C REV 9/88 - AML - for conforming to OSW coding conventions 49 C 50 C USER INTERFACE TO WRC CALCULATIONS -- 51 C ASSIGN APPROPRIATE VALUES TO GENSKU...HISTPD. IN COMMON WCFCM1. 52 C PLACE FLOWS IN PKQ VECTOR, HISTORIC EVENTS FIRST.. Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 2 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 53 C CALL WCFAGB (PKQ ... IRC). (WCFAGB CALLS ALL OTHER ROUTINES USED.) 54 C UPON RETURN, IF IRC .LT. 3, THE RESULTS WILL BE IN COMMONS WCFCM1 55 C AND WCFCM2 AND IN VECTORS IN ARG LIST. OTHERWISE, IF IRC .EQ. 3, 56 C THE CALCULATIONS COULD NOT BE COMPLETED FOR REASONS EXPLAINED IN 57 C NUMBERED MESSAGES. THE CONTENTS OF THE COMMONS AND VECTORS OTHER 58 C THAN PKQ WILL BE MEANINGLESS. THE PKQ VECTOR IS NOT ALTERED BY 59 C WCFAGB AND THE PKLOG ELEMENTS ARE IN THE SAME ORDER AS THE PKQS. 60 C A SAMPLE MAIN PROGRAM, WCFZSD, IS INCLUDED IN THIS DISTRIBUTION. 61 C 62 C + + + DUMMY ARGUMENTS + + + 63 INTEGER NPK, IRC 64 REAL PKQ(NPK), PKLOG(NPK), WRCPP(NPK), SYSPP(NPK) 65 C 66 C + + + ARGUMENT DEFINITION + + + 67 C PKQ - peak flow input vector--first historic, then systematic 68 C PKLOG - peak flow LOG10 work-output vector 69 C WRCPP - WRC prob plot positions systematic and historic PKS. 70 C SYSPP - prob plot positions systematic record peaks 71 C NPK - number of observed peaks--first historic, then systematic 72 C IRC - return code--0, 1, 2, or 3: 73 C 0 - no error 74 C 3 - error, calculation aborted 75 C 76 C + + + PARAMETERS + + + 77 INCLUDE 'pmxint.inc' 78 C 79 C + + + COMMON BLOCKS + + + 80 INCLUDE 'cwcf0.inc' 81 INCLUDE 'cwcf1.inc' 82 INCLUDE 'cwcf2.inc' 83 C 84 C + + + LOCAL VARIABLES + + + 85 INTEGER NSYS1, ISYS 86 C 87 C + + + EXTERNALS + + + 88 EXTERNAL WCFASP, WCFDHH, WCFCSA, WCFDLO, WCFAPI, WCFFCA, 89 # WCFFCX, WCFEPP 90 C 91 C + + + END SPECIFICATIONS + + + 92 C 93 IRC=0 94 IF(MSL.GE.3)WRITE(MSG,001) 95 1 FORMAT(4X,'WCF001J-FLOOD FREQUENCY, BULLETIN 17-B.', 96 $ ' VER 2.6P (12/19/83)', 97 $ /,5X,' -PRELIMINARY MACHINE COMPUTATIONS. USER IS RE-', 98 $ /,5X,' -SPONSIBLE FOR ASSESSMENT AND INTERPRETATION. ') 99 C 100 CALL WCFAPI (PKQ,PKLOG,WRCPP,SYSPP,NPK,IRC) 101 IF(IRC.GE.3)GOTO95 102 ISYS=NHIST+1 103 NSYS1=NPK-NHIST 104 C Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 3 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 105 CALL WCFASP (PKLOG(ISYS),SYSPP(ISYS),NSYS1,IRC) 106 IF(IRC.GE.3)GOTO95 107 C 108 IF(WRCASK.GT.EPS2)GOTO50 109 IF(WRCASK.GE.EPS1)GOTO30 110 C 111 CALL WCFDLO (PKLOG(ISYS),NSYS1,IRC) 112 Ckmf IF(NLWOUT.GT.0 .AND. IRC.LT.3)CALL WCFCSA (4H17B1, IRC) 113 IF(NLWOUT.GT.0 .AND. IRC.LT.3)CALL WCFCSA ('17B1', IRC) 114 IF(IRC.GE.3)GOTO95 115 CALL WCFDHH (PKLOG,NPK, IRC) 116 Ckmf IF(NHISTN+NHIOUT.GT.0 .AND. IRC.LT.3) CALL WCFCSA (4H17B2,IRC) 117 IF(NHISTN+NHIOUT.GT.0 .AND. IRC.LT.3) CALL WCFCSA ('17B2',IRC) 118 IF(IRC.GE.3)GOTO95 119 GOTO70 120 C 121 30 CALL WCFDLO (PKLOG(ISYS),NSYS1,IRC) 122 IF(IRC.GE.3)GOTO95 123 CALL WCFDHH (PKLOG,NPK,IRC) 124 IF(IRC.GE.3)GOTO95 125 Ckmf IF(NHISTN+NHIOUT+NLWOUT.GT.0) CALL WCFCSA (4H17B3,IRC) 126 IF(NHISTN+NHIOUT+NLWOUT.GT.0) CALL WCFCSA ('17B3',IRC) 127 IF(IRC.GE.3)GOTO95 128 GOTO70 129 C 130 50 CALL WCFDHH (PKLOG,NPK,IRC) 131 Ckmf IF(NHIOUT+NHISTN.GT.0 .AND. IRC.LT.3) CALL WCFCSA (4H17B4,IRC) 132 IF(NHIOUT+NHISTN.GT.0 .AND. IRC.LT.3) CALL WCFCSA ('17B4',IRC) 133 IF(IRC.GE.3)GOTO95 134 CALL WCFDLO (PKLOG(ISYS),NSYS1,IRC) 135 Ckmf IF(NLWOUT.GT.0 .AND. IRC.LT.3) CALL WCFCSA (4H17B5,IRC) 136 IF(NLWOUT.GT.0 .AND. IRC.LT.3) CALL WCFCSA ('17B5',IRC) 137 IF(IRC.GE.3)GOTO95 138 C 139 70 CONTINUE 140 C 141 C IF(NOPPOS.NE.1) CALL WCFEPP (WRCPP, 4H17B ,NSYS+NHIST) 142 C CALL WCFFCA (WRCFC,4H17B ,IRC) 143 IF(NOPPOS.NE.1) CALL WCFEPP (WRCPP, '17B ', NSYS+NHIST) 144 CALL WCFFCA (WRCFC, '17B ', IRC ) 145 IF(IRC.GE.3)GOTO95 146 C 147 CALL WCFFCX (IRC) 148 IF(IRC.GE.3)GOTO95 149 C 150 IF(MSL.GE.3 .OR. (IRC.GE.2.AND.MSL.GE.2) ) WRITE(MSG,2)IRC 151 002 FORMAT(43H WCF002J-CALCS COMPLETED. RETURN CODE = ,I3) 152 RETURN 153 C 154 95 IF(MSL.GT.0)WRITE(MSG,003)IRC 155 003 FORMAT(43H ***WCF003E-CALCS ABORTED. RETURN CODE = ,I3) 156 RETURN Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 4 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 157 END Bytes of stack required for this program unit: 24. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- BIGLOG (REAL) scalar 316 (I1)4s (I1)12s CLIML (REAL) array 256 (I2)2s (I2)4s CLIMU (REAL) array 384 (I2)3s (I2)5s CLSIZE (REAL) scalar 40 (I1)3s (I1)11s CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s EPFC (REAL) array 128 (I2)2s (I2)4s EPS1 (REAL) scalar 320 (I1)5s (I1)13s 109u EPS2 (REAL) scalar 324 (I1)5s (I1)13s 108u GAGEB (REAL) scalar 12 (I3)2s (I3)10s GENSKU (REAL) scalar 0 (I3)2s (I3)10s GSKWGT (REAL) scalar 108 (I3)6s (I3)10s GSMAX (REAL) scalar 32 (I1)3s (I1)11s GSMIN (REAL) scalar 28 (I1)3s (I1)11s HISTPD (REAL) scalar 28 (I3)3s (I3)10s HISTPN (REAL) scalar 100 (I3)6s (I3)10s HISTWT (REAL) scalar 104 (I3)6s (I3)10s HRECWO (REAL) scalar 328 (I1)5s (I1)13s HRECWS (REAL) scalar 336 (I1)5s (I1)13s HRECWX (REAL) scalar 332 (I1)5s (I1)13s IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s INDXPT (INTEGER) array 56 (I1)4s (I1)9s IRC (INTEGER) scalar 5s 63s 93= 100r 100r 100r 101u 105r 105r 105r 106u 111r 111r 111r 113u 113r 113r 113r 114u 115r 115r 115r 117u 117r 117r 117r 118u 121r 121r 121r 122u 123r 123r 123r 124u 126r 126r 126r 127u 130r 130r 130r 132u 132r 132r 132r 133u 134r 134r 134r 136u 136r 136r 136r 137u 144r 144r 144r 145u 147r 147r 147r 148u 150u 150o 154o ISYS (INTEGER) scalar 85s 102= 105u 105u 111u 121u 134u IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s MSG (INTEGER) scalar 0 (I1)2s (I1)8s 94o 150o 154o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 94u 150u 150u 154u MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u (I2)4u (I2)4u (I2)4u (I2)5u (I2)5u NBGB (INTEGER) scalar 84 (I3)5s (I3)9s NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s 117u 126u 132u NHIST (INTEGER) scalar 24 (I3)2s (I3)9s 102u 103u 143u NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s 117u 126u 132u NINDX (INTEGER) scalar 52 (I1)4s (I1)9s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 5 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s 113u 126u 136u NMISS (INTEGER) scalar 76 (I3)5s (I3)9s NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s 143u NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s NPK (INTEGER) scalar 5s 63s 64u 64u 64u 64u 100r 100r 100r 103u 115r 115r 115r 123r 123r 123r 130r 130r 130r NSYS (INTEGER) scalar 80 (I3)5s (I3)9s 143u NSYS1 (INTEGER) scalar 85s 103= 105r 105r 105r 111r 111r 111r 121r 121r 121r 134r 134r 134r PKLOG (REAL) array 5s 64s 100r 100r 100r 105r 105r 105r 111r 111r 111r 115r 115r 115r 121r 121r 121r 123r 123r 123r 130r 130r 130r 134r 134r 134r PKQ (REAL) array 5s 64s 100r 100r 100r QHIOUT (REAL) scalar 20 (I3)2s (I3)11s QLWOUT (REAL) scalar 16 (I3)2s (I3)11s RMSDGS (REAL) scalar 340 (I1)6s (I1)14s RMSEGS (REAL) scalar 8 (I3)2s (I3)10s SIGHOT (REAL) scalar 360 (I1)7s (I1)15s SIGLOT (REAL) scalar 356 (I1)7s (I1)15s SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s SYSMAX (REAL) scalar 148 (I3)8s (I3)11s SYSMIN (REAL) scalar 144 (I3)8s (I3)12s SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSPP (REAL) array 5s 64s 100r 100r 100r 105r 105r 105r SYSRFC (REAL) array 512 (I2)3s (I2)5s SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSUAV (REAL) scalar 120 (I3)6s (I3)12s SYSUSD (REAL) scalar 124 (I3)7s (I3)12s TXPROB (REAL) array 184 (I1)4s (I1)12s WCFAGB SUBROUTINE 4s WCFAPI SUBROUTINE 88u 100u WCFASP SUBROUTINE 88u 105u WCFCSA SUBROUTINE 88u 113u 117u 126u 132u 136u WCFDHH SUBROUTINE 88u 115u 123u 130u WCFDLO SUBROUTINE 88u 111u 121u 134u WCFEPP SUBROUTINE 89u 143u WCFFCA SUBROUTINE 88u 144u WCFFCX SUBROUTINE 89u 147u WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s WORK (REAL) array 200 (I3)8s (I3)12s WRCAAV (REAL) scalar 52 (I3)4s (I3)13s WRCASD (REAL) scalar 56 (I3)4s (I3)13s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 6 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- WRCASK (REAL) scalar 60 (I3)4s (I3)13s 108u 109u WRCBAS (REAL) scalar 32 (I3)3s (I3)13s WRCFC (REAL) array 0 (I2)2s (I2)4s 144r 144r 144r WRCHHB (REAL) scalar 64 (I3)4s (I3)13s WRCHOT (REAL) scalar 68 (I3)4s (I3)13s WRCLOW (REAL) scalar 72 (I3)5s (I3)13s WRCPAB (REAL) scalar 36 (I3)3s (I3)14s WRCPP (REAL) array 5s 64s 100r 100r 100r 143r 143r 143r WRCSKW (REAL) scalar 48 (I3)4s (I3)14s WRCUAV (REAL) scalar 40 (I3)3s (I3)14s WRCUSD (REAL) scalar 44 (I3)4s (I3)14s WSKLAT (REAL) scalar 312 (I1)4s (I1)12s -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 30 109g 121s 50 108g 130s 70 119g 128g 139s 95 101g 106g 114g 118g 122g 124g 127g 133g 137g 145g 148g 154s 1 94f 95s 2 150f 151s 3 154f 155s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 7 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFAPI Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 158 C 159 C 160 C 161 SUBROUTINE WCFAPI 162 # (PKQ,PKLOG,WRCPP,SYSPP,NPK,IER) 163 C 164 C + + + PURPOSE + + + 165 C ANALYZE PARAMETERS AND INITIALIZE 166 C (WRC Bulletin-17 Flood Frequency Analysis) 167 C 168 C + + + DUMMY ARGUMENTS + + + 169 INTEGER NPK, IER 170 REAL PKQ(NPK), PKLOG(NPK), WRCPP(NPK), SYSPP(NPK) 171 C 172 C + + + ARGUMENT DEFINITIONS + + + 173 C PKQ - peak flow input vector--historic first, then systematic 174 C PKLOG - peak flow LOG10 work-output vector 175 C WRCPP - WRC prob plot positions observed PKS 176 C SYSPP - prob plot positions systematic record PKS 177 C NPK - number of observed peaks 178 C IER - error return code 0, 1, 2, or 3: 179 C 0 - no error 180 C 3 - error, calculation aborted 181 C 182 C + + + PARAMETERS + + + 183 INCLUDE 'pmxint.inc' 184 C 185 C + + + COMMON BLOCKS + + + 186 INCLUDE 'cwcf0.inc' 187 INCLUDE 'cwcf1.inc' 188 INCLUDE 'cwcf2.inc' 189 Cprh replace FC with explicit variable names in cwcf2.inc 190 Cprh REAL FC(31,5) 191 Cprh COMMON/ WCFCM2 /FC 192 C 193 C + + + LOCAL VARIABLES + + + 194 INTEGER J, I 195 REAL BIGQ, T, FINIT(30) 196 C 197 C + + + INTRINSICS + + + 198 INTRINSIC SIGN, ALOG10, MAX0, FLOAT, INT 199 C 200 C + + + EQUIVALENCE + + + 201 EQUIVALENCE (FINIT(1),WRCBAS) 202 C 203 C + + + END SPECIFICATIONS + + + 204 C 205 Cmeb BIGQ = 10.**BIGLOG meb02/92 206 BIGQ = 10.**INT(BIGLOG) 207 Cprh DO10J=1,5 208 Cprh DO10I=1,31 209 Cprh 10 FC(I,J)=SIGN(BIGLOG,FLOAT(I)-15.5) Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 8 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 210 Cprh init common variables in WCFCM2 (include file cwcf2.inc) 211 DO 10 I=1,MXINT 212 WRCFC(I) = SIGN(BIGLOG,FLOAT(I)-15.5) 213 EPFC(I) = SIGN(BIGLOG,FLOAT(I)-15.5) 214 CLIML(I) = SIGN(BIGLOG,FLOAT(I)-15.5) 215 CLIMU(I) = SIGN(BIGLOG,FLOAT(I)-15.5) 216 SYSRFC(I)= SIGN(BIGLOG,FLOAT(I)-15.5) 217 10 CONTINUE 218 DO20I=1,30 219 20 FINIT(I)=BIGQ 220 C 221 IF(MSL.GE.4) WRITE(MSG,101) GENSKU,IGSOPT,RMSEGS,GAGEB, 222 $ QLWOUT,QHIOUT,NHIST,HISTPD 223 101 FORMAT(4X,'WCF101L-INPUT PARAMS- GENSKU OPT STD-ERR GAGEB ', 224 $ 'QLWOUT QHIOUT NHIST HISTPD'/ 225 $ 26X,F6.3,1X,I2,1X,F6.3,1X,F6.1,1X,F7.1,1X,F9.1,1X,I4,1X,F7.1) 226 C 227 IF(NPK.GT.NHIST.AND.NHIST.GE.0)GOTO30 228 IER=3 229 IF(MSL.GT.0)WRITE(MSG,102)NPK,NHIST 230 102 FORMAT(44H ***WCF102E-INVALID PEAK COUNTS. NPK,NHIST = ,2I10) 231 RETURN 232 30 CONTINUE 233 C 234 IF(MSL.LT.4)GOTO40 235 IF(NOTRAN.NE.1)WRITE(MSG,103)NPK,PKQ 236 IF(NOTRAN.EQ.1)WRITE(MSG,104)NPK,PKLOG 237 103 FORMAT(50H WCF103L-INPUT PEAKS,HISTORIC FIRST. TOTAL NO = ,I6/ 238 $ (14X,5F12.1)) 239 104 FORMAT(50H WCF104L-INPUT LOG PEAKS,HIST FIRST. TOTAL NO = ,I6/ 240 $ (9X,10F7.3)) 241 40 CONTINUE 242 C 243 IF(GENSKU.GE.GSMIN.AND.GENSKU.LE.GSMAX)GOTO50 244 IER=MAX0(IER,1) 245 IF(MSL.GE.2)WRITE(MSG,107)GENSKU,GSMIN,GSMAX 246 107 FORMAT(49H *WCF107I-ACCEPTED GEN SKEW OUTSIDE MAP LIMITS.,3F8.3) 247 50 CONTINUE 248 C 249 IF(NOPPOS.EQ.1)GOTO70 250 WRCPP(1)=BIGQ 251 IF(NOSYS.EQ.1)GOTO70 252 J=NHIST+1 253 DO60I=1,J 254 60 SYSPP(I)=-1. 255 SYSPP(J)=BIGQ 256 70 CONTINUE 257 C 258 NMISS=0 259 IF(NOTRAN.EQ.1)GOTO110 260 DO90I=1,NPK 261 T=PKQ(I) Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 9 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 262 IF(T.LE.0.)GOTO80 263 PKLOG(I)=ALOG10(T) 264 GOTO90 265 80 IF(T.LT.0.)GOTO85 266 PKLOG(I)=-(BIGLOG+.001) 267 GOTO90 268 85 NMISS=NMISS+1 269 IF(I.LE.NHIST)IER=13 270 IF(T.GT.-BIGLOG-1.)T=-777777. 271 PKLOG(I)=T 272 IF(NOPPOS.EQ.1)GOTO90 273 J=NPK+1-NMISS 274 WRCPP(J)=T 275 IF(NOSYS.NE.1)SYSPP(I)=T 276 90 CONTINUE 277 C 278 C wording of following 2 messages changed by AML 8/93 after 279 C exchange of edoc's between Raymond Slade and Bill Kirby 280 IF(NMISS.EQ.0)GOTO100 281 IER=MAX0(IER,2) 282 IF(MSL.GT.0)WRITE(MSG,109)NMISS 283 109 FORMAT( ' **WCF109W-PEAKS WITH MINUS-FLAGGED DISCHARGES ', 284 $ 'WERE BYPASSED. ',I6) 285 IF(IER.LE.3)GOTO100 286 IER=3 287 IF(MSL.GT.0)WRITE(MSG,111) 288 111 FORMAT( ' ***WCF111E-HISTORIC PEAK HAD MINUS-FLAGGED ', 289 $ ' DISCHARGE') 290 100 CONTINUE 291 110 CONTINUE 292 C 293 NSYS=NPK-NHIST-NMISS 294 IF(NMISS.GT.0 .AND. MSL.GT.0)WRITE(MSG,113)NSYS 295 Ckmf revised format because of problem with Lahey compiler 296 C 113 FORMAT(56H **WCF113W-NUMBER OF SYSTEMATIC PEAKS HAS BEEN REDUCED 297 C $ ,10HTO NSYS = ,I4) 298 113 FORMAT( ' **WCF113W-NUMBER OF SYSTEMATIC PEAKS HAS BEEN ', 299 $ 'REDUCED TO NSYS = ',I4) 300 IF(NSYS.GT.0)GOTO130 301 IER=3 302 IF(MSL.GT.0)WRITE(MSG,117)NSYS,NPK,NHIST,NMISS 303 117 FORMAT(42H ***WCF117E-NO DATA IN SYSTEMATIC RECORD. , 304 $ 25H NSYS,NPK,NHIST,NMISS= / 30X,4I10) 305 GOTO140 306 130 IF(NSYS.GE.10)GOTO140 307 IER=MAX0(IER,2) 308 IF(MSL.GT.0)WRITE(MSG,118)NSYS 309 118 FORMAT( ' **WCF118W-SYSTEMATIC RECORD SHORTER THAN 17B SPEC.',I6) 310 140 CONTINUE 311 RETURN 312 END Bytes of stack required for this program unit: 192. Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 10 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- BIGLOG (REAL) scalar 316 (I1)4s (I1)12s 206r 212r 213r 214r 215r 216r 266u 270u BIGQ (REAL) scalar 195s 206= 219u 250u 255u CLIML (REAL) array 256 (I2)2s (I2)4s 214= CLIMU (REAL) array 384 (I2)3s (I2)5s 215= CLSIZE (REAL) scalar 40 (I1)3s (I1)11s CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s EPFC (REAL) array 128 (I2)2s (I2)4s 213= EPS1 (REAL) scalar 320 (I1)5s (I1)13s EPS2 (REAL) scalar 324 (I1)5s (I1)13s FINIT (REAL) array 32 195s 201u 219= GAGEB (REAL) scalar 12 (I3)2s (I3)10s 221o GENSKU (REAL) scalar 0 (I3)2s (I3)10s 221o 243u 243u 245o GSKWGT (REAL) scalar 108 (I3)6s (I3)10s GSMAX (REAL) scalar 32 (I1)3s (I1)11s 243u 245o GSMIN (REAL) scalar 28 (I1)3s (I1)11s 243u 245o HISTPD (REAL) scalar 28 (I3)3s (I3)10s 222o HISTPN (REAL) scalar 100 (I3)6s (I3)10s HISTWT (REAL) scalar 104 (I3)6s (I3)10s HRECWO (REAL) scalar 328 (I1)5s (I1)13s HRECWS (REAL) scalar 336 (I1)5s (I1)13s HRECWX (REAL) scalar 332 (I1)5s (I1)13s I (INTEGER) scalar 194s 212u 212r 213u 213r 214u 214r 215u 215r 216u 216r 219u 254u 261u 263u 266u 269u 271u 275u IER (INTEGER) scalar 162s 169s 228= 244r 244r 244= 269= 281r 281r 281= 285u 286= 301= 307r 307r 307= IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s 221o INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s INDXPT (INTEGER) array 56 (I1)4s (I1)9s IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s J (INTEGER) scalar 194s 252= 253u 255u 273= 274u MSG (INTEGER) scalar 0 (I1)2s (I1)8s 221o 229o 235o 236o 245o 282o 287o 294o 302o 308o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 221u 229u 234u 245u 282u 287u 294u 302u 308u MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u (I2)4u (I2)4u (I2)4u (I2)5u (I2)5u 211u NBGB (INTEGER) scalar 84 (I3)5s (I3)9s NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s NHIST (INTEGER) scalar 24 (I3)2s (I3)9s 222o 227u 227u 229o 252u 269u 293u 302o NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 11 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- NINDX (INTEGER) scalar 52 (I1)4s (I1)9s NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s NMISS (INTEGER) scalar 76 (I3)5s (I3)9s 258= 268u 268= 273u 280u 282o 293u 294u 302o NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s 249u 272u NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s 251u 275u NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s 235u 236u 259u NPK (INTEGER) scalar 162s 169s 170u 170u 170u 170u 227u 229o 235o 236o 260u 273u 293u 302o NSYS (INTEGER) scalar 80 (I3)5s (I3)9s 293= 294o 300u 302o 306u 308o PKLOG (REAL) array 162s 170s 236o 263= 266= 271= PKQ (REAL) array 162s 170s 235o 261u QHIOUT (REAL) scalar 20 (I3)2s (I3)11s 222o QLWOUT (REAL) scalar 16 (I3)2s (I3)11s 222o RMSDGS (REAL) scalar 340 (I1)6s (I1)14s RMSEGS (REAL) scalar 8 (I3)2s (I3)10s 221o SIGHOT (REAL) scalar 360 (I1)7s (I1)15s SIGLOT (REAL) scalar 356 (I1)7s (I1)15s SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s SYSMAX (REAL) scalar 148 (I3)8s (I3)11s SYSMIN (REAL) scalar 144 (I3)8s (I3)12s SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSPP (REAL) array 162s 170s 254= 255= 275= SYSRFC (REAL) array 512 (I2)3s (I2)5s 216= SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSUAV (REAL) scalar 120 (I3)6s (I3)12s SYSUSD (REAL) scalar 124 (I3)7s (I3)12s T (REAL) scalar 195s 261= 262u 263r 265u 270u 270= 271u 274u 275u TXPROB (REAL) array 184 (I1)4s (I1)12s WCFAPI SUBROUTINE 161s WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s WORK (REAL) array 200 (I3)8s (I3)12s WRCAAV (REAL) scalar 52 (I3)4s (I3)13s WRCASD (REAL) scalar 56 (I3)4s (I3)13s WRCASK (REAL) scalar 60 (I3)4s (I3)13s WRCBAS (REAL) scalar 32 (I3)3s (I3)13s 201u WRCFC (REAL) array 0 (I2)2s (I2)4s 212= WRCHHB (REAL) scalar 64 (I3)4s (I3)13s WRCHOT (REAL) scalar 68 (I3)4s (I3)13s WRCLOW (REAL) scalar 72 (I3)5s (I3)13s WRCPAB (REAL) scalar 36 (I3)3s (I3)14s WRCPP (REAL) array 162s 170s 250= 274= WRCSKW (REAL) scalar 48 (I3)4s (I3)14s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 12 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- WRCUAV (REAL) scalar 40 (I3)3s (I3)14s WRCUSD (REAL) scalar 44 (I3)4s (I3)14s WSKLAT (REAL) scalar 312 (I1)4s (I1)12s -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 10 211d 217s 20 218d 219s 30 227g 232s 40 234g 241s 50 243g 247s 60 253d 254s 70 249g 251g 256s 80 262g 265s 90 260d 264g 267g 272g 276s 85 265g 268s 100 280g 285g 290s 110 259g 291s 130 300g 306s 140 305g 306g 310s 101 221f 223s 111 287f 288s 102 229f 230s 103 235f 237s 113 294f 298s 104 236f 239s 107 245f 246s 117 302f 303s 118 308f 309s 109 282f 283s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 13 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFASP Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 313 C 314 C 315 C 316 SUBROUTINE WCFASP 317 # (SYSLOG,SYSPP,NSYS1,IER) 318 C 319 C + + + PURPOSE + + + 320 C ANALYZE SYSTEMATIC PEAKS 321 C (WRC Bulletin-17 Flood Frequency Analysis) 322 C REV 7/20/79 WK TO MAKE ADDENDS TO SUMS DOUBLE PRECISION 323 C 324 C + + + DUMMY ARGUMENTS + + + 325 INTEGER NSYS1,IER 326 REAL SYSLOG(NSYS1),SYSPP(NSYS1) 327 C 328 C + + + ARGUMENT DEFINITIONS + + + 329 C SYSLOG - systematic peak logarithms 330 C SYSPP - prob plot positions systematic record PKS 331 C NSYS1 - number of systematic peaks 332 C IER - error return code 333 C 334 C + + + PARAMETERS + + + 335 INCLUDE 'pmxint.inc' 336 C 337 C + + + COMMON BLOCKS + + + 338 INCLUDE 'cwcf0.inc' 339 INCLUDE 'cwcf1.inc' 340 INCLUDE 'cwcf2.inc' 341 C 342 C + + + LOCAL VARIABLES + + + 343 INTEGER NAGB, I 344 REAL T, SYSTAT(8), WSTATS(8) 345 DOUBLE PRECISION Z 346 C 347 C + + + EQUIVALENCE + + + 348 EQUIVALENCE (WSTATS(1),WRCBAS),(SYSTAT(1),SYSBAS) 349 C 350 C + + + INTRINSICS + + + 351 INTRINSIC ALOG10 352 C 353 C + + + EXTERNALS + + + 354 EXTERNAL WCFCSA, WCFEPP, WCFFCA 355 C 356 C + + + END SPECIFICATIONS + + + 357 C 358 NLWOUT=0 359 NHIOUT=0 360 NHISTN=0 361 HISTPN=NSYS 362 SYSMAX=-BIGLOG 363 SYSMIN=+BIGLOG 364 NAGB=0 Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 14 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 365 DO10I=1,3 366 SUMH(I)=0D0 367 10 SUMS(I)=0D0 368 C 369 WRCBAS=-BIGLOG 370 IF(GAGEB.GT.0.)WRCBAS=ALOG10(GAGEB) 371 SYSBAS=WRCBAS 372 DO30I=1,NSYS1 373 T=SYSLOG(I) 374 IF(T.LE.WRCBAS)GOTO30 375 NAGB=NAGB+1 376 IF(T.LT.SYSMIN)SYSMIN=T 377 IF(T.GT.SYSMAX)SYSMAX=T 378 Z = T 379 SUMS(1)=SUMS(1)+Z 380 SUMS(2)=SUMS(2)+Z**2 381 SUMS(3)=SUMS(3)+Z**3 382 30 CONTINUE 383 NBGB=NSYS-NAGB 384 C 385 IF(MSL.LT.2)GOTO40 386 T=10.**WRCBAS 387 IF(NBGB.GT.0)WRITE(MSG,133)NBGB,T 388 IF(NBGB.LE.0)WRITE(MSG,134)T 389 133 FORMAT(56H WCF133I-SYSTEMATIC PEAKS BELOW GAGE BASE WERE NOTED. 390 $ ,I10,1X,F8.1) 391 134 FORMAT(53H WCF134I-NO SYSTEMATIC PEAKS WERE BELOW GAGE BASE., 392 $ 14X,F8.1) 393 40 CONTINUE 394 C 395 Ckmf CALL WCFCSA (4HSYS , IER ) 396 CALL WCFCSA ( 'SYS ', IER ) 397 IF(IER.GE.3) RETURN 398 C 399 IF(NOSYS.EQ.1) RETURN 400 Ckmf IF(NOPPOS.NE.1) CALL WCFEPP (SYSPP, 4HSYS ,NSYS) 401 Ckmf CALL WCFFCA (SYSRFC,4HSYS , IER) 402 IF(NOPPOS.NE.1) CALL WCFEPP (SYSPP, 'SYS ', NSYS) 403 CALL WCFFCA (SYSRFC, 'SYS ', IER ) 404 C STATISTICS ARE REPORTED BY WCFFCA IF REQUIRED. 405 DO 50 I=1,8 406 50 SYSTAT(I)=WSTATS(I) 407 RETURN 408 END Bytes of stack required for this program unit: 40. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- BIGLOG (REAL) scalar 316 (I1)4s (I1)12s 362u 363u 363u 369u Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 15 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- CLIML (REAL) array 256 (I2)2s (I2)4s CLIMU (REAL) array 384 (I2)3s (I2)5s CLSIZE (REAL) scalar 40 (I1)3s (I1)11s CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s EPFC (REAL) array 128 (I2)2s (I2)4s EPS1 (REAL) scalar 320 (I1)5s (I1)13s EPS2 (REAL) scalar 324 (I1)5s (I1)13s GAGEB (REAL) scalar 12 (I3)2s (I3)10s 370u 370r GENSKU (REAL) scalar 0 (I3)2s (I3)10s GSKWGT (REAL) scalar 108 (I3)6s (I3)10s GSMAX (REAL) scalar 32 (I1)3s (I1)11s GSMIN (REAL) scalar 28 (I1)3s (I1)11s HISTPD (REAL) scalar 28 (I3)3s (I3)10s HISTPN (REAL) scalar 100 (I3)6s (I3)10s 361= HISTWT (REAL) scalar 104 (I3)6s (I3)10s HRECWO (REAL) scalar 328 (I1)5s (I1)13s HRECWS (REAL) scalar 336 (I1)5s (I1)13s HRECWX (REAL) scalar 332 (I1)5s (I1)13s I (INTEGER) scalar 343s 366u 367u 373u 406u 406u IER (INTEGER) scalar 317s 325s 396r 396r 396r 397u 403r 403r 403r IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s INDXPT (INTEGER) array 56 (I1)4s (I1)9s IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s MSG (INTEGER) scalar 0 (I1)2s (I1)8s 387o 388o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 385u MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u (I2)4u (I2)4u (I2)4u (I2)5u (I2)5u NAGB (INTEGER) scalar 343s 364= 375u 375= 383u NBGB (INTEGER) scalar 84 (I3)5s (I3)9s 383= 387u 387o 388u NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s 359= NHIST (INTEGER) scalar 24 (I3)2s (I3)9s NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s 360= NINDX (INTEGER) scalar 52 (I1)4s (I1)9s NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s 358= NMISS (INTEGER) scalar 76 (I3)5s (I3)9s NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s 402u NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s 399u NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s NSYS (INTEGER) scalar 80 (I3)5s (I3)9s 361u 383u 402r 402r 402r NSYS1 (INTEGER) scalar 317s 325s 326u 326u 372u QHIOUT (REAL) scalar 20 (I3)2s (I3)11s QLWOUT (REAL) scalar 16 (I3)2s (I3)11s RMSDGS (REAL) scalar 340 (I1)6s (I1)14s RMSEGS (REAL) scalar 8 (I3)2s (I3)10s SIGHOT (REAL) scalar 360 (I1)7s (I1)15s SIGLOT (REAL) scalar 356 (I1)7s (I1)15s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 16 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s 366= SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s 367= 379u 379= 380u 380= 381u 381= SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s 348u 371= SYSLOG (REAL) array 317s 326s 373u SYSMAX (REAL) scalar 148 (I3)8s (I3)11s 362= 377u 377= SYSMIN (REAL) scalar 144 (I3)8s (I3)12s 363= 376u 376= SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSPP (REAL) array 317s 326s 402r 402r 402r SYSRFC (REAL) array 512 (I2)3s (I2)5s 403r 403r 403r SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSTAT (REAL) array 112 344s 348u 406= SYSUAV (REAL) scalar 120 (I3)6s (I3)12s SYSUSD (REAL) scalar 124 (I3)7s (I3)12s T (REAL) scalar 344s 373= 374u 376u 376u 377u 377u 378u 386= 387o 388o TXPROB (REAL) array 184 (I1)4s (I1)12s WCFASP SUBROUTINE 316s WCFCSA SUBROUTINE 354u 396u WCFEPP SUBROUTINE 354u 402u WCFFCA SUBROUTINE 354u 403u WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s WORK (REAL) array 200 (I3)8s (I3)12s WRCAAV (REAL) scalar 52 (I3)4s (I3)13s WRCASD (REAL) scalar 56 (I3)4s (I3)13s WRCASK (REAL) scalar 60 (I3)4s (I3)13s WRCBAS (REAL) scalar 32 (I3)3s (I3)13s 348u 369= 370= 371u 374u 386u WRCFC (REAL) array 0 (I2)2s (I2)4s WRCHHB (REAL) scalar 64 (I3)4s (I3)13s WRCHOT (REAL) scalar 68 (I3)4s (I3)13s WRCLOW (REAL) scalar 72 (I3)5s (I3)13s WRCPAB (REAL) scalar 36 (I3)3s (I3)14s WRCSKW (REAL) scalar 48 (I3)4s (I3)14s WRCUAV (REAL) scalar 40 (I3)3s (I3)14s WRCUSD (REAL) scalar 44 (I3)4s (I3)14s WSKLAT (REAL) scalar 312 (I1)4s (I1)12s WSTATS (REAL) array 32 344s 348u 406u Z (REAL(kind 8)) scalar 345s 378= 379u 380u 381u -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 10 365d 367s 30 372d 374g 382s 40 385g 393s 50 405d 406s 133 387f 389s 134 388f 391s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 17 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 18 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFCSA Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 409 C 410 C 411 C 412 SUBROUTINE WCFCSA 413 # (LABEL,IER) 414 C 415 C + + + PURPOSE + + + 416 C CALCULATE STATISTICS OF ABOVE-BASE PEAKS 417 C (WRC Bulletin-17 Flood Frequency Analysis) 418 C INCLUDING WRC HISTORIC ADJ, WHEN HIST INFO EXISTS. 419 C NOTE -- THESE FORMULAS ARE OK EVEN IF NHIOUT.GT.0 AND HISTPN.EQ.NSYS. 420 C (I.E., EVEN IF WCFDHH HAS DETECTED HIGH OUTLIERS BUT NO HISTORIC 421 C INFO HAS BEEN GIVEN.) IN THIS CASE, HISTWT = 1. 422 C 423 C REV 5/79 WK TO SIMPLIFY FN CALC AND TO TRAP NEGATIVE VARIANCES. 424 C 425 C + + + HISTORY + + + 426 C kmf - nov 9, 2000 - changed label from Hollerith to character 427 C 428 C + + + DUMMY ARGUMENTS + + + 429 Ckmf INTEGER LABEL, IER 430 INTEGER IER 431 CHARACTER*4 LABEL 432 C 433 C + + + ARGUMENTS DEFINITIONS + + + 434 C LABEL - input identification label printed in error message, 435 C indicates point from which routine was called (used 436 C only in error messages) 437 C IER - error return code 438 C 439 C + + + PARAMETERS + + + 440 INCLUDE 'pmxint.inc' 441 C 442 C + + + COMMON BLOCKS + + + 443 INCLUDE 'cwcf0.inc' 444 INCLUDE 'cwcf1.inc' 445 INCLUDE 'cwcf2.inc' 446 C 447 C + + + LOCAL VARIABLES + + + 448 INTEGER I 449 REAL FN, Z, V 450 DOUBLE PRECISION S(3),DBLAAV 451 C 452 C + + + INTRINSICS + + + 453 INTRINSIC FLOAT, SQRT 454 C 455 C + + + END SPECIFICATIONS + + + 456 C 457 FN =NSYS-NHIOUT-NBGB-NLWOUT 458 IF (FN.LT.3.)GOTO97 459 Z=NHISTN+NHIOUT 460 HISTWT=(HISTPN-Z)/FLOAT(NSYS-NHIOUT) Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 19 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 461 FN = Z + HISTWT*FN 462 DO50I=1,3 463 50 S(I)=SUMS(I)*HISTWT + SUMH(I) 464 C 465 WRCPAB=FN/HISTPN 466 DBLAAV=S(1)/FN 467 WRCAAV=DBLAAV 468 V=(S(2)-S(1)**2/FN)/(FN-1.) 469 IF(V.GE.0.) GO TO 60 470 IER=3 471 IF(MSL.GT.0) WRITE(MSG, 143) LABEL, V 472 143 FORMAT(38H ***WCF143E-NEGATIVE VARIANCE OF LOGS ,1A4,F10.5) 473 RETURN 474 60 CONTINUE 475 WRCASD=SQRT(V) 476 WRCASK=(S(3)-3.*DBLAAV*S(2) + 2.*DBLAAV**2*S(1))*FN 477 $ / ((FN-1.)*(FN-2.)*WRCASD**3) 478 RETURN 479 C 480 97 IER=3 481 IF(MSL.GT.0)WRITE(MSG,141) LABEL, 482 $ FN,NSYS,NBGB,NLWOUT,NHIOUT,NHISTN,HISTPN 483 141 FORMAT(48H ***WCF141E-SAMPLE SIZE TOO SMALL TO CALC STATS. ,1A6, 484 $ F8.1/11X, 37HNSYS,NBGB,NLWOUT,NHIOUT,NHISTN,HISTPN / 485 $ 8X,5I6, F10.1) 486 RETURN 487 END Bytes of stack required for this program unit: 40. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- BIGLOG (REAL) scalar 316 (I1)4s (I1)12s CLIML (REAL) array 256 (I2)2s (I2)4s CLIMU (REAL) array 384 (I2)3s (I2)5s CLSIZE (REAL) scalar 40 (I1)3s (I1)11s CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s DBLAAV (REAL(kind 8)) scalar 450s 466= 467u 476u 476u EPFC (REAL) array 128 (I2)2s (I2)4s EPS1 (REAL) scalar 320 (I1)5s (I1)13s EPS2 (REAL) scalar 324 (I1)5s (I1)13s FN (REAL) scalar 449s 457= 458u 461u 461= 465u 466u 468u 468u 476u 477u 477u 482o GAGEB (REAL) scalar 12 (I3)2s (I3)10s GENSKU (REAL) scalar 0 (I3)2s (I3)10s GSKWGT (REAL) scalar 108 (I3)6s (I3)10s GSMAX (REAL) scalar 32 (I1)3s (I1)11s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 20 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- GSMIN (REAL) scalar 28 (I1)3s (I1)11s HISTPD (REAL) scalar 28 (I3)3s (I3)10s HISTPN (REAL) scalar 100 (I3)6s (I3)10s 460u 465u 482o HISTWT (REAL) scalar 104 (I3)6s (I3)10s 460= 461u 463u HRECWO (REAL) scalar 328 (I1)5s (I1)13s HRECWS (REAL) scalar 336 (I1)5s (I1)13s HRECWX (REAL) scalar 332 (I1)5s (I1)13s I (INTEGER) scalar 448s 463u 463u 463u IER (INTEGER) scalar 413s 430s 470= 480= IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s INDXPT (INTEGER) array 56 (I1)4s (I1)9s IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s LABEL (CHARACTER) scalar 413s 431s 471o 471o 481o 481o MSG (INTEGER) scalar 0 (I1)2s (I1)8s 471o 481o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 471u 481u MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u (I2)4u (I2)4u (I2)4u (I2)5u (I2)5u NBGB (INTEGER) scalar 84 (I3)5s (I3)9s 457u 482o NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s 457u 459u 460u 482o NHIST (INTEGER) scalar 24 (I3)2s (I3)9s NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s 459u 482o NINDX (INTEGER) scalar 52 (I1)4s (I1)9s NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s 457u 482o NMISS (INTEGER) scalar 76 (I3)5s (I3)9s NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s NSYS (INTEGER) scalar 80 (I3)5s (I3)9s 457u 460u 482o QHIOUT (REAL) scalar 20 (I3)2s (I3)11s QLWOUT (REAL) scalar 16 (I3)2s (I3)11s RMSDGS (REAL) scalar 340 (I1)6s (I1)14s RMSEGS (REAL) scalar 8 (I3)2s (I3)10s S (REAL(kind 8)) array 450s 463= 466u 468u 468u 476u 476u 476u SIGHOT (REAL) scalar 360 (I1)7s (I1)15s SIGLOT (REAL) scalar 356 (I1)7s (I1)15s SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s 463u SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s 463u SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s SYSMAX (REAL) scalar 148 (I3)8s (I3)11s SYSMIN (REAL) scalar 144 (I3)8s (I3)12s SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSRFC (REAL) array 512 (I2)3s (I2)5s SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSUAV (REAL) scalar 120 (I3)6s (I3)12s SYSUSD (REAL) scalar 124 (I3)7s (I3)12s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 21 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- TXPROB (REAL) array 184 (I1)4s (I1)12s V (REAL) scalar 449s 468= 469u 471o 475r WCFCSA SUBROUTINE 412s WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s WORK (REAL) array 200 (I3)8s (I3)12s WRCAAV (REAL) scalar 52 (I3)4s (I3)13s 467= WRCASD (REAL) scalar 56 (I3)4s (I3)13s 475= 477u WRCASK (REAL) scalar 60 (I3)4s (I3)13s 476= WRCBAS (REAL) scalar 32 (I3)3s (I3)13s WRCFC (REAL) array 0 (I2)2s (I2)4s WRCHHB (REAL) scalar 64 (I3)4s (I3)13s WRCHOT (REAL) scalar 68 (I3)4s (I3)13s WRCLOW (REAL) scalar 72 (I3)5s (I3)13s WRCPAB (REAL) scalar 36 (I3)3s (I3)14s 465= WRCSKW (REAL) scalar 48 (I3)4s (I3)14s WRCUAV (REAL) scalar 40 (I3)3s (I3)14s WRCUSD (REAL) scalar 44 (I3)4s (I3)14s WSKLAT (REAL) scalar 312 (I1)4s (I1)12s Z (REAL) scalar 449s 459= 460u 461u -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 50 462d 463s 60 469g 474s 97 458g 480s 141 481f 483s 143 471f 472s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 22 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFCWS Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 488 C 489 C 490 C 491 SUBROUTINE WCFCWS 492 # (IRC) 493 C 494 C + + + PURPOSE + + + 495 C CALCULATE WRC WEIGHTED SKEW 496 C (WRC Bulletin-17 Flood Frequency Analysis) 497 C REV 10/79 WK - REDEFINE 'ROUND' TO HANDLE NEGATIVE X. 498 C REV 1/81 WK - BULL 17-B WEIGHTED SKEW CALC BASED ON MEAN SQUARE ERROR 499 C 500 C + + + DUMMY ARGUMENTS + + + 501 INTEGER IRC 502 C 503 C + + + ARGUMENTS DEFINITIONS + + + 504 C IRC - return code, 0, 1, 2, or 3: 505 C 0 - no error 506 C 3 - error, calculation aborted 507 C 508 C + + + PARAMETERS + + + 509 INCLUDE 'pmxint.inc' 510 C 511 C + + + COMMON BLOCKS + + + 512 INCLUDE 'cwcf0.inc' 513 INCLUDE 'cwcf1.inc' 514 INCLUDE 'cwcf2.inc' 515 C 516 C + + + LOCAL VARIABLES + + + 517 REAL ERL, AG, A, S, VSTA, VGS, T, HOLD, X 518 C 519 C + + + FUNCTIONS + + + 520 REAL ROUND 521 C 522 C + + + INTRINSICS + + + 523 INTRINSIC ALOG10, SIGN, MAX0, AINT, ABS 524 C 525 C + + + END SPECIFICATIONS + + + 526 C 527 ROUND(X) = SIGN( AINT(.5+ABS(X)/WSKLAT)*WSKLAT, X ) 528 C 529 ERL = NSYS 530 ERL = ERL + HRECWS*(HISTPN-ERL) 531 C 532 C MSE OF STATION SKEW 533 AG = ABS(WRCSKW) 534 A = -0.33 + 0.08*AG 535 IF(AG.GT.0.90) A = -0.52 + 0.30*AG 536 S = 0.94 - 0.26*AG 537 IF(AG.GT.1.50) S = 0.55 538 VSTA = 10.**(A - S*(ALOG10(ERL)-1.)) 539 C Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 23 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 540 C MSE OF GEN SKEW 541 VGS = RMSEGS**2 542 IF(RMSEGS.LE.0.) VGS = RMSDGS**2 543 C 544 C WEIGHT 545 GSKWGT = VSTA/(VSTA+VGS) 546 C 547 T = WRCSKW 548 WRCSKW = WRCSKW + GSKWGT*(GENSKU-WRCSKW) 549 IF(WSKLAT.GT.0.) WRCSKW = ROUND(WRCSKW) 550 IF(IGSOPT.EQ.0)RETURN 551 C 552 IRC=MAX0(IRC,1) 553 HOLD = WRCSKW 554 IF(IGSOPT.GT.0)GOTO10 555 GSKWGT=0. 556 WRCSKW = T 557 GOTO20 558 10 GSKWGT=1. 559 WRCSKW=GENSKU 560 20 CONTINUE 561 IF(WSKLAT.GT.0.) WRCSKW = ROUND(WRCSKW) 562 C 563 IF(MSL.GE.2) WRITE(MSG,151) HOLD, WRCSKW, IGSOPT 564 151 FORMAT( ' *WCF151I-17B WEIGHTED SKEW REPLACED BY USER OPTION.', 565 $ 2F10.3,I4) 566 RETURN 567 END Bytes of stack required for this program unit: 80. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- A (REAL) scalar 517s 534= 535= 538u AG (REAL) scalar 517s 533= 534u 535u 535u 536u 537u BIGLOG (REAL) scalar 316 (I1)4s (I1)12s CLIML (REAL) array 256 (I2)2s (I2)4s CLIMU (REAL) array 384 (I2)3s (I2)5s CLSIZE (REAL) scalar 40 (I1)3s (I1)11s CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s EPFC (REAL) array 128 (I2)2s (I2)4s EPS1 (REAL) scalar 320 (I1)5s (I1)13s EPS2 (REAL) scalar 324 (I1)5s (I1)13s ERL (REAL) scalar 517s 529= 530u 530u 530= 538r GAGEB (REAL) scalar 12 (I3)2s (I3)10s GENSKU (REAL) scalar 0 (I3)2s (I3)10s 548u 559u GSKWGT (REAL) scalar 108 (I3)6s (I3)10s 545= 548u 555= 558= GSMAX (REAL) scalar 32 (I1)3s (I1)11s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 24 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- GSMIN (REAL) scalar 28 (I1)3s (I1)11s HISTPD (REAL) scalar 28 (I3)3s (I3)10s HISTPN (REAL) scalar 100 (I3)6s (I3)10s 530u HISTWT (REAL) scalar 104 (I3)6s (I3)10s HOLD (REAL) scalar 517s 553= 563o HRECWO (REAL) scalar 328 (I1)5s (I1)13s HRECWS (REAL) scalar 336 (I1)5s (I1)13s 530u HRECWX (REAL) scalar 332 (I1)5s (I1)13s IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s 550u 554u 563o INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s INDXPT (INTEGER) array 56 (I1)4s (I1)9s IRC (INTEGER) scalar 492s 501s 552r 552r 552= IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s MSG (INTEGER) scalar 0 (I1)2s (I1)8s 563o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 563u MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u (I2)4u (I2)4u (I2)4u (I2)5u (I2)5u NBGB (INTEGER) scalar 84 (I3)5s (I3)9s NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s NHIST (INTEGER) scalar 24 (I3)2s (I3)9s NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s NINDX (INTEGER) scalar 52 (I1)4s (I1)9s NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s NMISS (INTEGER) scalar 76 (I3)5s (I3)9s NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s NSYS (INTEGER) scalar 80 (I3)5s (I3)9s 529u QHIOUT (REAL) scalar 20 (I3)2s (I3)11s QLWOUT (REAL) scalar 16 (I3)2s (I3)11s RMSDGS (REAL) scalar 340 (I1)6s (I1)14s 542u RMSEGS (REAL) scalar 8 (I3)2s (I3)10s 541u 542u ROUND FUNCTION 549u 561u ROUND (REAL) scalar 520s S (REAL) scalar 517s 536= 537= 538u SIGHOT (REAL) scalar 360 (I1)7s (I1)15s SIGLOT (REAL) scalar 356 (I1)7s (I1)15s SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s SYSMAX (REAL) scalar 148 (I3)8s (I3)11s SYSMIN (REAL) scalar 144 (I3)8s (I3)12s SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSRFC (REAL) array 512 (I2)3s (I2)5s SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSUAV (REAL) scalar 120 (I3)6s (I3)12s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 25 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- SYSUSD (REAL) scalar 124 (I3)7s (I3)12s T (REAL) scalar 517s 547= 556u TXPROB (REAL) array 184 (I1)4s (I1)12s VGS (REAL) scalar 517s 541= 542= 545u VSTA (REAL) scalar 517s 538= 545u 545u WCFCWS SUBROUTINE 491s WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s WORK (REAL) array 200 (I3)8s (I3)12s WRCAAV (REAL) scalar 52 (I3)4s (I3)13s WRCASD (REAL) scalar 56 (I3)4s (I3)13s WRCASK (REAL) scalar 60 (I3)4s (I3)13s WRCBAS (REAL) scalar 32 (I3)3s (I3)13s WRCFC (REAL) array 0 (I2)2s (I2)4s WRCHHB (REAL) scalar 64 (I3)4s (I3)13s WRCHOT (REAL) scalar 68 (I3)4s (I3)13s WRCLOW (REAL) scalar 72 (I3)5s (I3)13s WRCPAB (REAL) scalar 36 (I3)3s (I3)14s WRCSKW (REAL) scalar 48 (I3)4s (I3)14s 533r 547u 548u 548u 548= 549r 549r 549u 549= 553u 556= 559= 561r 561r 561u 561= 563o WRCUAV (REAL) scalar 40 (I3)3s (I3)14s WRCUSD (REAL) scalar 44 (I3)4s (I3)14s WSKLAT (REAL) scalar 312 (I1)4s (I1)12s 527u 527u 549u 561u -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 10 554g 558s 20 557g 560s 151 563f 564s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 26 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFDHH Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 568 C 569 C 570 C 571 SUBROUTINE WCFDHH 572 # (PKLOG,NDIM,IER) 573 C 574 C + + + PURPOSE + + + 575 C DETECT HISTORIC PEAKS AND HIGH OUTLIERS 576 C (WRC Bulletin-17 Flood Frequency Analysis) 577 C REV 7/20/79 WK TO MAKE ADDENDS TO SUMS DOUBLE PRECISION 578 C REV 1/81 WK - BULL 17-B REVISED OUTLIER TEST 579 C REV 4/28/81 WK - TO USE NUMBER OF PEAKS ABOVE FLOOD BASE IN HI-OUT 580 C TEST. ALSO PRINT HI-OUT TEST BEFORE SUPERSEDING WITH MIN HIST PEAK 581 C 582 C + + + DUMMY ARGUMENTS + + + 583 INTEGER NDIM, IER 584 REAL PKLOG(NDIM) 585 C 586 C + + + ARGUMENT DEFINITIONS + + + 587 C PKLOG - peak flow LOG10 work-input vector 588 C NDIM - size of PKLOG array 589 C IER - error return code 590 C 591 C + + + PARAMETERS + + + 592 INCLUDE 'pmxint.inc' 593 C 594 C + + + COMMON BLOCKS + + + 595 INCLUDE 'cwcf0.inc' 596 INCLUDE 'cwcf1.inc' 597 C 598 C + + + LOCAL VARIABLES + + + 599 INTEGER I,ISYS,J 600 REAL HPMIN,T,U 601 DOUBLE PRECISION Z 602 C 603 C + + + FUNCTIONS + + + 604 REAL OUTKGB 605 C 606 C + + + INTRINSICS + + + 607 INTRINSIC MAX0, ALOG10, AMIN1, FLOAT 608 C 609 C + + + EXTERNALS + + + 610 EXTERNAL OUTKGB 611 C 612 C + + + END SPECIFICATIONS + + + 613 C 614 C PRELIM CHECK NO. 153 IS NO LONGER NEEDED WITH NEW BULL 17-B 615 C HIGH OUTLIER TEST. SEE MSG 164 AND 167 FOR BAD HISTORIC INFO. 616 C 617 C 618 C NOTE HIST PKS 619 HPMIN = BIGLOG Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 27 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 620 IF(NHIST.LE.0)GOTO40 621 NHISTN=NHIST 622 DO 30 I=1,NHISTN 623 T=PKLOG(I) 624 IF(T.LT.HPMIN) HPMIN = T 625 Z = T 626 SUMH(1)=SUMH(1)+Z 627 SUMH(2)=SUMH(2)+Z**2 628 30 SUMH(3)=SUMH(3)+Z**3 629 40 CONTINUE 630 C 631 C SET HIGH-OUTLIER BASE 632 WRCHOT = WRCAAV + WRCASD*OUTKGB(SIGHOT, NSYS-NBGB-NLWOUT ) 633 IF(QHIOUT.LE.0.)GOTO60 634 IF(MSL.LT.2) GO TO 44 635 T = 10.**WRCHOT 636 WRITE(MSG,161) QHIOUT, T 637 161 FORMAT(' *WCF161I-USER HIGH OUTLIER CRITERION REPLACES 17B. ', 638 $ F10.1,1X,F10.1 ) 639 44 IER = MAX0(1, IER) 640 T=ALOG10(QHIOUT) 641 IF(T.LE.HPMIN) GO TO 50 642 IER=MAX0(IER,2) 643 U = 10.**HPMIN 644 IF(MSL.GT.0)WRITE(MSG,157)QHIOUT,U 645 157 FORMAT(45H **WCF157W-USER HIGH-OUTLIER CRIT LOWERED TO, 646 $ 13H MIN HIST PK., 2(1X,F9.1)) 647 T = HPMIN 648 50 WRCHHB=T 649 GO TO 66 650 60 CONTINUE 651 WRCHHB = WRCHOT 652 IF(WRCHHB.LE.HPMIN) GO TO 66 653 IF(MSL.LT.2) GO TO 65 654 T = 10.**WRCHOT 655 WRITE(MSG,156) T 656 156 FORMAT(4X,'WCF156I-17B HI-OUTLIER TEST SUPERSEDED BY MIN HIST PK', 657 $ 1X,F9.1) 658 65 WRCHHB = HPMIN 659 66 CONTINUE 660 IF(WRCHHB.GT.WRCBAS)GOTO70 661 IER=3 662 IF(MSL.LE.0)RETURN 663 U=10.**WRCHHB 664 T=10.**WRCBAS 665 WRITE(MSG,159)U,T 666 159 FORMAT(39H ***WCF159E-HIGH-OUT/HIST-PK BASE BELOW , 667 $ 19H LOW-OUT/GAGE BASE., 2(1X,F9.1)) 668 RETURN 669 70 CONTINUE 670 C 671 C DETECT HIGH OUTLIERS Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 28 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 672 IF(WRCHHB.GT.SYSMAX)GOTO100 673 ISYS=NHISTN+1 674 DO90I=ISYS,NDIM 675 T=PKLOG(I) 676 IF(T.LT.WRCHHB)GOTO90 677 NHIOUT=NHIOUT+1 678 Z = T 679 DO 85 J=1,3 680 SUMS(J)=SUMS(J)-Z 681 SUMH(J)=SUMH(J)+Z 682 85 Z = Z*T 683 90 CONTINUE 684 IF(NHIOUT.GT.0)GOTO100 685 IER=3 686 WRITE(MSG,197) SYSMAX, WRCHHB 687 197 FORMAT(/' ***WCF004*-INTERNAL PROGRAM LOGIC ERROR DHH-197 -- ', 688 $ 2F10.5) 689 100 CONTINUE 690 C 691 C REPORT NO HIGH-HIST 692 IF(NHIOUT+NHISTN.GT.0)GOTO110 693 IF(MSL.LT.2) GO TO 16301 694 U = 10.**WRCHHB 695 WRITE(MSG,163) U 696 16301 CONTINUE 697 163 FORMAT(' WCF163I-NO HIGH OUTLIERS OR HISTORIC PEAKS ', 698 $ 'EXCEEDED HHBASE. ',F10.1) 699 IF(HISTPD.LE.0.) RETURN 700 IER=MAX0(IER,2) 701 IF(MSL.LE.0) RETURN 702 WRITE(MSG,164) HISTPD 703 164 FORMAT(' **WCF164W-HISTORIC PERIOD IGNORED. ',F6.1 ) 704 RETURN 705 C 706 C REPORT HIGH-OUT-HIST COUNTS 707 110 CONTINUE 708 IF(HISTPD.LE.0) GO TO 210 709 HISTPN = HISTPD 710 T=10.**WRCHHB 711 IF(MSL.GE.2)WRITE(MSG,165)NHIOUT,NHISTN,T 712 165 FORMAT(44H WCF165I-HIGH OUTLIERS AND HISTORIC PEAKS, 713 $ 14H ABOVE HHBASE., 2(1X,I2),1X,F10.1) 714 IF(HISTPN.GT.FLOAT(NHISTN+NSYS))GOTO120 715 IER=3 716 IF(MSL.GT.0) WRITE(MSG,167) HISTPN, NSYS, NHISTN 717 167 FORMAT(54H ***WCF167E-HIST PERIOD NO LONGER THAN SYS + HIST PKS., 718 $ F9.1,2I5) 719 RETURN 720 120 T=AMIN1(300., 5.*FLOAT(NSYS)) 721 IF(HISTPN.LE.T)GOTO130 722 IER=MAX0(IER,1) 723 IF(MSL.GE.2)WRITE(MSG,169)HISTPN,T Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 29 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 724 169 FORMAT(48H *WCF169I-ACCEPTED HISTORIC PERIOD GTR THAN T.,2F9.1) 725 130 I=NHISTN+NHIOUT 726 IF(10*I.LE.NSYS .OR. I.LE.1) GO TO 140 727 IER=MAX0(IER,2) 728 IF(MSL.GT.0)WRITE(MSG,171)I,NSYS 729 171 FORMAT(42H **WCF171W-NUMBER HI-OUT/HIST PKS EXCEEDS, 730 $ 18H 10PCT OF SYS PKS. , 2I7) 731 140 CONTINUE 732 RETURN 733 210 CONTINUE 734 IF(NHISTN.LE.0) GO TO 215 735 IER = 3 736 WRITE(MSG,1601) NHIST, NHISTN, HISTPD 737 1601 FORMAT(/' ***WCF004*-INTERNAL PROGRAM LOGIC ERROR DHH-1601. ',2I6, 738 $ F7.1/) 739 215 CONTINUE 740 IF(MSL.LT.2) RETURN 741 T = 10.**WRCHHB 742 WRITE(MSG,162) NHIOUT, T 743 162 FORMAT(' WCF162I-SYSTEMATIC PEAKS EXCEEDED HIGH-', 744 $ 'OUTLIER CRITERION. ',I3,1X,F9.1) 745 RETURN 746 END Bytes of stack required for this program unit: 72. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- BIGLOG (REAL) scalar 316 (I1)4s (I1)12s 619u CLSIZE (REAL) scalar 40 (I1)3s (I1)11s CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s EPS1 (REAL) scalar 320 (I1)5s (I1)13s EPS2 (REAL) scalar 324 (I1)5s (I1)13s GAGEB (REAL) scalar 12 (I3)2s (I3)10s GENSKU (REAL) scalar 0 (I3)2s (I3)10s GSKWGT (REAL) scalar 108 (I3)6s (I3)10s GSMAX (REAL) scalar 32 (I1)3s (I1)11s GSMIN (REAL) scalar 28 (I1)3s (I1)11s HISTPD (REAL) scalar 28 (I3)3s (I3)10s 699u 702o 708u 709u 736o HISTPN (REAL) scalar 100 (I3)6s (I3)10s 709= 714u 716o 721u 723o HISTWT (REAL) scalar 104 (I3)6s (I3)10s HPMIN (REAL) scalar 600s 619= 624u 624= 641u 643u 647u 652u 658u HRECWO (REAL) scalar 328 (I1)5s (I1)13s HRECWS (REAL) scalar 336 (I1)5s (I1)13s HRECWX (REAL) scalar 332 (I1)5s (I1)13s I (INTEGER) scalar 599s 623u 675u 725= 726u 726u 728o IER (INTEGER) scalar 572s 583s 639r 639r 639= 642r 642r 642= 661= 685= 700r 700r 700= 715= 722r 722r 722= 727r Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 30 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 727r 727= 735= IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s INDXPT (INTEGER) array 56 (I1)4s (I1)9s ISYS (INTEGER) scalar 599s 673= 674u IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s J (INTEGER) scalar 599s 680u 680u 681u 681u MSG (INTEGER) scalar 0 (I1)2s (I1)8s 636o 644o 655o 665o 686o 695o 702o 711o 716o 723o 728o 736o 742o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 634u 644u 653u 662u 693u 701u 711u 716u 723u 728u 740u MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u NBGB (INTEGER) scalar 84 (I3)5s (I3)9s 632u NDIM (INTEGER) scalar 572s 583s 584u 674u NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s 677u 677= 684u 692u 711o 725u 742o NHIST (INTEGER) scalar 24 (I3)2s (I3)9s 620u 621u 736o NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s 621= 622u 673u 692u 711o 714u 716o 725u 734u 736o NINDX (INTEGER) scalar 52 (I1)4s (I1)9s NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s 632u NMISS (INTEGER) scalar 76 (I3)5s (I3)9s NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s NSYS (INTEGER) scalar 80 (I3)5s (I3)9s 632u 714u 716o 720r 726u 728o OUTKGB FUNCTION 604s 610u 632u PKLOG (REAL) array 572s 584s 623u 675u QHIOUT (REAL) scalar 20 (I3)2s (I3)11s 633u 636o 640r 644o QLWOUT (REAL) scalar 16 (I3)2s (I3)11s RMSDGS (REAL) scalar 340 (I1)6s (I1)14s RMSEGS (REAL) scalar 8 (I3)2s (I3)10s SIGHOT (REAL) scalar 360 (I1)7s (I1)15s 632r 632r 632r SIGLOT (REAL) scalar 356 (I1)7s (I1)15s SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s 626u 626= 627u 627= 628u 628= 681u 681= SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s 680u 680= SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s SYSMAX (REAL) scalar 148 (I3)8s (I3)11s 672u 686o SYSMIN (REAL) scalar 144 (I3)8s (I3)12s SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSUAV (REAL) scalar 120 (I3)6s (I3)12s SYSUSD (REAL) scalar 124 (I3)7s (I3)12s T (REAL) scalar 600s 623= 624u 624u 625u 635= 636o 640= 641u 647= 648u 654= 655o 664= 665o 675= 676u 678u 682u 710= 711o 720= 721u 723o 741= 742o Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 31 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- TXPROB (REAL) array 184 (I1)4s (I1)12s U (REAL) scalar 600s 643= 644o 663= 665o 694= 695o WCFDHH SUBROUTINE 571s WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s WORK (REAL) array 200 (I3)8s (I3)12s WRCAAV (REAL) scalar 52 (I3)4s (I3)13s 632u WRCASD (REAL) scalar 56 (I3)4s (I3)13s 632u WRCASK (REAL) scalar 60 (I3)4s (I3)13s WRCBAS (REAL) scalar 32 (I3)3s (I3)13s 660u 664u WRCHHB (REAL) scalar 64 (I3)4s (I3)13s 648= 651= 652u 658= 660u 663u 672u 676u 686o 694u 710u 741u WRCHOT (REAL) scalar 68 (I3)4s (I3)13s 632= 635u 651u 654u WRCLOW (REAL) scalar 72 (I3)5s (I3)13s WRCPAB (REAL) scalar 36 (I3)3s (I3)14s WRCSKW (REAL) scalar 48 (I3)4s (I3)14s WRCUAV (REAL) scalar 40 (I3)3s (I3)14s WRCUSD (REAL) scalar 44 (I3)4s (I3)14s WSKLAT (REAL) scalar 312 (I1)4s (I1)12s Z (REAL(kind 8)) scalar 601s 625= 626u 627u 628u 678= 680u 681u 682u 682= -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 30 622d 628s 40 620g 629s 50 641g 648s 60 633g 650s 70 660g 669s 90 674d 676g 683s 44 634g 639s 65 653g 658s 85 679d 682s 66 649g 652g 659s 100 672g 684g 689s 110 692g 707s 210 708g 733s 120 714g 720s 130 721g 725s 140 726g 731s 161 636f 637s 171 728f 729s 162 742f 743s 163 695f 697s 164 702f 703s 215 734g 739s 165 711f 712s 156 655f 656s 157 644f 645s 167 716f 717s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 32 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 197 686f 687s 159 665f 666s 169 723f 724s 16301 693g 696s 1601 736f 737s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 33 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFDLO Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 747 C 748 C 749 C 750 SUBROUTINE WCFDLO 751 # (SYSLOG,NSYS1,IER) 752 C 753 C + + + PURPOSE + + + 754 C DETECT LOW OUTLIERS 755 C REV 5/79 WK - MOVED WRC CHECK FOR TOO MANY BELOW BASE FROM FCA 756 C REV 7/20/79 WK TO MAKE ADDENDS TO SUMS DOUBLE PRECISION 757 C REV 1/81 WK - BULL 17-B REVISED OUTLIER TEST 758 C REV 4/28/81 WK - TO USE NUMBER OF PEAKS ABOVE GAGE BASE IN LOW-OUT 759 C TEST. USES HIST WT . 760 C 761 C + + + DUMMY ARGUMENTS + + + 762 INTEGER NSYS1,IER 763 REAL SYSLOG(NSYS1) 764 C 765 C + + + ARGUMENT DEFINITIONS + + + 766 C SYSLOG - systematic peak logarithms (input) 767 C NSYS1 - number of systematic peak logarithms 768 C IER - error return code 769 C 770 C + + + PARAMETERS + + + 771 INCLUDE 'pmxint.inc' 772 C 773 C + + + COMMON BLOCKS + + + 774 INCLUDE 'cwcf0.inc' 775 INCLUDE 'cwcf1.inc' 776 C 777 C + + + LOCAL VARIABLES + + + 778 INTEGER NBB,I,NBBMAX 779 REAL ERL,T,TEST 780 DOUBLE PRECISION Z 781 C 782 C + + + FUNCTIONS + + + 783 REAL OUTKGB 784 C 785 C + + + INTRINSICS + + + 786 INTRINSIC ALOG10, MAX0, FLOAT, INT 787 C 788 C + + + EXTERNALS + + + 789 EXTERNAL OUTKGB 790 C 791 C + + + END SPECIFICATIONS + + + 792 C 793 ERL = NSYS-NBGB 794 T = HISTPN - HISTWT*NBGB 795 ERL = ERL + HRECWO*(T-ERL) 796 T = OUTKGB(SIGLOT, INT(ERL+.5)) 797 WRCLOW=WRCAAV-WRCASD*T 798 TEST=WRCLOW Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 34 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 799 IF(QLWOUT.LE.0.)GOTO10 800 TEST=ALOG10(QLWOUT) 801 IER=MAX0(IER,1) 802 IF(MSL.LT.2)GOTO10 803 T=10.**WRCLOW 804 WRITE(MSG,191)QLWOUT,T 805 191 FORMAT( ' *WCF191I-USER LOW-OUTLIER CRITERION SUPERSEDES 17B.', 806 $ 4X,F8.1,1X,F8.1) 807 10 CONTINUE 808 C 809 IF(TEST.LT.WRCHHB)GOTO20 810 IER=3 811 IF(MSL.LE.0) RETURN 812 WRCHHB=10.**WRCHHB 813 TEST=10.**TEST 814 WRITE(MSG,193)TEST,WRCHHB 815 193 FORMAT(52H ***WCF193E-LOW-OUTLIER CRITERION EXCEEDS HIGH-HIST., 816 $ 2F11.1) 817 RETURN 818 C 819 20 NLWOUT=0 820 IF(TEST.GE.SYSMIN)GOTO30 821 IF(MSL.LT.2) GO TO 25 822 TEST=10.**TEST 823 WRITE(MSG,195)TEST 824 195 FORMAT(48H WCF195I-NO LOW OUTLIERS WERE DETECTED BELOW , 825 $ 10HCRITERION.,9X,F8.1) 826 25 IF(NBGB.GT.0) GO TO 65 827 RETURN 828 C 829 30 WRCBAS=TEST 830 DO 40 I=1,NSYS1 831 T=SYSLOG(I) 832 IF(T.GT.WRCBAS.OR.T.LE.SYSBAS)GOTO40 833 Z = T 834 SUMS(1)=SUMS(1)-Z 835 SUMS(2)=SUMS(2)-Z**2 836 SUMS(3)=SUMS(3)-Z**3 837 NLWOUT=NLWOUT+1 838 40 CONTINUE 839 IF(NLWOUT.GT.0)GOTO50 840 IER=3 841 IF(MSL.GT.0)WRITE(MSG,197)SYSMIN,WRCBAS,SYSBAS 842 197 FORMAT(/' ***WCF004*-INTERNAL PROGRAM LOGIC ERROR DLO-197. ', 843 $ 3F10.5/) 844 RETURN 845 C 846 50 IF(MSL.LT.2) GO TO 60 847 T=10.**WRCBAS 848 WRITE(MSG,198)NLWOUT,T 849 198 FORMAT(46H WCF198I-LOW OUTLIERS BELOW FLOOD BASE WERE, 850 $ 9H DROPPED., I11,1X,F8.1) Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 35 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 851 60 CONTINUE 852 C 853 C CHECK FOR TOO MANY BELOW-BASE PEAKS 854 65 CONTINUE 855 NBB = NBGB+NLWOUT 856 NBBMAX = 0.25*FLOAT(NSYS) 857 IF(NBB.LE.NBBMAX) GO TO 70 858 IER=MAX0(IER,2) 859 IF(MSL.LE.0) GO TO 70 860 T = 10.**WRCBAS 861 WRITE(MSG,199)NBB,T,NBBMAX 862 199 FORMAT( ' **WCF199W-NUMBER OF PEAKS BELOW FLOOD BASE', 863 $ ' EXCEEDS 17B SPEC.', 1X,I3,1X,F8.1,1X,I3) 864 70 CONTINUE 865 RETURN 866 END Bytes of stack required for this program unit: 40. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- BIGLOG (REAL) scalar 316 (I1)4s (I1)12s CLSIZE (REAL) scalar 40 (I1)3s (I1)11s CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s EPS1 (REAL) scalar 320 (I1)5s (I1)13s EPS2 (REAL) scalar 324 (I1)5s (I1)13s ERL (REAL) scalar 779s 793= 795u 795u 795= 796u GAGEB (REAL) scalar 12 (I3)2s (I3)10s GENSKU (REAL) scalar 0 (I3)2s (I3)10s GSKWGT (REAL) scalar 108 (I3)6s (I3)10s GSMAX (REAL) scalar 32 (I1)3s (I1)11s GSMIN (REAL) scalar 28 (I1)3s (I1)11s HISTPD (REAL) scalar 28 (I3)3s (I3)10s HISTPN (REAL) scalar 100 (I3)6s (I3)10s 794u HISTWT (REAL) scalar 104 (I3)6s (I3)10s 794u HRECWO (REAL) scalar 328 (I1)5s (I1)13s 795u HRECWS (REAL) scalar 336 (I1)5s (I1)13s HRECWX (REAL) scalar 332 (I1)5s (I1)13s I (INTEGER) scalar 778s 831u IER (INTEGER) scalar 751s 762s 801r 801r 801= 810= 840= 858r 858r 858= IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s INDXPT (INTEGER) array 56 (I1)4s (I1)9s IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s MSG (INTEGER) scalar 0 (I1)2s (I1)8s 804o 814o 823o 841o 848o 861o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 802u 811u 821u 841u 846u 859u Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 36 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u NBB (INTEGER) scalar 778s 855= 857u 861o NBBMAX (INTEGER) scalar 778s 856= 857u 861o NBGB (INTEGER) scalar 84 (I3)5s (I3)9s 793u 794u 826u 855u NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s NHIST (INTEGER) scalar 24 (I3)2s (I3)9s NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s NINDX (INTEGER) scalar 52 (I1)4s (I1)9s NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s 819= 837u 837= 839u 848o 855u NMISS (INTEGER) scalar 76 (I3)5s (I3)9s NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s NSYS (INTEGER) scalar 80 (I3)5s (I3)9s 793u 856r NSYS1 (INTEGER) scalar 751s 762s 763u 830u OUTKGB FUNCTION 783s 789u 796u QHIOUT (REAL) scalar 20 (I3)2s (I3)11s QLWOUT (REAL) scalar 16 (I3)2s (I3)11s 799u 800r 804o RMSDGS (REAL) scalar 340 (I1)6s (I1)14s RMSEGS (REAL) scalar 8 (I3)2s (I3)10s SIGHOT (REAL) scalar 360 (I1)7s (I1)15s SIGLOT (REAL) scalar 356 (I1)7s (I1)15s 796r 796r 796r SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s 834u 834= 835u 835= 836u 836= SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s 832u 841o SYSLOG (REAL) array 751s 763s 831u SYSMAX (REAL) scalar 148 (I3)8s (I3)11s SYSMIN (REAL) scalar 144 (I3)8s (I3)12s 820u 841o SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSUAV (REAL) scalar 120 (I3)6s (I3)12s SYSUSD (REAL) scalar 124 (I3)7s (I3)12s T (REAL) scalar 779s 794= 795u 796= 797u 803= 804o 831= 832u 832u 833u 847= 848o 860= 861o TEST (REAL) scalar 779s 798= 800= 809u 813u 813= 814o 820u 822u 822= 823o 829u TXPROB (REAL) array 184 (I1)4s (I1)12s WCFDLO SUBROUTINE 750s WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s WORK (REAL) array 200 (I3)8s (I3)12s WRCAAV (REAL) scalar 52 (I3)4s (I3)13s 797u WRCASD (REAL) scalar 56 (I3)4s (I3)13s 797u WRCASK (REAL) scalar 60 (I3)4s (I3)13s WRCBAS (REAL) scalar 32 (I3)3s (I3)13s 829= 832u 841o 847u 860u WRCHHB (REAL) scalar 64 (I3)4s (I3)13s 809u 812u 812= 814o WRCHOT (REAL) scalar 68 (I3)4s (I3)13s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 37 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- WRCLOW (REAL) scalar 72 (I3)5s (I3)13s 797= 798u 803u WRCPAB (REAL) scalar 36 (I3)3s (I3)14s WRCSKW (REAL) scalar 48 (I3)4s (I3)14s WRCUAV (REAL) scalar 40 (I3)3s (I3)14s WRCUSD (REAL) scalar 44 (I3)4s (I3)14s WSKLAT (REAL) scalar 312 (I1)4s (I1)12s Z (REAL(kind 8)) scalar 780s 833= 834u 835u 836u -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 10 799g 802g 807s 20 809g 819s 30 820g 829s 40 830d 832g 838s 50 839g 846s 60 846g 851s 70 857g 859g 864s 25 821g 826s 65 826g 854s 191 804f 805s 193 814f 815s 195 823f 824s 197 841f 842s 198 848f 849s 199 861f 862s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 38 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFEPP Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 867 C 868 C 869 C 870 SUBROUTINE WCFEPP 871 # (PPV, LABEL, NDIM) 872 C 873 C + + + PURPOSE + + + 874 C EMPIRICAL PROBABILITY PLOTTING POSITIONS, INCL WRC HIST ADJ AS APPRO 875 C 876 C + + + HISTORY + + + 877 C kmf - nov 9, 2000 - changed label from Hollerith to character 878 C 879 C + + + DUMMY ARGUMENTS + + + 880 C INTEGER LABEL,NDIM 881 INTEGER NDIM 882 REAL PPV(NDIM) 883 CHARACTER*4 LABEL 884 C 885 C + + + ARGUMENT DEFINITIONS + + + 886 C PPV - output vector for empirical plotting probabilities 887 C LABEL - identification label for printing message 888 C NDIM - dimension of PPV 889 C 890 C + + + PARAMETERS + + + 891 INCLUDE 'pmxint.inc' 892 C 893 C + + + COMMON BLOCKS + + + 894 INCLUDE 'cwcf0.inc' 895 INCLUDE 'cwcf1.inc' 896 C 897 C + + + LOCAL VARIABLES + + + 898 INTEGER NPK,NZ,I 899 REAL FAC,A 900 C 901 C + + + INTRINSICS + + + 902 INTRINSIC FLOAT,MIN0 903 C 904 C + + + END SPECIFICATIONS + + + 905 C 906 NPK=NSYS+NHISTN 907 NZ=NHIOUT+NHISTN 908 FAC=1./(HISTPN+1.-2.*WEIBA) 909 IF(NZ.LE.0)NZ=NPK 910 DO20I=1,NZ 911 20 PPV(I)=(FLOAT(I)-WEIBA)*FAC 912 IF(NZ.GE.NPK) GO TO 70 913 A=(HISTWT-1.)*(FLOAT(NZ)+.5)+WEIBA 914 NZ=NZ+1 915 DO40I=NZ,NPK 916 40 PPV(I) = (FLOAT(I)*HISTWT - A)*FAC 917 70 IF(MSL.LT.3)GOTO90 918 NZ=MIN0(10,NPK) Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 39 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 919 WRITE(MSG,203)LABEL,(PPV(I),I=1,NZ) 920 203 FORMAT(48H WCF203J-PLOTTING POSITIONS OF TOP TEN PEAKS. , 921 $ A6/(9X,10F7.4)) 922 90 CONTINUE 923 RETURN 924 END Bytes of stack required for this program unit: 32. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- A (REAL) scalar 899s 913= 916u BIGLOG (REAL) scalar 316 (I1)4s (I1)12s CLSIZE (REAL) scalar 40 (I1)3s (I1)11s CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s EPS1 (REAL) scalar 320 (I1)5s (I1)13s EPS2 (REAL) scalar 324 (I1)5s (I1)13s FAC (REAL) scalar 899s 908= 911u 916u GAGEB (REAL) scalar 12 (I3)2s (I3)10s GENSKU (REAL) scalar 0 (I3)2s (I3)10s GSKWGT (REAL) scalar 108 (I3)6s (I3)10s GSMAX (REAL) scalar 32 (I1)3s (I1)11s GSMIN (REAL) scalar 28 (I1)3s (I1)11s HISTPD (REAL) scalar 28 (I3)3s (I3)10s HISTPN (REAL) scalar 100 (I3)6s (I3)10s 908u HISTWT (REAL) scalar 104 (I3)6s (I3)10s 913u 916u HRECWO (REAL) scalar 328 (I1)5s (I1)13s HRECWS (REAL) scalar 336 (I1)5s (I1)13s HRECWX (REAL) scalar 332 (I1)5s (I1)13s I (INTEGER) scalar 898s 911u 911r 916u 916r 919= 919u 919u IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s INDXPT (INTEGER) array 56 (I1)4s (I1)9s IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s LABEL (CHARACTER) scalar 871s 883s 919o 919o MSG (INTEGER) scalar 0 (I1)2s (I1)8s 919o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 917u MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u NBGB (INTEGER) scalar 84 (I3)5s (I3)9s NDIM (INTEGER) scalar 871s 881s 882u NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s 907u NHIST (INTEGER) scalar 24 (I3)2s (I3)9s NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s 906u 907u NINDX (INTEGER) scalar 52 (I1)4s (I1)9s NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s NMISS (INTEGER) scalar 76 (I3)5s (I3)9s NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 40 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s NPK (INTEGER) scalar 898s 906= 909u 912u 915u 918r 918r NSYS (INTEGER) scalar 80 (I3)5s (I3)9s 906u NZ (INTEGER) scalar 898s 907= 909u 909= 910u 912u 913r 914u 914= 915u 918= 919o PPV (REAL) array 871s 882s 911= 916= 919o QHIOUT (REAL) scalar 20 (I3)2s (I3)11s QLWOUT (REAL) scalar 16 (I3)2s (I3)11s RMSDGS (REAL) scalar 340 (I1)6s (I1)14s RMSEGS (REAL) scalar 8 (I3)2s (I3)10s SIGHOT (REAL) scalar 360 (I1)7s (I1)15s SIGLOT (REAL) scalar 356 (I1)7s (I1)15s SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s SYSMAX (REAL) scalar 148 (I3)8s (I3)11s SYSMIN (REAL) scalar 144 (I3)8s (I3)12s SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSUAV (REAL) scalar 120 (I3)6s (I3)12s SYSUSD (REAL) scalar 124 (I3)7s (I3)12s TXPROB (REAL) array 184 (I1)4s (I1)12s WCFEPP SUBROUTINE 870s WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s 908u 911u 913u WORK (REAL) array 200 (I3)8s (I3)12s WRCAAV (REAL) scalar 52 (I3)4s (I3)13s WRCASD (REAL) scalar 56 (I3)4s (I3)13s WRCASK (REAL) scalar 60 (I3)4s (I3)13s WRCBAS (REAL) scalar 32 (I3)3s (I3)13s WRCHHB (REAL) scalar 64 (I3)4s (I3)13s WRCHOT (REAL) scalar 68 (I3)4s (I3)13s WRCLOW (REAL) scalar 72 (I3)5s (I3)13s WRCPAB (REAL) scalar 36 (I3)3s (I3)14s WRCSKW (REAL) scalar 48 (I3)4s (I3)14s WRCUAV (REAL) scalar 40 (I3)3s (I3)14s WRCUSD (REAL) scalar 44 (I3)4s (I3)14s WSKLAT (REAL) scalar 312 (I1)4s (I1)12s -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 20 910d 911s 40 915d 916s 70 912g 917s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 41 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 90 917g 922s 203 919f 920s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 42 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFFCA Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 925C 926C 927C 928 SUBROUTINE WCFFCA 929 # (FCQ,LABEL,IER) 930C 931C + + + PURPOSE + + + 932C FREQUENCY CURVE PEARSON TYPE III ORDINATES, 933C INCL COND PROB ADJST (WRC APPX 4) IF REQUIRED. 934C (WRC Bulletin-17 Flood Frequency Analysis) 935C ALSO INCLUDES WEIGHTED-SKEW CALC IF LABEL.EQ.'WRC'. (1/81- B-17-B.) 936C 937C REV 5/79 WK - MOVED WRC CHECK FOR TOO MANY BELOW BASE TO DLO 938C AFTER COND PROB ADJ. 939C REV 1/81 WK - BULL 17-B -- THREE-POINT FIT PEARSON TYPE III 940C -- COMPUTE WEIGHTED SKEW AFTER COND PROB ADJ (NOT BEFORE) WHEN 941C LABEL.EQ.'WRC'. (BUT NOT IN SYST-REC CALCS.) 942C 943C REV 11/5/81 WK - WCFFCA CHECK *ABS-VALUE* OF WRCSKW IN TABLE RANGE 944C 945C + + + HISTORY + + + 946C kmf - Nov 9, 2000, converted LABEL and LCWRC Hollerith 947C to character 948C 949C + + + DUMMY ARGUMENTS + + + 950Ckmf INTEGER LABEL,IER 951 INTEGER IER 952 REAL FCQ(*) 953 CHARACTER*4 LABEL 954C 955C + + + ARGUMENT DEFINITIONS + + + 956C FCQ - output vector of frequency curve ordinates 957C LABEL - identifier printed in error message and used to 958C identify systematic-records vs WRC calculations 959C IER - error return code 960C 961C + + + PARAMETERS + + + 962 INCLUDE 'pmxint.inc' 963C 964C + + + COMMON BLOCKS + + + 965 INCLUDE 'cwcf0.inc' 966 INCLUDE 'cwcf1.inc' 967C 968C + + + LOCAL VARIABLES + + + 969Ckmf INTEGER I, LCWRC 970 INTEGER I 971 REAL PBB, TAB(3), T 972 CHARACTER*4 LCWRC 973C 974C + + + FUNCTIONS + + + 975 REAL HARTRG, HARTK 976C Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 43 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 977C + + + INTRINSICS + + + 978 INTRINSIC ABS,AMAX1 979C 980C + + + EXTERNAL + + + 981 EXTERNAL HARTIV, HARTRG, HARTK, WCFCWS 982C 983Ckmf DATA LCWRC / 4H17B / 984 DATA LCWRC / '17B ' / 985C 986C + + + END SPECIFICATIONS + + + 987C 988 IF(ABS(WRCASK).LE.9.) GO TO 20 989 IER=3 990 IF(MSL.GT.0) WRITE(MSG,215) LABEL, WRCASK 991215 FORMAT(' ***WCF215E-SKEW OUT OF TABLE RANGE. ',1A4,3F9.3) 992 RETURN 993 20 CALL HARTIV(WRCASK, MXINT, TXPROB, 994 O WORK) 995 write(99,2000) wrcask 996 2000 format('WCFFCA: After 20 HARTIV, WRCASK ',F8.5,' WORK:') 997 write(99,2010) work 998 2010 format(8(1X,F8.5)) 999C 1000C CONDITIONAL PROB ADJUSTMT. 1001 PBB = 1. - WRCPAB 1002 IF(WRCPAB.GT.0.9999)GOTO40 1003Cprh use CPAMP1,CPAMP2,CPAMP3 as declared in include file cwcf0.inc 1004Cprh instead of CPAMP(3) array originally declared in local common here 1005Cprh IF(WRCPAB.LE.CPAMP(1)) GO TO 30 1006Cprh DO 25 I=1,3 1007Cprh 25 TAB(I) = WRCAAV + WRCASD*HARTK(1.-CPAMP(I)/WRCPAB, WORK) 1008 IF(WRCPAB.LE.CPAMP1) GO TO 30 1009 write(99,2015) WRCPAB,CPAMP1,CPAMP2,CPAMP3 1010 2015 format('WCFFCA: Before TAB calcs, WRCPAB, CPAMP1-3',4(1x,f8.5)) 1011 TAB(1) = WRCAAV + WRCASD*HARTK(1.-CPAMP1/WRCPAB, WRCASK) 1012 TAB(2) = WRCAAV + WRCASD*HARTK(1.-CPAMP2/WRCPAB, WRCASK) 1013 TAB(3) = WRCAAV + WRCASD*HARTK(1.-CPAMP3/WRCPAB, WRCASK) 1014 write(99,2020) TAB(1),TAB(2),TAB(3) 1015 2020 format('WCFFCA: TAB 1-3 ',3(1X,F8.5)) 1016 WRCSKW = HARTRG((TAB(3)-TAB(2))/(TAB(2)-TAB(1))) 1017 IF(WRCSKW.LE.9.) GO TO 27 1018 IER = 3 1019 IF(MSL.GT.0) WRITE(MSG,215)LABEL,WRCASK,WRCSKW 1020 RETURN 1021 27 CALL HARTIV(WRCSKW, MXINT, TXPROB, 1022 O WORK) 1023 write(99,2030) wrcskw 1024 2030 format('WCFFCA: After 27 HARTIV, WRCSKW ',F8.5,' WORK:') 1025 write(99,2010) work 1026Cprh see note above regarding replacing CPAMP(3) array 1027Cprh WRCUSD = (TAB(3)-TAB(2))/(HARTK(1.-CPAMP(3),WORK) 1028Cprh $ -HARTK(1.-CPAMP(2),WORK)) Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 44 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1029Cprh WRCUAV = TAB(1) - WRCUSD*HARTK(1.-CPAMP(1),WORK) 1030 WRCUSD = (TAB(3)-TAB(2))/(HARTK(1.-CPAMP3,WRCSKW) 1031 $ -HARTK(1.-CPAMP2,WRCSKW)) 1032 WRCUAV = TAB(1) - WRCUSD*HARTK(1.-CPAMP1,WRCSKW) 1033 write(99,2040) wrcusd,wrcuav 1034 2040 format('WCFFCA: WRCUSD ',F8.5,' WRCUAV ',F8.5) 1035 GO TO 45 1036C 1037 30 CONTINUE 1038 IER=3 1039 IF(MSL.LE.0) RETURN 1040 WRITE(MSG,213)LABEL,PBB 1041213 FORMAT(46H ***WCF213E-COND PROB ADJUST FAILED--EXCESSIVE, 1042 $ 1X,1A4, 18H PROB BELOW BASE. ,F8.4) 1043 RETURN 1044C 1045 40 CONTINUE 1046 WRCUAV=WRCAAV 1047 WRCUSD=WRCASD 1048 WRCSKW = WRCASK 1049C 1050C WEIGHTED SKEW 1051 45 CONTINUE 1052 IF(LABEL.NE.LCWRC) GO TO 50 1053 T = WRCSKW 1054 CALL WCFCWS(IER) 1055 IF(ABS(WRCSKW-T).LE.0.) GO TO 50 1056 IF(ABS(WRCSKW).LE.9) GO TO 48 1057 IER = 3 1058 IF(MSL.GT.0) WRITE(MSG,215) LABEL,WRCASK,WRCSKW,GENSKU 1059 RETURN 1060 48 CALL HARTIV(WRCSKW, MXINT, TXPROB, 1061 O WORK) 1062 write(99,2050) wrcskw 1063 2050 format('WCFFCA: After 48 HARTIV, WRCSKW ',F8.5,' WORK:') 1064 write(99,2010) work 1065C 1066 50 IF(MSL.GE.4)WRITE(MSG,217)LABEL,WRCPAB,WRCUAV,WRCUSD,WRCSKW, 1067 $ WRCAAV,WRCASD,WRCASK 1068217 FORMAT(' WCF217L-FREQUENCY CURVE PARAMS --', 1069 $ 1X,1A4,4F8.4/(50X,3F8.4)) 1070Cprh change to parameter MXINT, added CONTINUE to loop 1071Cprh DO 55 I=1,31 1072 DO 55 I=1,MXINT 1073 WORK(I) = AMAX1(WRCBAS, WRCUAV+WORK(I)*WRCUSD) 1074 55 CONTINUE 1075 DO 60 I=INDX1,INDX2 1076Cprh IF(TXPROB(I).GT.WRCPAB) GO TO 60 1077 IF(TXPROB(I).LE.WRCPAB) THEN 1078 FCQ(I) = WORK(I) 1079C *** removed so don't get break in plot, ie don't substitute 1080C *** the flood base under this condition Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 45 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1081C IF(ABS(TXPROB(I)-WRCPAB) .LE. 0.) FCQ(I) = WRCBAS 1082C *** 1083 END IF 1084 60 CONTINUE 1085C 1086Cprh replace RETURN with IF...END IF block 1087Cprh IF(MSL.LT.3) RETURN 1088 IF(MSL.GE.3) THEN 1089Cprh added continue to 70 loop 1090 DO 70 I=1,3 1091 TAB(I) = 10.**WORK(INDXPT(I)) 1092 70 CONTINUE 1093 WRITE(MSG,219)LABEL,TAB 1094219 FORMAT(33H WCF219J-FREQ CURVE ORDINATES ,1A4,4X, 1095 $ 38H 2-YR (.50) 10-YR (.10) 100-YR (.01) /39X,3F12.1) 1096 END IF 1097 1098 RETURN 1099 END Bytes of stack required for this program unit: 80. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- BIGLOG (REAL) scalar 316 (I1)4s (I1)12s CLSIZE (REAL) scalar 40 (I1)3s (I1)11s CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s 1008u 1009o 1011u 1032u CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s 1009o 1012u 1031u CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s 1009o 1013u 1030u EPS1 (REAL) scalar 320 (I1)5s (I1)13s EPS2 (REAL) scalar 324 (I1)5s (I1)13s FCQ (REAL) array 929s 952s 1078= GAGEB (REAL) scalar 12 (I3)2s (I3)10s GENSKU (REAL) scalar 0 (I3)2s (I3)10s 1058o GSKWGT (REAL) scalar 108 (I3)6s (I3)10s GSMAX (REAL) scalar 32 (I1)3s (I1)11s GSMIN (REAL) scalar 28 (I1)3s (I1)11s HARTIV SUBROUTINE 981u 993u 1021u 1060u HARTK FUNCTION 975s 981u 1011u 1012u 1013u 1030u 1031u 1032u HARTRG FUNCTION 975s 981u 1016u HISTPD (REAL) scalar 28 (I3)3s (I3)10s HISTPN (REAL) scalar 100 (I3)6s (I3)10s HISTWT (REAL) scalar 104 (I3)6s (I3)10s HRECWO (REAL) scalar 328 (I1)5s (I1)13s HRECWS (REAL) scalar 336 (I1)5s (I1)13s HRECWX (REAL) scalar 332 (I1)5s (I1)13s I (INTEGER) scalar 970s 1073u 1073u 1077u 1078u 1078u 1091u 1091u IER (INTEGER) scalar 929s 951s 989= 1018= 1038= 1054r 1054r 1054r 1057= IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 46 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s 1075u INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s 1075u INDXPT (INTEGER) array 56 (I1)4s (I1)9s 1091u IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s LABEL (CHARACTER) scalar 929s 953s 990o 990o 1019o 1019o 1040o 1040o 1052u 1058o 1058o 1066o 1066o 1093o 1093o 1052i LCWRC (CHARACTER) scalar 972s 984/ 1052u 1052i MSG (INTEGER) scalar 0 (I1)2s (I1)8s 990o 1019o 1040o 1058o 1066o 1093o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 990u 1019u 1039u 1058u 1066u 1088u MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u 993u 1021u 1060u 1072u NBGB (INTEGER) scalar 84 (I3)5s (I3)9s NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s NHIST (INTEGER) scalar 24 (I3)2s (I3)9s NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s NINDX (INTEGER) scalar 52 (I1)4s (I1)9s NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s NMISS (INTEGER) scalar 76 (I3)5s (I3)9s NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s NSYS (INTEGER) scalar 80 (I3)5s (I3)9s PBB (REAL) scalar 971s 1001= 1040o QHIOUT (REAL) scalar 20 (I3)2s (I3)11s QLWOUT (REAL) scalar 16 (I3)2s (I3)11s RMSDGS (REAL) scalar 340 (I1)6s (I1)14s RMSEGS (REAL) scalar 8 (I3)2s (I3)10s SIGHOT (REAL) scalar 360 (I1)7s (I1)15s SIGLOT (REAL) scalar 356 (I1)7s (I1)15s SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s SYSMAX (REAL) scalar 148 (I3)8s (I3)11s SYSMIN (REAL) scalar 144 (I3)8s (I3)12s SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSUAV (REAL) scalar 120 (I3)6s (I3)12s SYSUSD (REAL) scalar 124 (I3)7s (I3)12s T (REAL) scalar 971s 1053= 1055u TAB (REAL) array 971s 1011= 1012= 1013= 1014o 1014o 1014o 1016u 1016u 1016u 1016u 1030u 1030u 1032u 1091= 1093o TXPROB (REAL) array 184 (I1)4s (I1)12s 993r 993r 993r 1021r 1021r 1021r 1060r 1060r 1060r 1077u WCFCWS SUBROUTINE 981u 1054u WCFFCA SUBROUTINE 928s WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 47 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- WORK (REAL) array 200 (I3)8s (I3)12s 994r 993r 993r 997o 1022r 1021r 1021r 1025o 1061r 1060r 1060r 1064o 1073u 1073= 1078u 1091u WRCAAV (REAL) scalar 52 (I3)4s (I3)13s 1011u 1012u 1013u 1046u 1067o WRCASD (REAL) scalar 56 (I3)4s (I3)13s 1011u 1012u 1013u 1047u 1067o WRCASK (REAL) scalar 60 (I3)4s (I3)13s 988r 990o 993r 993r 993r 995o 1011r 1011r 1011r 1012r 1012r 1012r 1013r 1013r 1013r 1019o 1048u 1058o 1067o WRCBAS (REAL) scalar 32 (I3)3s (I3)13s 1073r 1073r WRCHHB (REAL) scalar 64 (I3)4s (I3)13s WRCHOT (REAL) scalar 68 (I3)4s (I3)13s WRCLOW (REAL) scalar 72 (I3)5s (I3)13s WRCPAB (REAL) scalar 36 (I3)3s (I3)14s 1001u 1002u 1008u 1009o 1011u 1012u 1013u 1066o 1077u WRCSKW (REAL) scalar 48 (I3)4s (I3)14s 1016= 1017u 1019o 1021r 1021r 1021r 1023o 1030r 1030r 1030r 1031r 1031r 1031r 1032r 1032r 1032r 1048= 1053u 1055u 1056r 1058o 1060r 1060r 1060r 1062o 1066o WRCUAV (REAL) scalar 40 (I3)3s (I3)14s 1032= 1033o 1046= 1066o 1073u WRCUSD (REAL) scalar 44 (I3)4s (I3)14s 1030= 1032u 1033o 1047= 1066o 1073u WSKLAT (REAL) scalar 312 (I1)4s (I1)12s -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 20 988g 993s 30 1008g 1037s 40 1002g 1045s 50 1052g 1055g 1066s 60 1075d 1084s 70 1090d 1092s 45 1035g 1051s 55 1072d 1074s 27 1017g 1021s 48 1056g 1060s 213 1040f 1041s 215 990f 991s 1019f 1058f 217 1066f 1068s 219 1093f 1094s 2000 995f 996s 2010 997f 998s 1025f 1064f 2020 1014f 1015s 2030 1023f 1024s 2040 1033f 1034s 2050 1062f 1063s 2015 1009f 1010s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 48 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. SUBROUTINE WCFFCX Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1100C 1101C 1102C 1103 SUBROUTINE WCFFCX 1104 # ( IER ) 1105C 1106C + + + PURPOSE + + + 1107C FREQ CURVE EXTENSIONS -- CONFIDENCE LIMS AND EXPECTED PROBS 1108C (WRC Bulletin-17 Flood Frequency Analysis) 1109C REV 10/79 WK - TO USE STUTX FOR T-QUANTILES IN EXPECTED PROB. 1110C 1111C + + + DUMMY ARGUMENTS + + + 1112 INTEGER IER 1113C 1114C + + + ARGUMENT DEFINITIONS + + + 1115C IER - error return code 1116C 1117C + + + PARAMETERS + + + 1118 INCLUDE 'pmxint.inc' 1119C 1120C + + + COMMON BLOCKS + + + 1121 INCLUDE 'cwcf0.inc' 1122 INCLUDE 'cwcf1.inc' 1123 INCLUDE 'cwcf2.inc' 1124C 1125C + + + LOCAL VARIABLES + + + 1126 INTEGER I,IIK,N233,NDF 1127 REAL A,ABT,C,EP,ERL,EXPFAC,RAD,SDA,SIGMA,T,TAB(3),TAB2(3),Z,ZC2 1128C 1129C + + + FUNCTIONS + + + 1130 REAL HARTK, GAUSCF, GAUSAB, STUTX 1131C 1132C + + + INTRINSICS + + + 1133 INTRINSIC MAX0,SQRT,ABS 1134C 1135C + + + EXTERNALS + + + 1136 EXTERNAL HARTK, GAUSCF, GAUSAB, STUTX 1137C 1138C + + + END SPECIFICATIONS + + + 1139C 1140 ERL = NSYS 1141 ERL = ERL + HRECWX*(HISTPN-ERL) 1142C 1143C EXPECTED PROBABILITY ADJUSTMENT, APPROX FOR SMALL SKEW 1144 IF(NOEPFC.EQ.1)GOTO70 1145 EXPFAC = SQRT(1. + 1./ERL) 1146 NDF = ERL - 1. 1147 N233 = 0 1148 DO60IIK=1,NINDX 1149 I=INDXPT(IIK) 1150 IF(TXPROB(I).GT.WRCPAB) GO TO 60 1151 T = STUTX(1.-TXPROB(I), NDF) Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 49 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1152 IF(ABS(T).LE.1E4) GO TO 30 1153 IER=3 1154 IF(MSL.GT.0) WRITE(MSG,237) NDF,TXPROB(I), T 1155237 FORMAT(/' ***WCF004*-INTERNAL PROGRAM LOGIC ERROR-STUTX. ', 1156 $ I6,2F10.5) 1157 RETURN 1158 30 EP=GAUSCF(T*EXPFAC) 1159 write(99,*) 'WCFFCX: B4 HARTK, I,EP,WRCSKW',I,EP,WRCSKW,' WORK:' 1160 write(99,2010) work 1161 2010 format(8(1X,F8.5)) 1162 Z = HARTK(EP, WRCSKW)* wrcusd + wrcuav 1163 write(99,*) 'WCFFCX: after HARTK, Z(',I,') is ',Z 1164 IF( Z.LT. BIGLOG) GO TO 50 1165 Z = BIGLOG 1166 N233 = N233 + 1 1167 IF(N233.GT.1) GO TO 50 1168 IER = MAX0(2, IER) 1169 EP=1.-EP 1170 IF(MSL.GT.0) WRITE(MSG,233) EP, TXPROB(I) 1171233 FORMAT(51H **WCF233W-EXPECTED PROB OUT OF RANGE AT TAB PROB., 1172 $ 2F10.5) 1173 50 IF(Z.LT.-BIGLOG) Z=-BIGLOG 1174 EPFC(I) = Z 1175 60 CONTINUE 1176 IF(MSL.LT.3)GOTO70 1177 DO65I=1,3 1178 65 TAB(I)=10.**EPFC(INDXPT(I)) 1179 WRITE(MSG,238)TAB 1180238 FORMAT( ' WCF238J-FREQ CURVE 17B-EXPECT-PROB.',3F12.1) 1181 70 CONTINUE 1182C 1183C 95 PCT CONFIDENCE LIMITS APPROX FOR SMALL SKEW 1184 IF(NOCLIM.EQ.1)GOTO90 1185 ZC2=GAUSAB(CLSIZE)**2 1186 A=1.-ZC2/(2.*(ERL-1.)) 1187 C=1.-A 1188 ABT = A*ZC2/ERL 1189 SDA=WRCUSD/A 1190 SIGMA=WRCUSD 1191 IF(ABS(SIGMA).LE.0.)SIGMA=1. 1192Cprh DO80IIK=1,NINDX 1193 DO 80 I = INDX1,INDX2 1194Cprh I=INDXPT(IIK) 1195 IF(TXPROB(I).GT.WRCPAB) GO TO 80 1196 Z = (WRCFC(I)-WRCUAV)/SIGMA 1197 RAD = SQRT(C*Z**2+ABT) 1198 CLIML(I) = WRCUAV + SDA*(Z-RAD) 1199 CLIMU(I) = WRCUAV + SDA*(Z+RAD) 1200 80 CONTINUE 1201 IF(MSL.LT.3)GOTO90 1202 DO85I=1,3 1203 TAB(I)= 10.**CLIML(INDXPT(I)) Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 50 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1204 85 TAB2(I)=10.**CLIMU(INDXPT(I)) 1205 Z=100.*CLSIZE 1206 WRITE(MSG,239) Z, TAB2, TAB 1207239 FORMAT( ' WCF239J-FREQ CURVE CONF LIMS B17B', F5.1,F9.1,2F12.1, 1208 $ /,39X, 3F12.1 ) 1209 90 CONTINUE 1210 RETURN 1211 END Bytes of stack required for this program unit: 184. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- A (REAL) scalar 1127s 1186= 1187u 1188u 1189u ABT (REAL) scalar 1127s 1188= 1197u BIGLOG (REAL) scalar 316 (I1)4s (I1)12s 1164u 1165u 1173u 1173u C (REAL) scalar 1127s 1187= 1197u CLIML (REAL) array 256 (I2)2s (I2)4s 1198= 1203u CLIMU (REAL) array 384 (I2)3s (I2)5s 1199= 1204u CLSIZE (REAL) scalar 40 (I1)3s (I1)11s 1185r 1185r 1185r 1205u CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s EP (REAL) scalar 1127s 1158= 1159o 1162r 1162r 1162r 1169u 1169= 1170o EPFC (REAL) array 128 (I2)2s (I2)4s 1174= 1178u EPS1 (REAL) scalar 320 (I1)5s (I1)13s EPS2 (REAL) scalar 324 (I1)5s (I1)13s ERL (REAL) scalar 1127s 1140= 1141u 1141u 1141= 1145u 1146u 1186u 1188u EXPFAC (REAL) scalar 1127s 1145= 1158u GAGEB (REAL) scalar 12 (I3)2s (I3)10s GAUSAB FUNCTION 1130s 1136u 1185u GAUSCF FUNCTION 1130s 1136u 1158u GENSKU (REAL) scalar 0 (I3)2s (I3)10s GSKWGT (REAL) scalar 108 (I3)6s (I3)10s GSMAX (REAL) scalar 32 (I1)3s (I1)11s GSMIN (REAL) scalar 28 (I1)3s (I1)11s HARTK FUNCTION 1130s 1136u 1162u HISTPD (REAL) scalar 28 (I3)3s (I3)10s HISTPN (REAL) scalar 100 (I3)6s (I3)10s 1141u HISTWT (REAL) scalar 104 (I3)6s (I3)10s HRECWO (REAL) scalar 328 (I1)5s (I1)13s HRECWS (REAL) scalar 336 (I1)5s (I1)13s HRECWX (REAL) scalar 332 (I1)5s (I1)13s 1141u I (INTEGER) scalar 1126s 1149= 1150u 1151u 1154u 1159o 1163o 1170u 1174u 1178u 1178u 1195u 1196u 1198u 1199u 1203u 1203u 1204u 1204u IER (INTEGER) scalar 1104s 1112s 1153= 1168r 1168r 1168= IGSOPT (INTEGER) scalar 4 (I3)2s (I3)9s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 51 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- IIK (INTEGER) scalar 1126s 1149u INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s 1193u INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s 1193u INDXPT (INTEGER) array 56 (I1)4s (I1)9s 1149u 1178u 1203u 1204u IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s MSG (INTEGER) scalar 0 (I1)2s (I1)8s 1154o 1170o 1179o 1206o MSL (INTEGER) scalar 4 (I1)2s (I1)8s 1154u 1170u 1176u 1201u MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u (I3)12u (I2)4u (I2)4u (I2)4u (I2)5u (I2)5u N233 (INTEGER) scalar 1126s 1147= 1166u 1166= 1167u NBGB (INTEGER) scalar 84 (I3)5s (I3)9s NDF (INTEGER) scalar 1126s 1146= 1151r 1151r 1151r 1154o NHIOUT (INTEGER) scalar 92 (I3)5s (I3)9s NHIST (INTEGER) scalar 24 (I3)2s (I3)9s NHISTN (INTEGER) scalar 96 (I3)5s (I3)9s NINDX (INTEGER) scalar 52 (I1)4s (I1)9s 1148u NLWOUT (INTEGER) scalar 88 (I3)5s (I3)9s NMISS (INTEGER) scalar 76 (I3)5s (I3)9s NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s 1184u NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s 1144u NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s NSYS (INTEGER) scalar 80 (I3)5s (I3)9s 1140u QHIOUT (REAL) scalar 20 (I3)2s (I3)11s QLWOUT (REAL) scalar 16 (I3)2s (I3)11s RAD (REAL) scalar 1127s 1197= 1198u 1199u RMSDGS (REAL) scalar 340 (I1)6s (I1)14s RMSEGS (REAL) scalar 8 (I3)2s (I3)10s SDA (REAL) scalar 1127s 1189= 1198u 1199u SIGHOT (REAL) scalar 360 (I1)7s (I1)15s SIGLOT (REAL) scalar 356 (I1)7s (I1)15s SIGMA (REAL) scalar 1127s 1190= 1191r 1191= 1196u STUTX FUNCTION 1130s 1136u 1151u SUMH (REAL(kind 8)) array 176 (I3)8s (I3)15s SUMS (REAL(kind 8)) array 152 (I3)8s (I3)15s SYSAAV (REAL) scalar 132 (I3)7s (I3)11s SYSASD (REAL) scalar 136 (I3)7s (I3)11s SYSASK (REAL) scalar 140 (I3)7s (I3)11s SYSBAS (REAL) scalar 112 (I3)6s (I3)11s SYSMAX (REAL) scalar 148 (I3)8s (I3)11s SYSMIN (REAL) scalar 144 (I3)8s (I3)12s SYSPAB (REAL) scalar 116 (I3)6s (I3)12s SYSRFC (REAL) array 512 (I2)3s (I2)5s SYSSKW (REAL) scalar 128 (I3)7s (I3)12s SYSUAV (REAL) scalar 120 (I3)6s (I3)12s SYSUSD (REAL) scalar 124 (I3)7s (I3)12s T (REAL) scalar 1127s 1151= 1152r 1154o 1158u TAB (REAL) array 1127s 1178= 1179o 1203= 1206o TAB2 (REAL) array 1127s 1204= 1206o TXPROB (REAL) array 184 (I1)4s (I1)12s 1150u 1151u 1154o 1170o 1195u WCFFCX SUBROUTINE 1103s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 52 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- WCXAUX (REAL) array 364 (I1)7s (I1)15s WEIBA (REAL) scalar 36 (I1)3s (I1)11s WORK (REAL) array 200 (I3)8s (I3)12s 1160o WRCAAV (REAL) scalar 52 (I3)4s (I3)13s WRCASD (REAL) scalar 56 (I3)4s (I3)13s WRCASK (REAL) scalar 60 (I3)4s (I3)13s WRCBAS (REAL) scalar 32 (I3)3s (I3)13s WRCFC (REAL) array 0 (I2)2s (I2)4s 1196u WRCHHB (REAL) scalar 64 (I3)4s (I3)13s WRCHOT (REAL) scalar 68 (I3)4s (I3)13s WRCLOW (REAL) scalar 72 (I3)5s (I3)13s WRCPAB (REAL) scalar 36 (I3)3s (I3)14s 1150u 1195u WRCSKW (REAL) scalar 48 (I3)4s (I3)14s 1159o 1162r 1162r 1162r WRCUAV (REAL) scalar 40 (I3)3s (I3)14s 1162u 1196u 1198u 1199u WRCUSD (REAL) scalar 44 (I3)4s (I3)14s 1162u 1189u 1190u WSKLAT (REAL) scalar 312 (I1)4s (I1)12s Z (REAL) scalar 1127s 1162= 1163o 1164u 1165= 1173u 1173= 1174u 1196= 1197u 1198u 1199u 1205= 1206o ZC2 (REAL) scalar 1127s 1185= 1186u 1188u -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 30 1152g 1158s 50 1164g 1167g 1173s 60 1148d 1150g 1175s 70 1144g 1176g 1181s 80 1193d 1195g 1200s 90 1184g 1201g 1209s 65 1177d 1178s 85 1202d 1204s 233 1170f 1171s 237 1154f 1155s 238 1179f 1180s 239 1206f 1207s 2010 1160f 1161s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 53 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. PROGRAM _LF_BLOCK_DATA Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 54 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. PROGRAM _LF_BLOCK_DATA Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1212C 1213C 1214C 1215 BLOCKDATA -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- BIGLOG (REAL) scalar 316 (I1)4s (I1)12s 1240/ CLSIZE (REAL) scalar 40 (I1)3s (I1)11s 1240/ CPAMP1 (REAL) scalar 344 (I1)6s (I1)14s 1261/ CPAMP2 (REAL) scalar 348 (I1)6s (I1)14s 1261/ CPAMP3 (REAL) scalar 352 (I1)6s (I1)14s 1261/ EPS1 (REAL) scalar 320 (I1)5s (I1)13s 1261/ EPS2 (REAL) scalar 324 (I1)5s (I1)13s 1261/ GSMAX (REAL) scalar 32 (I1)3s (I1)11s 1240/ GSMIN (REAL) scalar 28 (I1)3s (I1)11s 1240/ HRECWO (REAL) scalar 328 (I1)5s (I1)13s 1261/ HRECWS (REAL) scalar 336 (I1)5s (I1)13s 1261/ HRECWX (REAL) scalar 332 (I1)5s (I1)13s 1261/ INDX1 (INTEGER) scalar 44 (I1)4s (I1)9s 1240/ INDX2 (INTEGER) scalar 48 (I1)4s (I1)9s 1240/ INDXPT (INTEGER) array 56 (I1)4s (I1)9s 1246/ IWXMOD (INTEGER) scalar 372 (I1)7s (I1)10s 1265/ MSG (INTEGER) scalar 0 (I1)2s (I1)8s 1235/ MSL (INTEGER) scalar 4 (I1)2s (I1)8s 1235/ MXINT (INTEGER) scalar (I4)1s (I4)2s (I1)9u (I1)12u NINDX (INTEGER) scalar 52 (I1)4s (I1)9s 1246/ NOCLIM (INTEGER) scalar 24 (I1)2s (I1)8s 1235/ NOEPFC (INTEGER) scalar 20 (I1)2s (I1)8s 1235/ NOPPOS (INTEGER) scalar 8 (I1)2s (I1)8s 1235/ NOSYS (INTEGER) scalar 16 (I1)2s (I1)8s 1235/ NOTRAN (INTEGER) scalar 12 (I1)2s (I1)8s 1235/ RMSDGS (REAL) scalar 340 (I1)6s (I1)14s 1261/ SIGHOT (REAL) scalar 360 (I1)7s (I1)15s 1263/ SIGLOT (REAL) scalar 356 (I1)7s (I1)15s 1263/ TXPROB (REAL) array 184 (I1)4s (I1)12s 1255/ WCXAUX (REAL) array 364 (I1)7s (I1)15s 1265/ WEIBA (REAL) scalar 36 (I1)3s (I1)11s 1240/ WSKLAT (REAL) scalar 312 (I1)4s (I1)12s 1240/ _LF_BLOCK_DATA Procedure 1215s Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 55 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. FUNCTION WCFGSM Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1216C 1217C + + + PURPOSE + + + 1218C INITIALIZE COMMON WCFCM0 1219C REV 7/77 WK TO ENLARGE WCFCM0 FOR USE BY J407-X. 1220C REV 5/79 WK - SET IWXMOD TO 2 IF ANY WRC PARAMETERS ARE CHANGED. 1221C REV 11/03 PRH, AQUA TERRA Consultants, for batch version of PEAKFQ 1222C 1223C + + + PARAMETERS + + + 1224 INCLUDE 'pmxint.inc' 1225C 1226C + + + COMMON BLOCKS + + + 1227C WCFBDI is set here, but not used anywhere else 1228Cprh COMMON/ WCFBDI /WCFBDI 1229Cprh INTEGER WCFBDI(1) 1230C 1231 INCLUDE 'cwcf0.inc' 1232C 1233C + + + DATA INITIALIZATIONS + + + 1234Cprh DATA WCFBDI /0/ 1235 DATA MSG, MSL, NOPPOS, NOTRAN, NOSYS, NOEPFC, NOCLIM 1236 $ / 6, 2, 0, 0, 0, 0, 0 / 1237C 1238Cmeb: BIGLOG changed from 38 to 37 on Data General AViiON 1239Caml: and further changed to 29 for 5/94 compiler 1240 DATA GSMIN, GSMAX, WEIBA, CLSIZE, INDX1, INDX2, WSKLAT, BIGLOG 1241 $ /-.4001, .8001, 0., .95, 5, 29, 0., 29. / 1242C 1243c DATA NINDX/14/, INDXPT/16,21,26,5,6,10,11,12, 17,20,23,25,27,28, 1244c * 17*0 / 1245Cprh: expanded index array to include 1.5 yr (.6667) and 2.33 yr (.4292) floods 1246 DATA NINDX/15/, INDXPT/17,22,27,5,6,10,11,12,14,18,21,24,26,28,29, 1247 * 17*0 / 1248 1249c DATA TXPROB / 1250c $ 0.9999, 0.9995, 0.9990, 0.9980, 0.9950, 0.9900, 0.9800, 0.9750, 1251c $ 0.9600, 0.9500, 0.9000, 0.8000, 0.7000, 0.6000, 0.5704, 0.5000, 1252c $ 0.4296, 0.4000, 0.3000, 0.2000, 0.1000, 0.0500, 0.0400, 0.0250, 1253c $ 0.0200, 0.0100, 0.0050, 0.0020, 0.0010, 0.0005, 0.0001/ 1254Cprh: added 1.5 yr (.6667) flood to probability array 1255 DATA TXPROB / 1256 $ 0.9999, 0.9995, 0.9990, 0.9980, 0.9950, 0.9900, 0.9800, 0.9750, 1257 $ 0.9600, 0.9500, 0.9000, 0.8000, 0.7000, 0.6667, 0.6000, 0.5704, 1258 $ 0.5000, 0.4292, 0.4000, 0.3000, 0.2000, 0.1000, 0.0500, 0.0400, 1259 $ 0.0250, 0.0200, 0.0100, 0.0050, 0.0020, 0.0010, 0.0005, 0.0001/ 1260C 1261 DATA EPS1,EPS2, HRECWO,HRECWX,HRECWS, RMSDGS, CPAMP1,CPAMP2,CPAMP3 1262 $ / -0.4,+0.4, 1.0, 0.0, 1.0, 0.55, 0.50, 0.10, 0.01/ 1263 DATA SIGLOT, SIGHOT / 0.10, 0.10/ 1264C 1265 DATA WCXAUX / 2*0. /, IWXMOD / 0 / 1266C 1267C + + + END SPECIFICATIONS + + + Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 56 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1268C 1269 END 1270C 1271C 1272C 1273 REAL FUNCTION WCFGSM 1274 # ( FLATA, FLONGA ) 1275C 1276C + + + PURPOSE + + + 1277C WRC BULL. 17-A GENERALIZED SKEW MAP.C 1278C FOR LATITUDE AND LONGITUDE IN DEGREES AND FRACTIONS OF A DEGREE. 1279C (E.G., 45 DEG, 30 MIN = 45.5 DEG.) 1280C POINTS OUTSIDE U.S. ARE ASSIGNED LARGE NEGATIVE VALUE. 1281C 1282C DATA STATEMENT USES 2 CARDS PER MERIDIAN OF LONGITUDE, STARTING AT 1283C 67 DEG AND ENDING AT 125 DEG WEST LONGITUDE. SKEW VALUES ARE READ 1284C FROM SOUTH (25 DEG) TO NORTH (50 DEG) ALONG EACH MERIDIAN. 1285C SKEW VALUES TABULATED ARE IN UNITS OF 0.01 SKEW UNIT. POINTS 1286C OUTSIDE CONTINENTAL U.S. ARE SET TO -222. 1287C ALASKA AND HAWAII ARE NOT TABULATED BY LAT-LONG BUT ARE HANDLED 1288C CORRECTLY BY PROGRAMMED TESTS. ALSO PUERTO RICO. 1289C 1290C REV 8/6/79 WK TO IMPROVE THE REPRESENTATION OF THE GEN SKEW RIDGE 1291C AT GEN SKEW VALUE OF 0.2 RUNNING THROUGH EAST CALIFORNIA 1292C AND WEST OREGON AND WASHINGTON. 1293C 1294C REV 8/9/79 WK REVISED VARIOUS MAP SKEWS. 1295C 1296C REV 11/5/81 WK -- ADD PUERTO RICO GEN SKEW = 0.0 (12/17/86-WK) 1297C 1298C REV 4/86 WK TO CIRCUMVENT PRIME F77 REV19.4.4 COMPILER BUG(S) 1299C RELATING TO STATEMENT FUNCTIONS. 1300C 1301C 1302C + + + DUMMY ARGUMENTS + + + 1303 REAL FLATA, FLONGA 1304C 1305C + + + ARGUMENT DEFINITIONS + + + 1306C FLATA - Latitude North, degrees and fraction 1307C FLONGA - Longitude West, degrees and fraction 1308C 1309C + + + LOCAL VARIABLES + + + 1310 INTEGER*2 ISK(26,59) 1311 INTEGER I, J, LAT, LONG 1312 REAL FLAT, FLONG, GSK, S, WX, WY 1313 REAL FLATMN, FLATMX, FLONMN, FLONMX 1314 LOGICAL LLCELL, LLC 1315C 1316C + + + SAVE VARIABLES + + + 1317 INTEGER*2 ISK1A(26,5),ISK1B(26,4),ISK2A(26,5),ISK2B(26,4), 1318 1 ISK3A(26,5),ISK3B(26,4),ISK4A(26,5),ISK4B(26,4), 1319 2 ISK5A(26,5),ISK5B(26,4),ISK6A(26,5),ISK6B(26,4), Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 57 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1320 3 ISK7A(26,5) 1321 SAVE ISK1A,ISK1B,ISK2A,ISK2B, 1322 1 ISK3A,ISK3B,ISK4A,ISK4B, 1323 2 ISK5A,ISK5B,ISK6A,ISK6B, 1324 3 ISK7A 1325C 1326C + + + EQUIVALENCES + + + 1327 EQUIVALENCE (ISK1A(1,1),ISK(1, 1)),(ISK1B(1,1),ISK(1, 6)), 1328 1 (ISK2A(1,1),ISK(1,10)),(ISK2B(1,1),ISK(1,15)), 1329 2 (ISK3A(1,1),ISK(1,19)),(ISK3B(1,1),ISK(1,24)), 1330 3 (ISK4A(1,1),ISK(1,28)),(ISK4B(1,1),ISK(1,33)), 1331 4 (ISK5A(1,1),ISK(1,37)),(ISK5B(1,1),ISK(1,42)), 1332 5 (ISK6A(1,1),ISK(1,46)),(ISK6B(1,1),ISK(1,51)), 1333 6 (ISK7A(1,1),ISK(1,55)) 1334C 1335C + + + DATA INITIALIZATIONS + + + 1336 DATA ISK1A / 1337 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1338 $ -222,-222,-222,-222,-222,-222, 55, 38, 12, -8, -15,-222,-222, 1339 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1340 $ -222,-222,-222,-222,-222,-222, 50, 33, 10, -8, -15,-222,-222, 1341 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1342 $ -222,-222,-222, 70, 70, 70, 47, 34, 10, -8, -20,-222,-222, 1343 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1344 $ -222,-222,-222, 70, 70, 60, 44, 31, 5, -10, -20,-222,-222, 1345 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1346 $ -222,-222,-222, 70, 70, 57, 40, 28, 0, -10,-222,-222,-222/ 1347 DATA ISK1B / 1348 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1349 $ -222,-222, 70, 70, 70, 48, 35, 18, -2,-222,-222,-222,-222, 1350 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1351 $ -222,-222, 70, 70, 60, 42, 30, 8, -8,-222,-222,-222,-222, 1352 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1353 $ -222, 70, 70, 70, 48, 36, 20, 0, -12,-222,-222,-222,-222, 1354 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 45, 70, 70, 1355 $ 70, 70, 70, 56, 40, 28, 10, -5, -16,-222,-222,-222,-222/ 1356 DATA ISK2A / 1357 $ -222,-222,-222,-222,-222,-222,-222,-222,-222, 20, 40, 70, 70, 1358 $ 70, 70, 65, 44, 32, 18, 0, -10,-222,-222,-222,-222,-222, 1359 $ -222,-222,-222,-222,-222,-222,-222,-222,-222, 17, 38, 60, 70, 1360 $ 70, 70, 50, 33, 25, 5, -8, -17,-222,-222,-222,-222,-222, 1361 $ -222,-222,-222,-222,-222,-222,-222,-222, -3, 12, 32, 54, 68, 1362 $ 65, 52, 38, 28, 10, -5, -15,-222,-222,-222,-222,-222,-222, 1363 $ -222,-222,-222,-222,-222,-222,-222, -10, -10, 2, 30, 50, 60, 1364 $ 55, 42, 30, 12, -2, -13, -24,-222,-222,-222,-222,-222,-222, 1365 $ -10, -10, -10, -10, -10, -10, -10, -10, -10, -7, 25, 45, 51, 1366 $ 40, 28, 12, -1, -12, -24, -34,-222,-222,-222,-222,-222,-222/ 1367 DATA ISK2B / 1368 $ -3, -7, -10, -10, -10, -10, -10, -10, -10, -10, 15, 40, 40, 1369 $ 30, 15, 0, -10, -22, -33,-222,-222,-222,-222,-222,-222,-222, 1370 $ 0, -2, -5, -6, -10, -10, -10, -10, -10, -10, 10, 32, 29, 1371 $ 18, 0, -10, -20, -30, -40, -40,-222,-222,-222,-222,-222,-222, Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 58 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1372 $ -222, 0, 0, 0, -4, -7, -10, -10, -10, -10, 5, 18, 15, 1373 $ 0, -10, -20, -30, -40, -40, -40, -40, -33,-222,-222,-222,-222, 1374 $ -222,-222,-222,-222, 5, 0, -3, -8, -8, -3, 5, 8, 0, 1375 $ -9, -19, -30, -40, -40, -40, -40, -40, -25, -10,-222,-222,-222/ 1376 DATA ISK3A / 1377 $ -222,-222,-222,-222, 17, 13, 8, 3, 2, 2, 5, 2, -7, 1378 $ -15, -28, -40, -40, -40, -40, -40, -33, -20, -10,-222,-222,-222, 1379 $ -222,-222,-222,-222, 20, 20, 20, 17, 12, 7, 5, -2, -14, 1380 $ -24, -37, -40, -40, -40, -40, -40, -28, -10, -10,-222,-222,-222, 1381 $ -222,-222,-222,-222,-222, 20, 20, 20, 16, 8, 0, -12, -19, 1382 $ -32, -40, -40, -40, -40, -40, -38, -23, -10, -10, -10,-222,-222, 1383 $ -222,-222,-222,-222,-222, 20, 20, 20, 11, 3, -8, -17, -27, 1384 $ -38, -40, -40, -40, -40, -40, -35, -22, -10, -10, -10, -10, -10, 1385 $ -222,-222,-222,-222, 20, 20, 19, 12, 3, -5, -14, -24, -30, 1386 $ -38, -40, -40, -40, -40, -40, -36, -23, -10, -10, -10, -10, -10/ 1387 DATA ISK3B / 1388 $ -222,-222,-222,-222, 18, 13, 12, 2, -6, -15, -23, -30, -32, 1389 $ -38, -40, -40, -40, -40, -40, -38, -27, -14, -10, -10, -10, -10, 1390 $ -222,-222,-222,-222, 9, 4, 2, -9, -19, -24, -30, -30, -30, 1391 $ -35, -39, -40, -40, -40, -40, -40, -30, -20, -10, -10, -10, -10, 1392 $ -222,-222,-222,-222, 0, -7, -13, -20, -25, -30, -30, -30, -30, 1393 $ -34, -38, -40, -40, -40, -40, -40, -36, -26, -13, -10, -10, -10, 1394 $ -222,-222,-222,-222, -10, -15, -21, -27, -30, -30, -30, -30, -30, 1395 $ -31, -35, -38, -40, -40, -40, -40, -40, -32, -22, -15, -10, -10/ 1396 DATA ISK4A / 1397 $ -222,-222,-222,-222, -18, -24, -29, -30, -30, -30, -30, -30, -30, 1398 $ -30, -31, -34, -37, -40, -40, -40, -40, -39, -27, -22, -17, -10, 1399 $ -222,-222,-222, -23, -26, -30, -30, -30, -29, -29, -28, -28, -28, 1400 $ -28, -29, -30, -32, -38, -40, -40, -40, -40, -33, -28, -20, -15, 1401 $ -222,-222,-222, -29, -30, -29, -29, -28, -27, -26, -25, -25, -24, 1402 $ -24, -24, -25, -25, -30, -38, -40, -40, -40, -40, -34, -28, -20, 1403 $ -30, -30, -30, -29, -28, -28, -26, -25, -24, -23, -23, -22, -22, 1404 $ -20, -19, -17, -15, -14, -22, -40, -40, -40, -40, -40, -32, -27, 1405 $ -30, -29, -28, -28, -27, -25, -24, -23, -22, -21, -20, -19, -18, 1406 $ -15, -13, -9, -2, 30, 2, -25, -40, -40, -40, -40, -37, -32/ 1407 DATA ISK4B / 1408 $ -222, -27, -27, -26, -24, -23, -22, -21, -20, -19, -18, -16, -14, 1409 $ -13, -9, -5, 35, 60, 45, -5, -30, -40, -40, -40, -40, -35, 1410 $ -222, -25, -24, -23, -22, -21, -20, -18, -17, -16, -15, -13, -13, 1411 $ -10, -7, -2, 60, 60, 60, 10, -20, -40, -40, -40, -40, -40, 1412 $ -222,-222,-222, -21, -20, -18, -16, -15, -14, -14, -12, -11, -10, 1413 $ -10, -6, -1, 60, 60, 60, 23, -11, -33, -40, -40, -40, -40, 1414 $ -222,-222,-222,-222, -13, -13, -12, -12, -11, -11, -10, -10, -10, 1415 $ -10, -7, -2, 50, 60, 60, 32, -6, -27, -40, -40, -40, -40/ 1416 DATA ISK5A / 1417 $ -222,-222,-222,-222, 0, -4, -5, -7, -8, -9, -10, -10, -10, 1418 $ -10, -8, -4, 20, 60, 60, 30, -5, -24, -40, -40, -40, -40, 1419 $ -222,-222,-222,-222, -9, -5, 0, 0, -4, -6, -8, -10, -10, 1420 $ -10, -10, -8, -2, 40, 50, 12, -3, -18, -40, -40, -40, -40, 1421 $ -222,-222,-222,-222, -15, -13, -8, 0, 0, -3, -5, -8, -10, 1422 $ -10, -12, -13, -8, -1, 8, 0, -1, -10, -38, -40, -40, -40, 1423 $ -222,-222,-222,-222,-222, -20, -18, -12, -5, 0, -3, -6, -8, Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 59 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1424 $ -10, -23, -30, -25, -8, 0, 0, 0, -7, -27, -40, -40, -40, 1425 $ -222,-222,-222,-222,-222,-222, -20, -20, -14, -5, 0, -5, -8, 1426 $ -12, -30, -30, -30, -15, -4, 0, 10, 0, -17, -38, -40, -40/ 1427 DATA ISK5B / 1428 $ -222,-222,-222,-222,-222,-222, -20, -20, -20, -12, 0, -4, -7, 1429 $ -14, -30, -30, -30, -30, -8, 3, 12, 11, -5, -26, -40, -40, 1430 $ -222,-222,-222,-222,-222,-222, -20, -20, -20, -17, -7, 0, -7, 1431 $ -12, -26, -30, -30, -30, -20, -2, 15, 17, 8, -14, -30, -40, 1432 $ -222,-222,-222,-222,-222,-222, -20, -20, -20, -17, -10, 0, -4, 1433 $ -9, -22, -30, -30, -30, -30, -13, 12, 21, 17, 2, -18, -35, 1434 $ -222,-222,-222,-222,-222,-222, -13, -20, -20, -17, -11, -4, -3, 1435 $ -8, -18, -27, -30, -30, -30, -30, -5, 19, 24, 20, 0, -25/ 1436 DATA ISK6A / 1437 $ -222,-222,-222,-222,-222,-222, 8, -5, -11, -13, -10, -4, 0, 1438 $ -4, -13, -20, -28, -30, -30, -30, -25, 5, 28, 38, 28, -10, 1439 $ -222,-222,-222,-222,-222,-222, 18, 9, -3, -8, -7, -3, 0, 1440 $ 0, -8, -14, -23, -28, -30, -30, -30, -20, 18, 60, 60, 20, 1441 $ -222,-222,-222,-222,-222,-222, 0, 20, 8, -1, -3, 0, 0, 1442 $ 0, 0, -9, -17, -23, -29, -30, -30, -30, -20, 30, 60, 60, 1443 $ -222,-222,-222,-222,-222,-222,-222, 3, 20, 8, 0, 0, 0, 1444 $ 0, 0, -2, -10, -17, -23, -30, -30, -30, -30, -20, 30, 60, 1445 $ -222,-222,-222,-222,-222,-222,-222, -28, 5, 20, 10, 0, 0, 1446 $ 0, 0, 0, -3, -11, -16, -23, -28, -30, -30, -30, -10, 45/ 1447 DATA ISK6B / 1448 $ -222,-222,-222,-222,-222,-222,-222, -30, -22, 10, 20, 8, 0, 1449 $ 0, 0, 0, 0, -3, -8, -15, -20, -28, -30, -30, -30, -10, 1450 $ -222,-222,-222,-222,-222,-222,-222, -30, -30, -17, 7, 20, 10, 1451 $ 2, 0, 0, 0, 0, 0, -6, -12, -17, -25, -30, -30, -30, 1452 $ -222,-222,-222,-222,-222,-222,-222,-222, -30, -30, -15, 5, 20, 1453 $ 12, 4, 0, 0, 0, 0, 0, -2, -5, -12, -18, -22, -30, 1454 $ -222,-222,-222,-222,-222,-222,-222,-222,-222, -30, -30, -17, 2, 1455 $ 20, 20, 8, 1, 0, 0, 0, 0, 0, 0, -3, -7, -10/ 1456 DATA ISK7A / 1457 $ -222,-222,-222,-222,-222,-222,-222,-222,-222, -30, -30, -30, -23, 1458 $ -6, 10, 20, 17, 8, 0, 0, 0, 0, 0, 0, 0, 0, 1459 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222, -30, -30, -30, 1460 $ -30, -15, 0, 15, 20, 15, 10, 1, 0, 0, 0, 0, 0, 1461 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, -30, -30, 1462 $ -30, -30, -30, -18, -10, 15, 20, 20, 15, 3, 0, 0, 0, 1463 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1464 $ -30, -30, -30, -30, -30, -30, -18, -5, 15, 20, 20, 5, 0, 1465 $ -222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222,-222, 1466 $ -222,-222, -30, -30, -30, -30, -30, -30, -30, -20, 0, 0, 20/ 1467C 1468C + + + END SPECIFICATIONS + + + 1469C 1470 LLCELL(FLATMN,FLATMX,FLONMN,FLONMX) = FLONG.GE.FLONMN 1471 $ .AND. FLONG.LT.FLONMX .AND. FLAT.GE.FLATMN .AND. FLAT.LT.FLATMX 1472C 1473Caml exponent changed from 37 to 29 for 5/94 compiler 1474 WCFGSM = -1E29 1475 FLAT = FLATA Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 60 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- 1476 FLONG = FLONGA 1477C 1478 LLC = LLCELL(25., 50.01, 67., 125.01 ) 1479 IF(LLC) GO TO 100 1480C 1481C CHECK ALASKA, HAWAII, PUERTO RICO 1482C 1483 LLC = LLCELL(54., 61.,129., 140.) 1484 IF(LLC) WCFGSM = 0.33 1485C 1486 LLC = LLCELL(58., 72.,140., 170.) 1487 IF(LLC) WCFGSM = 0.70 1488C 1489 LLC = LLCELL(50., 58.,150., 360.) 1490 IF(LLC) WCFGSM = 0.70 1491C 1492 LLC = LLCELL(18., 30.,154., 162.) 1493 IF(LLC) WCFGSM = -0.05 1494C 1495 LLC = LLCELL(17., 19., 63., 68.) 1496 IF(LLC) WCFGSM = 0.0 1497C 1498 RETURN 1499C 1500C CONTINENTAL U.S. SKEW MAP 1501 100 CONTINUE 1502 LAT = FLAT 1503 LONG = FLONG 1504 WY = FLAT - LAT 1505 WX = FLONG - LONG 1506 GSK = 0. 1507 DO 150 I =1, 2 1508 WY = 1. - WY 1509 DO 150 J = 1, 2 1510 WX = 1. - WX 1511 S = ISK( LAT-24+I-1, LONG-66+J-1 ) 1512 IF(S.GT.-220.) GO TO 150 1513 IF(WX.LT.0.05) WX = 0. 1514 IF( WY.LT.0.05) WY = 0. 1515Caml 37 changed to 29 fro 5/94 compiler 1516 S = -1E29 1517 150 GSK = GSK + 0.01*S*WX*WY 1518 IF(GSK.GT.-1E29) WCFGSM = GSK 1519 RETURN 1520 END Bytes of stack required for this program unit: 32. -------------------------------------------------- Symbol Cross Reference --------------------------------------------------- Name (Type) Class Offset Context (d-DO, =-Assignment, /-Initialization, i-Input, o-Output, p-Pointer, r-Argument, s-Specification, u-Usage) ------------------------------------------------------------------------------------------------------------------------------- Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 61 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- FLAT (REAL) scalar 1312s 1471u 1471u 1475= 1502u 1504u FLATA (REAL) scalar 1274s 1303s 1475u FLONG (REAL) scalar 1312s 1470u 1471u 1476= 1503u 1505u FLONGA (REAL) scalar 1274s 1303s 1476u GSK (REAL) scalar 1312s 1506= 1517u 1517= 1518u 1518u I (INTEGER) scalar 1311s 1511u ISK (INTEGER(kind 2)) array 0 1310s 1327u 1327u 1328u 1328u 1329u 1329u 1330u 1330u 1331u 1331u 1332u 1332u 1333u 1511u ISK1A (INTEGER(kind 2)) array 0 1317s 1321s 1327u 1336/ ISK1B (INTEGER(kind 2)) array 260 1317s 1321s 1327u 1347/ ISK2A (INTEGER(kind 2)) array 468 1317s 1321s 1328u 1356/ ISK2B (INTEGER(kind 2)) array 728 1317s 1321s 1328u 1367/ ISK3A (INTEGER(kind 2)) array 936 1318s 1322s 1329u 1376/ ISK3B (INTEGER(kind 2)) array 1196 1318s 1322s 1329u 1387/ ISK4A (INTEGER(kind 2)) array 1404 1318s 1322s 1330u 1396/ ISK4B (INTEGER(kind 2)) array 1664 1318s 1322s 1330u 1407/ ISK5A (INTEGER(kind 2)) array 1872 1319s 1323s 1331u 1416/ ISK5B (INTEGER(kind 2)) array 2132 1319s 1323s 1331u 1427/ ISK6A (INTEGER(kind 2)) array 2340 1319s 1323s 1332u 1436/ ISK6B (INTEGER(kind 2)) array 2600 1319s 1323s 1332u 1447/ ISK7A (INTEGER(kind 2)) array 2808 1320s 1324s 1333u 1456/ J (INTEGER) scalar 1311s 1511u LAT (INTEGER) scalar 1311s 1502= 1504u 1511u LLC (LOGICAL) scalar 1314s 1478= 1479u 1483= 1484u 1486= 1487u 1489= 1490u 1492= 1493u 1495= 1496u LLCELL FUNCTION 1478u 1483u 1486u 1489u 1492u 1495u LLCELL (LOGICAL) scalar 1314s LONG (INTEGER) scalar 1311s 1503= 1505u 1511u S (REAL) scalar 1312s 1511= 1512u 1516= 1517u WCFGSM (REAL) scalar 1273s 1474= 1484= 1487= 1490= 1493= 1496= 1518= WCFGSM FUNCTION 1273s WX (REAL) scalar 1312s 1505= 1510u 1510= 1513u 1513= 1517u WY (REAL) scalar 1312s 1504= 1508u 1508= 1514u 1514= 1517u Lahey Fortran 90 Compiler Release 4.50b Tue Dec 07 16:12:17 2004 Page: 62 Copyright (C) 1994-1998 Lahey Computer Systems. All rights reserved. Compiling Options: -chk -ndal -ndbl -nf90 -fix -ng -nhed -nin -inln -lst -pca -stchk -nsyn -trap -nvax -w -nwo -xref ------------------------------------------------ Source Listing ------------------------------------------------------------- -------------------------------------------------- Label Cross Reference ---------------------------------------------------- Label Context (d-DO, =-ASSIGN, f-FORMAT, g-GOTO, i-IF statement, s-Specification, r-Argument) ------------------------------------------------------------------------------------------------------------------------------- 100 1479g 1501s 150 1507d 1509d 1512g 1517s