我正在使用cdo(气候数据运算符)以便将具有Lambert保形圆锥投影(分辨率为1km)的NetCDF文件重新映射到具有相同投影的较粗网格(2km)。
我的lcc投影具有球面基准(“ + proj = lcc + lat_1 = 30 + lat_2 = 60 + lat_0 = 26 + lon_0 = 13 + x_0 = 0 + y_0 = 0 + a = 6370000 + b = 6370000 + units = m + no_defs”)。
在cdo中,我定义目标网格描述如下
cat > $Database_outdir/mygrid.txt <<EOF
gridtype = projection
gridsize = 38809
xsize = 197
ysize = 197
xname = x
xlongname = "x coordinate of projection"
xunits = "m"
yname = y
ylongname = "y coordinate of projection"
yunits = "m"
xfirst = $ORIGX
xinc = 2000
yfirst = $ORIGY
yinc = 2000
grid_mapping = lambert_conformal_conic
grid_mapping_name = lambert_conformal_conic
longitude_of_central_meridian = 13
false_easting = 0.
false_northing = 0.
latitude_of_projection_origin = 26
standard_parallel = 30. 60.
longitude_of_prime_meridian = 0.
earth_radius=6370000.
longitudeOfFirstGridPointInDegrees = $orig_lon
latitudeOfFirstGridPointInDegrees =$orig_lat
EOF
当我尝试使用cdo重新映射
cdo remapcon, mygrid.txt grid1km.nc out_grid2km.nc
我收到以下错误cdo remapcon(Abort):项目错误:主轴或半径= 0或未给出。
对此有任何帮助吗?它与自定义球面基准有关吗?我如何在cdo中定义它?