[ncl-talk] vinth2p help

[Dennis Shea]

The documentation:
     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> 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") ; 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") ; 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
> 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/20170325/d8b441c1/attachment.html