gc_latlon球面距离的类似函数

kathy-

pre_field= where(((gc_latlon(tylat,tylon,pre_field&lat,pre_field&lon,2,4)).ge.radius),_fillvalue,pre_field)tylat,tylon是一个台风的经纬度
pre_field是一个73*144的全球降水观测场

想把一个循环简化成上述第一行的代码，但是gc_latlon函数要求pre_field&lat,pre_field&lon是同等维度，想问有没有函数是可以求这样的一个确定的维度和一张确定的网格点之间的距离？

不想编循环的原因，是因为循环跑起来实在太慢了