CTITLESAFAN3 -- MAP A 2-D GRID USING UNIRAS SOFTWARE 00010000 C*********************************************************************** C COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991 * C*********************************************************************** CA AUTHOR R. E. MCMILLAN 00020000 CA DESIGNER R. E. MCMILLAN 00030000 CA LANGUAGE FORTRAN 00040000 CA SYSTEM IBM AND CRAY 00050000 CA WRITTEN 8-08-86 00060000 C REVISED 11-25-86 ESN. FOR CONVERSION TO THE CRAY. 00070000 C REVISED 08-12-90 CLJ CORRECT AN ERROR IN THE REVERSING OF C THE RIGHT HALF OF THE ARRAY CA 00080000 CA 00090000 CA CALL SAFAN3 (ZIN, NZX, NZY) 00100000 CA 00110000 CA IN/OUT ARGUMENT TYPE DESCRIPTION 00120000 CA 00130000 CA IN/OUT ZIN R4 2 DIMENSIONAL ARRAY TO BE REARRANGED 00140000 CA IN NZX I4 INSIDE DIMENSION OF ZIN 00150000 CA IN NZY I4 OUTSIDE DIMENSION OF ZIN 00160000 CA 00170000 CA 00180000 CA 00190000 CA THIS SUBROUTINE WILL REARRANGE THE INPUT ARRAY ZIN SO THAT F/K DIP00200000 CA WILL APPEAR IN THE EXPECTED WAY WHEN PLOTTED. 00210000 CAEND 00220000 C EJECT 00230000 C 00240000 SUBROUTINE SAFAN3 ( ZIN, NZX, NZY) 00250000 C 00260000 REAL ZIN(NZX,NZY) 00270000 C 00280000 C 00290000 NZY2 = NZY / 2 00300000 NZY4 = NZY / 4 00310000 KKB = NZY2 + 1 00320000 KKE = KKB + NZY4 - 1 00330000 C 00340000 DO 300 JJ = 1, NZX 00350000 C 00360000 K = NZY2 00370000 C 00380000 DO 100 II = 1, NZY4 00390000 A = ZIN(JJ,K) 00400000 ZIN(JJ,K) = ZIN(JJ,II) 00410000 ZIN(JJ,II) = A 00420000 K = K - 1 00430000 100 CONTINUE 00440000 C 00450000 K = NZY 00460000 C 00470000 DO 200 II = KKB, KKE 00480000 A = ZIN(JJ,K) 00490000 ZIN(JJ,K) = ZIN(JJ,II) 00500000 ZIN(JJ,II) = A 00510000 K = K - 1 00520000 200 CONTINUE 00530000 C 00540000 300 CONTINUE 00550000 C 00560000 RETURN 00570000 END 00580000