求如何在没有undef的情况下写Ctl文件

tbag

大家好，

我现在想给一个netcdf文件写ctl，这个netcdf无法直接被grads打开，会显示“NetCDF: unknown data format”
我检查了该nc文件后得到如下信息:

netcdf CAT_PRECT_256x512.eam.h1 {
dimensions:
        lat = 256 ;
        lon = 512 ;
        nbnd = 2 ;
        time = UNLIMITED ; // (3936 currently)
variables:
        double lat(lat) ;
                lat:long_name = "Latitude of Grid Cell Centers" ;
                lat:standard_name = "latitude" ;
                lat:units = "degrees_north" ;
                lat:axis = "Y" ;
                lat:valid_min = -90. ;
                lat:valid_max = 90. ;
                lat:bounds = "lat_bnds" ;
        double lon(lon) ;
                lon:long_name = "Longitude of Grid Cell Centers" ;
                lon:standard_name = "longitude" ;
                lon:units = "degrees_east" ;
                lon:axis = "X" ;
                lon:valid_min = 0. ;
                lon:valid_max = 360. ;
                lon:bounds = "lon_bnds" ;
        double lat_bnds(lat, nbnd) ;
                lat_bnds:long_name = "Gridcell latitude interfaces" ;
        double lon_bnds(lon, nbnd) ;
                lon_bnds:long_name = "Gridcell longitude interfaces" ;
        double area(lat, lon) ;
                area:long_name = "Solid angle subtended by gridcell" ;
                area:standard_name = "solid_angle" ;
                area:units = "steradian" ;
                area:cell_mathods = "lat, lon: sum" ;
        float PRECT(time, lat, lon) ;
                PRECT:units = "m/s" ;
                PRECT:long_name = "Total (convective and large-scale) precipitation rate (liq + ice)" ;
                PRECT:cell_measures = "area: area" ;
                PRECT:cell_methods = "time: mean" ;
        double time(time) ;
                time:long_name = "time" ;
                time:units = "days since 2020-01-20 00:00:00" ;
                time:calendar = "noleap" ;
                time:bounds = "time_bnds" ;
                time:cell_methods = "time: mean" ;
        double time_bnds(time, nbnd) ;
                time_bnds:long_name = "time interval endpoints" ;
                time_bnds:cell_methods = "time: mean" ;



我写的ctl如下:

dset ^CAT_PRECT_256x512.eam.h1.nc
title CAT_PRECT
undef -9.99e8
xdef 512 linear 0 0.703125
ydef 256 linear -90 0.703125
zdef 128 linear 0 127
tdef 3936 linear 00Z20JAN2020 15mn
vars 1
precip=>PRECT  0  t,y,x  precipitation rate
endvars


打开之后显示的是乱码 （3e+38..）
我尝试过用cdo 添加missing values (cdo setmissval,-9.99e8 input.nc output.nc), 但结果没有差别。
求问该nc文件如何写ctl文件，谢谢！
我也试过在ctl中加入dtype netcdf，但加入之后文件都无法打开

tbag

自己鼓捣半天搞定了。。
我不知道为什么，我的这个netcdf不是classic netcdf，所以我先将这个文件转化为classic netcdf (nccopy -k 1 input.nc output.nc), 再通过cdo加入missing value，就可以用ctl读取了。