<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi, Jeff --&nbsp;<div>Perhaps the coordinates for RAINNC are different in m3 and one of the other members in ensemble set 2. Can you send the output from ncdump -h?&nbsp;</div><div>--Jennifer</div><div><br></div><div><br><div><br><div><div>On Apr 22, 2013, at 7:16 PM, Jeff Duda wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr"><div><div><div><div>Hello,<br>I have two WRF ensembles (both 7 members).&nbsp; In one ensemble, only the value of a constant in the dynamics section varies among the members (ensemble 1).&nbsp; In the other ensemble, the choice of microphysics, PBL, and land surface models vary (ensemble 2).&nbsp; All simulations are on the exact same grid with the same history dump interval.&nbsp; The control files I'm using to view the ensembles are shown below.<br>
<br>FOR ENSEMBLE 1<br>---------------------------------<br>dset /scratch/jdduda/ensemble_1/%e/wrfout_d01_%y4-%m2-%d2_%h2:00:00<br>dtype netcdf<br>undef -9999.0<br>options template<br>PDEF 599 549 lcc 25.7907 -108.7018 1 1 38.5 38.5 -97 4000 4000<br>
XDEF 599 linear -112.2245 0.050833<br>YDEF 549 linear 25.79018 0.036935<br>ZDEF 39 linear 1 1<br>TDEF 25 linear 12Z13APR2012 1hr<br>EDEF 7 names m1 m2 m3 m4 m5 m6 m7<br>VARS 5<br>P=&gt;p 39 t,z,y,x perturbation pressure (Pa)<br>
PB=&gt;pb 39 t,z,y,x base state pressure (Pa)<br>T=&gt;temp 39 t,z,y,x theta-290 K (K)<br>RT_TENDF_STOCH=&gt;utend 39 t,z,y,x temperature perturbation from SKEB (K)<br>RAINNC=&gt;rainnc 1 t,y,x precipitation (mm)<br>ENDVARS<br>
<br></div>FOR ENSEMBLE 2<br>--------------------------------<br>dset /scratch/jdduda/ensemble_2/%e/wrfout_d01_%y4-%m2-%d2_%h2:00:00<br>dtype netcdf<br>undef -9999.0<br>options template<br>PDEF 599 549 lcc 25.7907 -108.7018 1 1 38.5 38.5 -97 4000 4000<br>
XDEF 599 linear -112.2245 0.050833<br>YDEF 549 linear 25.79018 0.036935<br>ZDEF 39 linear 1 1<br>TDEF 25 linear 12Z13APR2012 1hr<br>EDEF 7 names m1 m2 m3 m4 m5 m6 m7<br>VARS 5<br>P=&gt;p 39 t,z,y,x perturbation pressure (Pa)<br>
PB=&gt;pb 39 t,z,y,x base state pressure (Pa)<br>T=&gt;temp 39 t,z,y,x theta-290 K (K)<br>RT_TENDF_STOCH=&gt;utend 39 t,z,y,x temperature perturbation from SKEB (K)<br>RAINNC=&gt;rainnc 1 t,y,x precipitation (mm)<br>ENDVARS<br>
<br></div>The difference between them is essentially one character in the dset entry ("1" vs. "2")<br><br></div>When I attempt to run a script to plot a time series of ensemble mean and spread of the variable RAINNC, it works fine for ensemble 1.&nbsp; However, for ensemble 2, I get this error:<br>
<br>NetCDF Error (gancrow, nc_get_vara_double): NetCDF: Index exceeds dimension bound<br><br></div>Looking further, I see that the history files indeed differ between the ensemble members of ensemble 2.&nbsp; They differ in the number and type of arrays contained.&nbsp; However, each member has each of the 5 variables indicated in the control file.&nbsp; The trouble seems to especially come from member m3.&nbsp; Despite the fact that member m3 differs from the other members in a similar way as the other members differ from each other, I cannot see why Grads would be giving me this problem.&nbsp; I can use ncview to view the RAINNC variable for m3 as well as for other members alone.&nbsp; Why am I having this problem when using the ensemble dimension to template in a control file?<br>
<br>Jeff Duda<br clear="all"><div><div><div><div><div><br>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
</div></div></div></div></div></div>
_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>http://gradsusr.org/mailman/listinfo/gradsusr<br></blockquote></div><br><div apple-content-edited="true">
<span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div>--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a></div><div><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></span></span></span>
</div>
<br></div></div></body></html>