From G.Monte at isac.cnr.it Tue Oct 3 03:19:52 2023 From: G.Monte at isac.cnr.it (G.Monte at isac.cnr.it) Date: Tue, 3 Oct 2023 11:19:52 +0200 Subject: [ncl-talk] Segmentation fault error Message-ID: Dears all, I get a segmentation fault error when i try to read a HDF5 file with the command: f=addfile("myfile.HDF5","r") I'm using NCL 6.6.2 installed with conda. Any suggestions, please? Thank you, Giulio Monte From setareh.rahimi at gmail.com Sat Oct 14 09:39:01 2023 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Sat, 14 Oct 2023 19:09:01 +0330 Subject: [ncl-talk] Area averaging Message-ID: Hi all, I have a question about area averaging in NCL. I have my WRF-Chem output file and would like to calculate average of aerosol optical depth (AOD) For an area at let's say 06:00 UTC. So I wonder if this would be possible through NCL? Thanks for your advise in advance. Best wishes, -------------- next part -------------- An HTML attachment was scrubbed... URL: From dave.allured at noaa.gov Sat Oct 14 12:46:10 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Sat, 14 Oct 2023 12:46:10 -0600 Subject: [ncl-talk] Area averaging In-Reply-To: References: Message-ID: Yes it is possible. If the grid has 1-D coordinates, then see the family of *wgt_areaave* functions. If 2-D, then more work is needed. Also, what kind of region definition do you need? Bounding box, polygon, shape file, list of grid points, or something else? On Sat, Oct 14, 2023 at 9:39?AM Setareh Rahimi via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Hi all, > > I have a question about area averaging in NCL. I have my WRF-Chem output > file and would like to calculate average of aerosol optical depth (AOD) > For an area at let's say 06:00 UTC. So I wonder if this would be possible > through NCL? > Thanks for your advise in advance. > > Best wishes, > -------------- next part -------------- An HTML attachment was scrubbed... URL: From leslie.linglai at gmail.com Sun Oct 15 04:19:06 2023 From: leslie.linglai at gmail.com (lzl lai) Date: Sun, 15 Oct 2023 18:19:06 +0800 Subject: [ncl-talk] Question about map and shapfile Message-ID: Hi I have a problem for deleting those thin lines (the red arrow pointed area) . Could you please help me to remove it? Thanks a lot! I use "ETOPO_2022_v1_60s_N90W180_bed.nc" ,for drawing topography, and "gadm41_TWN_0.shp" , for drawing coastline in Taiwan my scrip is below : ======================================== ;---Read data ; topo_file = "../topo/ETOPO_2022_v1_60s_N90W180_geoid.nc" topo_file = "../topo/ETOPO_2022_v1_60s_N90W180_bed.nc" f = addfile(topo_file,"r") elev = f->z ; ; Set all values below -100 to missing, hence removing all ; the ocean elevation values. The ocean will be filled in ; a light blue (see mpOceanFillColor below). ; elev = where(elev.lt.-100.,elev at _FillValue,elev) ;;file --- 2 shpdir = "/home/lzl/data2/lzl/CWBtrend/topo/TW/" filename ="gadm41_TWN_0.shp" wks = gsn_open_wks("pdf","./pic/TW_topo") ; send graphics to PNG file gsn_define_colormap(wks,"WhiteBlueGreenYellowRed") res = True res at gsnMaximize = True ; maximize plot in frame res at cnFillOn = True ; turn on contour fill ; res at cnFillPalette = cmap(2:,:) ; set color map res at cnLinesOn = False ; turn off contour lines res at cnLineLabelsOn = False ; turn off line labels res at cnInfoLabelOn = False ; turn off info label res at lbBoxLinesOn = False ; turn off labelbar box lines res at cnFillMode = "MeshFill" ; for faster draw ;---Pick "nice" contour levels ; mnmxint = nice_mnmxintvl( min(elev), max(elev), 18, False) mnmxint = nice_mnmxintvl( min(elev), max(elev), 5, False) printVarSummary(mnmxint) print(mnmxint(1)) print(mnmxint(2)) res at cnLevelSelectionMode = "ManualLevels" res at cnMinLevelValF = 100. res at cnMaxLevelValF = mnmxint(1) res at cnLevelSpacingF = mnmxint(2)/20. ; Increase the number of levels ; by choosing a smaller spacing. ;higher number mean higher resol res at gsnAddCyclic = False ; don't add longitude cyclic point ;---Zoom in on map res at mpLimitMode = "LatLon" res at mpMinLatF = 21.5 res at mpMaxLatF = 25.5 res at mpMinLonF = 119.5 res at mpMaxLonF = 122.5 ; res at mpCenterLonF = (res at mpMinLonF + res at mpMaxLonF) / 2. res at mpFillOn = True res at mpInlandWaterFillColor="Background" res at mpOceanFillColor = "transparent" res at mpLandFillColor = "transparent" res at mpFillDrawOrder = "PostDraw" res at mpDataBaseVersion = "MediumRes" res at mpOutlineOn = False res at pmTickMarkDisplayMode = "Always" res at mpGeophysicalLineThicknessF = 2 plot = gsn_csm_contour_map(wks,elev,res) nres = True lnres at gsLineColor = "black" lnres at gsLineThicknessF = 3.0 dum = gsn_add_shapefile_polylines(wks,plot,shpdir+filename,lnres) draw(plot) ; This will draw the map and the shapefile outlines. frame(wks) -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Q.docx Type: application/vnd.openxmlformats-officedocument.wordprocessingml.document Size: 111995 bytes Desc: not available URL: From dave.allured at noaa.gov Sun Oct 15 06:50:01 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Sun, 15 Oct 2023 06:50:01 -0600 Subject: [ncl-talk] Question about map and shapfile In-Reply-To: References: Message-ID: It looks like NCL example topo_4 does what you want. Try carefully changing this example for your region of interest. Those thin lines are the "fill value perimeter line" at the -100 meter elevation. You can experiment with resources that control or modify this feature. Look at the cnMissingValPerim* resources. Try res at cnMissingValPerimOn = False. Also look at your color map and color fill settings. You might try setting the ocean fill colors to white or light blue, instead of transparent. On Sun, Oct 15, 2023 at 4:21?AM lzl lai via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Hi > > I have a problem for deleting those thin lines (the red arrow pointed area) > . Could you please help me to remove it? Thanks a lot! > > I use "ETOPO_2022_v1_60s_N90W180_bed.nc" ,for drawing topography, > and "gadm41_TWN_0.shp" , for drawing coastline in Taiwan > > my scrip is below : > ======================================== > ;---Read data > ; topo_file = "../topo/ETOPO_2022_v1_60s_N90W180_geoid.nc" > topo_file = "../topo/ETOPO_2022_v1_60s_N90W180_bed.nc" > f = addfile(topo_file,"r") > elev = f->z > ; > ; Set all values below -100 to missing, hence removing all > ; the ocean elevation values. The ocean will be filled in > ; a light blue (see mpOceanFillColor below). > ; > elev = where(elev.lt.-100.,elev at _FillValue,elev) > > ;;file --- 2 > shpdir = "/home/lzl/data2/lzl/CWBtrend/topo/TW/" > filename ="gadm41_TWN_0.shp" > > wks = gsn_open_wks("pdf","./pic/TW_topo") ; send graphics to > PNG file > gsn_define_colormap(wks,"WhiteBlueGreenYellowRed") > > res = True > res at gsnMaximize = True ; maximize plot in frame > > res at cnFillOn = True ; turn on contour fill > ; res at cnFillPalette = cmap(2:,:) ; set color map > res at cnLinesOn = False ; turn off contour lines > res at cnLineLabelsOn = False ; turn off line labels > res at cnInfoLabelOn = False ; turn off info label > res at lbBoxLinesOn = False ; turn off labelbar box lines > > res at cnFillMode = "MeshFill" ; for faster draw > > ;---Pick "nice" contour levels > ; mnmxint = nice_mnmxintvl( min(elev), max(elev), 18, False) > mnmxint = nice_mnmxintvl( min(elev), max(elev), 5, False) > > printVarSummary(mnmxint) > print(mnmxint(1)) > print(mnmxint(2)) > res at cnLevelSelectionMode = "ManualLevels" > res at cnMinLevelValF = 100. > res at cnMaxLevelValF = mnmxint(1) > res at cnLevelSpacingF = mnmxint(2)/20. ; Increase the number of > levels > ; by choosing a smaller > spacing. > ;higher number mean higher > resol > > res at gsnAddCyclic = False ; don't add longitude cyclic > point > > ;---Zoom in on map > res at mpLimitMode = "LatLon" > res at mpMinLatF = 21.5 > res at mpMaxLatF = 25.5 > res at mpMinLonF = 119.5 > res at mpMaxLonF = 122.5 > ; res at mpCenterLonF = (res at mpMinLonF + res at mpMaxLonF) / 2. > > res at mpFillOn = True > > res at mpInlandWaterFillColor="Background" > res at mpOceanFillColor = "transparent" > res at mpLandFillColor = "transparent" > > res at mpFillDrawOrder = "PostDraw" > res at mpDataBaseVersion = "MediumRes" > res at mpOutlineOn = False > res at pmTickMarkDisplayMode = "Always" > res at mpGeophysicalLineThicknessF = 2 > plot = gsn_csm_contour_map(wks,elev,res) > > nres = True > lnres at gsLineColor = "black" > lnres at gsLineThicknessF = 3.0 > dum = gsn_add_shapefile_polylines(wks,plot,shpdir+filename,lnres) > > draw(plot) ; This will draw the map and the shapefile outlines. > frame(wks) > -------------- next part -------------- An HTML attachment was scrubbed... URL: From leslie.linglai at gmail.com Mon Oct 16 09:33:26 2023 From: leslie.linglai at gmail.com (lzl lai) Date: Mon, 16 Oct 2023 23:33:26 +0800 Subject: [ncl-talk] map question Message-ID: I have a problem deleting those thin lines in my map (the red arrow pointed area) . Could you please help me to remove them? Thanks a lot! I use "ETOPO_2022_v1_60s_N90W180_bed.nc" ,for drawing topography, and "gadm41_TWN_0.shp" , for drawing coastline in Taiwan and also set res at cnMissingValPerimOn = False my scrip is below : ======================================== ;---Read data ; topo_file = "../topo/ETOPO_2022_v1_60s_N90W180_geoid.nc" topo_file = "../topo/ETOPO_2022_v1_60s_N90W180_bed.nc" f = addfile(topo_file,"r") elev = f->z ; ; Set all values below -100 to missing, hence removing all ; the ocean elevation values. The ocean will be filled in ; a light blue (see mpOceanFillColor below). ; elev = where(elev.lt.-100.,elev at _FillValue,elev) ;;file --- 2 shpdir = "/home/lzl/data2/lzl/CWBtrend/topo/TW/" filename ="gadm41_TWN_0.shp" ;----plot----------------------------------------------------- wks = gsn_open_wks("pdf","./pic/TW_topo") ; send graphics to PNG file gsn_define_colormap(wks,"WhiteBlueGreenYellowRed") res = True res at gsnMaximize = True ; maximize plot in frame res at cnFillOn = True ; turn on contour fill res at cnLinesOn = False ; turn off contour lines res at cnLineLabelsOn = False ; turn off line labels res at cnInfoLabelOn = False ; turn off info label res at lbBoxLinesOn = False ; turn off labelbar box lines res at cnFillMode = "MeshFill" ; for faster draw res at cnLevelSelectionMode = "ManualLevels" res at cnMinLevelValF = 100. res at cnMaxLevelValF = 3600. res at cnLevelSpacingF = 200. res at cnMissingValPerimOn = False res at cnMissingValPerimGridBoundOn =False res at gsnAddCyclic = False ; don't add longitude cyclic point ;---color bar---- res at lbAutoManage = False ; we control label bar res at pmLabelBarDisplayMode = "Always" ; turns on label bar res at lbOrientation = "Horizontal" ; ncl default is vertical res at pmLabelBarSide = "Bottom" ; default is right res at lbLabelStride = 1 ; skip every other label res at pmLabelBarWidthF = 0.5 ; default is shorter res at pmLabelBarHeightF = 0.1 ; default is taller res at lbLabelFontHeightF = .018 ; default is HUGE res at lbPerimOn = False ; default has box res at lbBoxLinesOn = True ; turn on labelbar box lines res at lbBoxEndCapStyle = "TriangleBothEnds" ;---Zoom in on map res at mpLimitMode = "LatLon" res at mpMinLatF = 21.5 res at mpMaxLatF = 25.5 res at mpMinLonF = 119.5 res at mpMaxLonF = 122.5 res at mpFillOn = True res at mpGeophysicalLineColor = "transparent" ; color of cont. outlines res at mpOceanFillColor = "transparent" res at mpLandFillColor = "transparent" res at mpFillDrawOrder = "PostDraw" res at mpDataBaseVersion = "MediumRes" res at mpOutlineOn = False res at pmTickMarkDisplayMode = "Always" res at mpGeophysicalLineThicknessF = 2 plot = gsn_csm_contour_map(wks,elev,res) nres = True lnres at gsLineColor = "black" lnres at gsLineThicknessF = 2.0 dum = gsn_add_shapefile_polylines(wks,plot,shpdir+filename,lnres) ;---Generate random lat/lon locations and dummy text strings lont = (/121.513,120.684,120.197,120.746,121.613,121.154/) latt = (/25.049,24.12,23.07,22.013,24.05,22.765/) lons = (/121.513,120.684,120.197,120.746,121.613,121.154/) lats = (/25.037,24.145,22.995,22.003,23.975,22.752/) text = (/"Taipei","Taichung","Tainan","Hengchun","Hualian","Taitung"/) ;---Attach some dummy text strings txres = True txres at txFontHeightF = 0.02 txres at txJust = "BottomCenter" text_id = gsn_add_text(wks,plot,text,lont,latt,txres) ; mkres = True mkres at gsMarkerIndex = 16 ; Filled circle mkres at gsMarkerSizeF = 0.01 locations = gsn_add_polymarker(wks,plot,lons,lats,mkres) ;==================================================== draw(plot) ; This will draw the map and the shapefile outlines. frame(wks) -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: testmap.docx Type: application/vnd.openxmlformats-officedocument.wordprocessingml.document Size: 125523 bytes Desc: not available URL: From barry.h.lynn at gmail.com Tue Oct 17 05:46:38 2023 From: barry.h.lynn at gmail.com (Barry Lynn) Date: Tue, 17 Oct 2023 14:46:38 +0300 Subject: [ncl-talk] Maps flipped on their sides Message-ID: Hi: Is there some way to force NCL to plot these maps the correct way, with longitude on the x-axis? I am not doing anything different than I do with other map projections. This is a Lambert-Conformal map. plot = gsn_csm_contour_map(wks,precip,res) -- Barry H. Lynn, Ph.D Senior Scientist, Lecturer, The Institute of Earth Sciences, The Hebrew University of Jerusalem, Givat Ram, Jerusalem 91904, Israel Tel: 972 547 231 170 Fax: (972)-25662581 Weather It Is, LTD Weather and Climate Focus https://weather-it-is.com Jerusalem, Israel Local: 02 930 9525 Cell: 054 7 231 170 Int-IS: x972 2 930 9525 -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: mod_radar_map_2km_April_11_13_2013.pdf Type: application/pdf Size: 853329 bytes Desc: not available URL: From dave.allured at noaa.gov Tue Oct 17 10:42:36 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Tue, 17 Oct 2023 10:42:36 -0600 Subject: [ncl-talk] map question In-Reply-To: References: Message-ID: It looks like NCL example topo_4 does what you want. Here is my suggestion for debugging. Start with the original code for example topo_4. Change only one thing at a time, moving toward your code version, until you find something that causes the unwanted perimeter lines. Change that back to original. Then continue stepping through all the rest of your changes until you have a reduced set of changes that causes the perimeter lines. In this way, you will reduce the problem to something that is easier to understand. On Mon, Oct 16, 2023 at 9:35?AM lzl lai via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > I have a problem deleting those thin lines in my map (the red arrow > pointed area) > . Could you please help me to remove them? Thanks a lot! > > I use "ETOPO_2022_v1_60s_N90W180_bed.nc" ,for drawing topography, > and "gadm41_TWN_0.shp" , for drawing coastline in Taiwan > and also set res at cnMissingValPerimOn = False > > my scrip is below : > ======================================== > ;---Read data > ; topo_file = "../topo/ETOPO_2022_v1_60s_N90W180_geoid.nc" > topo_file = "../topo/ETOPO_2022_v1_60s_N90W180_bed.nc" > f = addfile(topo_file,"r") > elev = f->z > ; > ; Set all values below -100 to missing, hence removing all > ; the ocean elevation values. The ocean will be filled in > ; a light blue (see mpOceanFillColor below). > ; > elev = where(elev.lt.-100.,elev at _FillValue,elev) > > ;;file --- 2 > shpdir = "/home/lzl/data2/lzl/CWBtrend/topo/TW/" > filename ="gadm41_TWN_0.shp" > ;----plot----------------------------------------------------- > wks = gsn_open_wks("pdf","./pic/TW_topo") ; send graphics to > PNG file > gsn_define_colormap(wks,"WhiteBlueGreenYellowRed") > > res = True > res at gsnMaximize = True ; maximize plot in frame > > res at cnFillOn = True ; turn on contour fill > res at cnLinesOn = False ; turn off contour lines > res at cnLineLabelsOn = False ; turn off line labels > res at cnInfoLabelOn = False ; turn off info label > res at lbBoxLinesOn = False ; turn off labelbar box lines > res at cnFillMode = "MeshFill" ; for faster draw > > res at cnLevelSelectionMode = "ManualLevels" > res at cnMinLevelValF = 100. > res at cnMaxLevelValF = 3600. > res at cnLevelSpacingF = 200. > res at cnMissingValPerimOn = False > res at cnMissingValPerimGridBoundOn =False > res at gsnAddCyclic = False ; don't add longitude cyclic > point > > ;---color bar---- > res at lbAutoManage = False ; we control label bar > res at pmLabelBarDisplayMode = "Always" ; turns on label bar > res at lbOrientation = "Horizontal" ; ncl default is vertical > res at pmLabelBarSide = "Bottom" ; default is right > res at lbLabelStride = 1 ; skip every other label > res at pmLabelBarWidthF = 0.5 ; default is shorter > res at pmLabelBarHeightF = 0.1 ; default is taller > res at lbLabelFontHeightF = .018 ; default is HUGE > res at lbPerimOn = False ; default has box > res at lbBoxLinesOn = True ; turn on labelbar box lines > res at lbBoxEndCapStyle = "TriangleBothEnds" > > ;---Zoom in on map > res at mpLimitMode = "LatLon" > res at mpMinLatF = 21.5 > res at mpMaxLatF = 25.5 > res at mpMinLonF = 119.5 > res at mpMaxLonF = 122.5 > > res at mpFillOn = True > res at mpGeophysicalLineColor = "transparent" ; color of cont. > outlines > res at mpOceanFillColor = "transparent" > res at mpLandFillColor = "transparent" > > res at mpFillDrawOrder = "PostDraw" > res at mpDataBaseVersion = "MediumRes" > > res at mpOutlineOn = False > res at pmTickMarkDisplayMode = "Always" > res at mpGeophysicalLineThicknessF = 2 > plot = gsn_csm_contour_map(wks,elev,res) > > nres = True > lnres at gsLineColor = "black" > lnres at gsLineThicknessF = 2.0 > dum = gsn_add_shapefile_polylines(wks,plot,shpdir+filename,lnres) > ;---Generate random lat/lon locations and dummy text strings > lont = (/121.513,120.684,120.197,120.746,121.613,121.154/) > latt = (/25.049,24.12,23.07,22.013,24.05,22.765/) > lons = (/121.513,120.684,120.197,120.746,121.613,121.154/) > lats = (/25.037,24.145,22.995,22.003,23.975,22.752/) > text = (/"Taipei","Taichung","Tainan","Hengchun","Hualian","Taitung"/) > > > ;---Attach some dummy text strings > txres = True > txres at txFontHeightF = 0.02 > txres at txJust = "BottomCenter" > > text_id = gsn_add_text(wks,plot,text,lont,latt,txres) > ; > mkres = True > mkres at gsMarkerIndex = 16 ; Filled circle > mkres at gsMarkerSizeF = 0.01 > > locations = gsn_add_polymarker(wks,plot,lons,lats,mkres) > ;==================================================== > draw(plot) ; This will draw the map and the shapefile outlines. > frame(wks) > _______________________________________________ > -------------- next part -------------- An HTML attachment was scrubbed... URL: From jbuzan at purdue.edu Wed Oct 18 03:09:39 2023 From: jbuzan at purdue.edu (Jonathan Buzan) Date: Wed, 18 Oct 2023 09:09:39 +0000 Subject: [ncl-talk] Drawing Economic Zones derived from netcdf file Message-ID: Hi NCL-Talk, I am attempting to draw a the borders for economic zones on top of my temperature data. The Economic Zones map is the same resolution as the temperature file. The values for the economic zones are 1-16. I was looking through the map outlines page, but I don?t see anything that takes a netcdf file and converts that data into borders. Perhaps NCL users know a method? Thanks, -Jonathan -------------- next part -------------- A non-text attachment was scrubbed... Name: GTAP_FUND_Regions.All.nc Type: application/octet-stream Size: 225872 bytes Desc: GTAP_FUND_Regions.All.nc URL: From setareh.rahimi at gmail.com Wed Oct 18 06:23:05 2023 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Wed, 18 Oct 2023 15:53:05 +0330 Subject: [ncl-talk] Area averaging In-Reply-To: References: Message-ID: Hi Dave, Thanks for your email. Well, the area that I consider to have an area average is a bounding box. So how can I do this, please? any examples? Best wishes, On Sat, Oct 14, 2023 at 10:16?PM Dave Allured - NOAA Affiliate < dave.allured at noaa.gov> wrote: > Yes it is possible. If the grid has 1-D coordinates, then see the family > of *wgt_areaave* functions. If 2-D, then more work is needed. Also, > what kind of region definition do you need? Bounding box, polygon, shape > file, list of grid points, or something else? > > > On Sat, Oct 14, 2023 at 9:39?AM Setareh Rahimi via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> Hi all, >> >> I have a question about area averaging in NCL. I have my WRF-Chem output >> file and would like to calculate average of aerosol optical depth (AOD) >> For an area at let's say 06:00 UTC. So I wonder if this would be possible >> through NCL? >> Thanks for your advise in advance. >> >> Best wishes, >> > -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: From dave.allured at noaa.gov Wed Oct 18 06:25:51 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Wed, 18 Oct 2023 06:25:51 -0600 Subject: [ncl-talk] Drawing Economic Zones derived from netcdf file In-Reply-To: References: Message-ID: Converting a grid of classification data into outlines is possible with custom code, but it would produce blocky outlines. NCL's missing value perimeter lines use this outline algorithm. You may be able to force NCL to draw these perimeter lines as overlays, one at a time for each zone. You could also do this with fully custom code, stepping through all grid cells and erasing shared cell boundaries. Your zone map might be based on a more normally drawn original zone map. I suggest looking into the origins of the Netcdf zone map, to see whether there might be an outline shapefile available, or something like it. That would look good as an overlay on a temperature contour map. On Wed, Oct 18, 2023 at 3:10?AM Jonathan Buzan via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Hi NCL-Talk, > > I am attempting to draw a the borders for economic zones on top of my > temperature data. The Economic Zones map is the same resolution as the > temperature file. The values for the economic zones are 1-16. I was looking > through the map outlines page, but I don?t see anything that takes a netcdf > file and converts that data into borders. > > Perhaps NCL users know a method? > > Thanks, > -Jonathan > -------------- next part -------------- An HTML attachment was scrubbed... URL: From dave.allured at noaa.gov Wed Oct 18 07:06:42 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Wed, 18 Oct 2023 07:06:42 -0600 Subject: [ncl-talk] Area averaging In-Reply-To: References: Message-ID: Setareh, please see the examples page for Weighted Averages. The first two examples are both area averages. Example ave_2.ncl demonstrates coordinate subsetting on a grid with 1-D coordinates. Use that syntax on both the X and Y dimensions to define the bounding box. I suggest compute your own weights as the cosines of the latitudes. There are cosine weight functions available in the lat/lon function category, or just compute your own with cos(). On Wed, Oct 18, 2023 at 6:21?AM Setareh Rahimi wrote: > Hi Dave, > Thanks for your email. Well, the area that I consider to have an area > average is a bounding box. So how can I do this, please? any examples? > Best wishes, > > On Sat, Oct 14, 2023 at 10:16?PM Dave Allured - NOAA Affiliate < > dave.allured at noaa.gov> wrote: > >> Yes it is possible. If the grid has 1-D coordinates, then see the family >> of *wgt_areaave* functions. If 2-D, then more work is needed. Also, >> what kind of region definition do you need? Bounding box, polygon, shape >> file, list of grid points, or something else? >> >> >> On Sat, Oct 14, 2023 at 9:39?AM Setareh Rahimi via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >>> Hi all, >>> >>> I have a question about area averaging in NCL. I have my WRF-Chem output >>> file and would like to calculate average of aerosol optical depth (AOD) >>> For an area at let's say 06:00 UTC. So I wonder if this would be >>> possible through NCL? >>> Thanks for your advise in advance. >>> >>> Best wishes, >>> >> > > -- > S.Rahimi > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From jbuzan at purdue.edu Wed Oct 18 07:12:19 2023 From: jbuzan at purdue.edu (Jonathan Buzan) Date: Wed, 18 Oct 2023 13:12:19 +0000 Subject: [ncl-talk] Drawing Economic Zones derived from netcdf file In-Reply-To: References: Message-ID: Hi Dave, Thanks! I?ll look into this?I do have a higher resolution map. Cheers, -Jonathan On Oct 18, 2023, at 14:25, Dave Allured - NOAA Affiliate > wrote: ---- External Email: Use caution with attachments, links, or sharing data ---- Converting a grid of classification data into outlines is possible with custom code, but it would produce blocky outlines. NCL's missing value perimeter lines use this outline algorithm. You may be able to force NCL to draw these perimeter lines as overlays, one at a time for each zone. You could also do this with fully custom code, stepping through all grid cells and erasing shared cell boundaries. Your zone map might be based on a more normally drawn original zone map. I suggest looking into the origins of the Netcdf zone map, to see whether there might be an outline shapefile available, or something like it. That would look good as an overlay on a temperature contour map. On Wed, Oct 18, 2023 at 3:10?AM Jonathan Buzan via ncl-talk > wrote: Hi NCL-Talk, I am attempting to draw a the borders for economic zones on top of my temperature data. The Economic Zones map is the same resolution as the temperature file. The values for the economic zones are 1-16. I was looking through the map outlines page, but I don?t see anything that takes a netcdf file and converts that data into borders. Perhaps NCL users know a method? Thanks, -Jonathan -------------- next part -------------- An HTML attachment was scrubbed... URL: From leslie.linglai at gmail.com Wed Oct 18 08:43:57 2023 From: leslie.linglai at gmail.com (lzl lai) Date: Wed, 18 Oct 2023 22:43:57 +0800 Subject: [ncl-talk] adjust map longitude latitude font Message-ID: Dear all, I would like to ask how to change the longitude font from 119o30E, 120o30E,...121o30E...(please see the attached figure) . to 119.5E, 120.5E..., 121.5E.. when I try to use res at tmXBMode = "Explicit", and set tmXBValues ,tmXBLabels it shows "MapPlotSetValues: tm[XB|XT|YL|YR]Mode resources are not currently enabled for MapPlot tick marks" Thank you very much! Sincerely -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: TW_topo4.pdf Type: application/pdf Size: 264607 bytes Desc: not available URL: From dave.allured at noaa.gov Wed Oct 18 10:07:02 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Wed, 18 Oct 2023 10:07:02 -0600 Subject: [ncl-talk] Maps flipped on their sides In-Reply-To: References: Message-ID: Barry, I am trying to figure out if this is a side effect of the panel or overlay functions. What happens if you plot only the first map by itself, without any paneling? (Set res at gsnDraw and res at gsnFrame to True, right before gsn_csm_contour_map. Then exit immediately after.) On Tue, Oct 17, 2023 at 5:47?AM Barry Lynn via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Hi: > > Is there some way to force NCL to plot these maps the correct way, with > longitude on the x-axis? > > I am not doing anything different than I do with other map projections. > This is a Lambert-Conformal map. > > plot = gsn_csm_contour_map(wks,precip,res) > > Barry H. Lynn, Ph.D > Senior Scientist, Lecturer, > The Institute of Earth Sciences, > The Hebrew University of Jerusalem, > Givat Ram, Jerusalem 91904, Israel > Tel: 972 547 231 170 > Fax: (972)-25662581 > > Weather It Is, LTD > Weather and Climate Focus > https://weather-it-is.com > Jerusalem, Israel > Local: 02 930 9525 > Cell: 054 7 231 170 > Int-IS: x972 2 930 9525 > -------------- next part -------------- An HTML attachment was scrubbed... URL: From barry.h.lynn at gmail.com Wed Oct 18 10:23:48 2023 From: barry.h.lynn at gmail.com (Barry Lynn) Date: Wed, 18 Oct 2023 19:23:48 +0300 Subject: [ncl-talk] Maps flipped on their sides In-Reply-To: References: Message-ID: Hi Dave: I plotted the first (and the second), since I am not sure the code would have worked. They are the same. Barry On Wed, Oct 18, 2023 at 7:07?PM Dave Allured - NOAA Affiliate via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Barry, I am trying to figure out if this is a side effect of the panel or > overlay functions. What happens if you plot only the first map by itself, > without any paneling? (Set res at gsnDraw and res at gsnFrame to True, right > before gsn_csm_contour_map. Then exit immediately after.) > > > On Tue, Oct 17, 2023 at 5:47?AM Barry Lynn via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> Hi: >> >> Is there some way to force NCL to plot these maps the correct way, with >> longitude on the x-axis? >> >> I am not doing anything different than I do with other map projections. >> This is a Lambert-Conformal map. >> >> plot = gsn_csm_contour_map(wks,precip,res) >> >> Barry H. Lynn, Ph.D >> Senior Scientist, Lecturer, >> The Institute of Earth Sciences, >> The Hebrew University of Jerusalem, >> Givat Ram, Jerusalem 91904, Israel >> Tel: 972 547 231 170 >> Fax: (972)-25662581 >> >> Weather It Is, LTD >> Weather and Climate Focus >> https://weather-it-is.com >> Jerusalem, Israel >> Local: 02 930 9525 >> Cell: 054 7 231 170 >> Int-IS: x972 2 930 9525 >> > _______________________________________________ > ncl-talk mailing list > ncl-talk at mailman.ucar.edu > List instructions, subscriber options, unsubscribe: > https://mailman.ucar.edu/mailman/listinfo/ncl-talk > -- Barry H. Lynn, Ph.D Senior Scientist, Lecturer, The Institute of Earth Sciences, The Hebrew University of Jerusalem, Givat Ram, Jerusalem 91904, Israel Tel: 972 547 231 170 Fax: (972)-25662581 Weather It Is, LTD Weather and Climate Focus https://weather-it-is.com Jerusalem, Israel Local: 02 930 9525 Cell: 054 7 231 170 Int-IS: x972 2 930 9525 -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: mod_radar_map_2km_April_11_13_2013.pdf Type: application/pdf Size: 738750 bytes Desc: not available URL: From asphilli at ucar.edu Wed Oct 18 10:34:55 2023 From: asphilli at ucar.edu (Adam Phillips) Date: Wed, 18 Oct 2023 10:34:55 -0600 Subject: [ncl-talk] adjust map longitude latitude font In-Reply-To: References: Message-ID: The tm*Mode resources only work to alter the latitude/longitude labels on cylindrical equdistant projections. Please see example #9 here: https://www.ncl.ucar.edu/Applications/mptick.shtml#ex9 You will likely need to use the add_map_tickmarks function embedded in that example to get the labels you want. Hope that helps! Adam On Wed, Oct 18, 2023 at 8:44?AM lzl lai via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear all, > I would like to ask how to change the longitude font from > 119o30E, 120o30E,...121o30E...(please see the attached figure) . > to 119.5E, 120.5E..., 121.5E.. > > when I try to use res at tmXBMode = "Explicit", and set tmXBValues > ,tmXBLabels > it shows "MapPlotSetValues: tm[XB|XT|YL|YR]Mode resources are not > currently enabled for MapPlot tick marks" > > Thank you very much! > > Sincerely > > _______________________________________________ > ncl-talk mailing list > ncl-talk at mailman.ucar.edu > List instructions, subscriber options, unsubscribe: > https://mailman.ucar.edu/mailman/listinfo/ncl-talk > -- Adam Phillips Associate Scientist IV, Climate Analysis Section Climate and Global Dynamics Laboratory National Center for Atmospheric Research www.cgd.ucar.edu/staff/asphilli/ -------------- next part -------------- An HTML attachment was scrubbed... URL: From dave.allured at noaa.gov Wed Oct 18 10:36:32 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Wed, 18 Oct 2023 10:36:32 -0600 Subject: [ncl-talk] Maps flipped on their sides In-Reply-To: References: Message-ID: Yeah. I remember that Lambert mode can be a bit quirky. Please send your code as an attachment. You may send privately to my address, if you prefer. On Wed, Oct 18, 2023 at 10:24?AM Barry Lynn wrote: > Hi Dave: > > I plotted the first (and the second), since I am not sure the code would > have worked. They are the same. > > Barry > > On Wed, Oct 18, 2023 at 7:07?PM Dave Allured - NOAA Affiliate via ncl-talk > wrote: > >> Barry, I am trying to figure out if this is a side effect of the panel or >> overlay functions. What happens if you plot only the first map by itself, >> without any paneling? (Set res at gsnDraw and res at gsnFrame to True, right >> before gsn_csm_contour_map. Then exit immediately after.) >> >> >> On Tue, Oct 17, 2023 at 5:47?AM Barry Lynn via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >>> Hi: >>> >>> Is there some way to force NCL to plot these maps the correct way, with >>> longitude on the x-axis? >>> >>> I am not doing anything different than I do with other map projections. >>> This is a Lambert-Conformal map. >>> >>> plot = gsn_csm_contour_map(wks,precip,res) >>> >>> Barry H. Lynn, Ph.D >>> Senior Scientist, Lecturer, >>> The Institute of Earth Sciences, >>> The Hebrew University of Jerusalem, >>> Givat Ram, Jerusalem 91904, Israel >>> Tel: 972 547 231 170 >>> Fax: (972)-25662581 >>> >>> Weather It Is, LTD >>> Weather and Climate Focus >>> https://weather-it-is.com >>> Jerusalem, Israel >>> Local: 02 930 9525 >>> Cell: 054 7 231 170 >>> Int-IS: x972 2 930 9525 >>> >> -------------- next part -------------- An HTML attachment was scrubbed... URL: From barry.h.lynn at gmail.com Wed Oct 18 11:00:00 2023 From: barry.h.lynn at gmail.com (Barry Lynn) Date: Wed, 18 Oct 2023 20:00:00 +0300 Subject: [ncl-talk] Maps flipped on their sides In-Reply-To: References: Message-ID: Dear Dave: I was about to send you the code when I decided to add some comments. When I got to this section, I realized that the last setting in bold might be the problem. res at mpProjection = "LambertConformal" ; res at mpProjection = "Mercator" res at mpLambertParallel1F = 30 res at mpLambertParallel2F = 58 *res at mpLambertMeridianF = 30 ; was 258* The problem was that the Meridian was set for my previous study, or the central United States. Thank you for helping me solve this problem. Barry On Wed, Oct 18, 2023 at 7:23?PM Barry Lynn wrote: > Hi Dave: > > I plotted the first (and the second), since I am not sure the code would > have worked. They are the same. > > Barry > > On Wed, Oct 18, 2023 at 7:07?PM Dave Allured - NOAA Affiliate via ncl-talk > wrote: > >> Barry, I am trying to figure out if this is a side effect of the panel or >> overlay functions. What happens if you plot only the first map by itself, >> without any paneling? (Set res at gsnDraw and res at gsnFrame to True, right >> before gsn_csm_contour_map. Then exit immediately after.) >> >> >> On Tue, Oct 17, 2023 at 5:47?AM Barry Lynn via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >>> Hi: >>> >>> Is there some way to force NCL to plot these maps the correct way, with >>> longitude on the x-axis? >>> >>> I am not doing anything different than I do with other map projections. >>> This is a Lambert-Conformal map. >>> >>> plot = gsn_csm_contour_map(wks,precip,res) >>> >>> Barry H. Lynn, Ph.D >>> Senior Scientist, Lecturer, >>> The Institute of Earth Sciences, >>> The Hebrew University of Jerusalem, >>> Givat Ram, Jerusalem 91904, Israel >>> Tel: 972 547 231 170 >>> Fax: (972)-25662581 >>> >>> Weather It Is, LTD >>> Weather and Climate Focus >>> https://weather-it-is.com >>> Jerusalem, Israel >>> Local: 02 930 9525 >>> Cell: 054 7 231 170 >>> Int-IS: x972 2 930 9525 >>> >> _______________________________________________ >> ncl-talk mailing list >> ncl-talk at mailman.ucar.edu >> List instructions, subscriber options, unsubscribe: >> https://mailman.ucar.edu/mailman/listinfo/ncl-talk >> > > > -- > > Barry H. Lynn, Ph.D > Senior Scientist, Lecturer, > The Institute of Earth Sciences, > The Hebrew University of Jerusalem, > Givat Ram, Jerusalem 91904, Israel > Tel: 972 547 231 170 > Fax: (972)-25662581 > > Weather It Is, LTD > Weather and Climate Focus > https://weather-it-is.com > Jerusalem, Israel > Local: 02 930 9525 > Cell: 054 7 231 170 > Int-IS: x972 2 930 9525 > > -- Barry H. Lynn, Ph.D Senior Scientist, Lecturer, The Institute of Earth Sciences, The Hebrew University of Jerusalem, Givat Ram, Jerusalem 91904, Israel Tel: 972 547 231 170 Fax: (972)-25662581 Weather It Is, LTD Weather and Climate Focus https://weather-it-is.com Jerusalem, Israel Local: 02 930 9525 Cell: 054 7 231 170 Int-IS: x972 2 930 9525 -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: mod_radar_map_2km_April_11_13_2013.pdf Type: application/pdf Size: 736449 bytes Desc: not available URL: From barry.h.lynn at gmail.com Wed Oct 18 11:33:44 2023 From: barry.h.lynn at gmail.com (Barry Lynn) Date: Wed, 18 Oct 2023 20:33:44 +0300 Subject: [ncl-talk] Maps flipped on their sides In-Reply-To: References: Message-ID: Dave: Now that that is solved, I am perplexed why the first map would appear with land colored as grey, but not when I loop through again... Is there an obvious reason for this? Barry On Wed, Oct 18, 2023 at 7:37?PM Dave Allured - NOAA Affiliate via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Yeah. I remember that Lambert mode can be a bit quirky. Please send your > code as an attachment. You may send privately to my address, if you prefer. > > > On Wed, Oct 18, 2023 at 10:24?AM Barry Lynn > wrote: > >> Hi Dave: >> >> I plotted the first (and the second), since I am not sure the code would >> have worked. They are the same. >> >> Barry >> >> On Wed, Oct 18, 2023 at 7:07?PM Dave Allured - NOAA Affiliate via >> ncl-talk wrote: >> >>> Barry, I am trying to figure out if this is a side effect of the panel >>> or overlay functions. What happens if you plot only the first map by >>> itself, without any paneling? (Set res at gsnDraw and res at gsnFrame to >>> True, right before gsn_csm_contour_map. Then exit immediately after.) >>> >>> >>> On Tue, Oct 17, 2023 at 5:47?AM Barry Lynn via ncl-talk < >>> ncl-talk at mailman.ucar.edu> wrote: >>> >>>> Hi: >>>> >>>> Is there some way to force NCL to plot these maps the correct way, with >>>> longitude on the x-axis? >>>> >>>> I am not doing anything different than I do with other map >>>> projections. This is a Lambert-Conformal map. >>>> >>>> plot = gsn_csm_contour_map(wks,precip,res) >>>> >>>> Barry H. Lynn, Ph.D >>>> Senior Scientist, Lecturer, >>>> The Institute of Earth Sciences, >>>> The Hebrew University of Jerusalem, >>>> Givat Ram, Jerusalem 91904, Israel >>>> Tel: 972 547 231 170 >>>> Fax: (972)-25662581 >>>> >>>> Weather It Is, LTD >>>> Weather and Climate Focus >>>> https://weather-it-is.com >>>> Jerusalem, Israel >>>> Local: 02 930 9525 >>>> Cell: 054 7 231 170 >>>> Int-IS: x972 2 930 9525 >>>> >>> _______________________________________________ > ncl-talk mailing list > ncl-talk at mailman.ucar.edu > List instructions, subscriber options, unsubscribe: > https://mailman.ucar.edu/mailman/listinfo/ncl-talk > -- Barry H. Lynn, Ph.D Senior Scientist, Lecturer, The Institute of Earth Sciences, The Hebrew University of Jerusalem, Givat Ram, Jerusalem 91904, Israel Tel: 972 547 231 170 Fax: (972)-25662581 Weather It Is, LTD Weather and Climate Focus https://weather-it-is.com Jerusalem, Israel Local: 02 930 9525 Cell: 054 7 231 170 Int-IS: x972 2 930 9525 -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: mod_radar.ncl Type: application/octet-stream Size: 18698 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: mod_radar_map_2km_April_11_13_2013.pdf Type: application/pdf Size: 1384952 bytes Desc: not available URL: From setareh.rahimi at gmail.com Wed Oct 18 11:44:17 2023 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Wed, 18 Oct 2023 21:14:17 +0330 Subject: [ncl-talk] Area averaging In-Reply-To: References: Message-ID: Dear Dave, Thanks a lot for your help. Best wishes, On Wed, Oct 18, 2023 at 4:37?PM Dave Allured - NOAA Affiliate < dave.allured at noaa.gov> wrote: > Setareh, please see the examples page for Weighted Averages. The first > two examples are both area averages. Example ave_2.ncl demonstrates > coordinate subsetting on a grid with 1-D coordinates. Use that syntax on > both the X and Y dimensions to define the bounding box. I suggest compute > your own weights as the cosines of the latitudes. There are cosine weight > functions available in the lat/lon function category, or just compute your > own with cos(). > > > On Wed, Oct 18, 2023 at 6:21?AM Setareh Rahimi > wrote: > >> Hi Dave, >> Thanks for your email. Well, the area that I consider to have an area >> average is a bounding box. So how can I do this, please? any examples? >> Best wishes, >> >> On Sat, Oct 14, 2023 at 10:16?PM Dave Allured - NOAA Affiliate < >> dave.allured at noaa.gov> wrote: >> >>> Yes it is possible. If the grid has 1-D coordinates, then see the >>> family of *wgt_areaave* functions. If 2-D, then more work is needed. >>> Also, what kind of region definition do you need? Bounding box, polygon, >>> shape file, list of grid points, or something else? >>> >>> >>> On Sat, Oct 14, 2023 at 9:39?AM Setareh Rahimi via ncl-talk < >>> ncl-talk at mailman.ucar.edu> wrote: >>> >>>> Hi all, >>>> >>>> I have a question about area averaging in NCL. I have my WRF-Chem >>>> output file and would like to calculate average of aerosol optical depth >>>> (AOD) >>>> For an area at let's say 06:00 UTC. So I wonder if this would be >>>> possible through NCL? >>>> Thanks for your advise in advance. >>>> >>>> Best wishes, >>>> >>> >> >> -- >> S.Rahimi >> >> -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: From setareh.rahimi at gmail.com Wed Oct 18 11:50:25 2023 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Wed, 18 Oct 2023 21:20:25 +0330 Subject: [ncl-talk] Cross Sections Plot Message-ID: Dear all, I need to plot a cross-section of Aerosol Optical Depth (AOD) from the WRF-Chem model ( at say lat:20N, along lon : 0 to 23 E). I had a look at the below links which give an example of crosssection plotting. but this example provides a crosssection of an area and not a line. So How can I plot a crosssection as I described before, please? https://www.ncl.ucar.edu/Applications/wrfinterp.shtml https://www.ncl.ucar.edu/Applications/Scripts/wrf_interp_3.ncl Any suggestion would be appreciated. Kind regards. -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: From dave.allured at noaa.gov Wed Oct 18 13:03:41 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Wed, 18 Oct 2023 13:03:41 -0600 Subject: [ncl-talk] Cross Sections Plot In-Reply-To: References: Message-ID: I think start_lat = 20, end_lat = 20 should work here. On Wed, Oct 18, 2023 at 11:49?AM Setareh Rahimi via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear all, > > I need to plot a cross-section of Aerosol Optical Depth (AOD) from the > WRF-Chem model ( at say lat:20N, along lon : 0 to 23 E). > I had a look at the below links which give an example of crosssection > plotting. but this example provides a crosssection of an area and not a > line. So How can I plot a crosssection as I described before, please? > > https://www.ncl.ucar.edu/Applications/wrfinterp.shtml > https://www.ncl.ucar.edu/Applications/Scripts/wrf_interp_3.ncl > > Any suggestion would be appreciated. > > Kind regards. > -- > S.Rahimi > -------------- next part -------------- An HTML attachment was scrubbed... URL: From setareh.rahimi at gmail.com Sat Oct 21 06:13:37 2023 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Sat, 21 Oct 2023 15:43:37 +0330 Subject: [ncl-talk] Cross Sections Plot In-Reply-To: References: Message-ID: Hi Dave, Thanks for your advice. However, I calculated AOD, from 2 other variables namely: "*TAUAER2* "& "*TAUAER3*" ( attached). Therefore I cannot use the* "**wrf_user_vert_cross*" function to plot AOD in cross-section. So would you please help me plot AOD in cross-section please? Kind regards. -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: aod Type: application/octet-stream Size: 6655 bytes Desc: not available URL: From setareh.rahimi at gmail.com Sun Oct 22 04:59:36 2023 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Sun, 22 Oct 2023 14:29:36 +0330 Subject: [ncl-talk] Total daily precipitation from the WRF model Message-ID: Dear all, How can I calculate *total daily precipitation* from the WRF model, please ( for example total daily precipitation for 2021-01-014)? I have got hourly out from the model, (*history_interval=60, **frames_per_outfile=72*). Please advise me on this issue. Thanks in advance, -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: From zmumba at gmail.com Mon Oct 23 07:17:15 2023 From: zmumba at gmail.com (Zilore Mumba) Date: Mon, 23 Oct 2023 15:17:15 +0200 Subject: [ncl-talk] Total daily precipitation from the WRF model In-Reply-To: References: Message-ID: Setareh, Perhaps you can try the attached script, which I use to plot 24-hr rainfall up to 5 days ahead. Essentially if you get your WRF rainfall as RAINC = a[:]->RAINC(:,:,:) WRF rainfall is accumulated hourly, so to get 24hour total rainfall on day one from analysis, you would, e.g. plot = gsn_csm_contour_map(wks, RAIN((24,:,:)-RAIN((0,:,:), res) I hope it helps On Sun, Oct 22, 2023 at 12:58?PM Setareh Rahimi via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear all, > > How can I calculate *total daily precipitation* from the WRF model, > please ( for example total daily precipitation for 2021-01-014)? I have got > hourly out from the model, (*history_interval=60, **frames_per_outfile=72* > ). > Please advise me on this issue. > Thanks in advance, > -- > S.Rahimi > > _______________________________________________ > ncl-talk mailing list > ncl-talk at mailman.ucar.edu > List instructions, subscriber options, unsubscribe: > https://mailman.ucar.edu/mailman/listinfo/ncl-talk > -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 02_Plot24HrRainfall.ncl Type: application/octet-stream Size: 6888 bytes Desc: not available URL: From setareh.rahimi at gmail.com Sun Oct 29 11:12:02 2023 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Sun, 29 Oct 2023 20:42:02 +0330 Subject: [ncl-talk] Cross section plot for AOD Message-ID: Dear all, I have calculated Aerosol Optical Depth (AOD) from WRF-Chem output using "TAUAER2" and "TAUAER3" (script attached). I would like to plot AOD for a specific latitude, on height. So after calculating AOD, I used *"gsn_csm_pres_hgt" *to do so. However, I have got the below error: *fatal: Dimension (level) of (AOD550_3D) does not have an associated coordinate variable* So how can I remove this error, please? Thanks in advance. S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: aod Type: application/octet-stream Size: 6985 bytes Desc: not available URL: From dave.allured at noaa.gov Sun Oct 29 12:37:32 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Sun, 29 Oct 2023 12:37:32 -0600 Subject: [ncl-talk] Cross section plot for AOD In-Reply-To: References: Message-ID: Find the level coordinate on one of your input variables, and attach a copy to AOD after the AOD calculation. It should be something like "aod&level = input&level". Please refer to the section on coordinate variables in the user manual. If you have trouble, then please show printVarSummary for AOD, and for one of the input variables that has levels. On Sun, Oct 29, 2023 at 11:10?AM Setareh Rahimi via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear all, > I have calculated Aerosol Optical Depth (AOD) from WRF-Chem output using > "TAUAER2" and "TAUAER3" (script attached). I would like to plot AOD for a > specific latitude, on height. So after calculating AOD, I used *"gsn_csm_pres_hgt" > *to do so. However, I have got the below error: > > *fatal: Dimension (level) of (AOD550_3D) does not have an associated > coordinate variable* > > > > So how can I remove this error, please? > > Thanks in advance. > > > S.Rahimi > -------------- next part -------------- An HTML attachment was scrubbed... URL: From barry.h.lynn at gmail.com Mon Oct 30 08:07:09 2023 From: barry.h.lynn at gmail.com (Barry Lynn) Date: Mon, 30 Oct 2023 16:07:09 +0200 Subject: [ncl-talk] Strange problem related to linint1_Wrap Message-ID: Hello: I wanted to use linint1 to interpolate WRF data to height surfaces with constant spacing between levels. But, there are two mysteries. You can see the WRF data in a cross-section of radar in (a; the variable ff1_s below). When I interpolate the data (b, ff1) is compacted towards the bottom of the cross-sectional plot. Moreover, the height information on the y-axis will not plot. I understand why the data is compacted, but not why I can't eliminate the white space above. If I reduce the range of WRF data sampled in the z direction, I just have a truncated cross-section (it doesn't get any closer to the top of the cross-section). I don't see why using ff1 (from the interpolate) would make the labeling disappear. Variable: ff1_s Type: float Total Size: 33728 bytes 8432 values Number of Dimensions: 2 Dimensions and sizes: [height | 31] x [lon | 272] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 Variable: ff1 Type: float Total Size: 33728 bytes 8432 values Number of Dimensions: 2 Dimensions and sizes: [HEIGHT | 31] x [lon | 272] Coordinates: HEIGHT: [500..15500] Number Of Attributes: 1 _FillValue : 9.96921e+36 -- Barry H. Lynn, Ph.D Senior Scientist, Lecturer, The Institute of Earth Sciences, The Hebrew University of Jerusalem, Givat Ram, Jerusalem 91904, Israel Tel: 972 547 231 170 Fax: (972)-25662581 Weather It Is, LTD Weather and Climate Focus https://weather-it-is.com Jerusalem, Israel Local: 02 930 9525 Cell: 054 7 231 170 Int-IS: x972 2 930 9525 -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: radar_vert_crossection_10UTC.pdf Type: application/pdf Size: 83216 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: test2.ncl Type: application/octet-stream Size: 11462 bytes Desc: not available URL: From dave.allured at noaa.gov Mon Oct 30 09:16:45 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Mon, 30 Oct 2023 09:16:45 -0600 Subject: [ncl-talk] Strange problem related to linint1_Wrap In-Reply-To: References: Message-ID: Barry, please notice that there is a HEIGHT coordinate variable attached to ff1 after interpolation, but there is no such coordinate variable attached to ff1_s before interpolation. It looks like the attached HEIGHT coordinate variable is misleading the plot function in some way. Try delete ("ff1&HEIGHT"). I am not familiar with the cross section plotting function. But it looks like in the first plot, you are setting some plot resources to set the range, spacing, and units of the height axis. I suspect that these settings are in conflict with the attached HEIGHT coordinates which may have been generated by linint1_Wrap. The units are different. Instead of deleting the HEIGHT coordinate variable, you could try deleting or adjusting the conflicting resource settings. If you want to rescale the attached HEIGHT coordinates from meters to km, then simply "ff1&HEIGHT = ff1&HEIGHT / 1000". On Mon, Oct 30, 2023 at 8:08?AM Barry Lynn via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Hello: > > I wanted to use linint1 to interpolate WRF data to height surfaces with > constant spacing between levels. But, there are two mysteries. > > You can see the WRF data in a cross-section of radar in (a; the variable > ff1_s below). When I interpolate the data (b, ff1) is compacted towards > the bottom of the cross-sectional plot. Moreover, the height information > on the y-axis will not plot. > > I understand why the data is compacted, but not why I can't eliminate the > white space above. If I reduce the range of WRF data sampled in the z > direction, I just have a truncated cross-section (it doesn't get any closer > to the top of the cross-section). > > I don't see why using ff1 (from the interpolate) would make the labeling > disappear. > > Variable: ff1_s > > Type: float > > Total Size: 33728 bytes > > 8432 values > > Number of Dimensions: 2 > > Dimensions and sizes: [height | 31] x [lon | 272] > > Coordinates: > > Number Of Attributes: 1 > > _FillValue : 9.96921e+36 > > > Variable: ff1 > > Type: float > > Total Size: 33728 bytes > > 8432 values > > Number of Dimensions: 2 > > Dimensions and sizes: [HEIGHT | 31] x [lon | 272] > > Coordinates: > > HEIGHT: [500..15500] > > Number Of Attributes: 1 > > _FillValue : 9.96921e+36 > > -- > Barry H. Lynn, Ph.D > Senior Scientist, Lecturer, > The Institute of Earth Sciences, > The Hebrew University of Jerusalem, > Givat Ram, Jerusalem 91904, Israel > Tel: 972 547 231 170 > Fax: (972)-25662581 > > Weather It Is, LTD > Weather and Climate Focus > https://weather-it-is.com > Jerusalem, Israel > Local: 02 930 9525 > Cell: 054 7 231 170 > Int-IS: x972 2 930 9525 > -------------- next part -------------- An HTML attachment was scrubbed... URL: From dave.allured at noaa.gov Mon Oct 30 09:20:19 2023 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Mon, 30 Oct 2023 09:20:19 -0600 Subject: [ncl-talk] Strange problem related to linint1_Wrap In-Reply-To: References: Message-ID: Correction. Do not use quotes here: delete (ff1&HEIGHT) On Mon, Oct 30, 2023 at 9:16?AM Dave Allured - NOAA Affiliate < dave.allured at noaa.gov> wrote: > Barry, please notice that there is a HEIGHT coordinate variable attached > to ff1 after interpolation, but there is no such coordinate variable > attached to ff1_s before interpolation. It looks like the attached HEIGHT > coordinate variable is misleading the plot function in some way. Try > delete ("ff1&HEIGHT"). > > I am not familiar with the cross section plotting function. But it looks > like in the first plot, you are setting some plot resources to set the > range, spacing, and units of the height axis. I suspect that these > settings are in conflict with the attached HEIGHT coordinates which may > have been generated by linint1_Wrap. The units are different. Instead of > deleting the HEIGHT coordinate variable, you could try deleting or > adjusting the conflicting resource settings. > > If you want to rescale the attached HEIGHT coordinates from meters to km, > then simply "ff1&HEIGHT = ff1&HEIGHT / 1000". > > > On Mon, Oct 30, 2023 at 8:08?AM Barry Lynn via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> Hello: >> >> I wanted to use linint1 to interpolate WRF data to height surfaces with >> constant spacing between levels. But, there are two mysteries. >> >> You can see the WRF data in a cross-section of radar in (a; the variable >> ff1_s below). When I interpolate the data (b, ff1) is compacted towards >> the bottom of the cross-sectional plot. Moreover, the height information >> on the y-axis will not plot. >> >> I understand why the data is compacted, but not why I can't eliminate the >> white space above. If I reduce the range of WRF data sampled in the z >> direction, I just have a truncated cross-section (it doesn't get any closer >> to the top of the cross-section). >> >> I don't see why using ff1 (from the interpolate) would make the labeling >> disappear. >> >> Variable: ff1_s >> >> Type: float >> >> Total Size: 33728 bytes >> >> 8432 values >> >> Number of Dimensions: 2 >> >> Dimensions and sizes: [height | 31] x [lon | 272] >> >> Coordinates: >> >> Number Of Attributes: 1 >> >> _FillValue : 9.96921e+36 >> >> >> Variable: ff1 >> >> Type: float >> >> Total Size: 33728 bytes >> >> 8432 values >> >> Number of Dimensions: 2 >> >> Dimensions and sizes: [HEIGHT | 31] x [lon | 272] >> >> Coordinates: >> >> HEIGHT: [500..15500] >> >> Number Of Attributes: 1 >> >> _FillValue : 9.96921e+36 >> >> -- >> Barry H. Lynn, Ph.D >> Senior Scientist, Lecturer, >> The Institute of Earth Sciences, >> The Hebrew University of Jerusalem, >> Givat Ram, Jerusalem 91904, Israel >> Tel: 972 547 231 170 >> Fax: (972)-25662581 >> >> Weather It Is, LTD >> Weather and Climate Focus >> https://weather-it-is.com >> Jerusalem, Israel >> Local: 02 930 9525 >> Cell: 054 7 231 170 >> Int-IS: x972 2 930 9525 >> > -------------- next part -------------- An HTML attachment was scrubbed... URL: From barry.h.lynn at gmail.com Mon Oct 30 10:12:13 2023 From: barry.h.lynn at gmail.com (Barry Lynn) Date: Mon, 30 Oct 2023 18:12:13 +0200 Subject: [ncl-talk] Strange problem related to linint1_Wrap In-Reply-To: References: Message-ID: Hi Dave: I thought those quotes were "suspicious." The labeling is the same regardless of which variable is called for the plot, so that wasn't the problem. *Regardless*, your solution worked! (delete (ff1&HEIGHT)). Thank you! On Mon, Oct 30, 2023 at 5:20?PM Dave Allured - NOAA Affiliate < dave.allured at noaa.gov> wrote: > Correction. Do not use quotes here: delete (ff1&HEIGHT) > > > On Mon, Oct 30, 2023 at 9:16?AM Dave Allured - NOAA Affiliate < > dave.allured at noaa.gov> wrote: > >> Barry, please notice that there is a HEIGHT coordinate variable attached >> to ff1 after interpolation, but there is no such coordinate variable >> attached to ff1_s before interpolation. It looks like the attached HEIGHT >> coordinate variable is misleading the plot function in some way. Try >> delete ("ff1&HEIGHT"). >> >> I am not familiar with the cross section plotting function. But it looks >> like in the first plot, you are setting some plot resources to set the >> range, spacing, and units of the height axis. I suspect that these >> settings are in conflict with the attached HEIGHT coordinates which may >> have been generated by linint1_Wrap. The units are different. Instead of >> deleting the HEIGHT coordinate variable, you could try deleting or >> adjusting the conflicting resource settings. >> >> If you want to rescale the attached HEIGHT coordinates from meters to km, >> then simply "ff1&HEIGHT = ff1&HEIGHT / 1000". >> >> >> On Mon, Oct 30, 2023 at 8:08?AM Barry Lynn via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >>> Hello: >>> >>> I wanted to use linint1 to interpolate WRF data to height surfaces with >>> constant spacing between levels. But, there are two mysteries. >>> >>> You can see the WRF data in a cross-section of radar in (a; the variable >>> ff1_s below). When I interpolate the data (b, ff1) is compacted towards >>> the bottom of the cross-sectional plot. Moreover, the height information >>> on the y-axis will not plot. >>> >>> I understand why the data is compacted, but not why I can't eliminate >>> the white space above. If I reduce the range of WRF data sampled in the z >>> direction, I just have a truncated cross-section (it doesn't get any closer >>> to the top of the cross-section). >>> >>> I don't see why using ff1 (from the interpolate) would make the labeling >>> disappear. >>> >>> Variable: ff1_s >>> >>> Type: float >>> >>> Total Size: 33728 bytes >>> >>> 8432 values >>> >>> Number of Dimensions: 2 >>> >>> Dimensions and sizes: [height | 31] x [lon | 272] >>> >>> Coordinates: >>> >>> Number Of Attributes: 1 >>> >>> _FillValue : 9.96921e+36 >>> >>> >>> Variable: ff1 >>> >>> Type: float >>> >>> Total Size: 33728 bytes >>> >>> 8432 values >>> >>> Number of Dimensions: 2 >>> >>> Dimensions and sizes: [HEIGHT | 31] x [lon | 272] >>> >>> Coordinates: >>> >>> HEIGHT: [500..15500] >>> >>> Number Of Attributes: 1 >>> >>> _FillValue : 9.96921e+36 >>> >>> -- >>> Barry H. Lynn, Ph.D >>> Senior Scientist, Lecturer, >>> The Institute of Earth Sciences, >>> The Hebrew University of Jerusalem, >>> Givat Ram, Jerusalem 91904, Israel >>> Tel: 972 547 231 170 >>> Fax: (972)-25662581 >>> >>> Weather It Is, LTD >>> Weather and Climate Focus >>> https://weather-it-is.com >>> Jerusalem, Israel >>> Local: 02 930 9525 >>> Cell: 054 7 231 170 >>> Int-IS: x972 2 930 9525 >>> >> -- Barry H. Lynn, Ph.D Senior Scientist, Lecturer, The Institute of Earth Sciences, The Hebrew University of Jerusalem, Givat Ram, Jerusalem 91904, Israel Tel: 972 547 231 170 Fax: (972)-25662581 Weather It Is, LTD Weather and Climate Focus https://weather-it-is.com Jerusalem, Israel Local: 02 930 9525 Cell: 054 7 231 170 Int-IS: x972 2 930 9525 -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: radar_vert_crossection_10UTC.pdf Type: application/pdf Size: 86801 bytes Desc: not available URL: From barry.h.lynn at gmail.com Mon Oct 30 12:21:27 2023 From: barry.h.lynn at gmail.com (Barry Lynn) Date: Mon, 30 Oct 2023 20:21:27 +0200 Subject: [ncl-talk] Strange problem related to linint1_Wrap In-Reply-To: References: Message-ID: Dear All (& Dave): I still have the strange problem where there is a lot of white space above the mapped radar cross-section. If, however, I reduce the scope of the data (i_z is reduced), then the top of the contours get cut off. So, I have a feeling that the gsn plotting program is assuming some kind of non-linear height field, while my data is in height coordinates. Is this correct? If not, what might be the solution? Note, I'm not using any WRF functions and I converted the WRF data to equally spaced data levels. Thank you, Barry On Mon, Oct 30, 2023 at 6:12?PM Barry Lynn wrote: > Hi Dave: > > I thought those quotes were "suspicious." > > The labeling is the same regardless of which variable is called for the > plot, so that wasn't the problem. > > *Regardless*, your solution worked! (delete (ff1&HEIGHT)). > > Thank you! > > On Mon, Oct 30, 2023 at 5:20?PM Dave Allured - NOAA Affiliate < > dave.allured at noaa.gov> wrote: > >> Correction. Do not use quotes here: delete (ff1&HEIGHT) >> >> >> On Mon, Oct 30, 2023 at 9:16?AM Dave Allured - NOAA Affiliate < >> dave.allured at noaa.gov> wrote: >> >>> Barry, please notice that there is a HEIGHT coordinate variable attached >>> to ff1 after interpolation, but there is no such coordinate variable >>> attached to ff1_s before interpolation. It looks like the attached HEIGHT >>> coordinate variable is misleading the plot function in some way. Try >>> delete ("ff1&HEIGHT"). >>> >>> I am not familiar with the cross section plotting function. But it >>> looks like in the first plot, you are setting some plot resources to set >>> the range, spacing, and units of the height axis. I suspect that these >>> settings are in conflict with the attached HEIGHT coordinates which may >>> have been generated by linint1_Wrap. The units are different. Instead of >>> deleting the HEIGHT coordinate variable, you could try deleting or >>> adjusting the conflicting resource settings. >>> >>> If you want to rescale the attached HEIGHT coordinates from meters to >>> km, then simply "ff1&HEIGHT = ff1&HEIGHT / 1000". >>> >>> >>> On Mon, Oct 30, 2023 at 8:08?AM Barry Lynn via ncl-talk < >>> ncl-talk at mailman.ucar.edu> wrote: >>> >>>> Hello: >>>> >>>> I wanted to use linint1 to interpolate WRF data to height surfaces with >>>> constant spacing between levels. But, there are two mysteries. >>>> >>>> You can see the WRF data in a cross-section of radar in (a; the >>>> variable ff1_s below). When I interpolate the data (b, ff1) is compacted >>>> towards the bottom of the cross-sectional plot. Moreover, the height >>>> information on the y-axis will not plot. >>>> >>>> I understand why the data is compacted, but not why I can't eliminate >>>> the white space above. If I reduce the range of WRF data sampled in the z >>>> direction, I just have a truncated cross-section (it doesn't get any closer >>>> to the top of the cross-section). >>>> >>>> I don't see why using ff1 (from the interpolate) would make the >>>> labeling disappear. >>>> >>>> Variable: ff1_s >>>> >>>> Type: float >>>> >>>> Total Size: 33728 bytes >>>> >>>> 8432 values >>>> >>>> Number of Dimensions: 2 >>>> >>>> Dimensions and sizes: [height | 31] x [lon | 272] >>>> >>>> Coordinates: >>>> >>>> Number Of Attributes: 1 >>>> >>>> _FillValue : 9.96921e+36 >>>> >>>> >>>> Variable: ff1 >>>> >>>> Type: float >>>> >>>> Total Size: 33728 bytes >>>> >>>> 8432 values >>>> >>>> Number of Dimensions: 2 >>>> >>>> Dimensions and sizes: [HEIGHT | 31] x [lon | 272] >>>> >>>> Coordinates: >>>> >>>> HEIGHT: [500..15500] >>>> >>>> Number Of Attributes: 1 >>>> >>>> _FillValue : 9.96921e+36 >>>> >>>> -- >>>> Barry H. Lynn, Ph.D >>>> Senior Scientist, Lecturer, >>>> The Institute of Earth Sciences, >>>> The Hebrew University of Jerusalem, >>>> Givat Ram, Jerusalem 91904, Israel >>>> Tel: 972 547 231 170 >>>> Fax: (972)-25662581 >>>> >>>> Weather It Is, LTD >>>> Weather and Climate Focus >>>> https://weather-it-is.com >>>> Jerusalem, Israel >>>> Local: 02 930 9525 >>>> Cell: 054 7 231 170 >>>> Int-IS: x972 2 930 9525 >>>> >>> > > -- > > Barry H. Lynn, Ph.D > Senior Scientist, Lecturer, > The Institute of Earth Sciences, > The Hebrew University of Jerusalem, > Givat Ram, Jerusalem 91904, Israel > Tel: 972 547 231 170 > Fax: (972)-25662581 > > Weather It Is, LTD > Weather and Climate Focus > https://weather-it-is.com > Jerusalem, Israel > Local: 02 930 9525 > Cell: 054 7 231 170 > Int-IS: x972 2 930 9525 > > -- Barry H. Lynn, Ph.D Senior Scientist, Lecturer, The Institute of Earth Sciences, The Hebrew University of Jerusalem, Givat Ram, Jerusalem 91904, Israel Tel: 972 547 231 170 Fax: (972)-25662581 Weather It Is, LTD Weather and Climate Focus https://weather-it-is.com Jerusalem, Israel Local: 02 930 9525 Cell: 054 7 231 170 Int-IS: x972 2 930 9525 -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: radar_vert_crossection_10UTC.pdf Type: application/pdf Size: 58313 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: test2.ncl Type: application/octet-stream Size: 11759 bytes Desc: not available URL: From barry.h.lynn at gmail.com Tue Oct 31 00:36:02 2023 From: barry.h.lynn at gmail.com (Barry Lynn) Date: Tue, 31 Oct 2023 08:36:02 +0200 Subject: [ncl-talk] Strange problem related to linint1_Wrap In-Reply-To: References: Message-ID: Hello: Strangely, I was able to solve the problem with the extra white space this way (by limiting the vertical extent of the data passed into the gsm_contour_map program). *ff1_x = ff1(0:25,:)* plot = gsn_csm_contour(wks, ff1, res ) plots(i_file) = plot plot = gsn_csm_contour(wks, ff1_x, res ) plots(i_file+1) = plot I say strangely because reducing the vertical extent of the data read into the linear interpolation program just led to a cropping of the radar data. That makes no sense, as the map window clearly shows white space above (a). On Mon, Oct 30, 2023 at 8:21?PM Barry Lynn wrote: > Dear All (& Dave): > > I still have the strange problem where there is a lot of white space above > the mapped radar cross-section. > > If, however, I reduce the scope of the data (i_z is reduced), then the top > of the contours get cut off. > > So, I have a feeling that the gsn plotting program is assuming some kind > of non-linear height field, while my data is in height coordinates. > > Is this correct? If not, what might be the solution? Note, I'm not using > any WRF functions and I converted the WRF data to equally spaced data > levels. > > Thank you, > > Barry > > On Mon, Oct 30, 2023 at 6:12?PM Barry Lynn wrote: > >> Hi Dave: >> >> I thought those quotes were "suspicious." >> >> The labeling is the same regardless of which variable is called for the >> plot, so that wasn't the problem. >> >> *Regardless*, your solution worked! (delete (ff1&HEIGHT)). >> >> Thank you! >> >> On Mon, Oct 30, 2023 at 5:20?PM Dave Allured - NOAA Affiliate < >> dave.allured at noaa.gov> wrote: >> >>> Correction. Do not use quotes here: delete (ff1&HEIGHT) >>> >>> >>> On Mon, Oct 30, 2023 at 9:16?AM Dave Allured - NOAA Affiliate < >>> dave.allured at noaa.gov> wrote: >>> >>>> Barry, please notice that there is a HEIGHT coordinate variable >>>> attached to ff1 after interpolation, but there is no such coordinate >>>> variable attached to ff1_s before interpolation. It looks like the >>>> attached HEIGHT coordinate variable is misleading the plot function in some >>>> way. Try delete ("ff1&HEIGHT"). >>>> >>>> I am not familiar with the cross section plotting function. But it >>>> looks like in the first plot, you are setting some plot resources to set >>>> the range, spacing, and units of the height axis. I suspect that these >>>> settings are in conflict with the attached HEIGHT coordinates which may >>>> have been generated by linint1_Wrap. The units are different. Instead of >>>> deleting the HEIGHT coordinate variable, you could try deleting or >>>> adjusting the conflicting resource settings. >>>> >>>> If you want to rescale the attached HEIGHT coordinates from meters to >>>> km, then simply "ff1&HEIGHT = ff1&HEIGHT / 1000". >>>> >>>> >>>> On Mon, Oct 30, 2023 at 8:08?AM Barry Lynn via ncl-talk < >>>> ncl-talk at mailman.ucar.edu> wrote: >>>> >>>>> Hello: >>>>> >>>>> I wanted to use linint1 to interpolate WRF data to height surfaces >>>>> with constant spacing between levels. But, there are two mysteries. >>>>> >>>>> You can see the WRF data in a cross-section of radar in (a; the >>>>> variable ff1_s below). When I interpolate the data (b, ff1) is compacted >>>>> towards the bottom of the cross-sectional plot. Moreover, the height >>>>> information on the y-axis will not plot. >>>>> >>>>> I understand why the data is compacted, but not why I can't eliminate >>>>> the white space above. If I reduce the range of WRF data sampled in the z >>>>> direction, I just have a truncated cross-section (it doesn't get any closer >>>>> to the top of the cross-section). >>>>> >>>>> I don't see why using ff1 (from the interpolate) would make the >>>>> labeling disappear. >>>>> >>>>> Variable: ff1_s >>>>> >>>>> Type: float >>>>> >>>>> Total Size: 33728 bytes >>>>> >>>>> 8432 values >>>>> >>>>> Number of Dimensions: 2 >>>>> >>>>> Dimensions and sizes: [height | 31] x [lon | 272] >>>>> >>>>> Coordinates: >>>>> >>>>> Number Of Attributes: 1 >>>>> >>>>> _FillValue : 9.96921e+36 >>>>> >>>>> >>>>> Variable: ff1 >>>>> >>>>> Type: float >>>>> >>>>> Total Size: 33728 bytes >>>>> >>>>> 8432 values >>>>> >>>>> Number of Dimensions: 2 >>>>> >>>>> Dimensions and sizes: [HEIGHT | 31] x [lon | 272] >>>>> >>>>> Coordinates: >>>>> >>>>> HEIGHT: [500..15500] >>>>> >>>>> Number Of Attributes: 1 >>>>> >>>>> _FillValue : 9.96921e+36 >>>>> >>>>> -- >>>>> Barry H. Lynn, Ph.D >>>>> Senior Scientist, Lecturer, >>>>> The Institute of Earth Sciences, >>>>> The Hebrew University of Jerusalem, >>>>> Givat Ram, Jerusalem 91904, Israel >>>>> Tel: 972 547 231 170 >>>>> Fax: (972)-25662581 >>>>> >>>>> Weather It Is, LTD >>>>> Weather and Climate Focus >>>>> https://weather-it-is.com >>>>> Jerusalem, Israel >>>>> Local: 02 930 9525 >>>>> Cell: 054 7 231 170 >>>>> Int-IS: x972 2 930 9525 >>>>> >>>> >> >> -- >> >> Barry H. Lynn, Ph.D >> Senior Scientist, Lecturer, >> The Institute of Earth Sciences, >> The Hebrew University of Jerusalem, >> Givat Ram, Jerusalem 91904, Israel >> Tel: 972 547 231 170 >> Fax: (972)-25662581 >> >> Weather It Is, LTD >> Weather and Climate Focus >> https://weather-it-is.com >> Jerusalem, Israel >> Local: 02 930 9525 >> Cell: 054 7 231 170 >> Int-IS: x972 2 930 9525 >> >> > > -- > > Barry H. Lynn, Ph.D > Senior Scientist, Lecturer, > The Institute of Earth Sciences, > The Hebrew University of Jerusalem, > Givat Ram, Jerusalem 91904, Israel > Tel: 972 547 231 170 > Fax: (972)-25662581 > > Weather It Is, LTD > Weather and Climate Focus > https://weather-it-is.com > Jerusalem, Israel > Local: 02 930 9525 > Cell: 054 7 231 170 > Int-IS: x972 2 930 9525 > > -- Barry H. Lynn, Ph.D Senior Scientist, Lecturer, The Institute of Earth Sciences, The Hebrew University of Jerusalem, Givat Ram, Jerusalem 91904, Israel Tel: 972 547 231 170 Fax: (972)-25662581 Weather It Is, LTD Weather and Climate Focus https://weather-it-is.com Jerusalem, Israel Local: 02 930 9525 Cell: 054 7 231 170 Int-IS: x972 2 930 9525 -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: radar_vert_crossection_10UTC.pdf Type: application/pdf Size: 84958 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: test3.ncl Type: application/octet-stream Size: 11521 bytes Desc: not available URL: