grads中能否实现资料的极坐标化

ぁ砾

我想在这套矩阵资料中取一个网格点作为极点，然后将直角坐标系中的数据转换成极坐标系中的数据，原来用Fortran处理的想法是，先将直角坐标系中的格点转换成极坐标系中的格点，然后使用线性插值转换数据，现在遇到了问题是转换出来的极坐标系格点是不均匀的，不知道怎样插值才能避免较大的误差，我之前也查过资料，没有找到grads或者Fortran能用来直接转换坐标系的函数或者脚本程序什么的，想了很久，还是上论坛来请教下各位大牛，不知道能否得到一个Fortran编程转换的思想或者是我不知道的一个Fortran函数或者grads的命令脚本能帮我实现转化的。