<div>Dear all,&nbsp;</div><div><br></div><div>I'd like to create a netCDF file with the selected variables read off WRF output data, but I can't make it.</div><div><br></div><div>Here's my code:</div><div><br></div><div><div>&nbsp; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>&nbsp; load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div><br></div><div>begin</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>in = addfile("wrfout_d03_2011-06-13_18_00_00.nc","r")</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>out = addfile("eth850_1.nc","c")</div><div>&nbsp;&nbsp;</div><div>; First get the variables we will need &nbsp;</div><div><br></div><div>&nbsp; it = 1</div><div>&nbsp; lat = wrf_user_getvar(in,"lat",it) ; Equivalent PotentialTtemperature [K]&nbsp;</div><div>&nbsp; lon = wrf_user_getvar(in, "lon",it) ; total pressure</div><div>&nbsp;&nbsp;</div><div><br></div><div>&nbsp; eth = wrf_user_getvar(in,"eth",it) ; Equivalent PotentialTtemperature [K]&nbsp;</div><div>&nbsp; p = wrf_user_getvar(in, "pressure",it) ; total pressure</div><div><br></div><div>&nbsp; ; Horizontally interpolate to pressure coordinates ("p")</div><div>&nbsp; pressure = 850. &nbsp; ; 850 hPa</div><div>&nbsp; eth_plane = wrf_user_intrp3d(eth,p,"h",pressure,0.,False)&nbsp;</div><div><br></div><div><br></div><div>&nbsp;&nbsp;</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>out-&gt;lon = lon</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>out-&gt;lat = lat</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>out-&gt;eth850 = eth_plane</div><div><span class="Apple-tab-span" style="white-space:pre">        </span></div><div>end<span id="_editor_bookmark_start_0" style="display: none; line-height: 0px;">‍</span></div></div><div><br></div><div><br></div>