[ncl-talk] vinth2p help

[sidong]

Hi Dennis,

It works!  Thank you so much.

Best 

Dong Si


> On 26 Mar 2017, at 1:44 AM, Dennis Shea <shea at ucar.edu> wrote:
> 
> The documentation: 
>      https://www.ncl.ucar.edu/Document/Functions/Built-in/vinth2p.shtml <https://www.ncl.ucar.edu/Document/Functions/Built-in/vinth2p.shtml>
> 
> states:
> 
> hbcofa A one-dimensional array containing the hybrid A coefficients
> hbcofb A one-dimensional array containing the hybrid B coefficients
> 
> 
>   
>    hyam = f2[:]->hyam
>    hybm = f2[:]->hybm
>    printVarSummary(hyam)   ; <==== two dimensions 
>    p0mb = 1000.0
>    
>    Z3new =vinth2p (Z3_ens, hyam,hybm, pnew, PS_ens, 1, p0mb, 1, True)
>    
>    printVarSummary (Z3new)
> 
> =======
> Change 
> 
>    hyam = f2[:]->hyam
>    hybm = f2[:]->hybm
> 
> To:
> 
>    hyam = f2[0]->hyam          ; <=== one-dimension
>    hybm = f2[0]->hybm
>    printVarSummary(hyam)
> 
> 
> 
> 
> 
> 
> On Fri, Mar 24, 2017 at 11:11 PM, sidong <sidong at cma.gov.cn <mailto:sidong at cma.gov.cn>> wrote:
> Hi all,
> 
> I am trying to interpolate CCSM4 Decadal prediction dataset from hybrid coordinates to pressure coordinates.
> I used the vinth2p function and downloaded the PS file separately. However, I got the following error. 
> 
> Any suggestions will be highly appreciated.
> 
> Thank you
> 
> Dong
> 
> ———————————ERROR———————————————————
> fatal:Number of dimensions in parameter (1) of (vinth2p) is (2), (1) dimensions were expected 
> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 34 in file 
> ———————————  Script——————————————————————————————
> ;************************************
>   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/shea_util.ncl"
>   load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> ;************************************
>   begin
> ;-------------------------------------
>    diri1 = "//Volumes/Seagate/Z3/"   ; input directory   
>    fils  = systemfunc ("ls "+diri1+"b40.dp.hd-1955.20C.1d.*.cam2.h0.Z3.195501-196412.nc <http://cam2.h0.z3.195501-196412.nc/>") ; file paths     
>    f     = addfiles (fils, "r")   ; note the "s" of addfile
>    ListSetType (f, "join")   
>    Z3    = f[:]->Z3                ; read T from all files   
>    printVarSummary (Z3)
> ;-------------------------------------
>    diri2 = "//Volumes/Seagate/PS/"   ; input directory   
>    fils2 = systemfunc ("ls "+diri2+"b40.dp.hd-1955.20C.1d.*.cam2.h0.PS.195501-196412.nc <http://cam2.h0.ps.195501-196412.nc/>") ; file paths      
>    f2    = addfiles (fils2, "r")   ; note the "s" of addfile
>    ListSetType (f2, "join")  
>    PS    = f2[:]->PS                ; read T from all files  
>    printVarSummary (PS) 
> ;-------------------------------------
>    Z3_ens = dim_avg_n_Wrap(Z3,0)
>    PS_ens = dim_avg_n_Wrap(PS,0)
>    
>    printVarSummary (Z3_ens) 
>    printVarSummary (PS_ens) 
> ;-------------------------------------
>    pnew = (/950.0, 500.0, 300/) 
>    hyam = f2[:]->hyam
>    hybm = f2[:]->hybm
>    p0mb = 1000.0
>    
>    Z3new =vinth2p (Z3_ens, hyam,hybm, pnew, PS_ens, 1, p0mb, 1, True)
>    
>    printVarSummary (Z3new)
> 
> 
> ———————————VarSummary——————————————————————
> Variable: Z3
> Type: float
> Total Size: 6900940800 bytes
>             1725235200 values
> Number of Dimensions: 5
> Dimensions and sizes:	[ncl_join | 10] x [time | 120] x [lev | 26] x [lat | 192] x [lon | 288]
> Coordinates: 
>             time: [  31..3650]
>             lev: [3.54463800000001..992.5560999999998]
>             lat: [ -90..  90]
>             lon: [   0..358.75]
> Number Of Attributes: 4
>   mdims :	1
>   units :	m
>   long_name :	Geopotential Height (above sea level)
>   cell_methods :	time: mean
> 
> Variable: PS
> Type: float
> Total Size: 265420800 bytes
>             66355200 values
> Number of Dimensions: 4
> Dimensions and sizes:	[ncl_join | 10] x [time | 120] x [lat | 192] x [lon | 288]
> Coordinates: 
>             time: [  31..3650]
>             lat: [ -90..  90]
>             lon: [   0..358.75]
> Number Of Attributes: 3
>   units :	Pa
>   long_name :	Surface pressure
>   cell_methods :	time: mean
> 
> Variable: Z3_ens
> Type: float
> Total Size: 690094080 bytes
>             172523520 values
> Number of Dimensions: 4
> Dimensions and sizes:	[time | 120] x [lev | 26] x [lat | 192] x [lon | 288]
> Coordinates: 
>             time: [  31..3650]
>             lev: [3.54463800000001..992.5560999999998]
>             lat: [ -90..  90]
>             lon: [   0..358.75]
> Number Of Attributes: 5
>   cell_methods :	time: mean
>   long_name :	Geopotential Height (above sea level)
>   units :	m
>   mdims :	1
>   average_op_ncl :	dim_avg_n over dimension(s): ncl_join
> 
> Variable: PS_ens
> Type: float
> Total Size: 26542080 bytes
>             6635520 values
> Number of Dimensions: 3
> Dimensions and sizes:	[time | 120] x [lat | 192] x [lon | 288]
> Coordinates: 
>             time: [  31..3650]
>             lat: [ -90..  90]
>             lon: [   0..358.75]
> Number Of Attributes: 4
>   cell_methods :	time: mean
>   long_name :	Surface pressure
>   units :	Pa
>   average_op_ncl :	dim_avg_n over dimension(s): ncl_join
> 
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu <mailto:ncl-talk at ucar.edu>
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
> 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170326/db2923aa/attachment.html