[ncl-talk] A problem related to int2p_n_Wrap
sunmin park
mireiyue at gmail.com
Tue Sep 19 13:15:35 MDT 2017
Got it. Thank you!
sun-
<https://mailtrack.io/> Sent with Mailtrack
<https://mailtrack.io/install?source=signature&lang=en&referral=mireiyue@gmail.com&idSignature=22>
2017-09-18 22:06 GMT-07:00 Dennis Shea <shea at ucar.edu>:
>
> You are programming as if NCL is fortran. The reason why you are getting
> multiple repetitive messages is that each iteration call to 'int2p_n_Wrap'
> is working on only one vertical profile.
>
> "warning:int2p_n: *1 *input array(s) contained all missing data. No
> interpolation performed on these arrays"
>
> ---
> Interpreted languages like NCL, Matlab, IDL, Python, ,,, work on arrays
> and dynamically created arrays for the results
>
> I think your whole set of nested do loops could be replaced with something
> like:
>
> extnew = int2p_n_Wrap(pres,cal,po,1,1)
> printVarSummary(extnew)
> printMinMax(extnew, 0)
>
> ===
> Unlike (say) fortran where you have to pre-declare array space,
> Interpreted languages will create the array space dynamically. Hence, you
> do not need the following.
>
> extnew = new(dimsizes(cal),typeof(cal))
> extnew at _FillValue=9.96921e+36
>
>
>
> On Mon, Sep 18, 2017 at 10:23 PM, sunmin park <mireiyue at gmail.com> wrote:
>
>> Thank you RIck
>>
>> I tried and it doesn't give me the warming however, all the variable
>> contains all missing values...
>>
>> Sun-
>>
>> 2017-09-18 12:34 GMT-07:00 Rick Brownrigg <brownrig at ucar.edu>:
>>
>>> 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
>>>>
>>>>
>>>
>>
>>
>> --
>>
>> Sunmin Park,
>>
>> PhD Candidate, Department of Earth Science,
>>
>> University of California Riverside
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
--
Sunmin Park,
PhD Candidate, Department of Earth Science,
University of California Riverside
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170919/99107f99/attachment.html>
More information about the ncl-talk
mailing list