[ncl-talk] How to interpolate precipitation without losing much confidence?
Dennis Shea
shea at ucar.edu
Thu Oct 1 10:23:06 MDT 2020
There is no builtin function, As suggested by Dave:
**UNTESTED** you nay have to make sure this is correct.
nm = 3
nlat24 =
mlon24 =
prc8 = prc24(1:nlat24-2:mn,1:mlon24-2:mn)
prc8 = prc24 at _FillValue
printVarSummary(prc8)
dims_prc8 = dimsizes(prc8)
nlat8 = dims_prc8(0)
mlon8 = dims_prc8(1)
print(dims8_prc8)
; This does a straight arithmetic average.
do nl=1,nlat24-2,mm
do ml=1,mlon24-2,mn
prc8(nl/mn,ml/mn) = (prc24(nl+1,ml-1)+prc24(nl+1, ml
)+prc24(nl+1,ml+1) \
+prc24(nl ,ml-1)+prc24(nl ,
ml )+prc24(nl ,ml+1)\
+prc24(nl-1,ml-1)+prc24(nl-1, ml
)+prc24(nl-1,ml+1))/9.0
end do
end do
On Thu, Oct 1, 2020 at 9:59 AM Dave Allured - NOAA Affiliate via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:
> In NCL, this kind of operation falls under regridding, not interpolation.
> I recommend one of the area averaging functions, such
> as area_conserve_remap, area_hi2lores, or ESMF_regrid.
>
> None of those functions use the grid box calculation that you specified.
> If you want to perform that exact calculation over all grid boxes falling
> in a 1/8 degree box, you will need to write your own NCL code. I suggest
> that you use the dim_avg_n function inside a double loop over all output
> points on the 1/8 degree grid. "There is more than one way to do this."
>
>
> On Thu, Oct 1, 2020 at 8:58 AM Xueyan Zhang via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Hello NCL,
>>
>> I want to interpolate precipitation data from 1/24 degree to 1/8 degree.
>> The way I wanna do it is to: calculate the average of all grid boxes
>> falling in a 1/8 degree box? Is there any function in NCL?
>>
>> Thanks!
>>
>> Xueyan
>>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201001/f89e6645/attachment.html>
More information about the ncl-talk
mailing list