21 subroutine bilin(data,glamn,glomn,dla,dlo,
22 *nla,nlo,maxla,maxlo,xla,xlo,val)
37 implicit real*8 (a-h,o-z)
38 real*4 data(maxla,maxlo)
39 logical onedlat,onedlon
53 gla0 = glamn + (ila-1)*dla
59 glo0 = glomn + (ilo-1)*dlo
66 val = (1.d0-t)*(1.d0-u)*
data(ila ,ilo )
67 * + ( t)*(1.d0-u)*
data(ila ,ilo+1)
68 * + ( t)*( u)*
data(ila+1,ilo+1)
69 * + (1.d0-t)*( u)*
data(ila+1,ilo )
74 write(6,200)ila,ilo,gla0,glo0
75 write(6,199)xla-gla0,xlo-glo0
76 write(6,201)ila,ilo,gla0,glo0,
78 write(6,201)ila,ilo+1,gla0,glo0+dlo,
80 write(6,201)ila+1,ilo,gla0+dla,glo0,
82 write(6,201)ila+1,ilo+1,gla0+dla,glo0+dlo,
84 write(6,202)
'SW',(1.d0-t)*(1.d0-u)
85 write(6,202)
'SE',( t)*(1.d0-u)
86 write(6,202)
'NW',( t)*( u)
87 write(6,202)
'NE',(1.d0-t)*( u)
89 *
'Weight to ',a,
' corner = ',f20.15)
92 *6x,
'bilin.f: lat/lon = ',f10.6,1x,f10.6)
94 *6x,
'bilin.f: diff from SW corner of cell = ',f10.6,1x,f10.6)
96 *6x,
'bilin.f: SW corner: ',i6,1x,i6,1x,f10.6,1x,f10.6)
98 *6x,2(i6,1x),2(f10.6,1x),f20.15)
102 801
format(f15.8,1x,f15.8,1x,i6,1x,i6,1x,
103 *f15.8,1x,f15.8,1x,f8.6,1x,f8.6,1x,f20.10)
108 802
format(4(f15.8,1x))
subroutine bilin(data, glamn, glomn, dla, dlo, nla, nlo, maxla, maxlo, xla, xlo, val)
Subroutine to perform bilinear interpolation.