Я использую cdo (операторы климатических данных) для переназначения файла NetCDF с конформной конической проекцией Ламберта (с разрешением 1 км) в более грубую сетку (2 км) с такой же проекцией.
Моя 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 + единиц = 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 (Прервать): ошибка проекта: большая ось или радиус = 0 или не задано.
Помогите ли вы с этим?Связано ли это с пользовательским сферическим датумом?Как я мог определить это в cdo?