GetJK.m

function [J,K] = GetJK(Lat,Lon)

global LOv LAv jg kg

X1 = -2600000;
Y1 = 0;
X2 = 2300000;
Y2 = 3200000;
nx = 4901;
ny = 3201;

J = interp2(LOv,LAv,jg,Lon,Lat,'linear');
K = interp2(LOv,LAv,kg,Lon,Lat,'linear');