2012-11-28 01:57:45 +00:00
|
|
|
subroutine unpackgrid(ng,grid)
|
2012-05-22 17:09:48 +00:00
|
|
|
|
2012-11-28 01:57:45 +00:00
|
|
|
parameter (NGBASE=180*180)
|
|
|
|
|
character grid*4,grid6*6
|
2012-05-22 17:09:48 +00:00
|
|
|
|
2012-11-28 01:57:45 +00:00
|
|
|
grid=' '
|
|
|
|
|
if(ng.ge.32400) go to 10
|
|
|
|
|
dlat=mod(ng,180)-90
|
|
|
|
|
dlong=(ng/180)*2 - 180 + 2
|
|
|
|
|
call deg2grid(dlong,dlat,grid6)
|
|
|
|
|
grid=grid6(:4)
|
|
|
|
|
if(grid(1:2).eq.'KA') then
|
|
|
|
|
read(grid(3:4),*) n
|
|
|
|
|
n=n-50
|
|
|
|
|
write(grid,1001) n
|
|
|
|
|
1001 format(i3.2)
|
|
|
|
|
if(grid(1:1).eq.' ') grid(1:1)='+'
|
|
|
|
|
else if(grid(1:2).eq.'LA') then
|
|
|
|
|
read(grid(3:4),*) n
|
|
|
|
|
n=n-50
|
|
|
|
|
write(grid,1002) n
|
|
|
|
|
1002 format('R',i3.2)
|
|
|
|
|
if(grid(2:2).eq.' ') grid(2:2)='+'
|
|
|
|
|
endif
|
|
|
|
|
go to 900
|
2012-05-22 17:09:48 +00:00
|
|
|
|
2012-11-28 01:57:45 +00:00
|
|
|
10 n=ng-NGBASE-1
|
|
|
|
|
if(n.ge.1 .and.n.le.30) then
|
|
|
|
|
write(grid,1012) -n
|
|
|
|
|
1012 format(i3.2)
|
|
|
|
|
else if(n.ge.31 .and.n.le.60) then
|
|
|
|
|
n=n-30
|
|
|
|
|
write(grid,1022) -n
|
|
|
|
|
1022 format('R',i3.2)
|
|
|
|
|
else if(n.eq.61) then
|
|
|
|
|
grid='RO'
|
|
|
|
|
else if(n.eq.62) then
|
|
|
|
|
grid='RRR'
|
|
|
|
|
else if(n.eq.63) then
|
|
|
|
|
grid='73'
|
|
|
|
|
endif
|
2012-05-22 17:09:48 +00:00
|
|
|
|
2012-11-28 01:57:45 +00:00
|
|
|
900 return
|
|
|
|
|
end subroutine unpackgrid
|
2012-05-22 17:09:48 +00:00
|
|
|
|
