SUBROUTINE W3FI71 (IGRID, IGDS, IERR) C$$$ SUBPROGRAM DOCUMENTATION BLOCK C . . . . C SUBPROGRAM: W3FI71 MAKE ARRAY USED BY GRIB PACKER FOR GDS C PRGMMR: R.E.JONES ORG: W/NMC42 DATE: 93-03-26 C C ABSTRACT: W3FI71 MAKES A 18 OR 91 WORD INTEGER ARRAY USED BY W3FI72 C GRIB PACKER TO MAKE THE GRID DESCRIPTION SECTION (GDS) - SECTION 2. C C PROGRAM HISTORY LOG: C 92-02-21 R.E.JONES C 92-07-01 M. FARLEY ADDED REMARKS FOR 'IGDS' ARRAY ELEMENTS. C ADDED LAMBERT CONFORMAL GRIDS AND ENLARGED C IDGS ARRAY FROM 14 TO 18 WORDS. C 92-10-03 R.E.JONES ADDED CORRECTIONS TO AWIPS GRIB TABLES C 92-10-16 R.E.JONES ADD GAUSSIAN GRID 126 TO TABLES C 92-10-18 R.E.JONES CORRECTIONS TO LAMBERT CONFORMAL TABLES C AND OTHER TABLES C 92-10-19 R.E.JONES ADD GAUSSIAN GRID 98 TO TABLES C 93-01-25 R.E.JONES ADD ON84 GRIDS 87, 106, 107 TO TABLES C 93-03-10 R.E.JONES ADD ON84 GRIDS 1, 55, 56 TO TABLES C 93-03-26 R.E.JONES ADD GRIB GRIDS 2, 3 TO TABLES C 93-03-29 R.E.JONES ADD SAVE STATEMENT C 93-06-15 R.E.JONES ADD GRIB GRIDS 37 TO 44 TO TABLES C 93-09-29 R.E.JONES GAUSSIAN GRID DOCUMENT NOT CORRECT, C W3FI74 WILL BE CHANGED TO AGREE WITH C IT. GAUSSIAN GRID 98 TABLE HAS WRONG C VALUE. C 93-10-12 R.E.JONES CHANGES FOR ON388 REV. OCT 8,1993 FOR C GRID 204, 208. C 93-10-13 R.E.JONES CORRECTION FOR GRIDS 37-44, BYTES 7-8, C 24-25 SET TO ALL BITS 1 FOR MISSING. C 93-11-23 R.E.JONES ADD GRIDS 90-93 FOR ETA MODEL C ADD GRID 4 FOR 720*361 .5 DEG. GRID C 94-04-12 R.E.JONES CORRECTION FOR GRID 28 C 94-06-01 R.E.JONES ADD GRID 45, 288*145 1.25 DEG. GRID C 94-06-22 R.E.JONES ADD GRIDS 94, 95 FOR ETA MODEL C 95-04-11 R.E.JONES ADD GRIDS 96, 97 FOR ETA MODEL C 95-05-19 R.E.JONES ADD FROM 20 KM ETA MODEL AWIPS GRID 215 C 95-10-19 R.E.JONES ADD FROM 20 KM ETA MODEL ALASKA GRID 216 C 98-06-30 EBISUZAKI LINUX PORT C C USAGE: CALL W3FI71 (IGRID, IGDS, IERR) C INPUT ARGUMENT LIST: C IGRID - GRIB GRID NUMBER, OR OFFICE NOTE 84 GRID NUMBER C C OUTPUT ARGUMENT LIST: C IGDS - 18 OR 91 WORD INTEGER ARRAY WITH INFORMATION C TO MAKE A GRIB GRID DESCRIPTION SECTION. C IERR - 0 CORRECT EXIT C 1 GRID TYPE IN IGRID IS NOT IN TABLE C C REMARKS: C 1) OFFICE NOTE GRID TYPE 26 IS 6 IN GRIB, 26 IS AN C INTERNATIONAL EXCHANGE GRID. C C 2) VALUES RETURNED IN 18 OR 91 WORD INTEGER ARRAY IGDS C VARY DEPENDING ON GRID REPRESENTATION TYPE. C C LAT/LON GRID: C IGDS( 1) = NUMBER OF VERTICAL COORDINATES C IGDS( 2) = PV, PL OR 255 C IGDS( 3) = DATA REPRESENTATION TYPE (CODE TABLE 6) C IGDS( 4) = NO. OF POINTS ALONG A LATITUDE C IGDS( 5) = NO. OF POINTS ALONG A LONGITUDE MERIDIAN C IGDS( 6) = LATITUDE OF ORIGIN (SOUTH - IVE) C IGDS( 7) = LONGITUDE OF ORIGIN (WEST -IVE) C IGDS( 8) = RESOLUTION FLAG (CODE TABLE 7) C IGDS( 9) = LATITUDE OF EXTREME POINT (SOUTH - IVE) C IGDS(10) = LONGITUDE OF EXTREME POINT (WEST - IVE) C IGDS(11) = LATITUDE INCREMENT C IGDS(12) = LONGITUDE INCREMENT C IGDS(13) = SCANNING MODE FLAGS (CODE TABLE 8) C IGDS(14) = ... THROUGH ... C IGDS(18) = ... NOT USED FOR THIS GRID C IGDS(19) - IGDS(91) FOR GRIDS 37-44, NUMBER OF POINTS C IN EACH OF 73 ROWS. C C GAUSSIAN GRID: C IGDS( 1) = ... THROUGH ... C IGDS(10) = ... SAME AS LAT/LON GRID C IGDS(11) = NUMBER OF LATITUDE LINES BETWEEN A POLE C AND THE EQUATOR C IGDS(12) = LONGITUDE INCREMENT C IGDS(13) = SCANNING MODE FLAGS (CODE TABLE 8) C IGDS(14) = ... THROUGH ... C IGDS(18) = ... NOT USED FOR THIS GRID C C SPHERICAL HARMONICS: C IGDS( 1) = NUMBER OF VERTICAL COORDINATES C IGDS( 2) = PV, PL OR 255 C IGDS( 3) = DATA REPRESENTATION TYPE (CODE TABLE 6) C IGDS( 4) = J - PENTAGONAL RESOLUTION PARAMETER C IGDS( 5) = K - PENTAGONAL RESOLUTION PARAMETER C IGDS( 6) = M - PENTAGONAL RESOLUTION PARAMETER C IGDS( 7) = REPRESENTATION TYPE (CODE TABLE 9) C IGDS( 8) = REPRESENTATION MODE (CODE TABLE 10) C IGDS( 9) = ... THROUGH ... C IGDS(18) = ... NOT USED FOR THIS GRID C C POLAR STEREOGRAPHIC: C IGDS( 1) = NUMBER OF VERTICAL COORDINATES C IGDS( 2) = PV, PL OR 255 C IGDS( 3) = DATA REPRESENTATION TYPE (CODE TABLE 6) C IGDS( 4) = NO. OF POINTS ALONG X-AXIS C IGDS( 5) = NO. OF POINTS ALONG Y-AXIS C IGDS( 6) = LATITUDE OF ORIGIN (SOUTH -IVE) C IGDS( 7) = LONGITUTE OF ORIGIN (WEST -IVE) C IGDS( 8) = RESOLUTION FLAG (CODE TABLE 7) C IGDS( 9) = LONGITUDE OF MERIDIAN PARALLEL TO Y-AXIS C IGDS(10) = X-DIRECTION GRID LENGTH (INCREMENT) C IGDS(11) = Y-DIRECTION GRID LENGTH (INCREMENT) C IGDS(12) = PROJECTION CENTER FLAG (0=NORTH POLE ON PLANE, C 1=SOUTH POLE ON PLANE, C IGDS(13) = SCANNING MODE FLAGS (CODE TABLE 8) C IGDS(14) = ... THROUGH ... C IGDS(18) = .. NOT USED FOR THIS GRID C C MERCATOR: C IGDS( 1) = ... THROUGH ... C IGDS(12) = ... SAME AS LAT/LON GRID C IGDS(13) = LATITUDE AT WHICH PROJECTION CYLINDER C INTERSECTS EARTH C IGDS(14) = SCANNING MODE FLAGS C IGDS(15) = ... THROUGH ... C IGDS(18) = .. NOT USED FOR THIS GRID C C LAMBERT CONFORMAL: C IGDS( 1) = NUMBER OF VERTICAL COORDINATES C IGDS( 2) = PV, PL OR 255 C IGDS( 3) = DATA REPRESENTATION TYPE (CODE TABLE 6) C IGDS( 4) = NO. OF POINTS ALONG X-AXIS C IGDS( 5) = NO. OF POINTS ALONG Y-AXIS C IGDS( 6) = LATITUDE OF ORIGIN (SOUTH -IVE) C IGDS( 7) = LONGITUTE OF ORIGIN (WEST -IVE) C IGDS( 8) = RESOLUTION FLAG (CODE TABLE 7) C IGDS( 9) = LONGITUDE OF MERIDIAN PARALLEL TO Y-AXIS C IGDS(10) = X-DIRECTION GRID LENGTH (INCREMENT) C IGDS(11) = Y-DIRECTION GRID LENGTH (INCREMENT) C IGDS(12) = PROJECTION CENTER FLAG (0=NORTH POLE ON PLANE, C 1=SOUTH POLE ON PLANE, C IGDS(13) = SCANNING MODE FLAGS (CODE TABLE 8) C IGDS(14) = NOT USED C IGDS(15) = FIRST LATITUDE FROM THE POLE AT WHICH THE C SECANT CONE CUTS THE SPERICAL EARTH C IGDS(16) = SECOND LATITUDE ... C IGDS(17) = LATITUDE OF SOUTH POLE (MILLIDEGREES) C IGDS(18) = LONGITUDE OF SOUTH POLE (MILLIDEGREES) C C ARAKAWA SEMI-STAGGERED E-GRID ON ROTATED LAT/LON GRID C IGDS( 1) = NUMBER OF VERTICAL COORDINATES C IGDS( 2) = PV, PL OR 255 C IGDS( 3) = DATA REPRESENTATION TYPE (CODE TABLE 6) [201] C IGDS( 4) = NI - TOTAL NUMBER OF ACTUAL DATA POINTS C INCLUDED ON GRID C IGDS( 5) = NJ - DUMMY SECOND DIMENSION; SET=1 C IGDS( 6) = LA1 - LATITUDE OF FIRST GRID POINT C IGDS( 7) = LO1 - LONGITUDE OF FIRST GRID POINT C IGDS( 8) = RESOLUTION AND COMPONENT FLAG (CODE TABLE 7) C IGDS( 9) = LA2 - NUMBER OF MASS POINTS ALONG C SOUTHERNMOST ROW OF GRID C IGDS(10) = LO2 - NUMBER OF ROWS IN EACH COLUMN C IGDS(11) = DI - LONGITUDINAL DIRECTION INCREMENT C IGDS(12) = DJ - LATITUDINAL DIRECTION INCREMENT C IGDS(13) = SCANNING MODE FLAGS (CODE TABLE 8) C IGDS(14) = ... THROUGH ... C IGDS(18) = ... NOT USED FOR THIS GRID (SET TO ZERO) C C ARAKAWA FILLED E-GRID ON ROTATED LAT/LON GRID C IGDS( 1) = NUMBER OF VERTICAL COORDINATES C IGDS( 2) = PV, PL OR 255 C IGDS( 3) = DATA REPRESENTATION TYPE (CODE TABLE 6) [202] C IGDS( 4) = NI - TOTAL NUMBER OF ACTUAL DATA POINTS C INCLUDED ON GRID C IGDS( 5) = NJ - DUMMY SECOND DIMENTION; SET=1 C IGDS( 6) = LA1 - LATITUDE LATITUDE OF FIRST GRID POINT C IGDS( 7) = LO1 - LONGITUDE OF FIRST GRID POINT C IGDS( 8) = RESOLUTION AND COMPONENT FLAG (CODE TABLE 7) C IGDS( 9) = LA2 - NUMBER OF (ZONAL) POINTS IN EACH ROW C IGDS(10) = LO2 - NUMBER OF (MERIDIONAL) POINTS IN EACH C COLUMN C IGDS(11) = DI - LONGITUDINAL DIRECTION INCREMENT C IGDS(12) = DJ - LATITUDINAL DIRECTION INCREMENT C IGDS(13) = SCANNING MODE FLAGS (CODE TABLE 8) C IGDS(14) = ... THROUGH ... C IGDS(18) = ... NOT USED FOR THIS GRID CC C ATTRIBUTES: C LANGUAGE: SiliconGraphics 3.5 FORTRAN 77 C MACHINE: SiliconGraphics IRIS-4D/25, 35, INDIGO, Indy C LANGUAGE: IBM VS FORTRAN, CRAY CFT77 FORTRAN C MACHINE: HDS, CRAY C916-128, Y-MP8/864, CRAY Y-MP EL92/256 C C$$$ C INTEGER IGRID INTEGER IGDS (*) INTEGER GRD1 (18) INTEGER GRD2 (18) INTEGER GRD3 (18) INTEGER GRD4 (18) INTEGER GRD5 (18) INTEGER GRD6 (18) INTEGER GRD21 (18) INTEGER GRD22 (18) INTEGER GRD23 (18) INTEGER GRD24 (18) INTEGER GRD25 (18) INTEGER GRD26 (18) INTEGER GRD27 (18) INTEGER GRD28 (18) INTEGER GRD29 (18) INTEGER GRD30 (18) INTEGER GRD33 (18) INTEGER GRD34 (18) INTEGER GRD37 (91) INTEGER GRD38 (91) INTEGER GRD39 (91) INTEGER GRD40 (91) INTEGER GRD41 (91) INTEGER GRD42 (91) INTEGER GRD43 (91) INTEGER GRD44 (91) INTEGER GRD45 (18) C INTEGER GRD50 (18) INTEGER GRD55 (18) INTEGER GRD56 (18) INTEGER GRD61 (18) INTEGER GRD62 (18) INTEGER GRD63 (18) INTEGER GRD64 (18) INTEGER GRD85 (18) INTEGER GRD86 (18) INTEGER GRD87 (18) INTEGER GRD90 (18) INTEGER GRD91 (18) INTEGER GRD92 (18) INTEGER GRD93 (18) INTEGER GRD94 (18) INTEGER GRD95 (18) INTEGER GRD96 (18) INTEGER GRD97 (18) INTEGER GRD98 (18) INTEGER GRD100(18) INTEGER GRD101(18) INTEGER GRD103(18) INTEGER GRD104(18) INTEGER GRD105(18) INTEGER GRD106(18) INTEGER GRD107(18) INTEGER GRD126(18) INTEGER GRD201(18) INTEGER GRD202(18) INTEGER GRD203(18) INTEGER GRD204(18) INTEGER GRD205(18) INTEGER GRD206(18) INTEGER GRD207(18) INTEGER GRD208(18) INTEGER GRD209(18) INTEGER GRD210(18) INTEGER GRD211(18) INTEGER GRD212(18) INTEGER GRD213(18) INTEGER GRD214(18) INTEGER GRD215(18) INTEGER GRD216(18) C SAVE C DATA GRD1 / 0, 255, 1, 73, 23, -48090, 0, 128, 48090, & 0, 513669,513669, 22500, 64, 0, 0, 0, 0/ DATA GRD2 / 0, 255, 0, 144, 73, 90000, 0, 128, -90000, & -2500, 2500, 2500, 0, 0, 0, 0, 0, 0/ DATA GRD3 / 0, 255, 0, 360,181, 90000, 0, 128, -90000, & -1000, 1000, 1000, 0, 0, 0, 0, 0, 0/ DATA GRD4 / 0, 255, 0, 720,361, 90000, 0, 128, -90000, & -500, 500, 500, 0, 0, 0, 0, 0, 0/ DATA GRD5 / 0, 255, 5, 53, 57, 7647, -133443, 8, -105000, & 190500, 190500, 0, 64, 0, 0, 0, 0, 0/ DATA GRD6 / 0, 255, 5, 53, 45, 7647, -133443, 8, -105000, & 190500, 190500, 0, 64, 0, 0, 0, 0, 0/ DATA GRD21 / 0, 255, 0, 37, 37, 0, 0, 128, 90000, & 180000, 2500, 5000, 64, 0, 0, 0, 0, 0/ DATA GRD22 / 0, 255, 0, 37, 37, 0, -180000, 128, 90000, & 0, 2500, 5000, 64, 0, 0, 0, 0, 0/ DATA GRD23 / 0, 255, 0, 37, 37, -90000, 0, 128, 0, & 180000, 2500, 5000, 64, 0, 0, 0, 0, 0/ DATA GRD24 / 0, 255, 0, 37, 37, -90000, -180000, 128, 0, & 0, 2500, 5000, 64, 0, 0, 0, 0, 0/ DATA GRD25 / 0, 255, 0, 72, 19, 0, 0, 128, 90000, & 355000, 5000, 5000, 64, 0, 0, 0, 0, 0/ DATA GRD26 / 0, 255, 0, 72, 19, -90000, 0, 128, 0, & 355000, 5000, 5000, 64, 0, 0, 0, 0, 0/ DATA GRD27 / 0, 255, 5, 65, 65, -20286, -125000, 8, -80000, & 381000, 381000, 0, 64, 0, 0, 0, 0, 0/ DATA GRD28 / 0, 255, 5, 65, 65, 20286, 145000, 8, 100000, & 381000, 381000,128, 64, 0, 0, 0, 0, 0/ DATA GRD29 / 0, 255, 0, 145, 37, 0, 0, 128, 90000, & 360000, 2500, 2500, 64, 0, 0, 0, 0, 0/ DATA GRD30 / 0, 255, 0, 145, 37, 0, -90000, 128, 0, & 360000, 2500, 2500, 64, 0, 0, 0, 0, 0/ DATA GRD33 / 0, 255, 0, 181, 46, 0, 0, 128, 90000, & 360000, 2000, 2000, 64, 0, 0, 0, 0, 0/ DATA GRD34 / 0, 255, 0, 181, 46, -90000, 0, 128, 0, & 360000, 2000, 2000, 64, 0, 0, 0, 0, 0/ DATA GRD37 / 0, 33, 0,65535,73, 0, -30000, 128, 90000, & 60000, 1250,65535, 64, 0, 0, 0, 0, 0, & 73, 73, 73, 73, 73, 73, 73, 73, 72, 72, 72, 71, 71, 71, 70, & 70, 69, 69, 68, 67, 67, 66, 65, 65, 64, 63, 62, 61, 60, 60, & 59, 58, 57, 56, 55, 54, 52, 51, 50, 49, 48, 47, 45, 44, 43, & 42, 40, 39, 38, 36, 35, 33, 32, 30, 29, 28, 26, 25, 23, 22, & 20, 19, 17, 16, 14, 12, 11, 9, 8, 6, 5, 3, 2/ DATA GRD38 / 0, 33, 0,65535,73, 0, 60000, 128, 90000, & 150000, 1250,65535, 64, 0, 0, 0, 0, 0, & 73, 73, 73, 73, 73, 73, 73, 73, 72, 72, 72, 71, 71, 71, 70, & 70, 69, 69, 68, 67, 67, 66, 65, 65, 64, 63, 62, 61, 60, 60, & 59, 58, 57, 56, 55, 54, 52, 51, 50, 49, 48, 47, 45, 44, 43, & 42, 40, 39, 38, 36, 35, 33, 32, 30, 29, 28, 26, 25, 23, 22, & 20, 19, 17, 16, 14, 12, 11, 9, 8, 6, 5, 3, 2/ DATA GRD39 / 0, 33, 0,65535,73, 0, 150000, 128, 90000, & -120000, 1250,65535, 64, 0, 0, 0, 0, 0, & 73, 73, 73, 73, 73, 73, 73, 73, 72, 72, 72, 71, 71, 71, 70, & 70, 69, 69, 68, 67, 67, 66, 65, 65, 64, 63, 62, 61, 60, 60, & 59, 58, 57, 56, 55, 54, 52, 51, 50, 49, 48, 47, 45, 44, 43, & 42, 40, 39, 38, 36, 35, 33, 32, 30, 29, 28, 26, 25, 23, 22, & 20, 19, 17, 16, 14, 12, 11, 9, 8, 6, 5, 3, 2/ DATA GRD40 / 0, 33, 0,65535,73, 0, -120000, 128, 90000, & -30000, 1250,65535, 64, 0, 0, 0, 0, 0, & 73, 73, 73, 73, 73, 73, 73, 73, 72, 72, 72, 71, 71, 71, 70, & 70, 69, 69, 68, 67, 67, 66, 65, 65, 64, 63, 62, 61, 60, 60, & 59, 58, 57, 56, 55, 54, 52, 51, 50, 49, 48, 47, 45, 44, 43, & 42, 40, 39, 38, 36, 35, 33, 32, 30, 29, 28, 26, 25, 23, 22, & 20, 19, 17, 16, 14, 12, 11, 9, 8, 6, 5, 3, 2/ DATA GRD41 / 0, 33, 0,65535,73, -90000, -30000, 128, 0, & 60000, 1250,65535, 64, 0, 0, 0, 0, 0, & 2, 3, 5, 6, 8, 9, 11, 12, 14, 16, 17, 19, 20, 22, 23, & 25, 26, 28, 29, 30, 32, 33, 35, 36, 38, 39, 40, 42, 43, 44, & 45, 47, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 59, 60, 60, & 61, 62, 63, 64, 65, 65, 66, 67, 67, 68, 69, 69, 70, 70, 71, & 71, 71, 72, 72, 72, 73, 73, 73, 73, 73, 73, 73, 73/ DATA GRD42 / 0, 33, 0,65535,73, -90000, 60000, 128, 0, & 150000, 1250,65535, 64, 0, 0, 0, 0, 0, & 2, 3, 5, 6, 8, 9, 11, 12, 14, 16, 17, 19, 20, 22, 23, & 25, 26, 28, 29, 30, 32, 33, 35, 36, 38, 39, 40, 42, 43, 44, & 45, 47, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 59, 60, 60, & 61, 62, 63, 64, 65, 65, 66, 67, 67, 68, 69, 69, 70, 70, 71, & 71, 71, 72, 72, 72, 73, 73, 73, 73, 73, 73, 73, 73/ DATA GRD43 / 0, 33, 0,65535,73, -90000, 150000, 128, 0, & -120000, 1250,65535, 64, 0, 0, 0, 0, 0, & 2, 3, 5, 6, 8, 9, 11, 12, 14, 16, 17, 19, 20, 22, 23, & 25, 26, 28, 29, 30, 32, 33, 35, 36, 38, 39, 40, 42, 43, 44, & 45, 47, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 59, 60, 60, & 61, 62, 63, 64, 65, 65, 66, 67, 67, 68, 69, 69, 70, 70, 71, & 71, 71, 72, 72, 72, 73, 73, 73, 73, 73, 73, 73, 73/ DATA GRD44 / 0, 33, 0,65535,73, -90000, -120000, 128, 0, & -30000, 1250,65535, 64, 0, 0, 0, 0, 0, & 2, 3, 5, 6, 8, 9, 11, 12, 14, 16, 17, 19, 20, 22, 23, & 25, 26, 28, 29, 30, 32, 33, 35, 36, 38, 39, 40, 42, 43, 44, & 45, 47, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 59, 60, 60, & 61, 62, 63, 64, 65, 65, 66, 67, 67, 68, 69, 69, 70, 70, 71, & 71, 71, 72, 72, 72, 73, 73, 73, 73, 73, 73, 73, 73/ DATA GRD45 / 0, 255, 0, 288,145, 90000, 0, 128, -90000, & -1250, 1250, 1250, 0, 0, 0, 0, 0, 0/ DATA GRD55 / 0, 255, 5, 87, 71, -10947, -154289, 8, -105000, & 254000, 254000, 0, 64, 0, 0, 0, 0, 0/ DATA GRD56 / 0, 255, 5, 87, 71, 7647, -133443, 8, -105000, & 127000, 127000, 0, 64, 0, 0, 0, 0, 0/ DATA GRD61 / 0, 255, 0, 91, 46, 0, 0, 128, 90000, & 180000, 2000, 2000, 64, 0, 0, 0, 0, 0/ DATA GRD62 / 0, 255, 0, 91, 46, 0, -180000, 128, 90000, & 0, 2000, 2000, 64, 0, 0, 0, 0, 0/ DATA GRD63 / 0, 255, 0, 91, 46, 0, -90000, 128, 0, & 180000, 2000, 2000, 64, 0, 0, 0, 0, 0/ DATA GRD64 / 0, 255, 0, 91, 46, -90000, -180000, 128, 0, & 0, 2000, 2000, 64, 0, 0, 0, 0, 0/ DATA GRD85 / 0, 255, 0, 360, 90, 500, 500, 128, 89500, & 359500, 1000, 1000, 64, 0, 0, 0, 0, 0/ DATA GRD86 / 0, 255, 0, 360, 90, -89500, 500, 128, -500, & 359500, 1000, 1000, 64, 0, 0, 0, 0, 0/ DATA GRD87 / 0, 255, 5, 81, 62, 22876, -120491, 8, -105000, & 68153, 68153, 0, 64, 0, 0, 0, 0, 0/ DATA GRD90 / 0, 255,201,12902,1, 182, -149887, 136, 92, & 141, 577,538,64, 0, 0, 0, 0, 0/ DATA GRD91 / 0, 255,202,25803,1, 182, -149887, 136, 183, & 141, 577,538,64, 0, 0, 0, 0, 0/ DATA GRD92 / 0, 255,201,24162,1, 9678, -128826, 136, 127, & 191, 278,263,64, 0, 0, 0, 0, 0/ DATA GRD93 / 0, 255,202,48323,1, 9678, -128826, 136, 253, & 191, 278,263,64, 0, 0, 0, 0, 0/ DATA GRD94 / 0, 255,201,48916,1, 9678, -128826, 136, 181, & 271, 194,185,64, 0, 0, 0, 0, 0/ DATA GRD95 / 0, 255,202,97831,1, 9678, -128826, 136, 361, & 271, 194,185,64, 0, 0, 0, 0, 0/ DATA GRD96 / 0, 255,201,41630,1, -3441, -148799, 136, 160, & 261, 333,308,64, 0, 0, 0, 0, 0/ DATA GRD97 / 0, 255,202,83259,1, -3441, -148799, 136, 319, & 261, 333,308,64, 0, 0, 0, 0, 0/ DATA GRD98 / 0, 255, 4, 192, 94, 88542, 0, 128, -88542, & -938, 47,1875, 0, 0, 0, 0, 0, 0/ DATA GRD100/ 0, 255, 5, 83, 83, 17108, -129296, 8, -105000, & 91452, 91452, 0, 64, 0, 0, 0, 0, 0/ DATA GRD101/ 0, 255, 5, 113, 91, 10528, -137146, 8, -105000, & 91452, 91452, 0, 64, 0, 0, 0, 0, 0/ DATA GRD103/ 0, 255, 5, 65, 56, 22405, -121352, 8, -105000, & 91452, 91452, 0, 64, 0, 0, 0, 0, 0/ DATA GRD104/ 0, 255, 5, 147,110, -268, -139475, 8, -105000, & 90755, 90755, 0, 64, 0, 0, 0, 0, 0/ DATA GRD105/ 0, 255, 5, 83, 83, 17529, -129296, 8, -105000, & 90755, 90755, 0, 64, 0, 0, 0, 0, 0/ DATA GRD106/ 0, 255, 5, 165,117, 17533, -129296, 8, -105000, & 45373, 45373, 0, 64, 0, 0, 0, 0, 0/ DATA GRD107/ 0, 255, 5, 120, 92, 23438, -120168, 8, -105000, & 45373, 45373, 0, 64, 0, 0, 0, 0, 0/ DATA GRD126/ 0, 255, 4, 384,190, 89277, 0, 128, -89277, & -938, 95, 938, 0, 0, 0, 0, 0, 0/ DATA GRD201/ 0, 255, 5, 65, 65, -20826, -150000, 8, -105000, & 381000, 381000, 0, 64, 0, 0, 0, 0, 0/ DATA GRD202/ 0, 255, 5, 65, 43, 7838, -141028, 8, -105000, & 190500, 190500, 0, 64, 0, 0, 0, 0, 0/ DATA GRD203/ 0, 255, 5, 45, 39, 19132, -185837, 8, -150000, & 190500, 190500, 0, 64, 0, 0, 0, 0, 0/ DATA GRD204/ 0, 255, 1, 93, 68, -25000, 110000, 128, 60644, & -109129, 160000, 160000, 20000, 64, 0, 0, 0, 0/ DATA GRD205/ 0, 255, 5, 45, 39, 616, -84904, 8, -60000, & 190500, 190500, 0, 64, 0, 0, 0, 0, 0/ DATA GRD206/ 0, 255, 3, 51, 41, 22289, -117991, 8, - 95000, & 81271, 81271, 0, 64, 0, 25000, 25000, 0, 0/ DATA GRD207/ 0, 255, 5, 49, 35, 42085, -175641, 8, -150000, & 95250, 95250, 0, 64, 0, 0, 0, 0, 0/ DATA GRD208/ 0, 255, 1, 29, 27, 9343, -167315, 128, 28092, & -145878, 80000, 80000, 20000, 64, 0, 0, 0, 0/ DATA GRD209/ 0, 255, 3, 101, 81, 22289, -117991, 8, -95000, & 40635, 40635, 0, 64, 0, 25000, 25000, 0, 0/ DATA GRD210/ 0, 255, 1, 25, 25, 9000, -77000, 128, 26422, & -58625, 80000, 80000, 20000, 64, 0, 0, 0, 0/ DATA GRD211/ 0, 255, 3, 93, 65, 12190, -133459, 8, -95000, & 81271, 81271, 0, 64, 0, 25000, 25000, 0, 0/ DATA GRD212/ 0, 255, 3, 185,129, 12190, -133459, 8, -95000, & 40635, 40635, 0, 64, 0, 25000, 25000, 0, 0/ DATA GRD213/ 0, 255, 5, 129, 85, 7838, -141028, 8, -105000, & 95250, 95250, 0, 64, 0, 0, 0, 0, 0/ DATA GRD214/ 0, 255, 5, 97, 69, 42085, -175641, 8, -150000, & 47625, 47625, 0, 64, 0, 0, 0, 0, 0/ DATA GRD215/ 0, 255, 3, 369,257, 12190, -133459, 8, -95000, & 20318, 20318, 0, 64, 0, 25000, 25000, 0, 0/ DATA GRD216/ 0, 255, 5, 139,107, 30000, -173000, 8, -135000, & 45000, 45000, 0, 64, 0, 0, 0, 0, 0/ C IERR = 0 C DO 1 I = 1,18 IGDS(I) = 0 1 CONTINUE C IF (IGRID.GE.37.AND.IGRID.LE.44) THEN DO 2 I = 19,91 IGDS(I) = 0 2 CONTINUE END IF C IF (IGRID.EQ.1) THEN DO 3 I = 1,14 IGDS(I) = GRD1(I) 3 CONTINUE C ELSE IF (IGRID.EQ.2) THEN DO 4 I = 1,14 IGDS(I) = GRD2(I) 4 CONTINUE C ELSE IF (IGRID.EQ.3) THEN DO 5 I = 1,14 IGDS(I) = GRD3(I) 5 CONTINUE C ELSE IF (IGRID.EQ.4) THEN DO 6 I = 1,14 IGDS(I) = GRD4(I) 6 CONTINUE C ELSE IF (IGRID.EQ.5) THEN DO 10 I = 1,14 IGDS(I) = GRD5(I) 10 CONTINUE C ELSE IF (IGRID.EQ.6) THEN DO 20 I = 1,14 IGDS(I) = GRD6(I) 20 CONTINUE C ELSE IF (IGRID.EQ.21) THEN DO 30 I = 1,14 IGDS(I) = GRD21(I) 30 CONTINUE C ELSE IF (IGRID.EQ.22) THEN DO 40 I = 1,14 IGDS(I) = GRD22(I) 40 CONTINUE C ELSE IF (IGRID.EQ.23) THEN DO 50 I = 1,14 IGDS(I) = GRD23(I) 50 CONTINUE C ELSE IF (IGRID.EQ.24) THEN DO 60 I = 1,14 IGDS(I) = GRD24(I) 60 CONTINUE C ELSE IF (IGRID.EQ.25) THEN DO 70 I = 1,14 IGDS(I) = GRD25(I) 70 CONTINUE C ELSE IF (IGRID.EQ.26) THEN DO 80 I = 1,14 IGDS(I) = GRD26(I) 80 CONTINUE C ELSE IF (IGRID.EQ.27) THEN DO 90 I = 1,14 IGDS(I) = GRD27(I) 90 CONTINUE C ELSE IF (IGRID.EQ.28) THEN DO 100 I = 1,14 IGDS(I) = GRD28(I) 100 CONTINUE C ELSE IF (IGRID.EQ.29) THEN DO 110 I = 1,14 IGDS(I) = GRD29(I) 110 CONTINUE C ELSE IF (IGRID.EQ.30) THEN DO 120 I = 1,14 IGDS(I) = GRD30(I) 120 CONTINUE C ELSE IF (IGRID.EQ.33) THEN DO 130 I = 1,14 IGDS(I) = GRD33(I) 130 CONTINUE C ELSE IF (IGRID.EQ.34) THEN DO 140 I = 1,14 IGDS(I) = GRD34(I) 140 CONTINUE C ELSE IF (IGRID.EQ.37) THEN DO 141 I = 1,91 IGDS(I) = GRD37(I) 141 CONTINUE C ELSE IF (IGRID.EQ.38) THEN DO 142 I = 1,91 IGDS(I) = GRD38(I) 142 CONTINUE C ELSE IF (IGRID.EQ.39) THEN DO 143 I = 1,91 IGDS(I) = GRD39(I) 143 CONTINUE C ELSE IF (IGRID.EQ.40) THEN DO 144 I = 1,91 IGDS(I) = GRD40(I) 144 CONTINUE C ELSE IF (IGRID.EQ.41) THEN DO 145 I = 1,91 IGDS(I) = GRD41(I) 145 CONTINUE C ELSE IF (IGRID.EQ.42) THEN DO 146 I = 1,91 IGDS(I) = GRD42(I) 146 CONTINUE C ELSE IF (IGRID.EQ.43) THEN DO 147 I = 1,91 IGDS(I) = GRD43(I) 147 CONTINUE C ELSE IF (IGRID.EQ.44) THEN DO 148 I = 1,91 IGDS(I) = GRD44(I) 148 CONTINUE ELSE IF (IGRID.EQ.45) THEN DO 149 I = 1,14 IGDS(I) = GRD45(I) 149 CONTINUE C C ELSE IF (IGRID.EQ.50) THEN C DO 150 I = 1,14 C IGDS(I) = GRD50(I) C150 CONTINUE C ELSE IF (IGRID.EQ.55) THEN DO 152 I = 1,14 IGDS(I) = GRD55(I) 152 CONTINUE C ELSE IF (IGRID.EQ.56) THEN DO 154 I = 1,14 IGDS(I) = GRD56(I) 154 CONTINUE C ELSE IF (IGRID.EQ.61) THEN DO 160 I = 1,14 IGDS(I) = GRD61(I) 160 CONTINUE C ELSE IF (IGRID.EQ.62) THEN DO 170 I = 1,14 IGDS(I) = GRD62(I) 170 CONTINUE C ELSE IF (IGRID.EQ.63) THEN DO 180 I = 1,14 IGDS(I) = GRD63(I) 180 CONTINUE C ELSE IF (IGRID.EQ.64) THEN DO 190 I = 1,14 IGDS(I) = GRD64(I) 190 CONTINUE C ELSE IF (IGRID.EQ.85) THEN DO 192 I = 1,14 IGDS(I) = GRD85(I) 192 CONTINUE C ELSE IF (IGRID.EQ.86) THEN DO 194 I = 1,14 IGDS(I) = GRD86(I) 194 CONTINUE C ELSE IF (IGRID.EQ.87) THEN DO 195 I = 1,14 IGDS(I) = GRD87(I) 195 CONTINUE C ELSE IF (IGRID.EQ.90) THEN DO 196 I = 1,14 IGDS(I) = GRD90(I) 196 CONTINUE C ELSE IF (IGRID.EQ.91) THEN DO 197 I = 1,14 IGDS(I) = GRD91(I) 197 CONTINUE C ELSE IF (IGRID.EQ.92) THEN DO 198 I = 1,14 IGDS(I) = GRD92(I) 198 CONTINUE C ELSE IF (IGRID.EQ.93) THEN DO 199 I = 1,14 IGDS(I) = GRD93(I) 199 CONTINUE C ELSE IF (IGRID.EQ.94) THEN DO 200 I = 1,14 IGDS(I) = GRD94(I) 200 CONTINUE C ELSE IF (IGRID.EQ.95) THEN DO 201 I = 1,14 IGDS(I) = GRD95(I) 201 CONTINUE C ELSE IF (IGRID.EQ.96) THEN DO 202 I = 1,14 IGDS(I) = GRD96(I) 202 CONTINUE C ELSE IF (IGRID.EQ.97) THEN DO 203 I = 1,14 IGDS(I) = GRD97(I) 203 CONTINUE C ELSE IF (IGRID.EQ.98) THEN DO 204 I = 1,14 IGDS(I) = GRD98(I) 204 CONTINUE C ELSE IF (IGRID.EQ.100) THEN DO 205 I = 1,14 IGDS(I) = GRD100(I) 205 CONTINUE C ELSE IF (IGRID.EQ.101) THEN DO 210 I = 1,14 IGDS(I) = GRD101(I) 210 CONTINUE C ELSE IF (IGRID.EQ.103) THEN DO 220 I = 1,14 IGDS(I) = GRD103(I) 220 CONTINUE C ELSE IF (IGRID.EQ.104) THEN DO 230 I = 1,14 IGDS(I) = GRD104(I) 230 CONTINUE C ELSE IF (IGRID.EQ.105) THEN DO 240 I = 1,14 IGDS(I) = GRD105(I) 240 CONTINUE C ELSE IF (IGRID.EQ.106) THEN DO 242 I = 1,14 IGDS(I) = GRD106(I) 242 CONTINUE C ELSE IF (IGRID.EQ.107) THEN DO 244 I = 1,14 IGDS(I) = GRD107(I) 244 CONTINUE C ELSE IF (IGRID.EQ.126) THEN DO 245 I = 1,14 IGDS(I) = GRD126(I) 245 CONTINUE C ELSE IF (IGRID.EQ.201) THEN DO 250 I = 1,14 IGDS(I) = GRD201(I) 250 CONTINUE C ELSE IF (IGRID.EQ.202) THEN DO 260 I = 1,14 IGDS(I) = GRD202(I) 260 CONTINUE C ELSE IF (IGRID.EQ.203) THEN DO 270 I = 1,14 IGDS(I) = GRD203(I) 270 CONTINUE C ELSE IF (IGRID.EQ.204) THEN DO 280 I = 1,14 IGDS(I) = GRD204(I) 280 CONTINUE C ELSE IF (IGRID.EQ.205) THEN DO 290 I = 1,14 IGDS(I) = GRD205(I) 290 CONTINUE C ELSE IF (IGRID.EQ.206) THEN DO 300 I = 1,18 IGDS(I) = GRD206(I) 300 CONTINUE C ELSE IF (IGRID.EQ.207) THEN DO 310 I = 1,14 IGDS(I) = GRD207(I) 310 CONTINUE C ELSE IF (IGRID.EQ.208) THEN DO 320 I = 1,14 IGDS(I) = GRD208(I) 320 CONTINUE C ELSE IF (IGRID.EQ.209) THEN DO 330 I = 1,18 IGDS(I) = GRD209(I) 330 CONTINUE C ELSE IF (IGRID.EQ.210) THEN DO 340 I = 1,14 IGDS(I) = GRD210(I) 340 CONTINUE C ELSE IF (IGRID.EQ.211) THEN DO 350 I = 1,18 IGDS(I) = GRD211(I) 350 CONTINUE C ELSE IF (IGRID.EQ.212) THEN DO 360 I = 1,18 IGDS(I) = GRD212(I) 360 CONTINUE C ELSE IF (IGRID.EQ.213) THEN DO 370 I = 1,14 IGDS(I) = GRD213(I) 370 CONTINUE C ELSE IF (IGRID.EQ.214) THEN DO 380 I = 1,14 IGDS(I) = GRD214(I) 380 CONTINUE C ELSE IF (IGRID.EQ.215) THEN DO 390 I = 1,18 IGDS(I) = GRD215(I) 390 CONTINUE C ELSE IF (IGRID.EQ.216) THEN DO 400 I = 1,14 IGDS(I) = GRD216(I) 400 CONTINUE C ELSE IERR = 1 ENDIF C RETURN END