<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'><p class="ecxMsoNormal"><span style="" lang="EN-US">Good
morning to all!</span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US"> </span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US">Following the
earlier conversation (to which I respond) I think I have the obligation to
clarify some issues, since these conversations are available to all users, and
they are used to solve various kinds of problems. </span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US">After doing
some more tests, I found out that the otherwise elegant and technically
flawless script that was kindly prepared for me, does not illustrates correctly
the CTTH_HEIGHT field. The problem is related to the fact that the line/column
space of the raster data does not refer to a regular lat/lon space. In simple
words the pixels of the raster data are not equal in lat/lon dimensions.</span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US"> </span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US">Regarding
more technical issues: </span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US">XGEO_LOW_RIGHT,
XGEO_UP_LEFT, YGEO_LOW_RIGHT, YGEO_UP_LEFT are not given in degrees but in meters.</span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US"> </span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US">dims =
dimsizes(CT) , lat = fspan(minLat, maxLat, dims(0)) and lon = fspan(minLon,
maxLon, dims(1)) are producing <span style=""> </span>a regular
grid (which is not correct)</span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US"> </span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US">It would be
really helpful if a script could project correctly Meteosat fields (I think safnwc
is using the Meteosat grid) but I’m not sure that such a thing exists.</span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US"> </span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US">Kind
regards</span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US"> </span></p>
<p class="ecxMsoNormal"><span style="" lang="EN-US">Sakis</span></p><pre>Athanasios F. Karagiannidis<br>Physicist, PhD Meteorology and Climatology<br>E-mail: thankar@live.com</pre><br><br><div><hr id="stopSpelling">From: thankar@live.com<br>To: haley@ucar.edu<br>Date: Sat, 21 Feb 2015 08:29:02 +0000<br>CC: ncl-talk@ucar.edu<br>Subject: Re: [ncl-talk] Visualization of SAF NWC products<br><br>
<style><!--
.ExternalClass .ecxhmmessage P {
padding:0px;
}
.ExternalClass body.ecxhmmessage {
font-size:12pt;
font-family:Calibri;
}
--></style>
<div dir="ltr">Mary, thank you for your input!<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="ecxstopSpelling">Date: Fri, 20 Feb 2015 09:46:02 -0700<br>Subject: Re: [ncl-talk] Visualization of SAF NWC products<br>From: haley@ucar.edu<br>To: thankar@live.com<br>CC: huangwei@ucar.edu; ncl-talk@ucar.edu<br><br><div dir="ltr"><div class="ecxgmail_default" style="font-size:small;">The issue with "namedcolor2rgb" is not really an issue. The code should still work. The error message is just misleading, but can be ignored.</div><div class="ecxgmail_default" style="font-size:small;"><br></div><div class="ecxgmail_default" style="font-size:small;">We have fixed the incorrect error message for the V6.3.0 release coming next month.</div><div class="ecxgmail_default" style="font-size:small;"><br></div><div class="ecxgmail_default" style="font-size:small;">--Mary</div><div class="ecxgmail_default" style="font-size:small;"><br></div></div><div class="ecxgmail_extra"><br><div class="ecxgmail_quote">On Fri, Feb 20, 2015 at 1:46 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">Wei,<br><br>Your help is more than appreciated.<br><br>I will update my NCL version.<br><br>Kind regards<span><br><br>Sakis<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><br><br></span><div><hr>From: <a href="mailto:huangwei@ucar.edu" target="_blank">huangwei@ucar.edu</a><br>Date: Thu, 19 Feb 2015 09:23:14 -0700<div><div class="h5"><br>Subject: Re: [ncl-talk] Visualization of SAF NWC products<br>To: <a href="mailto:thankar@live.com" target="_blank">thankar@live.com</a><br>CC: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br><br><div dir="ltr"><div style="font-family:verdana,sans-serif;">Sakis,</div><div style="font-family:verdana,sans-serif;"><br></div><div style="font-family:verdana,sans-serif;">Attached is script which looks more like what you see with HDFview,</div><div style="font-family:verdana,sans-serif;">but used a whole user-defined color table.</div><div style="font-family:verdana,sans-serif;"><br></div><div style="font-family:verdana,sans-serif;">"gsdtol" is a standard color-table name, in NCL 6.2.1. Your NCL of 6.1.2</div><div style="font-family:verdana,sans-serif;">is a little bit old, and I strongly suggest you update your NCL to the newest version.</div><div style="font-family:verdana,sans-serif;"><br></div><div style="font-family:verdana,sans-serif;">Compare the newly attached script and the script before,</div><div style="font-family:verdana,sans-serif;">you'll be able to figure out that the difference is the projection applied.</div><div style="font-family:verdana,sans-serif;">So you can turn on/off those options based on your requirement.<br></div><div style="font-family:verdana,sans-serif;"><br></div><div style="font-family:verdana,sans-serif;">A new png file is attached as well.</div><div style="font-family:verdana,sans-serif;"><br></div><div style="font-family:verdana,sans-serif;">Regards,</div><div style="font-family:verdana,sans-serif;"><br></div><div style="font-family:verdana,sans-serif;">Wei</div><div style="font-family:verdana,sans-serif;"><br></div><div style="font-family:verdana,sans-serif;"><br></div></div><div><br clear="all"><div><div><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: <a target="_blank">303-497-8924</a></font></div></div></div></div>
<br><div>On Thu, Feb 19, 2015 at 4:07 AM, Athanasios Karagiannidis <span dir="ltr"><<a href="mailto:thankar@live.com" target="_blank">thankar@live.com</a>></span> wrote:<br><blockquote style="border-left:1px #ccc solid;padding-left:1ex;">
<div><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!<span><br><br>Sakis<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><br><br></span><div><hr>From: <a href="mailto:huangwei@ucar.edu" target="_blank">huangwei@ucar.edu</a><br>Date: Wed, 18 Feb 2015 13:26:21 -0700<br>Subject: Re: [ncl-talk] Visualization of SAF NWC products<br>To: <a href="mailto:thankar@live.com" target="_blank">thankar@live.com</a><br>CC: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><div><div><br><br><div dir="ltr"><div style="font-family:verdana,sans-serif;"><div style="font-size:12.8000001907349px;">Sakis,</div><div style="font-size:12.8000001907349px;"><br></div><div 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 style="font-size:12.8000001907349px;"><br></div><div style="font-size:12.8000001907349px;">Thanks,</div><div style="font-size:12.8000001907349px;"><br></div><div style="font-size:12.8000001907349px;">Wei</div><div><br></div><div><br></div></div></div><div><br clear="all"><div><div><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: <a target="_blank">303-497-8924</a></font></div></div></div></div>
<br><div>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 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></div> </div></div>
</blockquote></div><br></div></div></div></div> </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>
<br>_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk</div> </div></body>
</html>