<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hello,<br>
I hope you can see it now. Thank you.</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>Da:</b> Federico Buscemi - federico.buscemi@studio.unibo.it <federico.buscemi@studio.unibo.it><br>
<b>Inviato:</b> lunedì 5 luglio 2021 23:01<br>
<b>A:</b> ncl-talk@mailman.ucar.edu <ncl-talk@mailman.ucar.edu><br>
<b>Oggetto:</b> Wrong behaviour of lat/lon coordinates in an interpolated field</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<span style="color:rgb(0,0,0); font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt">Dear all,</span><br>
</div>
<div>
<div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
I am having some trouble with my script, copied and pasted below. I interpolated sea level pressure values (z), read from two of the attached text files, recorded in six locations around Venezia, Italy. I wanted to plot the interpolated field in the Venezia
Lagoon area. Even if values are well interpolated by means of "natgrid" function, it seems that they are not associated with latitude and longitude coordinates, since contours are always drawn in the same way, although I vary<span style="color:rgb(0,0,0); font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt"> lat/lon
coordinates of the map. Since I don't get errors in my output, I don't know how to solve this issue. I tried in several ways, but still can't figure it out. Any suggestion would be appreciated. Thank you.<br>
<br>
Federico Buscemi </span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<div>begin</div>
<div><br>
</div>
<div> fname1 = "filetxt/vento_direzione_201911.TXT"</div>
<div> fname2 = "filetxt/vento_velocita_201911.TXT"</div>
<div> fname3 = "filetxt/Vento_Direzione_ISPRA.txt"</div>
<div> fname4 = "filetxt/Vento_Velocità_ISPRA.txt"</div>
<div> fname5 = "filetxt/pressione_201911.TXT"</div>
<div> fname6 = "filetxt/Pressioni_ISPRA.txt"</div>
<div> lines1 = asciiread(fname1,-1,"string")</div>
<div> lines2 = asciiread(fname2,-1,"string")</div>
<div> lines3 = asciiread(fname3,-1,"string")</div>
<div> lines4 = asciiread(fname4,-1,"string")</div>
<div> lines5 = asciiread(fname5,-1,"string")</div>
<div> lines6 = asciiread(fname6,-1,"string")</div>
<div><br>
</div>
<div> delim = ";"</div>
<div> ndata = str_fields_count(lines1, delim)</div>
<div> </div>
<div> date = str_get_field(lines1, 1, delim)</div>
<div> time = str_get_field(lines1, 2, delim)</div>
<div><br>
</div>
<div> dpiatt = stringtoint(str_get_field(lines1, 3, delim))</div>
<div> dmalpor = stringtoint(str_get_field(lines1, 6, delim))</div>
<div> dchipor = stringtoint(str_get_field(lines1, 7, delim))</div>
<div> dsang = stringtoint(str_get_field(lines1, 9, delim))</div>
<div> dmeab = stringtoint(str_get_field(lines1, 11, delim))</div>
<div> dpame = stringtoint(str_get_field(lines3, 3, delim))</div>
<div> dfopo = stringtoint(str_get_field(lines3, 5, delim))</div>
<div><br>
</div>
<div> vpiatt = stringtofloat(str_get_field(lines2, 3, delim))</div>
<div> vmalpor = stringtofloat(str_get_field(lines2, 6, delim))</div>
<div> vchipor = stringtofloat(str_get_field(lines2, 7, delim))</div>
<div> vsang = stringtofloat(str_get_field(lines2, 9, delim))</div>
<div> vmeab = stringtofloat(str_get_field(lines2, 11, delim))</div>
<div> vpame = stringtofloat(str_get_field(lines4, 3, delim))</div>
<div> vfopo = stringtofloat(str_get_field(lines4, 5, delim))</div>
<div><br>
</div>
<div> ppiatt = stringtofloat(str_get_field(lines5, 3, delim))</div>
<div> ppaca = stringtofloat(str_get_field(lines5, 4, delim))</div>
<div> pmeab = stringtofloat(str_get_field(lines5, 5, delim))</div>
<div> ppame = stringtofloat(str_get_field(lines6, 3, delim))</div>
<div> pfopo = stringtofloat(str_get_field(lines6, 4, delim))</div>
<div> plime = stringtofloat(str_get_field(lines6, 5, delim))</div>
<div><br>
</div>
<div> wpiatt = wind_component(vpiatt,dpiatt,0)</div>
<div> wmalpor = wind_component(vmalpor,dmalpor,0)</div>
<div> wchipor = wind_component(vchipor,dchipor,0)</div>
<div> wsang = wind_component(vsang,dsang,0)</div>
<div> wmeab = wind_component(vmeab,dmeab,0)</div>
<div> wpame = wind_component(vpame,dpame,0)</div>
<div> wfopo = wind_component(vfopo,dfopo,0)</div>
<div><br>
</div>
<div> latpiatt = 45.3142417</div>
<div> lonpiatt = 12.5083139</div>
<div> latmalpor = 45.3398</div>
<div> lonmalpor = 12.2919667</div>
<div> latchipor = 45.2325389</div>
<div> lonchipor = 12.2805972</div>
<div> latsang = 45.4284083</div>
<div> lonsang = 12.3462639</div>
<div> latmeab = 45.235</div>
<div> lonmeab = 12.7716667</div>
<div> latpame = 45.4027694</div>
<div> lonpame = 11.8585417</div>
<div> latfopo = 44.8438333</div>
<div> lonfopo = 12.4624028</div>
<div> latpaca = 45.4364324</div>
<div> lonpaca = 12.3335352</div>
<div> latlime = 45.43015</div>
<div> lonlime = 12.38295</div>
<div><br>
</div>
<div> x = (/lonpame, lonfopo, lonpiatt, lonmeab, lonpaca, lonlime/)</div>
<div> y = (/latpame, latfopo, latpiatt, latmeab, latpaca, latlime/)</div>
<div> z = (/ppame(263), pfopo(263), ppiatt(263), pmeab(263), ppaca(263), plime(263)/)</div>
<div><br>
</div>
<div> xo = fspan(11.8000, 12.8000, 4000)</div>
<div> yo = fspan(44.8000, 45.4000, 2400)</div>
<div><br>
</div>
<div> yo@units = "degrees_north"</div>
<div> xo@units = "degrees_east"</div>
<div><br>
</div>
<div> yo!0 = "lat"</div>
<div> xo!0 = "lon"</div>
<div> </div>
<div> p_int = natgrid(x, y, z, xo, yo) </div>
<div> p_int_real = transpose(p_int)</div>
<div><br>
</div>
<div> p_int_real!0 = "lat"</div>
<div> p_int_real!1 = "lon"</div>
<div><br>
</div>
<div> p_int_real&lat = yo</div>
<div> p_int_real&lon = xo</div>
<div><br>
</div>
<div> p_int_real&lat@units = "degrees_north"</div>
<div> p_int_real&lon@units = "degrees_east"</div>
<div><br>
</div>
<div><br>
</div>
<div> wks = gsn_open_wks("pdf","station") ; send graphics to PNG file</div>
<div><br>
</div>
<div> mpres = True</div>
<div> mpres@gsnFrame = False</div>
<div> mpres@mpLimitMode = "LatLon"</div>
<div> mpres@mpMinLatF = 40.5</div>
<div> mpres@mpMinLonF = 10</div>
<div> mpres@mpMaxLatF = 48.8</div>
<div> mpres@mpMaxLonF = 15.8</div>
<div> mpres@mpDataBaseVersion = "HighRes"</div>
<div> mpres@mpOutlineBoundarySets = "AllBoundaries"</div>
<div> mpres@gsnMaximize = True</div>
<div><br>
</div>
<div> map = gsn_csm_map(wks,mpres)</div>
<div> </div>
<div> cnres = True</div>
<div> cnres@gsnAddCyclic = False</div>
<div> cnres@gsnFrame = False</div>
<div> cnres@cnFillOn = False</div>
<div> cnres@cnLinesOn = True</div>
<div> cnres@cnInfoLabelOn = False</div>
<div> cnres@cnLevelSelectionMode = "ManualLevels"</div>
<div> cnres@cnMinLevelValF = 980.0</div>
<div> cnres@cnMaxLevelValF = 1010.0</div>
<div> cnres@cnLevelSpacing = .5</div>
<div> cnres@cnFillDrawOrder = "PreDraw"</div>
<div> cnres@cnLineLabelsOn = True</div>
<div> cnres@cnLineLabelPlacementMode = "Constant"</div>
<div>; cnres@sfXArray = xo</div>
<div>; cnres@sfYArray = yo</div>
<div><br>
</div>
<div> cnres@tmXBBorderOn = False</div>
<div> cnres@tmXBOn = False</div>
<div> cnres@tmXTBorderOn = False</div>
<div> cnres@tmXTOn = False</div>
<div> cnres@tmYLBorderOn = False</div>
<div> cnres@tmYLOn = False</div>
<div> cnres@tmYRBorderOn = False</div>
<div> cnres@tmYROn = False</div>
<div><br>
</div>
<div><br>
</div>
<div> wmsetp("wbs", .045)</div>
<div><br>
</div>
<div> wmbarbmap(wks, latpiatt, lonpiatt, wpiatt(0,263), wpiatt(1,263))</div>
<div> wmbarbmap(wks, latsang, lonsang, wsang(0,263), wsang(1,263))</div>
<div> wmbarbmap(wks, latmalpor, lonmalpor, wmalpor(0,263), wmalpor(1,263))</div>
<div> wmbarbmap(wks, latchipor, lonchipor, wchipor(0,263), wchipor(1,263))</div>
<div> wmbarbmap(wks, latmeab, lonmeab, wmeab(0,263), wmeab(1,263))</div>
<div> wmbarbmap(wks, latpame, lonpame, wpame(0,263), wpame(1,263))</div>
<div> wmbarbmap(wks, latfopo, lonfopo, wfopo(0,263), wfopo(1,263))</div>
<div><br>
</div>
<div> cnp = gsn_csm_contour(wks,p_int_real,cnres)</div>
<div><br>
</div>
<div> overlay(map,cnp)</div>
<div><br>
</div>
<span>end </span><br>
</div>
</div>
</div>
</div>
</body>
</html>