<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=""><div class="">Because you’re using ****_Wrap() the function tries to take information about the dimension it is operating on and attach that information to the new variable. </div><div class=""><br class=""></div><div class="">e.g. if you averaged over a dimension named time: the new variable will have the below attribute. </div><div class="">average_op_ncl :<span class="Apple-tab-span" style="white-space: pre;"> </span>dim_avg_n over dimension(s): time</div><div class=""><br class=""></div><div class="">You could define T_plane (and others) to have a coordinate name for the left most dimension but it doesn’t matter to the result. </div><div class=""><br class=""></div><div class=""><blockquote type="cite" class=""><div class=""><div class="">T_plane = new((/ens_size,dimy,dimx/),"float”)</div></div></blockquote>T_plane!0 = “ensembles”</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Alan</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><br class=""><div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">##############################<br class="">Alan Brammer,<div class="">Post-Doc Researcher</div><div class=""><br class=""></div><div class="">Department of Atmospheric and Environmental Sciences,<br class="">University at Albany, State University of New York, Albany, NY, 12222<div class=""><a href="mailto:abrammer@albany.edu" class="">abrammer@albany.edu</a><br class="">##############################</div></div></div></span></div></div>
</div>
<br class=""><div><blockquote type="cite" class=""><div class="">On 9 Feb 2017, at 21:56, qiuxiaobin.tj <<a href="mailto:qiuxiaobin.tj@qq.com" class="">qiuxiaobin.tj@qq.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Hi All</div><div class=""><br class=""></div><div class="">I always met a warning:</div><div class="">Dimension (0) has not been defined </div><div class="">when use dim_stddev_n_Wrap or dim_variance_n_Wrap to compute ensemble std or var on given height. I'm not sure about the result calculate from dim_stddev_n_Wrap. Does anyone has met this warning and give me some advise? Thanks a lot. </div><div class=""><br class=""></div><div class="">-------------------------------------------------------------------------------</div><div class="">Here is the output</div><div class=""><br class=""></div><div class=""><div class=""> Copyright (C) 1995-2015 - All Rights Reserved</div><div class=""> University Corporation for Atmospheric Research</div><div class=""> NCAR Command Language Version 6.3.0</div><div class=""> The use of this software is governed by a License Agreement.</div><div class=""> See <a href="http://www.ncl.ucar.edu/" class="">http://www.ncl.ucar.edu/</a> for more details.</div><div class="">(0)<span class="Apple-tab-span" style="white-space:pre"> </span>================================================================</div><div class="">(0)<span class="Apple-tab-span" style="white-space:pre"> </span>Time 2013-05-31_21:00:00 0</div><div class="">(0)<span class="Apple-tab-span" style="white-space:pre"> </span>Read in prior</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Variable: fils</div><div class="">Type: string</div><div class="">Total Size: 16 bytes</div><div class=""> 2 values</div><div class="">Number of Dimensions: 1</div><div class="">Dimensions and sizes:<span class="Apple-tab-span" style="white-space:pre"> </span>[2]</div><div class="">Coordinates: </div><div class="">(0)<span class="Apple-tab-span" style="white-space:pre"> </span>/home/qiuxb/WRF_DART/Data/20130531/CONV_12/Plot/WRFOUT/wrffcst_d01_2013-05-31_21:00:00_1</div><div class="">(1)<span class="Apple-tab-span" style="white-space:pre"> </span>/home/qiuxb/WRF_DART/Data/20130531/CONV_12/Plot/WRFOUT/wrffcst_d01_2013-05-31_21:00:00_2</div><div class="">warning:Dimension (0) has not been defined</div><div class="">(0)<span class="Apple-tab-span" style="white-space:pre"> </span>Read in post</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Variable: fils</div><div class="">Type: string</div><div class="">Total Size: 16 bytes</div><div class=""> 2 values</div><div class="">Number of Dimensions: 1</div><div class="">Dimensions and sizes:<span class="Apple-tab-span" style="white-space:pre"> </span>[2]</div><div class="">Coordinates: </div><div class="">(0)<span class="Apple-tab-span" style="white-space:pre"> </span>/home/qiuxb/WRF_DART/Data/20130531/CONV_12/Plot/WRFOUT/wrfout_d01_2013-05-31_21:00:00_1</div><div class="">(1)<span class="Apple-tab-span" style="white-space:pre"> </span>/home/qiuxb/WRF_DART/Data/20130531/CONV_12/Plot/WRFOUT/wrfout_d01_2013-05-31_21:00:00_2</div><div class="">warning:Dimension (0) has not been defined</div><div class="">(0)<span class="Apple-tab-span" style="white-space:pre"> </span>------------------temp------------------</div></div><div class=""><br class=""></div><div class="">-------------------------------------------------------------------------------</div><div class="">Here is the scripts</div><div class=""><br class=""></div><div class=""><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/wrf/WRFUserARW.ncl"</div><div class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div class=""><br class=""></div><div class="">;-----------------------------------------------------------------------------------------</div><div class=""><br class=""></div><div class="">begin</div><div class=""><br class=""></div><div class="">;-----------------------------------------------------------------------------------------</div><div class=""><br class=""></div><div class="">wrfout_dir = "/home/qiuxb/WRF_DART/Data/20130531/CONV_12/Plot/WRFOUT"</div><div class="">plot_dir = "/home/qiuxb/WRF_DART/Data/20130531/CONV_12/Plot/WRFOUT"</div><div class=""><br class=""></div><div class="">prior_post_name = (/ "fcst", "out" /)</div><div class=""><br class=""></div><div class="">dimx = 250</div><div class="">dimy = 200</div><div class="">dimz = 50</div><div class="">ens_size = 2</div><div class=""><br class=""></div><div class="">x_start = 80; 75</div><div class="">x_end = dimx - 75 ;dimx-60</div><div class="">y_start = 65 ;60</div><div class="">y_end = dimy-60 ;dimy-50</div><div class=""><br class=""></div><div class="">type = "png"</div><div class="">wks = gsn_open_wks(type, plot_dir + "/ens_spread" ) ; open a ps file</div><div class=""><br class=""></div><div class="">ntimes = 1</div><div class="">time_ini = "20130531210000"</div><div class="">time_intvl = 5 ; minute</div><div class=""><br class=""></div><div class="">z_intvl = (/ 1, 5, 10, 15, 20, 22, 24, 26, 30, 33 /)</div><div class="">height = 2000 ; msl </div><div class=""><br class=""></div><div class="">;gsn_define_colormap(wks,"BlGrYeOrReVi200")</div><div class="">gsn_define_colormap(wks,"WhBlGrYeRe")</div><div class=""><br class=""></div><div class="">z_mem = new((/dimz,dimy,dimx/),"float") ; grid point height</div><div class="">dbz_mem = new((/ens_size,dimz,dimy,dimx/),"float")</div><div class="">dbz_plane = new((/ens_size,dimy,dimx/),"float")</div><div class="">dbz_prior = new((/dimy,dimx/),"float")</div><div class="">dbz_post = new((/dimy,dimx/),"float")</div><div class="">T_mem = new((/ens_size,dimz,dimy,dimx/),"float")</div><div class="">T_plane = new((/ens_size,dimy,dimx/),"float")</div><div class="">T_prior = new((/dimy,dimx/),"float")</div><div class="">T_post = new((/dimy,dimx/),"float")</div><div class="">P_mem = new((/ens_size,dimz,dimy,dimx/),"float")</div><div class="">P_plane = new((/ens_size,dimy,dimx/),"float")</div><div class="">P_prior = new((/dimy,dimx/),"float")</div><div class="">P_post = new((/dimy,dimx/),"float")</div><div class="">Q_mem = new((/ens_size,dimz,dimy,dimx/),"float")</div><div class="">Q_plane = new((/ens_size,dimy,dimx/),"float")</div><div class="">Q_prior = new((/dimy,dimx/),"float")</div><div class="">Q_post = new((/dimy,dimx/),"float")</div><div class="">U_mem = new((/ens_size,dimz,dimy,dimx/),"float")</div><div class="">U_plane = new((/ens_size,dimy,dimx/),"float")</div><div class="">U_prior = new((/dimy,dimx/),"float")</div><div class="">U_post = new((/dimy,dimx/),"float")</div><div class="">V_mem = new((/ens_size,dimz,dimy,dimx/),"float")</div><div class="">V_plane = new((/ens_size,dimy,dimx/),"float")</div><div class="">V_prior = new((/dimy,dimx/),"float")</div><div class="">V_post = new((/dimy,dimx/),"float")</div><div class="">W_mem = new((/ens_size,dimz,dimy,dimx/),"float")</div><div class="">W_plane = new((/ens_size,dimy,dimx/),"float")</div><div class="">W_prior = new((/dimy,dimx/),"float")</div><div class="">W_post = new((/dimy,dimx/),"float")</div><div class="">;-----------------------------------------------------------------------------------------</div><div class="">time_wrf = systemfunc("da_advance_time.exe " + time_ini + " " + "0m -w")</div><div class="">fin = addfile(wrfout_dir + "/wrfout_d01_" + time_wrf + "_1" + ".nc","r")</div><div class="">z = wrf_user_getvar(fin, "z",0) ; grid point height</div><div class="">XLAT = wrf_user_getvar(fin, "XLAT",0)</div><div class="">XLONG = wrf_user_getvar(fin, "XLONG",0)</div><div class=""><br class=""></div><div class="">do tm = 0,ntimes-1</div><div class=""> time_wrf = systemfunc("da_advance_time.exe " + time_ini + " " + tostring(time_intvl * tm) + "m -w")</div><div class=""> print( "================================================================" )</div><div class=""> print( "Time " + time_wrf + " " + tm) </div><div class=""><br class=""></div><div class=""> print( "Read in prior" )</div><div class=""> fils = systemfunc("ls " + wrfout_dir + "/wrffcst_d01_" + time_wrf + "_*")</div><div class=""> print(fils)</div><div class=""> f = addfiles(fils + ".nc","r")</div><div class=""><br class=""></div><div class=""> do j=0, ens_size - 1</div><div class=""> z_mem(:,:,:) = wrf_user_getvar(f[j], "z",0)</div><div class=""> dbz_mem(j,:,:,:) = wrf_user_getvar(f[j],"REFL_10CM",0)</div><div class=""> T_mem(j,:,:,:) = wrf_user_getvar(f[j],"T",0)</div><div class=""> P_mem(j,:,:,:) = wrf_user_getvar(f[j],"P",0)</div><div class=""> Q_mem(j,:,:,:) = wrf_user_getvar(f[j],"QVAPOR",0)</div><div class=""> U_mem(j,:,:,:) = wrf_user_getvar(f[j],"ua",0)</div><div class=""> V_mem(j,:,:,:) = wrf_user_getvar(f[j],"va",0)</div><div class=""> W_mem(j,:,:,:) = wrf_user_getvar(f[j],"wa",0)</div><div class=""><br class=""></div><div class=""> dbz_plane(j,:,:) = wrf_user_intrp3d(dbz_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> T_plane(j,:,:) = wrf_user_intrp3d(T_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> P_plane(j,:,:) = wrf_user_intrp3d(P_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> Q_plane(j,:,:) = wrf_user_intrp3d(Q_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> U_plane(j,:,:) = wrf_user_intrp3d(U_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> V_plane(j,:,:) = wrf_user_intrp3d(V_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> W_plane(j,:,:) = wrf_user_intrp3d(W_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> end do</div><div class="">; dbz_prior(:,:) = dim_variance_n_Wrap(dbz_plane, 0)</div><div class=""> T_prior = dim_stddev_n_Wrap(T_plane, 0)</div><div class="">; P_prior(:,:) = dim_variance_n_Wrap(P_plane, 0)</div><div class="">; Q_prior(:,:) = dim_variance_n_Wrap(Q_plane, 0)</div><div class="">; U_prior(:,:) = dim_variance_n_Wrap(U_plane, 0)</div><div class="">; V_prior(:,:) = dim_variance_n_Wrap(V_plane, 0)</div><div class="">; W_prior(:,:) = dim_variance_n_Wrap(W_plane, 0)</div><div class=""><br class=""></div><div class=""> delete(fils)</div><div class=""> delete(f)</div><div class=""> z_mem = 0.0</div><div class=""> dbz_mem = 0.0</div><div class=""> T_mem = 0.0</div><div class=""> P_mem = 0.0</div><div class=""> Q_mem = 0.0</div><div class=""> U_mem = 0.0</div><div class=""> V_mem = 0.0</div><div class=""> W_mem = 0.0</div><div class=""> dbz_plane = 0.0</div><div class=""> T_plane = 0.0</div><div class=""> P_plane = 0.0</div><div class=""> Q_plane = 0.0</div><div class=""> U_plane = 0.0</div><div class=""> V_plane = 0.0</div><div class=""> W_plane = 0.0</div><div class=""><br class=""></div><div class=""> print( "Read in post" )</div><div class=""> fils = systemfunc("ls " + wrfout_dir + "/wrfout_d01_" + time_wrf + "_*")</div><div class=""> print(fils)</div><div class=""> f = addfiles(fils + ".nc","r")</div><div class=""><br class=""></div><div class=""> do j=0, ens_size - 1</div><div class=""> z_mem(:,:,:) = wrf_user_getvar(f[j], "z",0)</div><div class=""> dbz_mem(j,:,:,:) = wrf_user_getvar(f[j],"REFL_10CM",0)</div><div class=""> T_mem(j,:,:,:) = wrf_user_getvar(f[j],"T",0)</div><div class=""> P_mem(j,:,:,:) = wrf_user_getvar(f[j],"P",0)</div><div class=""> Q_mem(j,:,:,:) = wrf_user_getvar(f[j],"QVAPOR",0)</div><div class=""> U_mem(j,:,:,:) = wrf_user_getvar(f[j],"ua",0)</div><div class=""> V_mem(j,:,:,:) = wrf_user_getvar(f[j],"va",0)</div><div class=""> W_mem(j,:,:,:) = wrf_user_getvar(f[j],"wa",0)</div><div class=""><br class=""></div><div class=""> dbz_plane(j,:,:) = wrf_user_intrp3d(dbz_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> T_plane(j,:,:) = wrf_user_intrp3d(T_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> P_plane(j,:,:) = wrf_user_intrp3d(P_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> Q_plane(j,:,:) = wrf_user_intrp3d(Q_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> U_plane(j,:,:) = wrf_user_intrp3d(U_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> V_plane(j,:,:) = wrf_user_intrp3d(V_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> W_plane(j,:,:) = wrf_user_intrp3d(W_mem(j,:,:,:),z_mem(:,:,:),"h",height,0.,False)</div><div class=""> end do</div><div class="">; dbz_post (:,:) = dim_variance_n_Wrap(dbz_plane, 0)</div><div class=""> T_post = dim_stddev_n_Wrap(T_plane, 0)</div><div class="">; P_post (:,:) = dim_variance_n_Wrap(P_plane, 0)</div><div class="">; Q_post (:,:) = dim_variance_n_Wrap(Q_plane, 0)</div><div class="">; U_post (:,:) = dim_variance_n_Wrap(U_plane, 0)</div><div class="">; V_post (:,:) = dim_variance_n_Wrap(V_plane, 0)</div><div class="">; W_post (:,:) = dim_variance_n_Wrap(W_plane, 0)</div></div><div class=""><br class=""></div><div class="">end do</div><div class="">end</div><div class=""><br class=""></div><div class="">Best,</div><div class="">Xiaobin</div>_______________________________________________<br class="">ncl-talk mailing list<br class=""><a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a><br class="">List instructions, subscriber options, unsubscribe:<br class="">http://mailman.ucar.edu/mailman/listinfo/ncl-talk<br class=""></div></blockquote></div><br class=""></body></html>