<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, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi Will</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I have seen this page. But the NSIDC sea-ice motion is in binary format, not nc. I was hoping if I can get any help on dealing with those.</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thanks</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Sourav<br>
</div>
<div>
<hr style="display:inline-block;width:98%" tabindex="-1">
</div>
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> ncl-talk <ncl-talk-bounces@ucar.edu> on behalf of ncl-talk-request@ucar.edu <ncl-talk-request@ucar.edu><br>
<b>Sent:</b> Thursday, July 26, 2018 6:00 PM<br>
<b>To:</b> ncl-talk@ucar.edu<br>
<b>Subject:</b> ncl-talk Digest, Vol 176, Issue 41</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">Send ncl-talk mailing list submissions to<br>
ncl-talk@ucar.edu<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
or, via email, send a message with subject or body 'help' to<br>
ncl-talk-request@ucar.edu<br>
<br>
You can reach the person managing the list at<br>
ncl-talk-owner@ucar.edu<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of ncl-talk digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: Reading NSIDC sea ice motion binary (.bin format) data<br>
(Will Hobbs)<br>
2. ploting SLP from a regrid (SARA CORNEJO BUENO)<br>
3. Re: problem with ESFM mapping (Michele Petrini - CITG)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 26 Jul 2018 06:07:46 +0000<br>
From: Will Hobbs <will.hobbs@utas.edu.au><br>
To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu><br>
Subject: Re: [ncl-talk] Reading NSIDC sea ice motion binary (.bin<br>
format) data<br>
Message-ID: <4C384360-A016-4130-AA13-9BC917FF72FF@utas.edu.au><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Sourav<br>
<br>
There?s an entire example page on the NCL help files on how to plot the NSIDC EASE grid.<br>
<br>
<a href="http://www.ncl.ucar.edu/Applications/ease.shtml">http://www.ncl.ucar.edu/Applications/ease.shtml</a><br>
<br>
This should resolve your cluelessness?.<br>
<br>
Will<br>
<br>
From: ncl-talk <ncl-talk-bounces@ucar.edu> on behalf of Sourav Chatterjee <srvsxc@outlook.com><br>
Date: Thursday, 26 July 2018 at 4:03 PM<br>
To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu><br>
Subject: [ncl-talk] Reading NSIDC sea ice motion binary (.bin format) data<br>
<br>
Dear NCL users<br>
<br>
Can anyone help me out with reading the NSIDC EASE-Grid Sea Ice Motion Vectors data (.bin format) in NCL.<br>
<br>
I have not used this format earlier in NCL, so clueless about it. Any help is appreciated.<br>
<br>
A sample data is attached. The data format can be found in <a href="http://nsidc.org/data/nsidc-0116/versions/3<http://nsidc.org/data/nsidc-0116/versions/3">
http://nsidc.org/data/nsidc-0116/versions/3<http://nsidc.org/data/nsidc-0116/versions/3</a>><br>
Polar Pathfinder Daily 25 km EASE-Grid Sea Ice Motion Vectors, Version 3 | National Snow and Ice Data Center<<a href="http://nsidc.org/data/nsidc-0116/versions/3">http://nsidc.org/data/nsidc-0116/versions/3</a>><br>
Polar Pathfinder Daily 25 km EASE-Grid Sea Ice Motion Vectors, Version 3<br>
nsidc.org<<a href="http://nsidc.org">http://nsidc.org</a>><br>
<br>
<br>
<br>
Thanks in advance<br>
Sourav<br>
<br>
<br>
University of Tasmania Electronic Communications Policy (December, 2014).<br>
This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained
in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180726/5682ecf3/attachment-0001.html">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180726/5682ecf3/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 26 Jul 2018 10:48:07 +0200<br>
From: SARA CORNEJO BUENO <sarcorne@ucm.es><br>
To: ncl-talk@ucar.edu<br>
Subject: [ncl-talk] ploting SLP from a regrid<br>
Message-ID:<br>
<CANq32LVv_5BDRMQDycqr1RjPcsXs_BTaZO9VdaVP9hNVhpbF8g@mail.gmail.com><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hello,<br>
<br>
I'm trying to make a sea level pressure plot from a Regrid; this is part of<br>
the code<br>
<br>
;;;;;; Set basic resources ;;;;;;<br>
res = True<br>
res@MainTitle = "Sea Level Pressure_811"<br>
time = wrf_user_getvar(a,"time",-1) ; get all time<br>
in the file<br>
ntime = dimsizes(time)<br>
do it = 0,ntime-1,1 ; TIME LOOP<br>
res@TimeLabel = time(it) ; Set Valid time<br>
to use on plots<br>
;exit()<br>
; Open workstation<br>
wks_type = "png"<br>
wks_type@wkWidth = 2000<br>
wks_type@wkHeight = 2000<br>
wks = gsn_open_wks(wks_type, "/Volumes/SimuPlots/slp_811_"<br>
+time(it)+ "")<br>
<br>
gsn_define_colormap(wks,"BlAqGrYeOrRe")<br>
;exit()<br>
<br>
;;;;;; Set plot and map resources ;;;;;;<br>
pltres = True<br>
mpres = True<br>
mpres@mpOutlineBoundarySets = "National"<br>
mpres@mpDataBaseVersion = "MediumRes"<br>
mpres@mpDataSetName = "Earth..4"<br>
mpres@mpGeophysicalLineColor = "Black"<br>
mpres@mpGeophysicalLineThicknessF = 4.0<br>
;exit()<br>
LON =a->longitude<br>
;print(lon)<br>
LAT=a->latitude<br>
;print(lat)<br>
LON = LON-360 ;as? son longitudes W<br>
;print(LON)<br>
LAT = LAT(::-1) ; overwrite: reverse using NCL syntax (para recolocar<br>
la latitud (tiene que estar S-N))<br>
;exit()<br>
if (isvar("zoomCoords")) then<br>
; lat2d = wrf_user_Agetvar(a,"latitude",0)<br>
;lon2d = wrf_user_getvar(a,"longitude",0)<br>
;exit()<br>
; Coordenadas de la zona del mapa que quiero dibujar<br>
mpres@mpLimitMode = "Corners"<br>
mpres@mpLeftCornerLatF = zoomCoords(0)<br>
mpres@mpLeftCornerLonF = zoomCoords(1)<br>
mpres@mpRightCornerLatF = zoomCoords(2)<br>
mpres@mpRightCornerLonF = zoomCoords(3)<br>
end if<br>
;exit()<br>
; Get the variables<br>
<br>
slp = wrf_user_getvar(a,"SLP",it)<br>
slp@units ="hPa"<br>
<br>
<br>
All goes right untill the following part:<br>
<br>
the problem appears at the end, on this part of the code:<br>
<br>
;;;;;; Plotting options for slp ;;;;;<br>
opts2 = res<br>
opts2@cnLevelSelectionMode = "ExplicitLevels"<br>
opts2@gsnContourLineThicknessesScale = 7.0<br>
opts2@cnLineColor = "Blue"<br>
opts2@cnLineLabelInterval = 1 ; default = 2<br>
opts2@cnLineLabelPlacementMode = "constant"<br>
opts2@cnLineLabelDensityF = 2 ;densidad de info-boxes en una linea de<br>
contorno<br>
; opts2@cnLevels = (/ 980, 984, 988, 990, 994, 998, 1000, 1002, 1004,<br>
1008, 1010, 1012, 1018, 1020, 1022 /)<br>
opts2@ContourParameters = (/ 980, 1022, 2 /)<br>
opts2@cnFillOn = False<br>
contour_slp = wrf_contour(a, wks, slp, opts2)<br>
delete(opts2)<br>
<br>
It appears these warnings:<br>
<br>
warning:Attempt to reference attribute (TITLE) which is undefined<br>
<br>
warning:Attempt to reference attribute (DX) which is undefined<br>
<br>
warning:Attempt to reference attribute (WEST-EAST_GRID_DIMENSION) which is<br>
undefined<br>
<br>
warning:Attempt to reference attribute (SOUTH-NORTH_GRID_DIMENSION) which<br>
is undefined<br>
<br>
warning:Attempt to reference attribute (BOTTOM-TOP_GRID_DIMENSION) which is<br>
undefined<br>
<br>
<br>
The code continues like this<br>
<br>
; MAKE PLOTS<br>
if (isvar("zoomCoords")) then<br>
pltres@LatLonOverlay = True<br>
end if<br>
plot = wrf_map_overlays(a, wks, (/contour_slp/), pltres, mpres)<br>
<br>
end do ; TIME LOOP<br>
;end do ; DAY LOOP<br>
end<br>
And then the error says:<br>
<br>
(0) wrf_map: Error: no MAP_PROJ attribute in input file<br>
<br>
fatal:overlay: missing value as input, can't continue<br>
<br>
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 83 in file<br>
Plot_SLP_simus.ncl,<br>
<br>
<br>
line 83 corresponds to:<br>
<br>
<br>
plot = wrf_map_overlays(a, wks, (/contour_slp/), pltres, mpres)<br>
<br>
<br>
What I could understand is that, now, I don't have LAT and LONG<br>
coordinates, but grid points. But I don't know how to call them.<br>
<br>
Thank you very much in advance.<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180726/3d2947da/attachment-0001.html">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180726/3d2947da/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 26 Jul 2018 14:00:31 +0000<br>
From: Michele Petrini - CITG <M.Petrini@tudelft.nl><br>
To: Dennis Shea <shea@ucar.edu><br>
Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu><br>
Subject: Re: [ncl-talk] problem with ESFM mapping<br>
Message-ID: <910ABA19-B9CC-4FE6-B461-40F8AB8FDDD8@tudelft.nl><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear Dennis,<br>
thanks for the reply.<br>
<br>
A fist immediate difference I saw was the output of gcc ?version. By combining alias and export path, I managed to obtain<br>
<br>
laptop-petrini:ESFM_regrid michele$ gcc --version<br>
gcc-7 (Homebrew GCC 7.3.0) 7.3.0<br>
Copyright (C) 2017 Free Software Foundation, Inc.<br>
This is free software; see the source for copying conditions. There is NO<br>
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.<br>
<br>
laptop-petrini:ESFM_regrid michele$ uname -a<br>
Darwin laptop-petrini.grs.vrlab.tudelft.nl<http://laptop-petrini.grs.vrlab.tudelft.nl> 17.5.0 Darwin Kernel Version 17.5.0: Mon Mar 5 22:24:32 PST 2018; root:xnu-4570.51.1~1/RELEASE_X86_64 x86_64<br>
<br>
<br>
But I still get the same error - I really don?t understand why.<br>
<br>
Right now, the only differences between my output and your output is the Kernel version (I have 17.5, you have 17.7) and the gcc version (I have 7.3.0, you have 6.2.0).<br>
<br>
Please note that I installed gcc-7 because with ncl 6.5.0 I was getting an error re the fortran library libgfortran.4.dylib, which is not included in gcc-6 - that?s another difference.<br>
<br>
any thoughts on that?<br>
<br>
thanks,<br>
Michele<br>
<br>
<br>
******************************************************************<br>
Michele Petrini, PhD<br>
<br>
Department of Geoscience and Remote Sensing, TUDelft<br>
Delft, The Netherlands<br>
<br>
M.Petrini@tudelft.nl<mailto:M.Petrini@tudelft.nl><br>
+39 3398367372<br>
michele.petrins (Skype alias)<br>
<br>
<br>
<br>
On 26 Jul 2018, at 03:19, Dennis Shea <shea@ucar.edu<mailto:shea@ucar.edu>> wrote:<br>
<br>
Hi Michele,<br>
<br>
I just ran the example on my MAC with no problems using 6.5.0 and had no problem.<br>
<br>
THX for including the output from your mchine (appended at the end of this email.)<br>
Perhaps there is a setup issue. Unfortunately, I am not knowledgeable enough about<br>
these types of thing. Hopefully, someone will chime in.<br>
<br>
%> ncl ESMF_regrid_1.ncl<br>
<br>
-rw-r--r-- 1 shea 2200 6840 Jul 25 20:56 ESMF_regrid_1.ncl<<a href="http://www.ncl.ucar.edu/Applications/Scripts/ESMF_regrid_1.ncl">http://www.ncl.ucar.edu/Applications/Scripts/ESMF_regrid_1.ncl</a>><br>
-rw-r--r-- 1 shea 2200 50104 Jul 25 20:58 sst.nc<<a href="http://sst.nc/">http://sst.nc/</a>><br>
-rw-r--r-- 1 shea 2200 973768 Jul 25 21:00 source_grid_file.nc<<a href="http://source_grid_file.nc/">http://source_grid_file.nc/</a>><br>
-rw-r--r-- 1 shea 2200 151908 Jul 25 21:00 destination_grid_file.nc<http://destination_grid_file.nc/><br>
-rw-r--r-- 1 shea 2200 231 Jul 25 21:00 PET0.RegridWeightGen.Log<br>
-rw-r--r-- 1 shea 2200 1644652 Jul 25 21:00 NCEP_2_Rect.nc<br>
-rw-r--r-- 1 shea 2200 9424 Jul 25 21:00 sst_regrid.nc<<a href="http://sst_regrid.nc/">http://sst_regrid.nc/</a>><br>
-rw-r--r-- 1 shea 2200 215451 Jul 25 21:00 ESMF_regrid.png<br>
<br>
I am attaching to output from<br>
<br>
%> ncdump -h sst.nc<<a href="http://sst.nc/">http://sst.nc/</a>> >! out.ncdump_sst.txt<br>
<br>
and the ESMF_regrid..png file<br>
<br>
Your MAC:<br>
<br>
laptop-petrini:Downloads michele$ uname -a<br>
Darwin Micheles-MacBook-Pro.local 17.5.0 Darwin Kernel Version 17.5.0: Mon Mar 5 22:24:32 PST 2018; root:xnu-4570.51.1~1/RELEASE_X86_64 x86_64<br>
<br>
Configured with: --prefix=/Library/Developer/CommandLineTools/usr --with-gxx-include-dir=/usr/include/c++/4.2.1<br>
Apple LLVM version 10.0.0 (clang-1000.10.25.5)<br>
Target: x86_64-apple-darwin17.5.0<br>
Thread model: posix<br>
InstalledDir: /Library/Developer/CommandLineTools/usr/bin<br>
<br>
<br>
My MAC:<br>
<br>
%> uname -a<br>
Darwin sherkin 17.7.0 Darwin Kernel Version 17.7.0: Thu Jun 21 22:53:14 PDT 2018; root:xnu-4570.71.2~1/RELEASE_X86_64 x86_64<br>
<br>
%> gcc --version<br>
gcc (GCC) 6.2.0<br>
<br>
<br>
<br>
<br>
On Wed, Jul 25, 2018 at 7:45 AM, Michele Petrini - CITG <M.Petrini@tudelft.nl<mailto:M.Petrini@tudelft.nl>> wrote:<br>
Dear ncl helpdesk,<br>
<br>
I am trying to remap ocean data from a global rectlinear grid to the CISM Greenland grid using ncl.<br>
<br>
As a fist step, I am trying to run ESMF_regrid_1.ncl<<a href="https://www.ncl.ucar.edu/Applications/Scripts/ESMF_regrid_1.ncl">https://www.ncl.ucar.edu/Applications/Scripts/ESMF_regrid_1.ncl</a>> from
<a href="https://www.ncl.ucar.edu/Applications/ESMF.shtml">https://www.ncl.ucar.edu/Applications/ESMF.shtml</a>.<br>
<br>
I get the following error<br>
<br>
sh: line 1: 33907 Abort trap: 6 ESMF_RegridWeightGen --help > /dev/null 2>&1<br>
(0) ESMF_regrid_gen_weights: could not find ESMF_RegridWeightGen executable.<br>
<br>
and I don?t understand why.<br>
<br>
I am trying to run the script using ncl version 6.5.0 (installed gcc-8) - but I had same problem with older ncl version (for which I had to load di ESFM_regrid module).<br>
<br>
What I am doing wrong?<br>
<br>
thanks,<br>
M<br>
<br>
******************************************************************<br>
Michele Petrini, PhD<br>
<br>
Department of Geoscience and Remote Sensing, TUDelft<br>
Delft, The Netherlands<br>
<br>
M.Petrini@tudelft.nl<mailto:M.Petrini@tudelft.nl><br>
+39 3398367372<br>
michele.petrins (Skype alias)<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
ncl-talk mailing list<br>
ncl-talk@ucar.edu<<a href="mailto:ncl-talk@ucar.edu">mailto:ncl-talk@ucar.edu</a>><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
<br>
<out.ncdump_sst.txt><ESMF_regrid.png><br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180726/69cf8ae1/attachment-0001.html">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180726/69cf8ae1/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
ncl-talk mailing list<br>
ncl-talk@ucar.edu<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
<br>
------------------------------<br>
<br>
End of ncl-talk Digest, Vol 176, Issue 41<br>
*****************************************<br>
</div>
</span></font></div>
</body>
</html>