<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'><div>Dear NCL help,</div><div>I am trying to plot the first domain of wrf with topography (also demarcating ocean and land area) and other two domains boundry embedded over it with the below script but i am unable to draw the boundry and also to get ocean land demarcation. I look forward for the support from ncl help in this regard. I also tried&nbsp;</div><div><a href="http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/EXPERIMENTAL/wrf_overlay_doms.htm" target="_blank">http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/EXPERIMENTAL/wrf_overlay_doms.htm</a>&nbsp;this link but it always gives segmentation fault error.</div><div><div><br></div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"&nbsp;</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"&nbsp;</div><div><br></div><div>&nbsp;begin&nbsp;</div><div>&nbsp; m1 = addfile("/mnt/gpfs1/data/private/fgvy024/Build_WRF/WPS/met_em.d01.2015-11-29_06:00:00.nc", "r")</div><div>&nbsp; m2 = addfile("/mnt/gpfs1/data/private/fgvy024/Build_WRF/WPS/met_em.d02.2015-11-29_06:00:00.nc", "r")</div><div>&nbsp; m3 = addfile("/mnt/gpfs1/data/private/fgvy024/Build_WRF/WPS/met_em.d03.2015-11-29_06:00:00.nc", "r")</div><div>&nbsp; m1t = wrf_user_getvar(m1,"HGT_M",0)</div><div>&nbsp; m2t = wrf_user_getvar(m2,"HGT_M",0) &nbsp; &nbsp; ; Finer resolution</div><div>&nbsp; m3t = wrf_user_getvar(m3,"HGT_M",0) ; Finest resolution</div><div>; Read lat/lon coordinates for each dataset.</div><div>&nbsp; m1lat = wrf_user_getvar(m1,"XLAT",0)&nbsp;</div><div>&nbsp; m1lon = wrf_user_getvar(m1,"XLONG",0) &nbsp;&nbsp;</div><div>&nbsp; m2lat = wrf_user_getvar(m2,"XLAT",0)&nbsp;</div><div>&nbsp; m2lon = wrf_user_getvar(m2,"XLONG",0)</div><div>&nbsp; m3lat = wrf_user_getvar(m3,"XLAT",0)</div><div>&nbsp; m3lon = wrf_user_getvar(m3,"XLONG",0)</div><div>&nbsp; wks = gsn_open_wks("png", "ccl") &nbsp;; send graphics to PNG file</div><div><br></div><div>&nbsp; res = True&nbsp;</div><div>&nbsp; res@InitTime = False &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Do not plot time or footers</div><div>&nbsp; res@Footer = False&nbsp;</div><div>&nbsp; res@gsnDraw &nbsp; &nbsp; &nbsp; &nbsp; = False &nbsp; ; Don't draw plot (will do later)</div><div>&nbsp; res@gsnFrame &nbsp; &nbsp; &nbsp; &nbsp;= False &nbsp; ; Don't advance framce &nbsp;(will do later)</div><div>;y axis levels</div><div>&nbsp; res@tiYAxisString = "Latitude"</div><div>&nbsp; res@tiYAxisFontHeightF = 0.015</div><div>&nbsp; res@tiYAxisFontColor = "black"</div><div>&nbsp; res@tiYAxisSide = "left"</div><div>&nbsp; res@tiYAxisAngleF = 90</div><div>; x asis label</div><div>&nbsp; res@tiXAxisString = "Longitude"</div><div>&nbsp; res@tiXAxisFontHeightF = 0.015</div><div>&nbsp; res@tiXAxisFontColor = "black"&nbsp;</div><div>&nbsp; pltres = True</div><div>&nbsp; ;pltres@PanelPlot &nbsp;= True</div><div>&nbsp; mpres1 = True</div><div>&nbsp; mpres1@mpGeophysicalLineColor = "Black"</div><div>&nbsp; mpres1@mpNationalLineColor &nbsp; &nbsp;= "Black"</div><div>&nbsp; mpres1@mpGridLineColor &nbsp; &nbsp; &nbsp; &nbsp;= "Black"</div><div>&nbsp; mpres1@mpLimbLineColor &nbsp; &nbsp; &nbsp; &nbsp;= "Black"</div><div>&nbsp; mpres1@mpPerimLineColor &nbsp; &nbsp; &nbsp; = "Black"</div><div>&nbsp; mpres1@mpFillOn = False</div><div>&nbsp; mpres1@mpGridAndLimbOn &nbsp; &nbsp; &nbsp; = False &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Turn off lat/lon lines</div><div>&nbsp; mpres1@pmTickMarkDisplayMode = "Always"</div><div>&nbsp; mpres1@mpProjection &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= "LambertConformal" &nbsp;</div><div>&nbsp; mpres1@mpDataBaseVersion &nbsp; &nbsp; = "MediumRes" &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Default is LowRes</div><div>&nbsp; mpres1@mpOutlineDrawOrder &nbsp; &nbsp;= "PostDraw" &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; Draw map outlines last &nbsp; &nbsp; &nbsp;&nbsp;</div><div>&nbsp; opts_r = res &nbsp;&nbsp;</div><div>&nbsp; opts_r@cnFillOn = True &nbsp;</div><div>&nbsp; contour = wrf_contour(m1,wks,m1t,opts_r)</div><div>&nbsp; plot = wrf_map_overlays(m1,wks,contour,pltres,mpres1) ; Plot field over map background</div><div>&nbsp; overlay(plot,contour)&nbsp;</div><div>&nbsp; lnres1 = True&nbsp;</div><div>&nbsp; lnres1@gsLineThicknessF = 1.8&nbsp;</div><div>&nbsp; lnres1@gsLineColor = "white"&nbsp;</div><div>&nbsp; lnres1@gsLineLabelFontColor = "red"&nbsp;</div><div>&nbsp; lnres1@sLineLabelFontHeightF = 0.0145&nbsp;</div><div>&nbsp; &nbsp;</div><div>; Add some boxes to the map, showing the two finer-resolution</div><div>; map areas.</div><div>&nbsp; dims = dimsizes(m2t)</div><div>&nbsp; r2 &nbsp; = dimsizes(dims)</div><div>&nbsp; r21 &nbsp;= dims(r2-1)</div><div>&nbsp; r22 &nbsp;= dims(r2-2)</div><div><br></div><div>&nbsp; xbox = (/m2lon(0,0),m2lon(0,r21-1),m2lon(r22-1,r21-1), \</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;m2lon(r22-1,0),m2lon(0,0)/)</div><div>&nbsp; ybox = (/m2lat(0,0),m2lat(0,r21-1),m2lat(r22-1,r21-1), \</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;m2lat(r22-1,0),m2lat(0,0)/)</div><div><br></div><div>&nbsp; x_out = new(dimsizes(xbox),typeof(xbox))</div><div>&nbsp; y_out = new(dimsizes(ybox),typeof(ybox))</div><div><br></div><div>; Can't use gsn_polyline here, because will get curved box lines.</div><div>&nbsp; datatondc(plot, xbox, ybox, x_out, y_out)</div><div>&nbsp; gsn_polyline_ndc(wks, x_out, y_out, lnres1)</div><div><br></div><div>&nbsp; dims = dimsizes(m3t)</div><div>&nbsp; r3 &nbsp; = dimsizes(dims)</div><div>&nbsp; r31 &nbsp;= dims(r3-1)</div><div>&nbsp; r32 &nbsp;= dims(r3-2)</div><div><br></div><div>&nbsp; xbox = (/m3lon(0,0),m3lon(0,r31-1),m3lon(r32-1,r31-1), \</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;m3lon(r32-1,0),m3lon(0,0)/)</div><div>&nbsp; ybox = (/m3lat(0,0),m3lat(0,r31-1),m3lat(r32-1,r31-1), \</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;m3lat(r32-1,0),m3lat(0,0)/)</div><div><br></div><div>&nbsp; datatondc(plot, xbox, ybox, x_out, y_out)</div><div>&nbsp; gsn_polyline_ndc(wks, x_out, y_out, lnres1)</div><div>&nbsp;draw(plot)</div><div>&nbsp;frame(wks)&nbsp;</div><div>end&nbsp;</div><div>&nbsp;</div></div><div><br></div><div><br></div><br><br><i>Regards<br>Ramchandra&nbsp;<br></i><br><br><div><hr id="stopSpelling">From: rammetro@hotmail.com<br>To: ncl-talk@ucar.edu<br>Subject: help ncl script for extraction of station data<br>Date: Tue, 19 Apr 2016 16:59:55 +0200<br><br>

<style><!--
.ExternalClass .ecxhmmessage P {
padding:0px;
}

.ExternalClass body.ecxhmmessage {
font-size:12pt;
font-family:Calibri;
}

--></style>
<div dir="ltr">Dear NCL group,<div>I tried to extract hourly air temperature data from wrf out and used the following scripts. would you figure out if there is any errors in it because when i run it it stops with showing time information and then pressing q several times it can execute the script but it is too time consuming and difficult to work because i have to press q for the size of times in file. I look foward for kind response</div><div><br></div><div><br><div>&nbsp; &nbsp; ; -------------- &nbsp;LOAD FUNCTIONS AND PROCEDURES ----------------</div><div><br></div><div>&nbsp; &nbsp; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>&nbsp; &nbsp; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div>&nbsp; &nbsp; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div>&nbsp; &nbsp; load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div><br></div><div>&nbsp; &nbsp; ; -------------- &nbsp;BEGINING OF NCL SCRIPT ----------------</div><div>begin</div><div>DATADir = "/mnt/gpfs1/data/private/fgvy024/Build_WRF/WRFV3/run/"</div><div>FILES = systemfunc (" ls -1 " + DATADir + "wrfout_d03* ")</div><div>numFILES = dimsizes(FILES)</div><div>a = addfiles(FILES+".nc","r")</div><div>;times = wrf_user_list_times(a) &nbsp; &nbsp; ; get times in the file</div><div>ntimes = 4320 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; number of times in the file</div><div>temp = new(ntimes,float)&nbsp;</div><div>do it = 0, ntimes-1 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ;Loop for the time:&nbsp;</div><div>time = it</div><div>;************************************************</div><div>&nbsp;; &nbsp;- Select lon &amp; lat of the point of interest -</div><div>&nbsp;;************************************************</div><div>&nbsp; ; - The function wrf_user_ll_to_ij find the nearest grid point for a specific lat and lon</div><div><br></div><div>Latitude = 27.959167 &nbsp;;phakding</div><div>Longitude = 86.81278</div><div>res = True &nbsp; &nbsp; &nbsp;</div><div>res@returnInt = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;</div><div>point = wrf_user_ll_to_ij(a,Longitude,Latitude,res)&nbsp;</div><div>x = point(0)-1</div><div>y = point(1)-1</div><div>;tc2 = wrf_user_getvar(a,"HGT",it) &nbsp; &nbsp; ; heigth extraction</div><div>tc2 = wrf_user_getvar(a,"T2",it) &nbsp; &nbsp; &nbsp; ; temp extraction</div><div>tc2 = tc2 -273.16</div><div><br></div><div>temp(it) = tc2(y,x)&nbsp;</div><div>end do</div><div><br></div><div>&nbsp; &nbsp; npts=ntimes</div><div>&nbsp; &nbsp; fName = "pyramid_1km.txt"</div><div>&nbsp; &nbsp; data &nbsp;= new( npts, "string")</div><div>&nbsp; &nbsp; &nbsp; &nbsp;</div><div>&nbsp; &nbsp; &nbsp;</div><div>&nbsp; &nbsp; print(" &nbsp;Time &nbsp; &nbsp; &nbsp;temp &nbsp; &nbsp; &nbsp; ")</div><div>&nbsp; &nbsp; &nbsp;</div><div>&nbsp; &nbsp; do it = 0, ntimes-1 ;ntimes-1</div><div><br></div><div>&nbsp; &nbsp; &nbsp; &nbsp; print (sprintf("%5.0f",it) &nbsp; &nbsp;+" " \</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; +sprintf("%21.2f", temp(it)) +" &nbsp;" )</div><div><br></div><div>&nbsp; &nbsp; end do &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; end of time loop</div><div><br></div><div>&nbsp; &nbsp; do it = 0,ntimes-1 ;ntimes-1</div><div><br></div><div>&nbsp; &nbsp; &nbsp; &nbsp;data (it)= (sprintf("%5.0f",it) &nbsp; &nbsp;+" " \</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; +sprintf("%21.2f", temp(it)) +" &nbsp;" &nbsp;)</div><div><br></div><div>&nbsp;end do &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; end of time loop</div><div>&nbsp;asciiwrite (fName , data)</div><div>end</div><i><div><i><br></i></div>Regards<br>Ramchandra Karki</i><div><i>PhD Student, &nbsp;</i></div><div><i>Institute of Geography,&nbsp;</i><i style="font-size:12pt;">University of Hamburg</i></div><div><i style="font-size:12pt;">Bundesstrase 55, Hamburg Germany</i></div><div><i><br></i></div></div>                                               </div></div>                                               </div></body>
</html>