subroutine GetOrog (xsize, ysize, orog, optn) dimension orog(37,73) integer topo(72, 36) integer xsize, ysize, xsiz, ysiz, optn external GTDA xsiz= xsize-1 ysiz= ysize-1 call GTDA(xsiz, ysiz, topo, optn) c NOTE: This code is hardcoded to handle 5 degree resolution only. c GTDA is capable of 10 minute resolution but GetOrog will need to be c modified to support it (via modification of topo()). c Now switch the data around and fill in the gaps since the calling c routines array is 37x73 and GTDA only does 5 degree increments in c 72x36. do lat= 1, ysiz do long= 1, xsiz orog(lat, long)= topo(long, lat) enddo enddo orog(37,73)= topo(72, 36) return end