[ncl-talk] A problem related to int2p_n_Wrap
sunmin park
mireiyue at gmail.com
Mon Sep 18 13:12:16 MDT 2017
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>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170918/8f754333/attachment.html>
More information about the ncl-talk
mailing list