<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'>Dear Wei<br><br>I'm more than thankful for your script! It is excellent in correctly reading my data and more or less what i wanted regarding the visualization.<br>I have just 3 comments:<br>1. I would really like to see the drawn field (in example the CT) even over land. How can i do this?<br>2. When i run the script in my machine i get the following warnings:<br><br><b><i>(0) Warning: namedcolor2rgb: 'gsdtol' is not a valid named color.<br>(0) Will return missing values for this color.<br></i></b><br>Is 'gsdtol' a custom made color? In this case maybe i should replace it.<br><br>3. In the example HDFView visualization i provided the grid was not "tilted" to the left (a projection issue i guess!). How should i make your visualization, so it looks like the HDFView one?<br><br>Thank you in advance!<br><br>Sakis<br><br><pre>Athanasios F. Karagiannidis<br>Physicist, PhD Meteorology and Climatology<br>E-mail: thankar@live.com</pre><br><br><div><hr id="stopSpelling">From: huangwei@ucar.edu<br>Date: Wed, 18 Feb 2015 13:26:21 -0700<br>Subject: Re: [ncl-talk] Visualization of SAF NWC products<br>To: thankar@live.com<br>CC: ncl-talk@ucar.edu<br><br><div dir="ltr"><div class="ecxgmail_default" style="font-family:verdana,sans-serif;"><div class="ecxgmail_default" style="font-size:12.8000001907349px;">Sakis,</div><div class="ecxgmail_default" style="font-size:12.8000001907349px;"><br></div><div class="ecxgmail_default" style="font-size:12.8000001907349px;">You may take a look of the attached script and pdf file to see if these are what you wanted.</div><div class="ecxgmail_default" style="font-size:12.8000001907349px;"><br></div><div class="ecxgmail_default" style="font-size:12.8000001907349px;">Thanks,</div><div class="ecxgmail_default" style="font-size:12.8000001907349px;"><br></div><div class="ecxgmail_default" style="font-size:12.8000001907349px;">Wei</div><div><br></div><div><br></div></div></div><div class="ecxgmail_extra"><br clear="all"><div><div class="ecxgmail_signature"><div dir="ltr"><font face="tahoma, sans-serif">================================================</font><div><font face="tahoma, sans-serif">1850 Table Mesa Dr.</font></div><div><font face="tahoma, sans-serif">Boulder, CO 80307</font></div><div><font face="tahoma, sans-serif">Phone: 303-497-8924</font></div></div></div></div>
<br><div class="ecxgmail_quote">On Wed, Feb 18, 2015 at 1:38 AM, Athanasios Karagiannidis <span dir="ltr"><<a href="mailto:thankar@live.com" target="_blank">thankar@live.com</a>></span> wrote:<br><blockquote class="ecxgmail_quote" style="border-left:1px #ccc solid;padding-left:1ex;">
<div><div dir="ltr">Dear sirs/madams<br><br>I want to visualize SAF NWC output data using ncl. They are included in an h5 file.<br>These data are not in a simple lat-lon grid.<br>After an earlier contact with you i managed to visualize them by creating such a grid, but i feel that the image is not accurate.<br>What i want to help me with, is to tell me weather there is a standard way or an already existing script to visualize correctly that kind of data.<br>To help you I provide the following:<br><br>1. A sample datafile.<br>2. The printVarSummary output<br>3. A screenshot of the HDFView visualization of the variable where the Variable info are shown (more or less the same info with printVarSummary).<br>4. The script i used (although, i do not think that you need that!). The lat/lon limits that i have n my script are not acurate. I used them to make the visualization.<br><br><br>---------- printVarSummary------------------------------<br><br>ncl 2> printVarSummary(CT)<br><br>Variable: CT<br>Type: ubyte<br>Total Size: 370688 bytes<br> 370688 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [DIM_000 | 512] x [DIM_001 | 724]<br>Coordinates: <br>Number Of Attributes: 11<br> CLASS : IMAGE_INDEXED<br> ID : CTTH_HEIGHTED<br> IMAGE_COLORMODEL : RGBH_HEIGHTED<br> IMAGE_SUBCLASS : IMAGE_INDEXED<br> IMAGE_VERSION : 1.0GE_INDEXED<br> N_COLS : 724<br> N_LINES : 512<br> OFFSET : -2000<br> PALETTE : ��8�3<br> PRODUCT : CTTHE_INDEXED<br> SCALING_FACTOR : 200<br><br>---------------------------------------------------------------------------------<br><br>---------------- My Script ---------------------------------------------------<br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" <br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" <br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br><br>begin<br>cdf_file = addfile("SAFNWC_MSG2_CTTH_200812121100_Greece______.h5","r")<br>CT = cdf_file->CTTH_HEIGHT(:,:)<br><br> res = True<br><br> lat = fspan(50,20,512) ; A.K. -- reverse the matrix<br> lon = fspan(10,30,724)<br> CT!0 = "lat" ;-- set name of CT dimension 0<br> CT!1 = "lon" ;-- set name of CT dimension 1<br> CT&lat = lat ;-- assign lat array to CT<br> CT&lon = lon ;-- assign lon array to CT<br> CT&lat@units = "degrees_north" ;-- set the units of the named coordinate lat<br> CT&lon@units = "degrees_east" ;-- set the units of the named coordinate lon<br><br>res@gsnAddCyclic = False<br>res@tiMainString = "Default Color" ; main title<br>res@cnFillOn = True ; turn on color fill<br> res@mpMaxLatF = 50 ; choose subregion <br> res@mpMinLatF = 20<br> res@mpMaxLonF = 30<br> res@mpMinLonF = 10<br><br>xwks = gsn_open_wks("pdf","CloudTops")<br>plot = gsn_csm_contour_map(xwks,CT,res)<br>end<br><br>--------------------------------------------------------------------------------<br><br>Thank you in advance.<br><br><br>Sakis<br><br><br><br><br><pre>Athanasios F. Karagiannidis<br>Physicist, PhD Meteorology and Climatology<br>E-mail: <a href="mailto:thankar@live.com" target="_blank">thankar@live.com</a></pre> </div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div> </div></body>
</html>