[ncl-talk] A problem related to int2p_n_Wrap
Rick Brownrigg
brownrig at ucar.edu
Mon Sep 18 13:34:17 MDT 2017
Hi,
I don't know for certain what's going on, but my best guess is that out of
the 84 monthly timesteps, there 9 months that either don't have values for
the variables "pres" or "cal". The messages are just warnings -- does this
materially affect your calculations. If you need to detect these case and
do something different, you might try something like:
if (any(ismissing(pi)) .or. any(ismissing(xi))) then
; do something different.....
; perhaps something like: extnew(t,:,j,i) = extnew@
_FillValue
else
extnew(t,:,j,i) = int2p_n_Wrap (pi,xi,po,1,0)
Hope that helps...
Rick
On Mon, Sep 18, 2017 at 1:12 PM, sunmin park <mireiyue at gmail.com> wrote:
> Hello
>
> I face a problem using int2p_n_Warp.
> It gives me with following message
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays
>
> Here is my code. I appreciate any helps comments
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>
> begin
>
> diri = "/Users/spark/03_cmfprec550aer/03_output/"
> fili = "calipsodaynight.nc"
> filenc = addfile (diri+fili+".nc", "r")
> varname = "dnext"
> cali = filenc->$varname$
> printVarSummary(cali)
> cali at _FillValue = -9999
>
> cal = cali(month|:,lev|:,lat|:,lon|:) ;84*208*85*72
> printVarSummary(cal)
> printMinMax(cal,0)
>
> diri1 = "/Users/spark/03_cmfprec550aer/03_output/"
> fili1 = "pressure_hpa.nc"
> filenc1 = addfile (diri1+fili1+".nc", "r")
> varname1 = "Pre"
> pres = filenc1->$varname1$ ;84*208*85*72
> ;pres = apres(:,::-1,:,:)
> printVarSummary(pres)
> printMinMax(pres,0)
>
> x=fspan(1013.24,88.2926,208) ;max(p)=1013.24 min(P)=88.2926
> x at units = "hPa"
> printVarSummary(x)
>
> extnew = new(dimsizes(cal),typeof(cal))
> extnew at _FillValue=9.96921e+36
> printVarSummary(extnew)
> printMinMax(extnew,0)
>
> po=x
> po!0 = "p"
> po at units = "hPa"
>
> do t=0,83
> do j=0,84
> do i=0,71
> pi=pres(t,:,j,i)
> xi=cal(t,:,j,i)
> extnew(t,:,j,i) = int2p_n_Wrap (pi,xi,po,1,0)
> delete([/xi,pi/])
> end do
> end do
> end do
> extnew at _FillValue=9.96921e+36
>
> printVarSummary(extnew)
> printMinMax(extnew,0)
>
> end
>
> ;;;results
> Variable: cal
> Type: float
> Total Size: 427714560 bytes
> 106928640 values
> Number of Dimensions: 4
> Dimensions and sizes: [month | 84] x [lev | 208] x [lat | 85] x [lon | 72]
> Coordinates:
> month: [0..83]
> lev: [-0.5..12]
> lat: [-85..85]
> lon: [-180..180]
> Number Of Attributes: 1
> _FillValue : -9999
> (0) min=-5.20219 max=13.6829
>
> Variable: pres
> Type: float
> Total Size: 427714560 bytes
> 106928640 values
> Number of Dimensions: 4
> Dimensions and sizes: [month | 84] x [lev | 208] x [lat | 85] x [lon | 72]
> Coordinates:
> month: [0..83]
> lev: [12..-0.5]
> lat: [-85..85]
> lon: [-180..180]
> Number Of Attributes: 3
> _FillValue : -999
> long_name : pressure; U.S. Standard Atmosphere
> units : hPa
> (0) pressure; U.S. Standard Atmosphere: min=88.2926 max=1013.24
>
> Variable: x
> Type: float
> Total Size: 832 bytes
> 208 values
> Number of Dimensions: 1
> Dimensions and sizes: [208]
> Coordinates:
> Number Of Attributes: 1
> units : hPa
>
> Variable: extnew
> Type: float
> Total Size: 427714560 bytes
> 106928640 values
> Number of Dimensions: 4
> Dimensions and sizes: [84] x [208] x [85] x [72]
> Coordinates:
> Number Of Attributes: 1
> _FillValue : 9.96921e+36
> (0) min=9.96921e+36 max=9.96921e+36
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays
> warning:int2p_n: 1 input array(s) contained all missing data. No
> interpolation performed on these arrays...
> continuous
>
> --
>
> Sunmin Park,
>
> PhD Candidate, Department of Earth Science,
>
> University of California Riverside
>
>
>
>
>
>
>
>
> <https://mailtrack.io/> Sent with Mailtrack
> <https://mailtrack.io/install?source=signature&lang=en&referral=mireiyue@gmail.com&idSignature=22>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170918/d9ed76cc/attachment.html>
More information about the ncl-talk
mailing list