25 1 IE(J2.GI.JEND)J2=JEND POUT 24 



26 URIIE(6 f 102)(J,J>Jl,J2) POUT 25 



27 102 F0RMAI(/,3X,'I/J:',31I4) POUT 26 



28 WRIIE(6,103) POUT 27 

 2g 103 FORHAIUX.' ', POUT 28 



30 A ' ', POUT 29 



31 k ' ') POUT 30 



32 DO 2 1=111, 112, IVAL POUT 31 



33 DO 3 J=J1,J2 POUI 32 



34 RND=0.5 POUI 33 



35 IF(DUM1(I,J).LI.O.O)RND=-0.5 POUI 34 



36 3 IX(J)=INT(EACT*DUM1(I.J)+RND) POUT 35 



37 2 WRIIE(6, 104)1, (IX(J),J=J1,J2) POUI 36 



38 104 E0RHAI(1X,I3,2X,':'.31I4) POUT 37 



39 J1=J1+NC POUT 38 



40 J2=J2+NC POUI 39 



41 IE(J1.LE.JEND)G0 10 1 POUI 40 



42 WRIIE(6,101! POUT 41 



43 101 FORMAK//) POUT 42 



44 RETURN POUT 43 



45 END POUT 44 



1 C*A**AA*AAAAAAAA*Aikjk*AAAA*A*A***A*AA*AA*A****AAAAAAA*A*AA*AA)i!**A*AA*AAA* GRID 2 



2 C GRID 3 



3 SUBROUTINE GRID GRID 4 



4 C GRID 5 



5 CAAAAAAAAAAAAA*AA*AAAAAAA*A*AAAAAAAAAAAAAAAAAAAAA**AAAAAAAAAAAAAAAAAAAAA GRID 6 



6 C GRID 7 



7 C THIS SUBROUTINE READS THE VARIABLE GRID HAPPING INFORHATION FROH GRID 8 



8 C FILE CODE 3 (X REGION FIRST, IHEN Y) AND GENERAIES THE VARIABLE GRID 9 



9 C GRID SCALE FACTORS AND DISTANCES TO THE GRID CENTERS GRID 10 



10 C GRID 11 



11 C WHEN GENERATING A VARIABLY SIZED GRID IT IS RECCQHENDED THAT IHE GRID 12 



12 C MAPPING BE DONE IN HAP INCEHES FROH SOHE BAIHYHETRIC CHART OF A GRID 13 



13 C KNOWN SCALE. IHE 'DX' AND 'DY' USED AS INPUT INTO IHE MODEL ARE GRID 14 



14 C IHEN SIHPLY THE NUMBER OF FEEI OR HETERS PER HAP INCH. IF MAP GRID 15 



15 C INCHES ARE NOT USED IN THE MAPPING, IHE VALUES OF 'XSCALE' AND GRID 16 



16 C 'YSCALE' CAN BE SEI 10 SOMETHING OTHER THAN 1.0 (IN CONJUNCTION GRID 17 



17 C WITH THE CHOICE OF 'DX' AND 'DY') 10 CONVERT FROM THE MAPPED GRID 18 



18 C SPACE INTO REAL SPACE GRID 19 



19 C GRID 20 



20 CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA GRID 21 



21 C GRID 22 

 C ^DEFINITIONS OF VARIABLE GRID MAPPING INPUT REQUIRED BY THE MODEL GRID 23 



23 C GRID 24 



24 C NREG(X.Y) - IS IHE NUMBER OF MAPPING REGIONS IN IHE X OR Y GRID 25 

 C DIRECTION GRID 26 

 C NA(X.Y) - IS THE NUMBER OF GRID CELLS IN A PARTICULAR MAPPING GRID 27 

 C REGION IN THE X OR Y DIRECTION GRID 28 

 C A(X,Y),B(X,Y),C(X,Y) - ARE IHE MAPPING COEFFICIENTS FOR A GRID 29 



29 C PARTICULAR REGION IN THE X OR Y DIRECTION GRID 30 



30 C GRID 31 



31 CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA GRID 32 



32 C PARAM 2 



33 PARAHEIER(IQ=95,JQ=95) PARAM 3 



34 C PARAM 4 



35 COHHON/CONSI/ G,PI,PI2,RAD,HC0NVR,SC0NVR,DX,DY,DX2,DY2,I,0MEG, GRID 34 



36 AM,N,M1,M2 ( NH1,NM2,IWET,IDRY,IUETP1,CNTRANG,HFACI,DLEVEL GRID 35 



37 COHHON/C0NS2/IGRID,IREF,ITAMX,ITHHX,IDIFF GRID 36 



38 COHMON/MUSC/XHUC(IQ),XMUS(IQ),YHUC(JQ).YMUS(JQ) GRID 37 



F18 



