<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 </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> 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" </div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" </div><div><br></div><div> begin </div><div> m1 = addfile("/mnt/gpfs1/data/private/fgvy024/Build_WRF/WPS/met_em.d01.2015-11-29_06:00:00.nc", "r")</div><div> m2 = addfile("/mnt/gpfs1/data/private/fgvy024/Build_WRF/WPS/met_em.d02.2015-11-29_06:00:00.nc", "r")</div><div> m3 = addfile("/mnt/gpfs1/data/private/fgvy024/Build_WRF/WPS/met_em.d03.2015-11-29_06:00:00.nc", "r")</div><div> m1t = wrf_user_getvar(m1,"HGT_M",0)</div><div> m2t = wrf_user_getvar(m2,"HGT_M",0) ; Finer resolution</div><div> m3t = wrf_user_getvar(m3,"HGT_M",0) ; Finest resolution</div><div>; Read lat/lon coordinates for each dataset.</div><div> m1lat = wrf_user_getvar(m1,"XLAT",0) </div><div> m1lon = wrf_user_getvar(m1,"XLONG",0) </div><div> m2lat = wrf_user_getvar(m2,"XLAT",0) </div><div> m2lon = wrf_user_getvar(m2,"XLONG",0)</div><div> m3lat = wrf_user_getvar(m3,"XLAT",0)</div><div> m3lon = wrf_user_getvar(m3,"XLONG",0)</div><div> wks = gsn_open_wks("png", "ccl") ; send graphics to PNG file</div><div><br></div><div> res = True </div><div> res@InitTime = False ; Do not plot time or footers</div><div> res@Footer = False </div><div> res@gsnDraw = False ; Don't draw plot (will do later)</div><div> res@gsnFrame = False ; Don't advance framce (will do later)</div><div>;y axis levels</div><div> res@tiYAxisString = "Latitude"</div><div> res@tiYAxisFontHeightF = 0.015</div><div> res@tiYAxisFontColor = "black"</div><div> res@tiYAxisSide = "left"</div><div> res@tiYAxisAngleF = 90</div><div>; x asis label</div><div> res@tiXAxisString = "Longitude"</div><div> res@tiXAxisFontHeightF = 0.015</div><div> res@tiXAxisFontColor = "black" </div><div> pltres = True</div><div> ;pltres@PanelPlot = True</div><div> mpres1 = True</div><div> mpres1@mpGeophysicalLineColor = "Black"</div><div> mpres1@mpNationalLineColor = "Black"</div><div> mpres1@mpGridLineColor = "Black"</div><div> mpres1@mpLimbLineColor = "Black"</div><div> mpres1@mpPerimLineColor = "Black"</div><div> mpres1@mpFillOn = False</div><div> mpres1@mpGridAndLimbOn = False ; Turn off lat/lon lines</div><div> mpres1@pmTickMarkDisplayMode = "Always"</div><div> mpres1@mpProjection = "LambertConformal" </div><div> mpres1@mpDataBaseVersion = "MediumRes" ; Default is LowRes</div><div> mpres1@mpOutlineDrawOrder = "PostDraw" ; Draw map outlines last </div><div> opts_r = res </div><div> opts_r@cnFillOn = True </div><div> contour = wrf_contour(m1,wks,m1t,opts_r)</div><div> plot = wrf_map_overlays(m1,wks,contour,pltres,mpres1) ; Plot field over map background</div><div> overlay(plot,contour) </div><div> lnres1 = True </div><div> lnres1@gsLineThicknessF = 1.8 </div><div> lnres1@gsLineColor = "white" </div><div> lnres1@gsLineLabelFontColor = "red" </div><div> lnres1@sLineLabelFontHeightF = 0.0145 </div><div> </div><div>; Add some boxes to the map, showing the two finer-resolution</div><div>; map areas.</div><div> dims = dimsizes(m2t)</div><div> r2 = dimsizes(dims)</div><div> r21 = dims(r2-1)</div><div> r22 = dims(r2-2)</div><div><br></div><div> xbox = (/m2lon(0,0),m2lon(0,r21-1),m2lon(r22-1,r21-1), \</div><div> m2lon(r22-1,0),m2lon(0,0)/)</div><div> ybox = (/m2lat(0,0),m2lat(0,r21-1),m2lat(r22-1,r21-1), \</div><div> m2lat(r22-1,0),m2lat(0,0)/)</div><div><br></div><div> x_out = new(dimsizes(xbox),typeof(xbox))</div><div> 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> datatondc(plot, xbox, ybox, x_out, y_out)</div><div> gsn_polyline_ndc(wks, x_out, y_out, lnres1)</div><div><br></div><div> dims = dimsizes(m3t)</div><div> r3 = dimsizes(dims)</div><div> r31 = dims(r3-1)</div><div> r32 = dims(r3-2)</div><div><br></div><div> xbox = (/m3lon(0,0),m3lon(0,r31-1),m3lon(r32-1,r31-1), \</div><div> m3lon(r32-1,0),m3lon(0,0)/)</div><div> ybox = (/m3lat(0,0),m3lat(0,r31-1),m3lat(r32-1,r31-1), \</div><div> m3lat(r32-1,0),m3lat(0,0)/)</div><div><br></div><div> datatondc(plot, xbox, ybox, x_out, y_out)</div><div> gsn_polyline_ndc(wks, x_out, y_out, lnres1)</div><div> draw(plot)</div><div> frame(wks) </div><div>end </div><div> </div></div><div><br></div><div><br></div><br><br><i>Regards<br>Ramchandra <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> ; -------------- LOAD FUNCTIONS AND PROCEDURES ----------------</div><div><br></div><div> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div><br></div><div> ; -------------- 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) ; get times in the file</div><div>ntimes = 4320 ; number of times in the file</div><div>temp = new(ntimes,float) </div><div>do it = 0, ntimes-1 ;Loop for the time: </div><div>time = it</div><div>;************************************************</div><div> ; - Select lon & lat of the point of interest -</div><div> ;************************************************</div><div> ; - 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 ;phakding</div><div>Longitude = 86.81278</div><div>res = True </div><div>res@returnInt = True </div><div>point = wrf_user_ll_to_ij(a,Longitude,Latitude,res) </div><div>x = point(0)-1</div><div>y = point(1)-1</div><div>;tc2 = wrf_user_getvar(a,"HGT",it) ; heigth extraction</div><div>tc2 = wrf_user_getvar(a,"T2",it) ; temp extraction</div><div>tc2 = tc2 -273.16</div><div><br></div><div>temp(it) = tc2(y,x) </div><div>end do</div><div><br></div><div> npts=ntimes</div><div> fName = "pyramid_1km.txt"</div><div> data = new( npts, "string")</div><div> </div><div> </div><div> print(" Time temp ")</div><div> </div><div> do it = 0, ntimes-1 ;ntimes-1</div><div><br></div><div> print (sprintf("%5.0f",it) +" " \</div><div> +sprintf("%21.2f", temp(it)) +" " )</div><div><br></div><div> end do ; end of time loop</div><div><br></div><div> do it = 0,ntimes-1 ;ntimes-1</div><div><br></div><div> data (it)= (sprintf("%5.0f",it) +" " \</div><div> +sprintf("%21.2f", temp(it)) +" " )</div><div><br></div><div> end do ; end of time loop</div><div> asciiwrite (fName , data)</div><div>end</div><i><div><i><br></i></div>Regards<br>Ramchandra Karki</i><div><i>PhD Student, </i></div><div><i>Institute of Geography, </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>