【已解决】grd用cdo转成nc格式，但是grd有pdef参数

心如止水的饭团

我把grd用cdo转成nc格式的，但是grd的ctl有pdef参数，因为cdo不支出，所以我把pdef参数删了。现在用ncl画图，发现这个图的数据有偏移，我觉得应该是投影的问题，可是大家还有什么解决方法吗？就是想用ncl画这个grd，不用grads。
这个是grd的ctl
dset C:\Users\fan\Downloads\datagrid\wrfd01.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.5 MODEL
pdef  180 170 lcc  24.000  106.000   91.500   86.500  44.00000   4.00000  106.00000  45000.000  45000.000
xdef  180 linear   47.12270   0.20270270
ydef  170 linear  -16.35261   0.20270270
tdef    1 linear 00Z01MAY2014      60MN      
zdef    1  levels 1 1
tdef    1 linear 00Z01MAY2014      60MN      
VARS    4
LU_INDEX       1  0  LAND USE CATEGORY (-)
HGT            1  0  Terrain Height (m)
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
ENDVARS

我画的图，这个地形高度往左边偏了

图片

解决方法在5楼

心如止水的饭团

过了几天我终于弄出来了，所以来这里自问自答一下。ctl里面的xdef、ydef其实是数据的格点数，并不是真实的经纬度，我所用的资料是1*1的，并且是兰博托正形投影，所以直接把xdef和ydef改成
xdef  180 linear   1  1.
ydef  170 linear   1  1.
而变量xlat和xlong才是真正的经纬度，在ncl的官网中有一个例子就是讲的只知道本地格点的例子，我参照官网例子画的图，其中
mpLambertParallel1F =
mpLambertParallel2F =
mpLambertMeridianF =
这三个参数很重要，两个标准纬度，一个标准经度，要找好这三个数的正确取值，画出来的图就不会偏移啦，官网上也有教我们如何有逻辑的去猜这三个数字的方法
这个是本地格点的链接http://www.ncl.ucar.edu/Applications/lcnative.shtml
这个是大胆猜想标准经纬度的脚本http://www.ncl.ucar.edu/Applications/Scripts/lcnative_3.ncl

心如止水的饭团

输入    cdo -f nc import_binary inputfile.ctl outputfile.nc   
inputfile.ctl 是你需要转格式的grd 的ctl，outputfile.nc是输出的转了格式之后的文件

/xin村儿/

求分享 wrfd01.dat， ctl, 使用例子

心如止水的饭团

/xin村儿/ 发表于 2015-7-31 17:27
求分享 wrfd01.dat， ctl, 使用例子

额，什么使用例子呀？我这里都没弄出来

Figo

可以考虑用cdo先差值改一下分辨率

心如止水的饭团

C:\Users\fan\Downloads\lcnative.pdf最后画的图

kanyu

心如止水的饭团 发表于 2015-8-21 13:07
过了几天我终于弄出来了，所以来这里自问自答一下。ctl里面的xdef、ydef其实是数据的格点数，并不是真实的 ...

你好，请问用cdo将grd数据转为nc是怎么做的啊？可以详细说一下吗？谢谢啦

大森林

谢谢分享，好东东

东风急流

楼主，看了这个帖子有些困惑。
1、如果你的原始数据是1度乘以1度的，为什么ctl不直接写成
xdef  180 linear   1  1.
ydef  170 linear   1  1.
请问你改了ctl以后还可以正常用grads绘图吗？

2、我觉得原始数据确实是在兰伯特投影上的，但应该是等距离的，所以有pdef 180 170 lcc 24.000 106.000 91.500 86.500 44.00000 4.00000 106.00000 45000.000 45000.000这一行。
pedf里面写的是数据是等距离45km格距的，而1度大概是110km啊，这个资料是从哪里下载的呢？像是模式转出来的啊。