From ehsantaghizadeh at yahoo.com Tue Apr 1 23:50:25 2025 From: ehsantaghizadeh at yahoo.com (Ehsan Taghizadeh) Date: Wed, 2 Apr 2025 05:50:25 +0000 (UTC) Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming References: <1920293968.370628.1743573026010.ref@mail.yahoo.com> Message-ID: <1920293968.370628.1743573026010@mail.yahoo.com> Dear NCL team, I would like to plot 10 m wind vectors, 2 m temperature as contours, and topography as shaded from a WRF output file. Additionally, I need to zoom in on a domain approximately 1 ? 1 degree. For instance, I would like to display one wind vector for every five grid points. I reviewed the provided page (https://www.ncl.ucar.edu/Applications/wrfzoom.shtml) and tried the attached script. While it is relatively clean, it has some issues?particularly with the labels on the axes, as shown. Any assistance in creating such a script would be greatly appreciated. ----------------------------------------------------------------------------- Best regards,Ehsan TaghizadehPh.D. Student in Natural Resources & Environmental Studies (NRES),????University of Northern British Columbia (UNBC) -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: wind10_temp1.png Type: image/png Size: 177758 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: wind10_temp2.png Type: image/png Size: 150126 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: WRF_wind10_tempV2.ncl Type: application/octet-stream Size: 9598 bytes Desc: not available URL: From kemal.gurer at arb.ca.gov Wed Apr 2 14:33:40 2025 From: kemal.gurer at arb.ca.gov (Gurer, Kemal@ARB) Date: Wed, 2 Apr 2025 20:33:40 +0000 Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming In-Reply-To: <1920293968.370628.1743573026010@mail.yahoo.com> References: <1920293968.370628.1743573026010.ref@mail.yahoo.com> <1920293968.370628.1743573026010@mail.yahoo.com> Message-ID: Hi Ehsan, Try the attached ncl script to plot horizontal cross sections of several variables from your WRF output. Just read the command lines and posted comments to understand how the script operates and disregard the warnings that ncl issues. You can plot other variables by looking at the examples given in the script. Kemal. From: ncl-talk On Behalf Of Ehsan Taghizadeh via ncl-talk Sent: Tuesday, April 1, 2025 10:50 PM To: Ncl-talk Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Dear NCL team, I would like to plot 10 m wind vectors, 2 m temperature as contours, and topography as shaded from a WRF output file. Additionally, I need to zoom in on a domain approximately 1 ? 1 degree. For instance, I would like to display one wind vector for every five grid points. I reviewed the provided page (https://www.ncl.ucar.edu/Applications/wrfzoom.shtml) and tried the attached script. While it is relatively clean, it has some issues-particularly with the labels on the axes, as shown. Any assistance in creating such a script would be greatly appreciated. ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh Ph.D. Student in Natural Resources & Environmental Studies (NRES), University of Northern British Columbia (UNBC) -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: wrf.hrz.cross.ncl Type: application/octet-stream Size: 48454 bytes Desc: wrf.hrz.cross.ncl URL: From ehsantaghizadeh at yahoo.com Sun Apr 6 20:42:02 2025 From: ehsantaghizadeh at yahoo.com (Ehsan Taghizadeh) Date: Mon, 7 Apr 2025 02:42:02 +0000 (UTC) Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming In-Reply-To: References: <1920293968.370628.1743573026010.ref@mail.yahoo.com> <1920293968.370628.1743573026010@mail.yahoo.com> Message-ID: <1411930447.1375295.1743993722489@mail.yahoo.com> Hi Kemal, Thank you so much for sharing your script! It looks great and very well-commented, making it easier to understand its functionality. I truly appreciate the effort you've put into creating and providing this. I'm planning to use your script to plot WRF data for my research. Just to clarify, my main challenge involves zooming in on a specific region, superimposing three variables, and plotting the wind vector for a specific interval of grid points. I hope it's okay to ask for a bit more guidance?specifically on how to modify the dataonmap_zoom_10.ncl script to meet these requirements. If anyone has any insights on achieving this or can suggest adjustments, that would be incredibly helpful. Additionally, if any GeoCAT or Python scripts might suit these needs, I'd greatly appreciate your recommendations. Any further assistance would mean a lot to me. ----------------------------------------------------------------------------- Best regards,Ehsan Taghizadeh On Wednesday, April 2, 2025 at 08:06:35 PM PDT, Gurer, Kemal at ARB wrote: Hi Ehsan, ? Try the attached ncl script to plot horizontal cross sections of several variables from your WRF output. Just read the command lines and posted comments to understand how the script operates and disregard the warnings that ncl issues. You can plot other variables by looking at the examples given in the script. ? Kemal. ? From: ncl-talk On Behalf Of Ehsan Taghizadeh via ncl-talk Sent: Tuesday, April 1, 2025 10:50 PM To: Ncl-talk Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming ? CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Dear NCL team, ? I would like to plot 10 m wind vectors, 2 m temperature as contours, and topography as shaded from a WRF output file. Additionally, I need to zoom in on a domain approximately 1 ? 1 degree. For instance, I would like to display one wind vector for every five grid points. I reviewed the provided page (https://www.ncl.ucar.edu/Applications/wrfzoom.shtml) and tried the attached script. While it is relatively clean, it has some issues?particularly with the labels on the axes, as shown. Any assistance in creating such a script would be greatly appreciated. ? ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh Ph.D. Student in Natural Resources & Environmental Studies (NRES), ????University of Northern British Columbia (UNBC) -------------- next part -------------- An HTML attachment was scrubbed... URL: From kemal.gurer at arb.ca.gov Mon Apr 7 11:26:16 2025 From: kemal.gurer at arb.ca.gov (Gurer, Kemal@ARB) Date: Mon, 7 Apr 2025 17:26:16 +0000 Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming In-Reply-To: <1411930447.1375295.1743993722489@mail.yahoo.com> References: <1920293968.370628.1743573026010.ref@mail.yahoo.com> <1920293968.370628.1743573026010@mail.yahoo.com> <1411930447.1375295.1743993722489@mail.yahoo.com> Message-ID: Hi Ehsan, There is a transition from ncl to python usage to plot variables from WRF outputs and, for example, you can get the routine from here: https://github.com/NCAR/wrf-python I have certainly experimented with python to plot WRF outputs. However, since I have been using ncl for quite some time, I feel comfortable tweaking the ncl scripts to get the plots that I want. In the ncl script that I gave you, you set the lat/lon values of the corners of your plot area using the following variables for all variables: vecres1 at mpLeftCornerLatF = xlatmin vecres1 at mpRightCornerLatF = xlatmax vecres1 at mpLeftCornerLonF = xlonmin vecres1 at mpRightCornerLonF = xlonmax to zoom into the area that you want to focus on. Experiment with the script and the rest will be a breeze for you. Kemal. From: Ehsan Taghizadeh Sent: Sunday, April 6, 2025 7:42 PM To: Ncl-talk ; Gurer, Kemal at ARB Subject: Re: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Hi Kemal, Thank you so much for sharing your script! It looks great and very well-commented, making it easier to understand its functionality. I truly appreciate the effort you've put into creating and providing this. I'm planning to use your script to plot WRF data for my research. Just to clarify, my main challenge involves zooming in on a specific region, superimposing three variables, and plotting the wind vector for a specific interval of grid points. I hope it's okay to ask for a bit more guidance-specifically on how to modify the dataonmap_zoom_10.ncl script to meet these requirements. If anyone has any insights on achieving this or can suggest adjustments, that would be incredibly helpful. Additionally, if any GeoCAT or Python scripts might suit these needs, I'd greatly appreciate your recommendations. Any further assistance would mean a lot to me. ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh On Wednesday, April 2, 2025 at 08:06:35 PM PDT, Gurer, Kemal at ARB > wrote: Hi Ehsan, Try the attached ncl script to plot horizontal cross sections of several variables from your WRF output. Just read the command lines and posted comments to understand how the script operates and disregard the warnings that ncl issues. You can plot other variables by looking at the examples given in the script. Kemal. From: ncl-talk > On Behalf Of Ehsan Taghizadeh via ncl-talk Sent: Tuesday, April 1, 2025 10:50 PM To: Ncl-talk > Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Dear NCL team, I would like to plot 10 m wind vectors, 2 m temperature as contours, and topography as shaded from a WRF output file. Additionally, I need to zoom in on a domain approximately 1 ? 1 degree. For instance, I would like to display one wind vector for every five grid points. I reviewed the provided page (https://www.ncl.ucar.edu/Applications/wrfzoom.shtml) and tried the attached script. While it is relatively clean, it has some issues-particularly with the labels on the axes, as shown. Any assistance in creating such a script would be greatly appreciated. ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh Ph.D. Student in Natural Resources & Environmental Studies (NRES), University of Northern British Columbia (UNBC) -------------- next part -------------- An HTML attachment was scrubbed... URL: From kemal.gurer at arb.ca.gov Mon Apr 7 11:37:34 2025 From: kemal.gurer at arb.ca.gov (Gurer, Kemal@ARB) Date: Mon, 7 Apr 2025 17:37:34 +0000 Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming In-Reply-To: References: <1920293968.370628.1743573026010.ref@mail.yahoo.com> <1920293968.370628.1743573026010@mail.yahoo.com> <1411930447.1375295.1743993722489@mail.yahoo.com> Message-ID: I forgot to address one specific question that you had asked, which is you want to plot 10m wind vectors at every 5th grid points. In ncl, you can specify that via, for example, U = u10(it,::5,::5) where the dimensions of u are: [Time x [south_north] x [west_east] and 5 is the incremental value that you want to plot every 5th grid points. Kemal. From: ncl-talk On Behalf Of Gurer, Kemal at ARB via ncl-talk Sent: Monday, April 7, 2025 10:26 AM To: Ehsan Taghizadeh ; Ncl-talk Subject: Re: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Hi Ehsan, There is a transition from ncl to python usage to plot variables from WRF outputs and, for example, you can get the routine from here: https://github.com/NCAR/wrf-python I have certainly experimented with python to plot WRF outputs. However, since I have been using ncl for quite some time, I feel comfortable tweaking the ncl scripts to get the plots that I want. In the ncl script that I gave you, you set the lat/lon values of the corners of your plot area using the following variables for all variables: vecres1 at mpLeftCornerLatF = xlatmin vecres1 at mpRightCornerLatF = xlatmax vecres1 at mpLeftCornerLonF = xlonmin vecres1 at mpRightCornerLonF = xlonmax to zoom into the area that you want to focus on. Experiment with the script and the rest will be a breeze for you. Kemal. From: Ehsan Taghizadeh > Sent: Sunday, April 6, 2025 7:42 PM To: Ncl-talk >; Gurer, Kemal at ARB > Subject: Re: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Hi Kemal, Thank you so much for sharing your script! It looks great and very well-commented, making it easier to understand its functionality. I truly appreciate the effort you've put into creating and providing this. I'm planning to use your script to plot WRF data for my research. Just to clarify, my main challenge involves zooming in on a specific region, superimposing three variables, and plotting the wind vector for a specific interval of grid points. I hope it's okay to ask for a bit more guidance-specifically on how to modify the dataonmap_zoom_10.ncl script to meet these requirements. If anyone has any insights on achieving this or can suggest adjustments, that would be incredibly helpful. Additionally, if any GeoCAT or Python scripts might suit these needs, I'd greatly appreciate your recommendations. Any further assistance would mean a lot to me. ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh On Wednesday, April 2, 2025 at 08:06:35 PM PDT, Gurer, Kemal at ARB > wrote: Hi Ehsan, Try the attached ncl script to plot horizontal cross sections of several variables from your WRF output. Just read the command lines and posted comments to understand how the script operates and disregard the warnings that ncl issues. You can plot other variables by looking at the examples given in the script. Kemal. From: ncl-talk > On Behalf Of Ehsan Taghizadeh via ncl-talk Sent: Tuesday, April 1, 2025 10:50 PM To: Ncl-talk > Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Dear NCL team, I would like to plot 10 m wind vectors, 2 m temperature as contours, and topography as shaded from a WRF output file. Additionally, I need to zoom in on a domain approximately 1 ? 1 degree. For instance, I would like to display one wind vector for every five grid points. I reviewed the provided page (https://www.ncl.ucar.edu/Applications/wrfzoom.shtml) and tried the attached script. While it is relatively clean, it has some issues-particularly with the labels on the axes, as shown. Any assistance in creating such a script would be greatly appreciated. ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh Ph.D. Student in Natural Resources & Environmental Studies (NRES), University of Northern British Columbia (UNBC) -------------- next part -------------- An HTML attachment was scrubbed... URL: From ehsantaghizadeh at yahoo.com Mon Apr 7 18:08:59 2025 From: ehsantaghizadeh at yahoo.com (Ehsan Taghizadeh) Date: Tue, 8 Apr 2025 00:08:59 +0000 (UTC) Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming In-Reply-To: References: <1920293968.370628.1743573026010.ref@mail.yahoo.com> <1920293968.370628.1743573026010@mail.yahoo.com> <1411930447.1375295.1743993722489@mail.yahoo.com> Message-ID: <1186412743.1614622.1744070939984@mail.yahoo.com> Hi Kemal (and all), Thank you so much for your kind help. Your guidance has been incredibly helpful to me. For more clarification, I?d like to elaborate on a couple of points. Regarding zooming in and plotting wind vectors at every 5th grid point, my goal is to obtain a better resolution to analyze wind direction in more detail. By selecting every 5th grid point as a sample, I aim to control the density of wind vectors on the plot to achieve the most effective visualization. Additionally, I?d like to plot wind vectors for every 5th grid point as an example, while still including all grid points of the WRF domain on the plot. This approach would help balance clarity and detail in the visualization. When I refer to "zooming in," I don't mean simply isolating a portion of the map or plot. Rather, I?m looking to enhance resolution to reveal finer details, increase the grid points displayed, and present more specific variable information. I tried to illustrate more clearly in the following hypothetical image what I mean by plotting every grid point while plotting wind vectors for every 5th grid point.Any assistance would be appreciated. ----------------------------------------------------------------------------- Best regards,Ehsan Taghizadeh On Monday, April 7, 2025 at 11:11:28 AM PDT, Gurer, Kemal at ARB wrote: I forgot to address one specific question that you had asked, which is you want to plot 10m wind vectors at every 5th grid points. In ncl, you can specify that via, for example, ? ? U = u10(it,::5,::5) ? where the dimensions of u are: ? ? [Time x [south_north] ?x [west_east] ? and 5 is the incremental value that you want to plot every 5th grid points. ? Kemal. ? From: ncl-talk On Behalf Of Gurer, Kemal at ARB via ncl-talk Sent: Monday, April 7, 2025 10:26 AM To: Ehsan Taghizadeh ; Ncl-talk Subject: Re: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming ? CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Hi Ehsan, ? There is a transition from ncl to python usage to plot variables from WRF outputs and, for example, you can get the routine from here: ? https://github.com/NCAR/wrf-python ? I have certainly experimented with python to plot WRF outputs. However, since I have been using ncl for quite some time, I feel comfortable tweaking the ncl scripts to get the plots that I want. In the ncl script that I gave you, you set the lat/lon values of the corners of your plot area using the following variables for all variables: ? ????? vecres1 at mpLeftCornerLatF? = xlatmin ????? vecres1 at mpRightCornerLatF = xlatmax ????? vecres1 at mpLeftCornerLonF? = xlonmin ????? vecres1 at mpRightCornerLonF = xlonmax ? to zoom into the area that you want to focus on. Experiment with the script and the rest will be a breeze for you. ? Kemal. ? From: Ehsan Taghizadeh Sent: Sunday, April 6, 2025 7:42 PM To: Ncl-talk ; Gurer, Kemal at ARB Subject: Re: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming ? CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Hi Kemal, ? Thank you so much for sharing your script! It looks great and very well-commented, making it easier to understand its functionality. I truly appreciate the effort you've put into creating and providing this. I'm planning to use your script to plot WRF data for my research. ? Just to clarify, my main challenge involves zooming in on a specific region, superimposing three variables, and plotting the wind vector for a specific interval of grid points. I hope it's okay to ask for a bit more guidance?specifically on how to modify thedataonmap_zoom_10.ncl script to meet these requirements. If anyone has any insights on achieving this or can suggest adjustments, that would be incredibly helpful. Additionally, if any GeoCAT or Python scripts might suit these needs, I'd greatly appreciate your recommendations. ? Any further assistance would mean a lot to me. ? ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh ? ? On Wednesday, April 2, 2025 at 08:06:35 PM PDT, Gurer, Kemal at ARB wrote: ? ? Hi Ehsan, ? Try the attached ncl script to plot horizontal cross sections of several variables from your WRF output. Just read the command lines and posted comments to understand how the script operates and disregard the warnings that ncl issues. You can plot other variables by looking at the examples given in the script. ? Kemal. ? From: ncl-talk On Behalf Of Ehsan Taghizadeh via ncl-talk Sent: Tuesday, April 1, 2025 10:50 PM To: Ncl-talk Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming ? CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Dear NCL team, ? I would like to plot 10 m wind vectors, 2 m temperature as contours, and topography as shaded from a WRF output file. Additionally, I need to zoom in on a domain approximately 1 ? 1 degree. For instance, I would like to display one wind vector for every five grid points. I reviewed the provided page (https://www.ncl.ucar.edu/Applications/wrfzoom.shtml) and tried the attached script. While it is relatively clean, it has some issues?particularly with the labels on the axes, as shown. Any assistance in creating such a script would be greatly appreciated. ? ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh Ph.D. Student in Natural Resources & Environmental Studies (NRES), ????University of Northern British Columbia (UNBC) -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 1744070672197blob.jpg Type: image/png Size: 9916 bytes Desc: not available URL: From kemal.gurer at arb.ca.gov Tue Apr 8 08:14:48 2025 From: kemal.gurer at arb.ca.gov (Gurer, Kemal@ARB) Date: Tue, 8 Apr 2025 14:14:48 +0000 Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming In-Reply-To: <1186412743.1614622.1744070939984@mail.yahoo.com> References: <1920293968.370628.1743573026010.ref@mail.yahoo.com> <1920293968.370628.1743573026010@mail.yahoo.com> <1411930447.1375295.1743993722489@mail.yahoo.com> <1186412743.1614622.1744070939984@mail.yahoo.com> Message-ID: Hi Ehsan, In this case, you need to overlay (or underlay grid points on the top of wind vectors. Here is an example script that plots grid points on the top of contours: https://www.ncl.ucar.edu/Applications/Scripts/shapefiles_16.ncl from the main location of other similar scripts that plot contours, data points as well as grid points, grid boxes, etc: https://www.ncl.ucar.edu/Applications/datagrid.shtml >From the script, you can take the portion below to add onto the script that I gave you: ;---Add dots at the lat/lon grid locations mkres = True mkres at gsMarkerIndex = 16 ; Filled dots mkres at gsMarkerSizeF = 0.001 ; Make them small mkres at gsMarkerColor = "darkorchid4" mkres at gsnCoordsAttach = True gsn_coordinates(wks,plot_orig_coarse,data_coarse,mkres) gsn_coordinates(wks,plot_mask_coarse,data_mask_coarse,mkres) gsn_coordinates(wks,plot_orig_dense,data_dense,mkres) gsn_coordinates(wks,plot_mask_dense,data_mask_dense,mkres) Kemal. From: Ehsan Taghizadeh Sent: Monday, April 7, 2025 5:09 PM To: Gurer, Kemal at ARB ; Ncl-talk Subject: Re: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Hi Kemal (and all), Thank you so much for your kind help. Your guidance has been incredibly helpful to me. For more clarification, I'd like to elaborate on a couple of points. Regarding zooming in and plotting wind vectors at every 5th grid point, my goal is to obtain a better resolution to analyze wind direction in more detail. By selecting every 5th grid point as a sample, I aim to control the density of wind vectors on the plot to achieve the most effective visualization. Additionally, I'd like to plot wind vectors for every 5th grid point as an example, while still including all grid points of the WRF domain on the plot. This approach would help balance clarity and detail in the visualization. When I refer to "zooming in," I don't mean simply isolating a portion of the map or plot. Rather, I'm looking to enhance resolution to reveal finer details, increase the grid points displayed, and present more specific variable information. I tried to illustrate more clearly in the following hypothetical image what I mean by plotting every grid point while plotting wind vectors for every 5th grid point. Any assistance would be appreciated. [Inline image] ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh On Monday, April 7, 2025 at 11:11:28 AM PDT, Gurer, Kemal at ARB > wrote: I forgot to address one specific question that you had asked, which is you want to plot 10m wind vectors at every 5th grid points. In ncl, you can specify that via, for example, U = u10(it,::5,::5) where the dimensions of u are: [Time x [south_north] x [west_east] and 5 is the incremental value that you want to plot every 5th grid points. Kemal. From: ncl-talk > On Behalf Of Gurer, Kemal at ARB via ncl-talk Sent: Monday, April 7, 2025 10:26 AM To: Ehsan Taghizadeh >; Ncl-talk > Subject: Re: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Hi Ehsan, There is a transition from ncl to python usage to plot variables from WRF outputs and, for example, you can get the routine from here: https://github.com/NCAR/wrf-python I have certainly experimented with python to plot WRF outputs. However, since I have been using ncl for quite some time, I feel comfortable tweaking the ncl scripts to get the plots that I want. In the ncl script that I gave you, you set the lat/lon values of the corners of your plot area using the following variables for all variables: vecres1 at mpLeftCornerLatF = xlatmin vecres1 at mpRightCornerLatF = xlatmax vecres1 at mpLeftCornerLonF = xlonmin vecres1 at mpRightCornerLonF = xlonmax to zoom into the area that you want to focus on. Experiment with the script and the rest will be a breeze for you. Kemal. From: Ehsan Taghizadeh > Sent: Sunday, April 6, 2025 7:42 PM To: Ncl-talk >; Gurer, Kemal at ARB > Subject: Re: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Hi Kemal, Thank you so much for sharing your script! It looks great and very well-commented, making it easier to understand its functionality. I truly appreciate the effort you've put into creating and providing this. I'm planning to use your script to plot WRF data for my research. Just to clarify, my main challenge involves zooming in on a specific region, superimposing three variables, and plotting the wind vector for a specific interval of grid points. I hope it's okay to ask for a bit more guidance-specifically on how to modify the dataonmap_zoom_10.ncl script to meet these requirements. If anyone has any insights on achieving this or can suggest adjustments, that would be incredibly helpful. Additionally, if any GeoCAT or Python scripts might suit these needs, I'd greatly appreciate your recommendations. Any further assistance would mean a lot to me. ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh On Wednesday, April 2, 2025 at 08:06:35 PM PDT, Gurer, Kemal at ARB > wrote: Hi Ehsan, Try the attached ncl script to plot horizontal cross sections of several variables from your WRF output. Just read the command lines and posted comments to understand how the script operates and disregard the warnings that ncl issues. You can plot other variables by looking at the examples given in the script. Kemal. From: ncl-talk > On Behalf Of Ehsan Taghizadeh via ncl-talk Sent: Tuesday, April 1, 2025 10:50 PM To: Ncl-talk > Subject: [ncl-talk] Visualization of WRF Output: Superimposition and Zooming CAUTION: This email originated from outside of the organization. Do not click links or open attachments unless you recognize the sender and know the content is safe. Dear NCL team, I would like to plot 10 m wind vectors, 2 m temperature as contours, and topography as shaded from a WRF output file. Additionally, I need to zoom in on a domain approximately 1 ? 1 degree. For instance, I would like to display one wind vector for every five grid points. I reviewed the provided page (https://www.ncl.ucar.edu/Applications/wrfzoom.shtml) and tried the attached script. While it is relatively clean, it has some issues-particularly with the labels on the axes, as shown. Any assistance in creating such a script would be greatly appreciated. ----------------------------------------------------------------------------- Best regards, Ehsan Taghizadeh Ph.D. Student in Natural Resources & Environmental Studies (NRES), University of Northern British Columbia (UNBC) -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image001.png Type: image/png Size: 9916 bytes Desc: image001.png URL: From missquaddus at gmail.com Thu Apr 17 02:32:40 2025 From: missquaddus at gmail.com (qudsia zafar) Date: Thu, 17 Apr 2025 13:32:40 +0500 Subject: [ncl-talk] SPEI (drought index) calculation... Message-ID: Hi i am using cru datasets (ts 4.08) for precipitation and potential evapotranspiration. I have done (precip_minus_pet = prc - pet) and then I have used the function 'dim_gamfit_n' to estimate shape and scale parameters for the gamma distribution fitting. However in the next step of spei calculation I need to do the following: spei = gmm_transform(precip_minus_pet, shape, scale) I have now all the three data (precip_minus_pet, shape, scale) but i don't know if NCL has a gamma transformation function available? I have seen 'cdfgam_p' but I don't know if this corresponds to the gamma transform?. Kindly advise which function can be used in the final step of spei calculation. Many Thx.. qudsia -------------- next part -------------- An HTML attachment was scrubbed... URL: From shea at ucar.edu Thu Apr 17 19:40:52 2025 From: shea at ucar.edu (Dennis Shea) Date: Thu, 17 Apr 2025 19:40:52 -0600 Subject: [ncl-talk] SPEI (drought index) calculation... In-Reply-To: References: Message-ID: NCL's *cdfgam_p *documentation mentions 'gamma' functions from R, Matlab and Excel. What is your 'gmm_transform' supposed to return? On Thu, Apr 17, 2025 at 2:32?AM qudsia zafar via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Hi i am using cru datasets (ts 4.08) for precipitation and potential > evapotranspiration. I have done (precip_minus_pet = prc - pet) and then > I have used the function 'dim_gamfit_n' to estimate shape and scale > parameters for the gamma distribution fitting. > > However in the next step of spei calculation I need to do the following: > > spei = gmm_transform(precip_minus_pet, shape, scale) > > I have now all the three data (precip_minus_pet, shape, scale) but i > don't know if NCL has a gamma transformation function available? I have > seen 'cdfgam_p' but I don't know if this corresponds to the gamma > transform?. Kindly advise which function can be used in the final step of > spei calculation. Many Thx.. > > qudsia > > > > > > > _______________________________________________ > 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: From missquaddus at gmail.com Thu Apr 17 23:59:08 2025 From: missquaddus at gmail.com (qudsia zafar) Date: Fri, 18 Apr 2025 10:59:08 +0500 Subject: [ncl-talk] SPEI (drought index) calculation... In-Reply-To: References: Message-ID: Hi Dennis, As I couldn't find the direct function for spei calculation just like there is for spi in NCL therefore i was trying to do things manually. Till now I have used the 'dim_gamfit_n' function output (shape and scale parameters) as an input in the *cdfgam_p *function. Now the probabilities (gridded data) acquired have to be transformed into the 'standard normal distribution' resulting in the spei values with a mean of 0 and standard deviation of 1. An spei of 0 indicates near normal conditions, positive values indicate wetter than average conditions and negative values are drier than average conditions. I was wondering if there is a function that transforms from cdfgam_p to the standard normal distribution. As you have asked 'gmm_transform' is supposed to acquire the necessary transformation to 'standard normal distribution'. If you could have a look. Thx. qudsia - On Fri, 18 Apr 2025 at 06:41, Dennis Shea wrote: > NCL's *cdfgam_p > *documentation > mentions 'gamma' functions from R, Matlab and Excel. > > What is your 'gmm_transform' supposed to return? > > > > > > > > > > On Thu, Apr 17, 2025 at 2:32?AM qudsia zafar via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> Hi i am using cru datasets (ts 4.08) for precipitation and potential >> evapotranspiration. I have done (precip_minus_pet = prc - pet) and then >> I have used the function 'dim_gamfit_n' to estimate shape and scale >> parameters for the gamma distribution fitting. >> >> However in the next step of spei calculation I need to do the following: >> >> spei = gmm_transform(precip_minus_pet, shape, scale) >> >> I have now all the three data (precip_minus_pet, shape, scale) but i >> don't know if NCL has a gamma transformation function available? I have >> seen 'cdfgam_p' but I don't know if this corresponds to the gamma >> transform?. Kindly advise which function can be used in the final step of >> spei calculation. Many Thx.. >> >> qudsia >> >> >> >> >> >> >> _______________________________________________ >> 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: From shea at ucar.edu Fri Apr 25 17:36:57 2025 From: shea at ucar.edu (Dennis Shea) Date: Fri, 25 Apr 2025 17:36:57 -0600 Subject: [ncl-talk] SPEI (drought index) calculation... In-Reply-To: References:

Message-ID: Hello Qudsia re: "... if there is a function that transforms from cdfgam_p to the standard normal distribution. As you have asked 'gmm_transform' is supposed to acquire the necessary transformation to 'standard normal distribution'. If you could have a look." I could not find any function to transform from cdfgam_p to the standard normal distribution. Sorry. On Thu, Apr 17, 2025 at 11:59?PM qudsia zafar wrote: > Hi Dennis, > > As I couldn't find the direct function for spei calculation just like > there is for spi in NCL therefore i was trying to do things manually. > > Till now I have used the 'dim_gamfit_n' function output (shape and scale > parameters) as an input in the *cdfgam_p > *function. > Now the probabilities (gridded data) acquired have to be transformed into > the 'standard normal distribution' resulting in the spei values with a mean > of 0 and standard deviation of 1. > > An spei of 0 indicates near normal conditions, positive values indicate > wetter than average conditions and negative values are drier than average > conditions. > > I was wondering if there is a function that transforms from cdfgam_p to > the standard normal distribution. As you have asked 'gmm_transform' is > supposed to acquire the necessary transformation to 'standard normal > distribution'. If you could have a look. > > Thx. > > qudsia > > > > - > > > > On Fri, 18 Apr 2025 at 06:41, Dennis Shea wrote: > >> NCL's *cdfgam_p >> *documentation >> mentions 'gamma' functions from R, Matlab and Excel. >> >> What is your 'gmm_transform' supposed to return? >> >> >> >> >> >> >> >> >> >> On Thu, Apr 17, 2025 at 2:32?AM qudsia zafar via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >>> Hi i am using cru datasets (ts 4.08) for precipitation and potential >>> evapotranspiration. I have done (precip_minus_pet = prc - pet) and then >>> I have used the function 'dim_gamfit_n' to estimate shape and scale >>> parameters for the gamma distribution fitting. >>> >>> However in the next step of spei calculation I need to do the following: >>> >>> spei = gmm_transform(precip_minus_pet, shape, scale) >>> >>> I have now all the three data (precip_minus_pet, shape, scale) but i >>> don't know if NCL has a gamma transformation function available? I have >>> seen 'cdfgam_p' but I don't know if this corresponds to the gamma >>> transform?. Kindly advise which function can be used in the final step of >>> spei calculation. Many Thx.. >>> >>> qudsia >>> >>> >>> >>> >>> >>> >>> _______________________________________________ >>> 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: