<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><span style="font-size: 14px;" class="">Hi all,</span><div class=""><span style="font-size: 14px;" class=""><br class=""></span></div><div class=""><span style="font-size: 14px;" class="">I am trying to interpolate CCSM4 Decadal prediction dataset from hybrid coordinates to pressure coordinates.</span></div><div class=""><span style="font-size: 14px;" class="">I used the vinth2p function and downloaded the PS file separately. However, I got the following error. </span></div><div class=""><span style="font-size: 14px;" class=""><br class=""></span></div><div class=""><span style="font-size: 14px;" class="">Any suggestions will be highly appreciated.</span></div><div class=""><span style="font-size: 14px;" class=""><br class=""></span></div><div class=""><span style="font-size: 14px;" class="">Thank you</span></div><div class=""><span style="font-size: 14px;" class=""><br class=""></span></div><div class=""><span style="font-size: 14px;" class="">Dong</span></div><div class=""><br class=""></div><div class="">———————————<font color="#ff2600" class="">ERROR</font>———————————————————</div><div class="">fatal:Number of dimensions in parameter (1) of (vinth2p) is (2), (1) dimensions were expected </div><div class="">fatal:["Execute.c":8575]:Execute: Error occurred at or near line 34 in file </div><div class="">——————————— <font color="#0433ff" class=""><b class=""> Script</b></font>——————————————————————————————</div><div class=""><div class="">;************************************</div><div class=""> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div class=""> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div class=""> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"</div><div class=""> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div class="">;************************************</div><div class=""> begin</div><div class="">;-------------------------------------</div><div class=""> diri1 = "//Volumes/Seagate/Z3/" ; input directory </div><div class=""> fils = systemfunc ("ls "+diri1+"b40.dp.hd-1955.20C.1d.*.cam2.h0.Z3.195501-196412.nc") ; file paths </div><div class=""> f = addfiles (fils, "r") ; note the "s" of addfile</div><div class=""> ListSetType (f, "join") </div><div class=""> Z3 = f[:]->Z3 ; read T from all files </div><div class=""> printVarSummary (Z3)</div><div class="">;-------------------------------------</div><div class=""> diri2 = "//Volumes/Seagate/PS/" ; input directory </div><div class=""> fils2 = systemfunc ("ls "+diri2+"b40.dp.hd-1955.20C.1d.*.cam2.h0.PS.195501-196412.nc") ; file paths </div><div class=""> f2 = addfiles (fils2, "r") ; note the "s" of addfile</div><div class=""> ListSetType (f2, "join") </div><div class=""> PS = f2[:]->PS ; read T from all files </div><div class=""> printVarSummary (PS) </div><div class="">;-------------------------------------</div><div class=""> Z3_ens = dim_avg_n_Wrap(Z3,0)</div><div class=""> PS_ens = dim_avg_n_Wrap(PS,0)</div><div class=""> </div><div class=""> printVarSummary (Z3_ens) </div><div class=""> printVarSummary (PS_ens) </div><div class="">;-------------------------------------</div><div class=""> pnew = (/950.0, 500.0, 300/) </div><div class=""> hyam = f2[:]->hyam</div><div class=""> hybm = f2[:]->hybm</div><div class=""> p0mb = 1000.0</div><div class=""> </div><div class=""> Z3new =vinth2p (Z3_ens, hyam,hybm, pnew, PS_ens, 1, p0mb, 1, True)</div><div class=""> </div><div class=""> printVarSummary (Z3new)</div></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">———————————<font color="#ff40ff" class="">VarSummary</font>——————————————————————</div><div class=""><div class="">Variable: Z3</div><div class="">Type: float</div><div class="">Total Size: 6900940800 bytes</div><div class=""> 1725235200 values</div><div class="">Number of Dimensions: 5</div><div class="">Dimensions and sizes:<span class="Apple-tab-span" style="white-space:pre"> </span>[ncl_join | 10] x [time | 120] x [lev | 26] x [lat | 192] x [lon | 288]</div><div class="">Coordinates: </div><div class=""> time: [ 31..3650]</div><div class=""> lev: [3.54463800000001..992.5560999999998]</div><div class=""> lat: [ -90.. 90]</div><div class=""> lon: [ 0..358.75]</div><div class="">Number Of Attributes: 4</div><div class=""> mdims :<span class="Apple-tab-span" style="white-space:pre"> </span>1</div><div class=""> units :<span class="Apple-tab-span" style="white-space:pre"> </span>m</div><div class=""> long_name :<span class="Apple-tab-span" style="white-space:pre"> </span>Geopotential Height (above sea level)</div><div class=""> cell_methods :<span class="Apple-tab-span" style="white-space:pre"> </span>time: mean</div><div class=""><br class=""></div><div class="">Variable: PS</div><div class="">Type: float</div><div class="">Total Size: 265420800 bytes</div><div class=""> 66355200 values</div><div class="">Number of Dimensions: 4</div><div class="">Dimensions and sizes:<span class="Apple-tab-span" style="white-space:pre"> </span>[ncl_join | 10] x [time | 120] x [lat | 192] x [lon | 288]</div><div class="">Coordinates: </div><div class=""> time: [ 31..3650]</div><div class=""> lat: [ -90.. 90]</div><div class=""> lon: [ 0..358.75]</div><div class="">Number Of Attributes: 3</div><div class=""> units :<span class="Apple-tab-span" style="white-space:pre"> </span>Pa</div><div class=""> long_name :<span class="Apple-tab-span" style="white-space:pre"> </span>Surface pressure</div><div class=""> cell_methods :<span class="Apple-tab-span" style="white-space:pre"> </span>time: mean</div><div class=""><br class=""></div><div class="">Variable: Z3_ens</div><div class="">Type: float</div><div class="">Total Size: 690094080 bytes</div><div class=""> 172523520 values</div><div class="">Number of Dimensions: 4</div><div class="">Dimensions and sizes:<span class="Apple-tab-span" style="white-space:pre"> </span>[time | 120] x [lev | 26] x [lat | 192] x [lon | 288]</div><div class="">Coordinates: </div><div class=""> time: [ 31..3650]</div><div class=""> lev: [3.54463800000001..992.5560999999998]</div><div class=""> lat: [ -90.. 90]</div><div class=""> lon: [ 0..358.75]</div><div class="">Number Of Attributes: 5</div><div class=""> cell_methods :<span class="Apple-tab-span" style="white-space:pre"> </span>time: mean</div><div class=""> long_name :<span class="Apple-tab-span" style="white-space:pre"> </span>Geopotential Height (above sea level)</div><div class=""> units :<span class="Apple-tab-span" style="white-space:pre"> </span>m</div><div class=""> mdims :<span class="Apple-tab-span" style="white-space:pre"> </span>1</div><div class=""> average_op_ncl :<span class="Apple-tab-span" style="white-space:pre"> </span>dim_avg_n over dimension(s): ncl_join</div><div class=""><br class=""></div><div class="">Variable: PS_ens</div><div class="">Type: float</div><div class="">Total Size: 26542080 bytes</div><div class=""> 6635520 values</div><div class="">Number of Dimensions: 3</div><div class="">Dimensions and sizes:<span class="Apple-tab-span" style="white-space:pre"> </span>[time | 120] x [lat | 192] x [lon | 288]</div><div class="">Coordinates: </div><div class=""> time: [ 31..3650]</div><div class=""> lat: [ -90.. 90]</div><div class=""> lon: [ 0..358.75]</div><div class="">Number Of Attributes: 4</div><div class=""> cell_methods :<span class="Apple-tab-span" style="white-space:pre"> </span>time: mean</div><div class=""> long_name :<span class="Apple-tab-span" style="white-space:pre"> </span>Surface pressure</div><div class=""> units :<span class="Apple-tab-span" style="white-space:pre"> </span>Pa</div><div class=""> average_op_ncl :<span class="Apple-tab-span" style="white-space:pre"> </span>dim_avg_n over dimension(s): ncl_join</div></div></body></html>