关于grads读取grd数据中头文件ctl问题

海洋小卫士

向各位大佬请教一下：
grd数据格式是48(time) x 20(level) x lat(64) x lon(128)的逐日气温数据，根据我自己写的ctl文件不能合理绘制出气温空间分布图
望各位大佬指点迷津！！
图1是基于自己写的ctl文件绘制的出图
相关命令
open forc_winter.ctl
d tmp


正确的图应该如图2所示


我的ctl文件为
=====================================
* NCEP climatology interporated into T42 Gaussian grid
DSET ^forc_winter.grd
OPTIONS SEQUENTIAL
TITLE dumy
UNDEF -999.
XDEF 128 LINEAR 0. 2.81250
YDEF 64  LEVELS
-87.864 -85.097 -82.313 -79.526 -76.737 -73.948 -71.158 -68.368 -65.578
-62.787 -59.997 -57.207 -54.416 -51.626 -48.835 -46.045 -43.254 -40.464
-37.673 -34.883 -32.092 -29.301 -26.511 -23.720 -20.930 -18.139 -15.348
-12.558  -9.767  -6.976  -4.186  -1.395   1.395   4.186  6.976   9.767  
12.558  15.348  18.139  20.930  23.720  26.511  29.301 32.092  34.883  
37.673  40.464  43.254  46.045  48.835  51.626  54.416  57.207  59.997  
62.787  65.578  68.368  71.158  73.948  76.737  79.526  82.313  85.097  
87.864
ZDEF 20 LEVELS 0.99500 0.97999 0.94995 0.89988 0.82977 0.74468
0.64954 0.54946 0.45447 0.36948 0.29450 0.22953 0.17457 0.12440
0.0846830 0.0598005 0.0449337 0.0349146 0.0248800 0.00829901
TDEF 48 LINEAR 15jan0000 1dy
VARS 1
tmp  20 99 heating forcing [K s**-1]

ENDVARS

=============================================

注：1. 生成grd格式的命令是基于NCL中的fbinwrite()
2. 通过网盘分享的文件：forc_winter.grd
链接: https://pan.baidu.com/s/176Wl5wQQ4ANAAoSD2OAE1A 提取码: 1234

edwardli

ydef 为什么是levels？

海洋小卫士

edwardli 发表于 2025-11-8 22:43
ydef 为什么是levels？

我是参考线性斜压模式（LBM）里面的ctl文件写的