<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="">Hi Tabish,<div class=""><br class=""></div><div class="">please send the replies always to ncl-talk.</div><div class=""><br class=""></div><div class="">I took a look at your script and my question is do you really know what the parts are doing especially the&nbsp;</div><div class="">resources? The part of the polymarkers you mentioned in your first email aren’t included but different resources&nbsp;</div><div class="">been set without any use.</div><div class=""><br class=""></div><div class="">Here is the cleaned but untested script which includes the missing resource settings and where the&nbsp;</div><div class="">needless resources were deleted.</div><div class=""><div class=""><br class=""></div><div class=""><font face="Courier" color="#0042aa" class="">;+++++++++++++++++++++++++++++++++++++++++++++++</font></div></div><div class=""><div class=""><font face="Courier" color="#0042aa" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</font></div><div class=""><font face="Courier" color="#0042aa" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">begin</font></div><div class=""><font face="Courier" color="#0042aa" class="">; The WRF ARW input file. &nbsp;</font></div><div class=""><font face="Courier" color="#0042aa" class="">; This needs to have a ".nc" appended, so just do it.</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; a = addfile("/home/admin/IBM-Research_intern/Emissions_data/Anthro/SEAC4RS_North_Hemisphere.nc","r")</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; emiss = a-&gt;co &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ;emissions;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; x_start = 167</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; x_end &nbsp; = 175</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; y_start = 282</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; y_end &nbsp; = 289</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; emiss_zoom = emiss(y_start:y_end,x_start:x_end)</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; emiss_zoom = emiss_zoom * 1000 ; converting into Tonnes per 10km^2 per yr&nbsp;</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">;-- open a workstation</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; type = "x11"</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; wks &nbsp;= gsn_open_wks(type,"CO_SEAC4RS_Delhi") &nbsp; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; gsn_define_colormap(wks,"BlAqGrYeOrReVi200")</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">; Set some Basic Plot options</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= True&nbsp;</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@cnFillOn &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = True</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@gsnSpreadColors &nbsp; &nbsp; &nbsp;= True</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@cnLevelSelectionMode = "ManualLevels"</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@cnMinLevelValF &nbsp; &nbsp; &nbsp; = 0</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@cnMaxLevelValF &nbsp; &nbsp; &nbsp; = 30000</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@cnLevelSpacingF &nbsp; &nbsp; &nbsp;= 300</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@cnLinesOn &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= False ;turning off contour lines and labels</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@cnFillMode &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = "RasterFill"</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@gsnDraw &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= False</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; res@gsnFrame &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = False</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">;-- contour plot</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; plot = gsn_csm_contour(wks,emiss_zoom,res) &nbsp;&nbsp;</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; &nbsp; &nbsp;</font></div><div class=""><font face="Courier" color="#0042aa" class="">;-- shapefile</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; india_shp_name &nbsp;= "/home/admin/IBM-Research_intern/Shapefiles/India/IND_adm3.shp"</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; lnres &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= True</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; lnres@gsLineColor &nbsp; &nbsp; &nbsp;= "black"</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; lnres@gsLineThicknessF = 1.5</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; india_id = gsn_add_shapefile_polylines(wks,plot,india_shp_name,lnres)</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">;-- polymarker</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; pmres &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= True</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; pmres@gsMarkerColor &nbsp; &nbsp;= "Black"</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; pmres@gsMarkerIndex &nbsp; &nbsp;= 16</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; pmres@gsMarkerSizeF &nbsp; &nbsp;= 0.01</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp;&nbsp;</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; polym_id &nbsp;= &nbsp;gsn_add_polymarker(wks, plot,28.45, 77.00,pmres)</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp;</font></div><div class=""><font face="Courier" color="#0042aa" class="">;-- create the plot</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; draw(plot) &nbsp; &nbsp; &nbsp; ; This will draw the map and the shapefile outlines.</font></div><div class=""><font face="Courier" color="#0042aa" class="">&nbsp; frame(wks) &nbsp; &nbsp; &nbsp; ; Advance the frame</font></div><div class=""><font face="Courier" color="#0042aa" class=""><br class=""></font></div><div class=""><font face="Courier" color="#0042aa" class="">end</font></div></div><div class=""><font face="Courier" color="#0042aa" class="">;+++++++++++++++++++++++++++++++++++++++++++++++</font></div><div class=""><br class=""></div><div class="">Bye,</div><div class="">Karin</div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">Am 10.07.2015 um 14:09 schrieb Tabish Ansari &lt;<a href="mailto:tabishumaransari@gmail.com" class="">tabishumaransari@gmail.com</a>&gt;:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Hi Karin<br class=""><br class=""></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Thanks for your reply. I have done as you said, but it's still not working. Below is my entire script:<br class=""><br class=""><span style="color:rgb(255,0,0)" class="">;&nbsp;&nbsp; Example script to produce plots for a WRF real-data run,<br class="">;&nbsp;&nbsp; with the ARW coordinate dynamics option.<br class="">;&nbsp;&nbsp; The script shows how to zoom into a given area<br class=""><br class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br class=""><br class="">begin<br class="">; The WRF ARW input file.&nbsp; <br class="">; This needs to have a ".nc" appended, so just do it.<br class="">&nbsp; a = addfile("/home/admin/IBM-Research_intern/Emissions_data/Anthro/SEAC4RS_North_Hemisphere.nc","r")<br class=""><br class="">; We generate plots, but what kind do we prefer?<br class="">&nbsp; type = "x11"<br class="">&nbsp; wks = gsn_open_wks(type,"CO_SEAC4RS_Delhi");;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br class="">&nbsp; gsn_define_colormap(wks,"BlAqGrYeOrReVi200")<br class="">&nbsp; ;gsn_define_colormap(wks,"hotres")<br class="">; Set some Basic Plot options<br class="">&nbsp; res = True<br class="">&nbsp; ;res@MainTitle&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "REAL-TIME WRF"<br class=""><br class="">&nbsp; pltres = True<br class="">&nbsp; mpres = True<br class="">&nbsp; mpres1 = True<br class="">res = True <br class="">res@cnFillOn = True<br class="">res@gsnSpreadColors = True<br class="">res@cnLevelSelectionMode = "ManualLevels"<br class="">res@cnMinLevelValF = 0<br class="">res@cnMaxLevelValF = 30000<br class="">res@cnLevelSpacingF = 300<br class="">res@cnLinesOn = False ;turning off contour lines and labels<br class="">res@cnFillMode = "RasterFill"<br class="">res@gsnDraw=False<br class="">res@gsnFrame=False<br class="">&nbsp;&nbsp;&nbsp; emiss = a-&gt;co ;emissions;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br class="">&nbsp;&nbsp;&nbsp; x_start = 167<br class="">&nbsp;&nbsp;&nbsp; x_end&nbsp;&nbsp; = 175<br class="">&nbsp;&nbsp;&nbsp; y_start = 282<br class="">&nbsp;&nbsp;&nbsp; y_end&nbsp;&nbsp; = 289<br class=""><br class="">&nbsp;&nbsp;&nbsp; mpres1@ZoomIn = True<br class="">&nbsp;&nbsp;&nbsp; mpres1@Xstart = x_start<br class="">&nbsp;&nbsp;&nbsp; mpres1@Ystart = y_start<br class="">&nbsp;&nbsp;&nbsp; mpres1@Xend = x_end<br class="">&nbsp;&nbsp;&nbsp; mpres1@Yend = y_end<br class=""><br class="">&nbsp;&nbsp;&nbsp; emiss_zoom = emiss(y_start:y_end,x_start:x_end)<br class="">&nbsp;&nbsp;&nbsp; emiss_zoom = emiss_zoom * 1000 ; converting into Tonnes per 10km^2 per yr <br class="">&nbsp;&nbsp; ;emiss_zoom@description = emiss@description + " - Zoom over India"<br class="">&nbsp;&nbsp;&nbsp; opts = res&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br class="">&nbsp;&nbsp;&nbsp; opts@cnFillOn = True&nbsp; <br class="">&nbsp;&nbsp;&nbsp; plot = gsn_csm_contour(wks,emiss_zoom,res)&nbsp;&nbsp;&nbsp; <br class="">&nbsp;&nbsp;&nbsp; delete(opts)<br class="">&nbsp; india_shp_name&nbsp;&nbsp;&nbsp; = "/home/admin/IBM-Research_intern/Shapefiles/India/IND_adm3.shp"<br class="">&nbsp; ;lanka_shp_name&nbsp;&nbsp;&nbsp; = "/media/ResearchWork/Shapefiles/Sri-Lanka/LKA_adm0.shp"<br class="">&nbsp; lnres&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br class="">&nbsp; lnres@gsLineColor&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "black"<br class="">&nbsp; lnres@gsLineThicknessF = 1.5<br class="">&nbsp; india_id = gsn_add_shapefile_polylines(wks,plot,india_shp_name,lnres)<br class="">&nbsp; ;lanka_id = gsn_add_shapefile_polylines(wks,plot,lanka_shp_name,lnres)<br class=""><br class="">&nbsp; str&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; unique_string("mark")<br class="">&nbsp; plot@$str$ =&nbsp; gsn_add_polymarker(wks, plot,28.45, 77.00,pmres)<br class="">&nbsp;<br class="">&nbsp; draw(plot)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; This will draw the map and the shapefile outlines.<br class="">&nbsp; frame(wks)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; Advance the frame<br class=""><br class=""><br class="">end</span><br class=""><br class=""><br class=""><br class=""></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Thanks<br class=""><br class=""></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Tabish<br class=""></div></div><div class="gmail_extra"><br clear="all" class=""><div class=""><div class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div class=""><font color="#666666" face="trebuchet ms, sans-serif" size="2" class="">Tabish U Ansari<br class="">IBM Research India Laboratory<br class=""></font></div><div class=""><font color="#666666" face="trebuchet ms, sans-serif" size="2" class="">New Delhi<br class=""></font></div></div></div></div></div></div></div>
<br class=""><div class="gmail_quote">On 10 July 2015 at 17:33, Karin Meier-Fleischer <span dir="ltr" class="">&lt;<a href="mailto:meier-fleischer@dkrz.de" target="_blank" class="">meier-fleischer@dkrz.de</a>&gt;</span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div bgcolor="#FFFFFF" text="#000000" class="">
    Hi Tabish,<br class="">
    <br class="">
    if you created a map plot than you can use the function
    gsn_add_polymarker to draw marker <br class="">
    at specific lat/lon positions:<br class="">
    <br class="">
    &nbsp; str&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; unique_string("mark")<br class="">
    &nbsp; plot@$str$ =&nbsp; gsn_add_polymarker(wks, plot, 77.00, 28.45., pmres)<br class="">
    <br class="">
    Each polymarker has to have its own id so it is a good idea to use
    the unique_string function to generate a unique id name.<br class="">
    <br class="">
    Bye,<br class="">
    Karin<br class="">
    <br class="">
    <div class="">Am 10.07.15 um 13:55 schrieb Tabish
      Ansari:<br class="">
    </div><div class=""><div class="h5">
    <blockquote type="cite" class="">
      <div dir="ltr" class="">
        <div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Hi<br class="">
          <br class="">
        </div>
        <div class="gmail_default" style="font-family:trebuchet ms,sans-serif">I wish to add dots to demarcate specific
          locations of my interest in my map. I know the lat-lon of
          those locations. I am including the following lines in my
          script:<br class="">
          <b class=""><br class="">
            &nbsp; pmres = True<br class="">
            &nbsp; pmres@gsMarkerColor = "Black"<br class="">
            &nbsp; pmres@gsMarkerIndex = 16<br class="">
            &nbsp; pmres@gsMarkerSizeF = 0.01<br class="">
            &nbsp; gsn_polymarker(wks,plot,77.00,28.45,pmres)<br class="">
            &nbsp; draw(plot)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; This will draw the map and the
            shapefile outlines.<br class="">
            &nbsp; frame(wks)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; Advance the frame<br class="">
          </b><br class="">
        </div>
        <div class="gmail_default" style="font-family:trebuchet ms,sans-serif">However, no dot is visible in the map. Also I
          am not getting any errors or warning.<br class="">
          <br class="">
        </div>
        <div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Please help me out with this<br class="">
          <br class="">
        </div>
        <div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Thanks in advance!<br class="">
          <br class="">
        </div>
        <div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Tabish<br class="">
        </div>
        <div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br clear="all" class="">
        </div>
        <div class="">
          <div class="">
            <div dir="ltr" class="">
              <div class="">
                <div dir="ltr" class="">
                  <div class="">
                    <div class=""><font face="trebuchet ms,
                        sans-serif" color="#666666" size="2" class="">Tabish U Ansari<br class="">
                        IBM Research India Laboratory<br class="">
                      </font></div>
                    <div class=""><font face="trebuchet ms,
                        sans-serif" color="#666666" size="2" class="">New Delhi</font></div></div></div></div></div></div></div></div></blockquote></div></div></div></blockquote></div></div></div></blockquote></div><div apple-content-edited="true" class=""><br class="Apple-interchange-newline">
</div>
<br class=""></div></body></html>