subroutine geo2grid(ra,ro,ia,io,ca,co) C converts geographical coordinates of the center of 1x1 square C ( SD[DD].5 where S is + for EAST(NORTH) or - for WEST(SOUTH) C INTO C corresponding grid numbers (1 to 360 for LONG and 1 to 180 C for LAT) and C its character images (e.g., for outfile name construction) character*3 ca character*4 co C convert into characters C LATITUDE if(ra.lt.0.)then ix=int(-1.*ra) ca(3:3)='s' else ix=int(ra) ca(3:3)='n' endif if(ix.le.9)then ca(1:1)='0' write(ca(2:2),'(i1)') ix else write(ca(1:2),'(i2)') ix endif C LONGITUDE if(ro.lt.0.)then ix=int(-1.*ro) co(4:4)='w' else ix=int(ro) co(4:4)='e' endif if(ix.le.9)then co(1:2)='00' write(co(3:3),'(i1)') ix elseif(ix.le.99)then co(1:1)='0' write(co(2:3),'(i2)') ix else write(co(1:3),'(i3)') ix endif C convert into grid numbers ra=ra+90.5 ia=int(ra) if(ro.le.0.)then io=nint(ro+360.) else io=nint(ro) endif return end C************************************************************* subroutine ch12dep(cx) character*12 cx(33) cx(1) =' 0m' cx(2) =' 10m' cx(3) =' 20m' cx(4) =' 30m' cx(5) =' 50m' cx(6) =' 75m' cx(7) =' 100m' cx(8) =' 125m' cx(9) =' 150m' cx(10)=' 200m' cx(11)=' 250m' cx(12)=' 300m' cx(13)=' 400m' cx(14)=' 500m' cx(15)=' 600m' cx(16)=' 700m' cx(17)=' 800m' cx(18)=' 900m' cx(19)=' 1000m' cx(20)=' 1100m' cx(21)=' 1200m' cx(22)=' 1300m' cx(23)=' 1400m' cx(24)=' 1500m' cx(25)=' 1750m' cx(26)=' 2000m' cx(27)=' 2500m' cx(28)=' 3000m' cx(29)=' 3500m' cx(30)=' 4000m' cx(31)=' 4500m' cx(32)=' 5000m' cx(33)=' 5500m' return end