Generating weights and remapping MPAS-limited area outputs
Added by Kelly Núñez Ocasio about 2 years ago
Hello,
I am trying to generate weights to remap 3-km MPAS limited-area outputs on a 30km lat/lon grid. Unfortunately, I haven't been successful. Below are the commands I have used in the past to create weights and to remap global 15-km mesh MPAS outputs. Has anyone successfully regridded MPAS limited-area data before?
Generate weights:
cdo -P 1 --cellsearchmethod spherepart -gencon,r1440x720 -setgrid,mpas:${src_grid} -selgrid,1 ${ifile} ${ofile}
Remapping:
for ifile in diag.2006-09-*; do
ofile=post_process/deg.25_${ifile}
cdo -P 1 -f nc4 \
remap,r1440x720,$weights \
-setgrid,mpas:${src_grid} \
-selname,uzonal_isobaric,w_isobaric,umeridional_isobaric \
${ifile} ${ofile}
done
Kelly
Replies (2)
RE: Generating weights and remapping MPAS-limited area outputs - Added by Ralf Mueller about 2 years ago
hi Kelly!
your calls look ok and you don't provide information about what went wrong. So it's kind of hard to help you (at least for me).
You also mention limited data, but in your calls everything seems to be global. What exactly does not work for you?
cheers
ralf
RE: Generating weights and remapping MPAS-limited area outputs - Added by Kelly Núñez Ocasio about 2 years ago
Hi Ralf,
I provided code in my post that works for MPAS global outputs to make the point that the same lines of code did not work for limited-area outputs. I needed help translating this interpolation method for limited-area outputs.
Florian was able to help me with this. I am sharing below the line of code that he shared with me that was successful in generating the limited-area interpolation:
REMAP_EXTRAPOLATE=off cdo -sellonlatbox,-30,40,-40,40 -remapcon,t256grid -setgrid,mpas:mpas.grid.nc -selvar,t2m diag.2006-09-11_10.00.00.nc test_remap.nc
Note that using remap instead of remapcon is also viable and a little faster than using remapcon. For remap you must create the weights first.
Best,
Kelly