load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

fin = addfile("air.2m.gauss.2001.nc", "r")

LON = 260.625
LAT = fin->lat(35:37)
print(LAT)
print("")
print("---------- AT GRID CENTER --------")
do i = 0,2

  print("")
  lat1 = LAT(i)

  do delta = -3,3

    lat2 = lat1 + delta*0.0001

    temp = short2flt(fin->air(time|:, level|0, {lat|lat2}, {lon|LON}))

    print(sprintf("%7.4f",lat2)+" | "+temp(0)+"	"+temp(1)+"	"+temp(2)+"	"+temp(3))
  end do

end do


print("-------- HALFWAY BETWEEN ------")
lat1 = (19.9997 + 18.095)/2.0

do delta = -3,3
  lat2 = lat1 + delta*0.0001
  temp = short2flt(fin->air(time|:, level|0, {lat|lat2}, {lon|LON}))
  print(sprintf("%7.5f",lat2)+" | "+temp(0)+"	"+temp(1)+"	"+temp(2)+"	"+temp(3))
end do

print("")
lat1 = (19.9997 + 21.9044)/2.0

do delta = -3,3
  lat2 = lat1 + delta*0.0001  
  temp = short2flt(fin->air(time|:, level|0, {lat|lat2}, {lon|LON}))
  print(sprintf("%7.4f",lat2)+" | "+temp(0)+"	"+temp(1)+"	"+temp(2)+"	"+temp(3))
end do