From setareh.rahimi at gmail.com Wed Apr 3 04:32:32 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Wed, 3 Apr 2024 14:02:32 +0330 Subject: [ncl-talk] *******Position of shape-file move to the right by narrowing the lat/lon******* Message-ID: Dear NCL users, I tried to plot the total precipitation for 2021-03_12 (24 hours), using the attached script. However, once I tried to zoom in on the study area to better show the result by shortening the lat/lon range I noticed the position of the shapefile on the plot was not correct ( it moved to the right). How can I fix the shape file to the right position? ( please check the attached images). please help me to sole this issue. Thanks. -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 3.png Type: image/png Size: 116747 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 2.png Type: image/png Size: 172204 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: prc Type: application/octet-stream Size: 3444 bytes Desc: not available URL: From asphilli at ucar.edu Fri Apr 5 17:54:53 2024 From: asphilli at ucar.edu (Adam Phillips) Date: Fri, 5 Apr 2024 17:54:53 -0600 Subject: [ncl-talk] *******Position of shape-file move to the right by narrowing the lat/lon******* In-Reply-To: References: Message-ID: Hi Setareh, See example #2 here for how to zoom in on a plot when you are using wrf_map_resources: https://www.ncl.ucar.edu/Applications/wrfzoom.shtml#ex2 Best, Adam On Wed, Apr 3, 2024 at 4:26?AM Setareh Rahimi via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear NCL users, > > I tried to plot the total precipitation for 2021-03_12 (24 hours), using > the attached script. However, once I tried to zoom in on the study area to > better show the result by shortening the lat/lon range I noticed the > position of the shapefile on the plot was not correct ( it moved to the > right). How can I fix the shape file to the right position? ( please check > the attached images). > > please help me to sole this issue. > Thanks. > > -- > 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 > -- 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 setareh.rahimi at gmail.com Sat Apr 6 05:20:51 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Sat, 6 Apr 2024 14:50:51 +0330 Subject: [ncl-talk] *******Position of shape-file move to the right by narrowing the lat/lon******* In-Reply-To: References: Message-ID: Dear Adam, Thank you so much for your suggestion. I added the following lines (highlighted) to my script but I still have the same problem. res at mpLimitMode = "LatLon" ; choose region of map *res at ZoomIn = True* res at mpMaxLatF = 40 ; choose subregion res at mpMinLatF = 30 res at mpMaxLonF = 61 res at mpMinLonF = 42 r*es = wrf_map_resources(f,res) ** ; must call this again after setting ZoomIn/X/Ystart/end* please check the attached and please advise me where I made the mistakes. Thanks, On Sat, Apr 6, 2024 at 3:25?AM Adam Phillips wrote: > Hi Setareh, > See example #2 here for how to zoom in on a plot when you are > using wrf_map_resources: > https://www.ncl.ucar.edu/Applications/wrfzoom.shtml#ex2 > Best, > Adam > > On Wed, Apr 3, 2024 at 4:26?AM Setareh Rahimi via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> Dear NCL users, >> >> I tried to plot the total precipitation for 2021-03_12 (24 hours), using >> the attached script. However, once I tried to zoom in on the study area to >> better show the result by shortening the lat/lon range I noticed the >> position of the shapefile on the plot was not correct ( it moved to the >> right). How can I fix the shape file to the right position? ( please check >> the attached images). >> >> please help me to sole this issue. >> Thanks. >> >> -- >> 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 >> > > > -- > Adam Phillips > Associate Scientist IV, Climate Analysis Section > Climate and Global Dynamics Laboratory > National Center for Atmospheric Research > www.cgd.ucar.edu/staff/asphilli/ > > > -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- f = addfile ("wrfout_d01_2021-03-06_19:00:00","r") rainc = f->RAINC ; (Time, south_north, west_east) rainnc = f->RAINNC u10 = f->U10 ; (Time, south_north, west_east) v10 = f->V10 times = wrf_user_getvar(f,"times",-1) ntim = dimsizes(times) ; # time steps print(times) dims = dimsizes(rainc) print(dims) ; ; Use NCL operator > to make sure all values >=0.0 ; Sum components and assign attributes ; rainc = rainc > 0.0 rainnc = rainnc > 0.0 rainTot = rainc + rainnc rainTot at description = "Total Precipitation" rainTot at units = rainc at units wks = gsn_open_wks("png","WRF_lc") res = True ; plot mods desired res at gsnMaximize = True ; maximize size ;res at gsnScalarContour = True ; contours desired res at gsnLeftString = "Wind Vectors (m/s)" res at gsnRightString = "Total Precipitation (mm)" res at gsnDraw = False ; don't draw plot yet res at gsnFrame = False ; don't advance frame yet res at mpDataBaseVersion = "MediumRes" ; slightly better resolution res at cnFillOn = True ; color plot desired res at cnFillPalette = "BlAqGrYeOrReVi200" ; select color map res at cnLinesOn = False ; turn off contour lines res at cnLineLabelsOn = False ; turn off contour labels res at cnFillMode = "RasterFill" ; raster res at cnLevelSelectionMode = "ExplicitLevels" ; explicit [unequal] cn levels res at cnLevels = (/0,2,4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34,36,38,40,42,44,46,48,50,52,54,56,58,60/) res = wrf_map_resources(f,res) res at gsnAddCyclic = False ; regional data: not cyclic res at tfDoNDCOverlay = True res at mpFillOn = False res at mpGeophysicalLineColor = "black" ; wrf_map_resources uses "gray" res at mpUSStateLineColor = "black" res at mpGeophysicalLineThicknessF = 2.0 ; wrf_map_resources uses 0.5 res at mpUSStateLineThicknessF = 2.0 ;res at mpDataBaseVersion = "MediumRes" ; better map outlines res at mpLimitMode = "LatLon" ; choose region of map res at ZoomIn = True res at mpMaxLatF = 40 ; choose subregion res at mpMinLatF = 30 res at mpMaxLonF = 61 res at mpMinLonF = 42 res = wrf_map_resources(f,res) ; must call this again after setting ZoomIn/X/Ystart/end ; ; Plot one time and level for demo ; . create u and v on a common grid for visualization: nothing fancy ; ;res at tiMainString = times(nt) h=(rainTot(148,:,:) - rainTot(125,:,:)) ;plot = gsn_csm_vector_scalar_map(wks,u10(nt,:,:),v10(nt,:,:),(rainTot(nt,:,:) - rainTot(nt-1,:,:)),res) plot = gsn_csm_contour_map(wks,h,res) shapefile_name1 = "Ostan.shp" shapefile_name2 = "City-Poligon6.shp" lnres = True lnres at gsLineThicknessF = 4.0 lnres at gsLineColor = "Black" shape1 = gsn_add_shapefile_polylines(wks,plot,shapefile_name1,lnres) shape2 = gsn_add_shapefile_polylines(wks,plot,shapefile_name2,lnres) res at gsnDraw = False res at gsnFrame= False draw(plot) frame(wks) -------------- next part -------------- A non-text attachment was scrubbed... Name: no zoom.png Type: image/png Size: 122376 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: zoom.png Type: image/png Size: 107734 bytes Desc: not available URL: From asphilli at ucar.edu Tue Apr 9 10:27:16 2024 From: asphilli at ucar.edu (Adam Phillips) Date: Tue, 9 Apr 2024 10:27:16 -0600 Subject: [ncl-talk] *******Position of shape-file move to the right by narrowing the lat/lon******* In-Reply-To: References: Message-ID: Hi Setareh, Note that in the example I pointed you to that mpMinLatF/mpMaxLonF etc are not the resources that are used. You need to set the special Xstart/Xend/Ystart/Yend resources and pass them to wrf_map_resources: ;---Draw subsetted domain over zoomed in map res at ZoomIn = True ; Tell wrf_map_resources we want to zoom in. res at Xstart = 42 ; Set these four special WRF resources res at Xend = 61 ; required for zooming. res at Ystart = 30 res at Yend = 40 res = wrf_map_resources(a,res) ; must call this again after setting ZoomIn/X/Ystart/end plot = gsn_csm_contour_map(wks,h,res) shapefile_name1 = "Ostan.shp" shapefile_name2 = "City-Poligon6.shp" lnres = True lnres at gsLineThicknessF = 4.0 lnres at gsLineColor = "Black" shape1 = gsn_add_shapefile_polylines(wks,plot,shapefile_name1,lnres) shape2 = gsn_add_shapefile_polylines(wks,plot,shapefile_name2,lnres) draw(plot) frame(wks) Best, Adam On Sat, Apr 6, 2024 at 5:14?AM Setareh Rahimi wrote: > Dear Adam, > > Thank you so much for your suggestion. I added the following lines > (highlighted) to my script but I still have the same problem. > > res at mpLimitMode = "LatLon" ; choose region of map > > *res at ZoomIn = True* > > res at mpMaxLatF = 40 ; choose subregion > > res at mpMinLatF = 30 > > res at mpMaxLonF = 61 > > res at mpMinLonF = 42 > > > r*es = wrf_map_resources(f,res) ** ; must call this again after setting > ZoomIn/X/Ystart/end* > > please check the attached and please advise me where I made the mistakes. > Thanks, > > On Sat, Apr 6, 2024 at 3:25?AM Adam Phillips wrote: > >> Hi Setareh, >> See example #2 here for how to zoom in on a plot when you are >> using wrf_map_resources: >> https://www.ncl.ucar.edu/Applications/wrfzoom.shtml#ex2 >> Best, >> Adam >> >> On Wed, Apr 3, 2024 at 4:26?AM Setareh Rahimi via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >>> Dear NCL users, >>> >>> I tried to plot the total precipitation for 2021-03_12 (24 hours), using >>> the attached script. However, once I tried to zoom in on the study area to >>> better show the result by shortening the lat/lon range I noticed the >>> position of the shapefile on the plot was not correct ( it moved to the >>> right). How can I fix the shape file to the right position? ( please check >>> the attached images). >>> >>> please help me to sole this issue. >>> Thanks. >>> >>> -- >>> 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 >>> >> >> >> -- >> Adam Phillips >> Associate Scientist IV, Climate Analysis Section >> Climate and Global Dynamics Laboratory >> National Center for Atmospheric Research >> www.cgd.ucar.edu/staff/asphilli/ >> >> >> > > > -- > S.Rahimi > > -- 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 setareh.rahimi at gmail.com Sat Apr 20 10:46:29 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Sat, 20 Apr 2024 20:16:29 +0330 Subject: [ncl-talk] ****Urgent help**** Message-ID: Dear NCL users, I need to plot Planet Boundry Height (PBL) like the attached one. But have no idea how to plot it. Would you please help me in this regard? Many thanks. -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Contour-plot-of-Planetary-Boundary-Layer-Height-PBLH-differences-m-between-a-WRF-data.png Type: image/png Size: 326343 bytes Desc: not available URL: From setareh.rahimi at gmail.com Sun Apr 21 01:15:50 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Sun, 21 Apr 2024 10:45:50 +0330 Subject: [ncl-talk] *******Position of shape-file move to the right by narrowing the lat/lon******* In-Reply-To: References: Message-ID: The issue solved, I used this example: https://www.ncl.ucar.edu/Applications/Scripts/dataonmap_zoom_10.ncl I have another question. Once I zoom over my study area since the model's spatial resolution was 12 km, the colors are not smooth and there are squares on the plot. Is there any way to plot it with a smooth texture, please? Check the attached plot. Thanks. Thanks, On Tue, Apr 9, 2024 at 7:57?PM Adam Phillips wrote: > Hi Setareh, > Note that in the example I pointed you to that mpMinLatF/mpMaxLonF etc are > not the resources that are used. You need to set the special > Xstart/Xend/Ystart/Yend resources and pass them to wrf_map_resources: > ;---Draw subsetted domain over zoomed in map > res at ZoomIn = True ; Tell wrf_map_resources we want to zoom in. > res at Xstart = 42 ; Set these four special WRF resources > res at Xend = 61 ; required for zooming. > res at Ystart = 30 > res at Yend = 40 > res = wrf_map_resources(a,res) ; must call this again after setting > ZoomIn/X/Ystart/end > plot = gsn_csm_contour_map(wks,h,res) > > shapefile_name1 = "Ostan.shp" > shapefile_name2 = "City-Poligon6.shp" > lnres = True > lnres at gsLineThicknessF = 4.0 > lnres at gsLineColor = "Black" > shape1 = gsn_add_shapefile_polylines(wks,plot,shapefile_name1,lnres) > shape2 = gsn_add_shapefile_polylines(wks,plot,shapefile_name2,lnres) > draw(plot) > frame(wks) > > Best, > Adam > > On Sat, Apr 6, 2024 at 5:14?AM Setareh Rahimi > wrote: > >> Dear Adam, >> >> Thank you so much for your suggestion. I added the following lines >> (highlighted) to my script but I still have the same problem. >> >> res at mpLimitMode = "LatLon" ; choose region of map >> >> *res at ZoomIn = True* >> >> res at mpMaxLatF = 40 ; choose subregion >> >> res at mpMinLatF = 30 >> >> res at mpMaxLonF = 61 >> >> res at mpMinLonF = 42 >> >> >> r*es = wrf_map_resources(f,res) ** ; must call this again after >> setting ZoomIn/X/Ystart/end* >> >> please check the attached and please advise me where I made the mistakes. >> Thanks, >> >> On Sat, Apr 6, 2024 at 3:25?AM Adam Phillips wrote: >> >>> Hi Setareh, >>> See example #2 here for how to zoom in on a plot when you are >>> using wrf_map_resources: >>> https://www.ncl.ucar.edu/Applications/wrfzoom.shtml#ex2 >>> Best, >>> Adam >>> >>> On Wed, Apr 3, 2024 at 4:26?AM Setareh Rahimi via ncl-talk < >>> ncl-talk at mailman.ucar.edu> wrote: >>> >>>> Dear NCL users, >>>> >>>> I tried to plot the total precipitation for 2021-03_12 (24 hours), >>>> using the attached script. However, once I tried to zoom in on the study >>>> area to better show the result by shortening the lat/lon range I noticed >>>> the position of the shapefile on the plot was not correct ( it moved to the >>>> right). How can I fix the shape file to the right position? ( please check >>>> the attached images). >>>> >>>> please help me to sole this issue. >>>> Thanks. >>>> >>>> -- >>>> 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 >>>> >>> >>> >>> -- >>> Adam Phillips >>> Associate Scientist IV, Climate Analysis Section >>> Climate and Global Dynamics Laboratory >>> National Center for Atmospheric Research >>> www.cgd.ucar.edu/staff/asphilli/ >>> >>> >>> >> >> >> -- >> S.Rahimi >> >> > > -- > Adam Phillips > Associate Scientist IV, Climate Analysis Section > Climate and Global Dynamics Laboratory > National Center for Atmospheric Research > www.cgd.ucar.edu/staff/asphilli/ > > > -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: test_lc.000002.png Type: image/png Size: 79402 bytes Desc: not available URL: From dave.allured at noaa.gov Sun Apr 21 10:41:34 2024 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Sun, 21 Apr 2024 10:41:34 -0600 Subject: [ncl-talk] *******Position of shape-file move to the right by narrowing the lat/lon******* In-Reply-To: References: Message-ID: Setareh, it looks like you have *res at cnFillMode* set to *RasterFill* or *CellFill*. Read about this under the Resources tab on the website. Either remove this resource and let it go to default, or else set it to *AreaFill*, to get visual smoothing. On Sun, Apr 21, 2024 at 1:08?AM Setareh Rahimi via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > The issue solved, I used this example: > > https://www.ncl.ucar.edu/Applications/Scripts/dataonmap_zoom_10.ncl > > > I have another question. > > Once I zoom over my study area since the model's spatial resolution was 12 > km, the colors are not smooth and there are squares on the plot. Is there > any way to plot it with a smooth texture, please? Check the attached plot. > > Thanks. > > Thanks, > -------------- next part -------------- An HTML attachment was scrubbed... URL: From setareh.rahimi at gmail.com Sun Apr 21 11:03:16 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Sun, 21 Apr 2024 20:33:16 +0330 Subject: [ncl-talk] *******Position of shape-file move to the right by narrowing the lat/lon******* In-Reply-To: References: Message-ID: Dear Dave, So many thanks for your advice. I did it and it worked. Best wishes, On Sun, Apr 21, 2024 at 8:12?PM Dave Allured - NOAA Affiliate < dave.allured at noaa.gov> wrote: > Setareh, it looks like you have *res at cnFillMode* set to *RasterFill* or > *CellFill*. Read about this under the Resources tab on the website. > > Either remove this resource and let it go to default, or else set it to > *AreaFill*, to get visual smoothing. > > > On Sun, Apr 21, 2024 at 1:08?AM Setareh Rahimi via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> The issue solved, I used this example: >> >> https://www.ncl.ucar.edu/Applications/Scripts/dataonmap_zoom_10.ncl >> >> >> I have another question. >> >> Once I zoom over my study area since the model's spatial resolution was >> 12 km, the colors are not smooth and there are squares on the plot. Is >> there any way to plot it with a smooth texture, please? Check the attached >> plot. >> >> Thanks. >> >> Thanks, >> > -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: From m.morichetti at pm.univpm.it Mon Apr 22 02:11:46 2024 From: m.morichetti at pm.univpm.it (MORICHETTI MAURO) Date: Mon, 22 Apr 2024 08:11:46 +0000 Subject: [ncl-talk] gsn_panel plot functions Message-ID: Hello everyone, I'm encountering an issue with my script. I'm using two sets of data to create a boxplot and an XY graph for 5 different case studies. The objective is to incorporate these two graphs into a panel plot similar to the example provided (rcp2p6_GPP.png). I understand that the script is not structured optimally, as I combined two separate scripts into one. The current problem arises when I attempt to include the individual plots in the panel using a final loop. It seems that the function is retaining the individual plots, resulting in the absence of the panel plots. I have attached the script and the data. Thank you in advance. Best regards, Mauro -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: rcp2p6_GPP.png Type: image/png Size: 587258 bytes Desc: rcp2p6_GPP.png URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: gsn_panel.zip Type: application/zip Size: 271085 bytes Desc: gsn_panel.zip URL: From m.morichetti at pm.univpm.it Mon Apr 22 02:16:33 2024 From: m.morichetti at pm.univpm.it (MORICHETTI MAURO) Date: Mon, 22 Apr 2024 08:16:33 +0000 Subject: [ncl-talk] gsn_panel plot functions Message-ID: Hello everyone, I'm encountering an issue with my script. I'm using two sets of data to create a boxplot and an XY graph for 5 different case studies. The objective is to incorporate these two graphs into a panel plot similar to the example provided (rcp2p6_GPP.png). I understand that the script is not structured optimally, as I combined two separate scripts into one. The current problem arises when I attempt to include the individual plots in the panel using a final loop. It seems that the function is retaining the individual plots, resulting in the absence of the panel plots. I have attached the script and the data. Thank you in advance. Best regards, Mauro -------------- next part -------------- An HTML attachment was scrubbed... URL: From asphilli at ucar.edu Mon Apr 22 15:34:41 2024 From: asphilli at ucar.edu (Adam Phillips) Date: Mon, 22 Apr 2024 15:34:41 -0600 Subject: [ncl-talk] gsn_panel plot functions In-Reply-To: References: Message-ID: Hi Mauro, Thank you for supplying an easy to run version of your script; that always makes it easier to see what is going on. The typical workflow for workstations/panel plotting is: 1 - Open the workstation. (wks = gsn_open_wks("ps","test") 2 - Create a graphic object plot array to hold the plots. (plot = new(15,graphic)) 3 - Fill in the array created in 2 and associate it with the workstation opened in 1. plot(cntr) = gsn_csm_xy(wks,.....) ) 4 - Call gsn_panel, specifying the workstation opened in 1) and the plotting array created and filled in 2 and 3. (gsn_panel(wks,plot,(/5,3/),pancreas) Your script contains numerous workstations (both x11 and png), you are copying/moving plots around after their creation, and you are calling gsn_panel while specifying a workstation different than that specified for the individual plots. While you can do some of this, it makes it rather confusing to keep track of everything. If you want to stick with the structure of this script, I would not use x11 workstations. Use .ps or .png instead, and only use 1 workstation if possible to keep things cleaner. I think your main issue is how you are handling the workstations. It will help to simplify that by shifting to a single workstation. On a more minor note, set gsnMaximize to False in the individual resource lists. When paneling only set this resource to True in the panel resource list. Hope that at least gets you going in the right direction! Best, Adam On Mon, Apr 22, 2024 at 2:12?AM MORICHETTI MAURO via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > > > Hello everyone, > > I'm encountering an issue with my script. > > I'm using two sets of data to create a boxplot and an XY graph for 5 > different case studies. The objective is to incorporate these two graphs > into a panel plot similar to the example provided (rcp2p6_GPP.png). > > I understand that the script is not structured optimally, as I combined > two separate scripts into one. The current problem arises when I attempt to > include the individual plots in the panel using a final loop. It seems that > the function is retaining the individual plots, resulting in the absence of > the panel plots. > > I have attached the script and the data. > > Thank you in advance. > > Best regards, > > Mauro > > > > > _______________________________________________ > 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 ehsantaghizadeh at yahoo.com Tue Apr 23 07:27:53 2024 From: ehsantaghizadeh at yahoo.com (Ehsan Taghizadeh) Date: Tue, 23 Apr 2024 13:27:53 +0000 (UTC) Subject: [ncl-talk] ****Urgent help**** In-Reply-To: References: Message-ID: <684005297.1454033.1713878874086@mail.yahoo.com> Dear Setareh, May I inquire if the attached figure is accurate? Specifically, are negative values for the Planetary Boundary Layer (PBL) Height reasonable? ----------------------------------------------------------------------------- Sincerely,Ehsan Taghizade On Saturday, April 20, 2024 at 08:09:27 PM GMT+3:30, Setareh Rahimi via ncl-talk wrote: Dear NCL users, I need to plot Planet Boundry?Height (PBL) like the?attached one. But have no idea how to plot it. Would you please help me in this regard?Many thanks.-- 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: Contour-plot-of-Planetary-Boundary-Layer-Height-PBLH-differences-m-between-a-WRF-data.png Type: image/png Size: 326343 bytes Desc: not available URL: From setareh.rahimi at gmail.com Thu Apr 25 02:35:46 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Thu, 25 Apr 2024 12:05:46 +0330 Subject: [ncl-talk] ****Urgent help**** In-Reply-To: <684005297.1454033.1713878874086@mail.yahoo.com> References: <684005297.1454033.1713878874086@mail.yahoo.com> Message-ID: Hi, It was just an example. It is 100% clear that PBL can not be negative! On Tue, Apr 23, 2024 at 4:57?PM Ehsan Taghizadeh wrote: > Dear Setareh, > > May I inquire if the attached figure is accurate? Specifically, are > negative values for the Planetary Boundary Layer (PBL) Height reasonable? > > > *-----------------------------------------------------------------------------* > *Sincerely,* > > *Ehsan Taghizade* > > On Saturday, April 20, 2024 at 08:09:27 PM GMT+3:30, Setareh Rahimi via > ncl-talk wrote: > > > > Dear NCL users, > > I need to plot Planet Boundry Height (PBL) like the attached one. But have > no idea how to plot it. Would you please help me in this regard? > Many thanks. > -- > 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 > -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: From setareh.rahimi at gmail.com Thu Apr 25 02:38:46 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Thu, 25 Apr 2024 12:08:46 +0330 Subject: [ncl-talk] =?utf-8?q?Indicating_=28=CE=BC=29_sign_in_the_title_s?= =?utf-8?q?ection_of_a_plot?= Message-ID: Hi all, How can I indicate the "micro" (?) sign in the title section of a plot, please? Thanks, -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: From meier-fleischer at dkrz.de Thu Apr 25 02:38:19 2024 From: meier-fleischer at dkrz.de (Karin Meier-Fleischer) Date: Thu, 25 Apr 2024 10:38:19 +0200 Subject: [ncl-talk] =?utf-8?q?Indicating_=28=CE=BC=29_sign_in_the_title_s?= =?utf-8?q?ection_of_a_plot?= In-Reply-To: References: Message-ID: <481ebd39-50a0-43b4-bdb5-c8fdc1aebbc5@dkrz.de> Hi Setareh, use "~F33~m~F21~" More about fonts and how to use it have a look at https://www.ncl.ucar.edu/Document/Graphics/font_tables.shtml Regards, Karin Am 25.04.24 um 10:38 schrieb Setareh Rahimi via ncl-talk: > Hi all, > > How can I indicate?the "micro" (?)?sign in the title section of a > plot, please? > > Thanks, > > -- > 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: From setareh.rahimi at gmail.com Thu Apr 25 02:50:42 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Thu, 25 Apr 2024 12:20:42 +0330 Subject: [ncl-talk] =?utf-8?q?Indicating_=28=CE=BC=29_sign_in_the_title_s?= =?utf-8?q?ection_of_a_plot?= In-Reply-To: <481ebd39-50a0-43b4-bdb5-c8fdc1aebbc5@dkrz.de> References: <481ebd39-50a0-43b4-bdb5-c8fdc1aebbc5@dkrz.de> Message-ID: Dear Karin, So many thanks for your help. Best wishes, On Thu, Apr 25, 2024 at 12:08?PM Karin Meier-Fleischer via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Hi Setareh, > > use "~F33~m~F21~" > > More about fonts and how to use it have a look at > https://www.ncl.ucar.edu/Document/Graphics/font_tables.shtml > > Regards, > Karin > > > Am 25.04.24 um 10:38 schrieb Setareh Rahimi via ncl-talk: > > Hi all, > > How can I indicate the "micro" (?) sign in the title section of a plot, > please? > > Thanks, > > -- > S.Rahimi > > > _______________________________________________ > ncl-talk mailing listncl-talk at mailman.ucar.edu > List instructions, subscriber options, unsubscribe:https://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > _______________________________________________ > ncl-talk mailing list > ncl-talk at mailman.ucar.edu > List instructions, subscriber options, unsubscribe: > https://mailman.ucar.edu/mailman/listinfo/ncl-talk > -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: From mughalali655 at gmail.com Mon Apr 29 00:49:36 2024 From: mughalali655 at gmail.com (ali mughal) Date: Mon, 29 Apr 2024 14:49:36 +0800 Subject: [ncl-talk] stat_desp code Message-ID: Hi My following little code is not providing value of the percentiles probably because of missing values dir_agcd = "/data2/WRF_ERA5_out/SWWA/AGCD_regrid_mon/means/" dir_wrf = "/data2/WRF_ERA5_out/SWWA/wrf_seasmean/" ens = "R3" yr_start = "1980" yr_end = "2019" nlon = 523 nlat = 346 opt = True system("mkdir -p " + diro) f_agcd_tmax = addfile(dir_agcd + "agcd_v1-0-1_" + "tmax_" + yr_start + "_" + yr_end + "_seasmean.nc", "r") agcd_tmax = f_agcd_tmax->tmax delete(f_agcd_tmax) f_wrf_tasmax = addfile(dir_wrf + "tasmax_" + ens + "_" + yr_start + "_" + yr_end + "_seasmean.nc", "r") wrf_tasmax = (f_wrf_tasmax->tasmax) wrf_tasmax = wrf_tasmax-273 delete(f_wrf_tasmax) do i = 0, nlat - 1, 1 do j = 0, nlon - 1, 1 agcd_tmax_out = agcd_tmax(:, i, j) if all(ismissing(agcd_tmax_out)) .eq. False then wrf_tasmax_out = wrf_tasmax(:, i, j) stat_agcd_tmax = stat_dispersion(agcd_tmax_out, opt) stat_wrf_tasmax = stat_dispersion(wrf_tasmax_out, opt) else delete(agcd_tmax_out) end if end do end do print(stat_agcd_tmax) print(stat_wrf_tasmax) Output Variable: stat_agcd_tmax Type: float Total Size: 120 bytes 30 values Number of Dimensions: 1 Dimensions and sizes: [30] Coordinates: Number Of Attributes: 2 long_name : Robust Dispersion Statistics _FillValue : 1e+20 (0) 31.07983 (1) 0.768555 (2) 29.00586 (3) 29.96567 (4) 30.07379 (5) 30.28866 (6) 30.6044 (7) 30.78744 (8) 31.18414 (9) 31.38833 (10) 31.58258 (11) 31.85179 (12) 31.96011 (13) 32.05429 (14) 32.88403 (15) 3.878166 (16) 5.046049 (17) 0.7660876 (18) 156 (19) 156 (20) 0 (21) 0 (22) 1e+20 (23) 1e+20 (24) 1e+20 (25) 1e+20 (26) 1e+20 (27) 1e+20 (28) -0.2533386 (29) -0.2925659 Variable: stat_wrf_tasmax Type: float Total Size: 120 bytes 30 values Number of Dimensions: 1 Dimensions and sizes: [30] Coordinates: Number Of Attributes: 3 units : K long_name : Robust Dispersion Statistics: Daily Maximum Near-Surface Air Temperature _FillValue : 1e+20 (0) 25.65089 (1) 0.8997155 (2) 23.26544 (3) 24.28622 (4) 24.36398 (5) 24.49762 (6) 25.16293 (7) 25.38058 (8) 25.80746 (9) 26.04553 (10) 26.27335 (11) 26.43805 (12) 26.66354 (13) 26.73364 (14) 28.01169 (15) 4.746246 (16) 5.275275 (17) 0.8968271 (18) 156 (19) 156 (20) 0 (21) 0 (22) 1e+20 (23) 1e+20 (24) 1e+20 (25) 1e+20 (26) 1e+20 (27) 1e+20 (28) -0.4017957 (29) -0.3656675 -------------- next part -------------- An HTML attachment was scrubbed... URL: From mughalali655 at gmail.com Mon Apr 29 00:49:36 2024 From: mughalali655 at gmail.com (ali mughal) Date: Mon, 29 Apr 2024 14:49:36 +0800 Subject: [ncl-talk] stat_desp code Message-ID: Hi My following little code is not providing value of the percentiles probably because of missing values dir_agcd = "/data2/WRF_ERA5_out/SWWA/AGCD_regrid_mon/means/" dir_wrf = "/data2/WRF_ERA5_out/SWWA/wrf_seasmean/" ens = "R3" yr_start = "1980" yr_end = "2019" nlon = 523 nlat = 346 opt = True system("mkdir -p " + diro) f_agcd_tmax = addfile(dir_agcd + "agcd_v1-0-1_" + "tmax_" + yr_start + "_" + yr_end + "_seasmean.nc", "r") agcd_tmax = f_agcd_tmax->tmax delete(f_agcd_tmax) f_wrf_tasmax = addfile(dir_wrf + "tasmax_" + ens + "_" + yr_start + "_" + yr_end + "_seasmean.nc", "r") wrf_tasmax = (f_wrf_tasmax->tasmax) wrf_tasmax = wrf_tasmax-273 delete(f_wrf_tasmax) do i = 0, nlat - 1, 1 do j = 0, nlon - 1, 1 agcd_tmax_out = agcd_tmax(:, i, j) if all(ismissing(agcd_tmax_out)) .eq. False then wrf_tasmax_out = wrf_tasmax(:, i, j) stat_agcd_tmax = stat_dispersion(agcd_tmax_out, opt) stat_wrf_tasmax = stat_dispersion(wrf_tasmax_out, opt) else delete(agcd_tmax_out) end if end do end do print(stat_agcd_tmax) print(stat_wrf_tasmax) Output Variable: stat_agcd_tmax Type: float Total Size: 120 bytes 30 values Number of Dimensions: 1 Dimensions and sizes: [30] Coordinates: Number Of Attributes: 2 long_name : Robust Dispersion Statistics _FillValue : 1e+20 (0) 31.07983 (1) 0.768555 (2) 29.00586 (3) 29.96567 (4) 30.07379 (5) 30.28866 (6) 30.6044 (7) 30.78744 (8) 31.18414 (9) 31.38833 (10) 31.58258 (11) 31.85179 (12) 31.96011 (13) 32.05429 (14) 32.88403 (15) 3.878166 (16) 5.046049 (17) 0.7660876 (18) 156 (19) 156 (20) 0 (21) 0 (22) 1e+20 (23) 1e+20 (24) 1e+20 (25) 1e+20 (26) 1e+20 (27) 1e+20 (28) -0.2533386 (29) -0.2925659 Variable: stat_wrf_tasmax Type: float Total Size: 120 bytes 30 values Number of Dimensions: 1 Dimensions and sizes: [30] Coordinates: Number Of Attributes: 3 units : K long_name : Robust Dispersion Statistics: Daily Maximum Near-Surface Air Temperature _FillValue : 1e+20 (0) 25.65089 (1) 0.8997155 (2) 23.26544 (3) 24.28622 (4) 24.36398 (5) 24.49762 (6) 25.16293 (7) 25.38058 (8) 25.80746 (9) 26.04553 (10) 26.27335 (11) 26.43805 (12) 26.66354 (13) 26.73364 (14) 28.01169 (15) 4.746246 (16) 5.275275 (17) 0.8968271 (18) 156 (19) 156 (20) 0 (21) 0 (22) 1e+20 (23) 1e+20 (24) 1e+20 (25) 1e+20 (26) 1e+20 (27) 1e+20 (28) -0.4017957 (29) -0.3656675 -------------- next part -------------- An HTML attachment was scrubbed... URL: From setareh.rahimi at gmail.com Mon Apr 29 07:02:49 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Mon, 29 Apr 2024 16:32:49 +0330 Subject: [ncl-talk] Liquid Water Path (LWP) Message-ID: Dear all, I would like to plot *liquid water path (LWP) *from WRF-Chem model using *"QCLOUD" *variable. So is there any example available in NCL? or any suggestion on how to calculate and plot it? Kind regards. -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: