<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 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: thankar@live.com</pre> </div></body>
</html>