From ehsantaghizadeh at yahoo.com Wed Jan 3 08:46:18 2024 From: ehsantaghizadeh at yahoo.com (Ehsan Taghizadeh) Date: Wed, 3 Jan 2024 15:46:18 +0000 (UTC) Subject: [ncl-talk] WRF output plot over surface References: <1862690046.4428736.1704296778773.ref@mail.yahoo.com> Message-ID: <1862690046.4428736.1704296778773@mail.yahoo.com> Dear NCL group,I'm not sure if this question belongs to the WRF forum or the NCL forum, but I would appreciate any help on this. I want to plot a variable (TOT_DUST in my case) on the surface level from a WRF-Chem output. I used an NCL script and found that the wrf_user_interp_level function in this language can interpolate a horizontal slice from a three-dimensional WRF-ARW field at the given vertical level(s). I'm not sure how to use this function (or maybe another function) to plot a 4D variable on the surface level.I look forward to hearing from you. -----------------------------------------------------------------------------Sincerely,Ehsan Taghizade -------------- next part -------------- An HTML attachment was scrubbed... URL: From dmjalal90 at gmail.com Fri Jan 5 06:15:31 2024 From: dmjalal90 at gmail.com (Md. Jalal Uddin) Date: Fri, 5 Jan 2024 22:15:31 +0900 Subject: [ncl-talk] Do not show significance level in a correlation map Message-ID: Dear all, I have tried to use the gsn_contour_shade function to show the significance level as a shaded pattern where p value is less than 0.05 and 0.01. However, the output map does not show the shaded areas (please see the attached map). Here are some codes: r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; ccr(lat,lon) copy_VarCoords(ano_sstSON(0,:,:), r_sst) r_sst at long_name = "Correlation: r_sst" p = rtest(r_sst,42,0) copy_VarCoords(r_sst, p) opt = True opt at gsnShadeHigh = 17 ; stipple opt at gsnShadeFillScaleF = 0.5 ; add extra density opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot plot2 = gsn_csm_contour(wks,p, res2) plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) overlay (plot, plot2) draw(plot) frame(wks) Please find the full code in the attached file. Any suggestions would be appreciated. -- *Regards,* Dr. Jalal *Postdoctoral Researcher,* Typhoon Research Center, Jeju National University, South Korea *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China *M.Sc.* in Applied Meteorology, NUIST, China *B.Sc.* in Disaster Management, PSTU, Bangladesh *Founder and Director* of Research Society *Website*: https://researchsociety20.org/founder-and-director/ *E-mails*: founder-and-director at researchsociety20.org jalal at jejunu.ac.kr 20205103002 at nuist.edu.cn -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: sst_SepNov_cor.png Type: image/png Size: 572694 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: sst_corr.ncl Type: application/octet-stream Size: 4484 bytes Desc: not available URL: From asphilli at ucar.edu Fri Jan 12 14:55:10 2024 From: asphilli at ucar.edu (Adam Phillips) Date: Fri, 12 Jan 2024 14:55:10 -0700 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Hi Jalal, As you are pattern filling, and the default is color filling, you need to set: opt at gsnShadeFillType = "pattern" I am also unclear what you mean by this statement: "...to show the significance level as a shaded pattern where p value is less than 0.05 and 0.01" Do you want to stipple areas less than .05? If so the following code should do that: opt = True opt at gsnShade*Low* = 17 ; stipple opt at gsnShadeFillScaleF = 0.5 ; add extra density plot2 = gsn_contour_shade(plot2,0.05,999.,opt) If instead you want to stipple areas less than .01, you will need to modify the contour levels in your res2 resource list to have a contour at .01: res2 at cnMinLevelValF = 0.00 ; set min contour level res2 at cnMaxLevelValF = .10 ; set max contour level res2 at cnLevelSpacingF = 0.01 ; set contour spacing opt = True opt at gsnShade*Low* = 17 ; stipple opt at gsnShadeFillScaleF = 0.5 ; add extra density plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ Hope that answers your question. As always please respond to the ncl-talk email list with further questions. Best, Adam On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear all, > > I have tried to use the gsn_contour_shade function to show the > significance level as a shaded pattern where p value is less than 0.05 and > 0.01. However, the output map does not show the shaded areas (please see > the attached map). > > Here are some codes: > r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; > ccr(lat,lon) > copy_VarCoords(ano_sstSON(0,:,:), r_sst) > r_sst at long_name = "Correlation: r_sst" > > p = rtest(r_sst,42,0) > copy_VarCoords(r_sst, p) > > opt = True > opt at gsnShadeHigh = 17 ; stipple > opt at gsnShadeFillScaleF = 0.5 ; add extra density > opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger > > plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot > > plot2 = gsn_csm_contour(wks,p, res2) > plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) > > overlay (plot, plot2) > > draw(plot) > frame(wks) > > Please find the full code in the attached file. > > Any suggestions would be appreciated. > > -- > > *Regards,* > > Dr. Jalal > > *Postdoctoral Researcher,* Typhoon Research Center, Jeju National > University, South Korea > > *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China > > *M.Sc.* in Applied Meteorology, NUIST, China > > *B.Sc.* in Disaster Management, PSTU, Bangladesh > > *Founder and Director* of Research Society > > *Website*: https://researchsociety20.org/founder-and-director/ > > *E-mails*: founder-and-director at researchsociety20.org > > jalal at jejunu.ac.kr > > 20205103002 at nuist.edu.cn > _______________________________________________ > 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 dmjalal90 at gmail.com Mon Jan 15 20:58:20 2024 From: dmjalal90 at gmail.com (Md. Jalal Uddin) Date: Tue, 16 Jan 2024 12:58:20 +0900 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Thank you very much, sir. I used the following resources: opt = True opt at gsnShadeFillType = "pattern" opt at gsnShadeLow = 17 ; stipple opt at gsnShadeFillScaleF = 0.5 ; add extra density plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot plot2 = gsn_csm_contour(wks,p, res2) plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ overlay (plot, plot2) draw(plot) frame(wks) However, I did not get output as stipple areas where p value is less than 0.05. I added the full code including Nino data in the attachment. Should I send the SST data (sst.mon.mean.nc) via ftp? On Sat, 13 Jan 2024 at 06:55, Adam Phillips wrote: > Hi Jalal, > As you are pattern filling, and the default is color filling, you need to > set: > opt at gsnShadeFillType = "pattern" > I am also unclear what you mean by this statement: > "...to show the significance level as a shaded pattern where p value is > less than 0.05 and 0.01" > > Do you want to stipple areas less than .05? If so the following code > should do that: > opt = True > opt at gsnShade*Low* = 17 ; stipple > opt at gsnShadeFillScaleF = 0.5 ; add extra density > plot2 = gsn_contour_shade(plot2,0.05,999.,opt) > > If instead you want to stipple areas less than .01, you will need to > modify the contour levels in your res2 resource list to have a contour at > .01: > res2 at cnMinLevelValF = 0.00 ; set min contour level > res2 at cnMaxLevelValF = .10 ; set max contour level > res2 at cnLevelSpacingF = 0.01 ; set contour spacing > > opt = True > opt at gsnShade*Low* = 17 ; stipple > opt at gsnShadeFillScaleF = 0.5 ; add extra density > plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ > > Hope that answers your question. As always please respond to the ncl-talk > email list with further questions. > Best, > Adam > > On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> Dear all, >> >> I have tried to use the gsn_contour_shade function to show the >> significance level as a shaded pattern where p value is less than 0.05 and >> 0.01. However, the output map does not show the shaded areas (please see >> the attached map). >> >> Here are some codes: >> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; >> ccr(lat,lon) >> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >> r_sst at long_name = "Correlation: r_sst" >> >> p = rtest(r_sst,42,0) >> copy_VarCoords(r_sst, p) >> >> opt = True >> opt at gsnShadeHigh = 17 ; stipple >> opt at gsnShadeFillScaleF = 0.5 ; add extra density >> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >> >> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >> >> plot2 = gsn_csm_contour(wks,p, res2) >> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >> >> overlay (plot, plot2) >> >> draw(plot) >> frame(wks) >> >> Please find the full code in the attached file. >> >> Any suggestions would be appreciated. >> >> -- >> >> *Regards,* >> >> Dr. Jalal >> >> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >> University, South Korea >> >> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >> >> *M.Sc.* in Applied Meteorology, NUIST, China >> >> *B.Sc.* in Disaster Management, PSTU, Bangladesh >> >> *Founder and Director* of Research Society >> >> *Website*: https://researchsociety20.org/founder-and-director/ >> >> *E-mails*: founder-and-director at researchsociety20.org >> >> jalal at jejunu.ac.kr >> >> 20205103002 at nuist.edu.cn >> _______________________________________________ >> 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/ > > > -- *Regards,* Dr. Jalal *Postdoctoral Researcher,* Typhoon Research Center, Jeju National University, South Korea *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China *M.Sc.* in Applied Meteorology, NUIST, China *B.Sc.* in Disaster Management, PSTU, Bangladesh *Founder and Director* of Research Society *Website*: https://researchsociety20.org/founder-and-director/ *E-mails*: founder-and-director at researchsociety20.org jalal at jejunu.ac.kr 20205103002 at nuist.edu.cn -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: sst_corr.ncl Type: application/octet-stream Size: 4471 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: AMM_DMI_NINO_PDO_ELI_SepNovAvg_1980_2021.csv Type: text/csv Size: 2586 bytes Desc: not available URL: From asphilli at ucar.edu Tue Jan 16 13:08:36 2024 From: asphilli at ucar.edu (Adam Phillips) Date: Tue, 16 Jan 2024 13:08:36 -0700 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Hi Jalal, I was able to plug in a monthly SST dataset that I have here, so I was able to run your script. So I think what is happening is that the stippling is so small/faint/sparse that it is hard to see it. I made two changes to your script: 1 - I was getting stippling over land as the rtest function was seemingly not recognizing the _FillValue of the input array. (I find this strange, but have not used rtest before.) To fix this, I put in a where statement and added the _FillValue attribute to your p array: p = rtest(r_sst,42,0) * p at _FillValue = 1.e20* copy_VarCoords(r_sst, p) * p = where(ismissing(r_sst),p at _FillValue,p)* 2 - I modified/added the following gsn_contour_shade options to make the stippling more visible: * opt at gsnShadeFillScaleF = 0.25 ; add extra density opt at gsnShadeFillDotSizeF = 0.0015* I have attached the modified script, along with the plot that I am getting. Hope that helps! Best, Adam On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin wrote: > Thank you very much, sir. I used the following resources: > > opt = True > opt at gsnShadeFillType = "pattern" > opt at gsnShadeLow = 17 ; stipple > opt at gsnShadeFillScaleF = 0.5 ; add extra density > > plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot > plot2 = gsn_csm_contour(wks,p, res2) > plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ > > overlay (plot, plot2) > > draw(plot) > frame(wks) > > However, I did not get output as stipple areas where p value is less than > 0.05. I added the full code including Nino data in the attachment. Should I > send the SST data (sst.mon.mean.nc) via ftp? > > On Sat, 13 Jan 2024 at 06:55, Adam Phillips wrote: > >> Hi Jalal, >> As you are pattern filling, and the default is color filling, you need to >> set: >> opt at gsnShadeFillType = "pattern" >> I am also unclear what you mean by this statement: >> "...to show the significance level as a shaded pattern where p value is >> less than 0.05 and 0.01" >> >> Do you want to stipple areas less than .05? If so the following code >> should do that: >> opt = True >> opt at gsnShade*Low* = 17 ; stipple >> opt at gsnShadeFillScaleF = 0.5 ; add extra density >> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >> >> If instead you want to stipple areas less than .01, you will need to >> modify the contour levels in your res2 resource list to have a contour at >> .01: >> res2 at cnMinLevelValF = 0.00 ; set min contour level >> res2 at cnMaxLevelValF = .10 ; set max contour level >> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >> >> opt = True >> opt at gsnShade*Low* = 17 ; stipple >> opt at gsnShadeFillScaleF = 0.5 ; add extra density >> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >> >> Hope that answers your question. As always please respond to the ncl-talk >> email list with further questions. >> Best, >> Adam >> >> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >>> Dear all, >>> >>> I have tried to use the gsn_contour_shade function to show the >>> significance level as a shaded pattern where p value is less than 0.05 and >>> 0.01. However, the output map does not show the shaded areas (please see >>> the attached map). >>> >>> Here are some codes: >>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; >>> ccr(lat,lon) >>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>> r_sst at long_name = "Correlation: r_sst" >>> >>> p = rtest(r_sst,42,0) >>> copy_VarCoords(r_sst, p) >>> >>> opt = True >>> opt at gsnShadeHigh = 17 ; stipple >>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>> >>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>> >>> plot2 = gsn_csm_contour(wks,p, res2) >>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>> >>> overlay (plot, plot2) >>> >>> draw(plot) >>> frame(wks) >>> >>> Please find the full code in the attached file. >>> >>> Any suggestions would be appreciated. >>> >>> -- >>> >>> *Regards,* >>> >>> Dr. Jalal >>> >>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>> University, South Korea >>> >>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>> >>> *M.Sc.* in Applied Meteorology, NUIST, China >>> >>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>> >>> *Founder and Director* of Research Society >>> >>> *Website*: https://researchsociety20.org/founder-and-director/ >>> >>> *E-mails*: founder-and-director at researchsociety20.org >>> >>> jalal at jejunu.ac.kr >>> >>> 20205103002 at nuist.edu.cn >>> _______________________________________________ >>> 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/ >> >> >> > > > -- > > *Regards,* > > Dr. Jalal > > *Postdoctoral Researcher,* Typhoon Research Center, Jeju National > University, South Korea > > *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China > > *M.Sc.* in Applied Meteorology, NUIST, China > > *B.Sc.* in Disaster Management, PSTU, Bangladesh > > *Founder and Director* of Research Society > > *Website*: https://researchsociety20.org/founder-and-director/ > > *E-mails*: founder-and-director at researchsociety20.org > > jalal at jejunu.ac.kr > > 20205103002 at nuist.edu.cn > -- 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: -------------- next part -------------- A non-text attachment was scrubbed... Name: sst_corr.ncl Type: application/octet-stream Size: 4613 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: sst_SepNov_cor.png Type: image/png Size: 545972 bytes Desc: not available URL: From dmjalal90 at gmail.com Thu Jan 18 23:32:04 2024 From: dmjalal90 at gmail.com (Md. Jalal Uddin) Date: Fri, 19 Jan 2024 15:32:04 +0900 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: This is great, thank you so much! On Wed, 17 Jan 2024 at 05:09, Adam Phillips wrote: > Hi Jalal, > I was able to plug in a monthly SST dataset that I have here, so I was > able to run your script. So I think what is happening is that the stippling > is so small/faint/sparse that it is hard to see it. I made two changes to > your script: > 1 - I was getting stippling over land as the rtest function was seemingly > not recognizing the _FillValue of the input array. (I find this strange, > but have not used rtest before.) To fix this, I put in a where statement > and added the _FillValue attribute to your p array: > p = rtest(r_sst,42,0) > * p at _FillValue = 1.e20* > copy_VarCoords(r_sst, p) > * p = where(ismissing(r_sst),p at _FillValue,p)* > > 2 - I modified/added the following gsn_contour_shade options to make the > stippling more visible: > > * opt at gsnShadeFillScaleF = 0.25 ; add extra density > opt at gsnShadeFillDotSizeF = 0.0015* > > I have attached the modified script, along with the plot that I am getting. > Hope that helps! > Best, > Adam > > > On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin > wrote: > >> Thank you very much, sir. I used the following resources: >> >> opt = True >> opt at gsnShadeFillType = "pattern" >> opt at gsnShadeLow = 17 ; stipple >> opt at gsnShadeFillScaleF = 0.5 ; add extra density >> >> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >> plot2 = gsn_csm_contour(wks,p, res2) >> plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ >> >> overlay (plot, plot2) >> >> draw(plot) >> frame(wks) >> >> However, I did not get output as stipple areas where p value is less than >> 0.05. I added the full code including Nino data in the attachment. Should I >> send the SST data (sst.mon.mean.nc) via ftp? >> >> On Sat, 13 Jan 2024 at 06:55, Adam Phillips wrote: >> >>> Hi Jalal, >>> As you are pattern filling, and the default is color filling, you need >>> to set: >>> opt at gsnShadeFillType = "pattern" >>> I am also unclear what you mean by this statement: >>> "...to show the significance level as a shaded pattern where p value is >>> less than 0.05 and 0.01" >>> >>> Do you want to stipple areas less than .05? If so the following code >>> should do that: >>> opt = True >>> opt at gsnShade*Low* = 17 ; stipple >>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >>> >>> If instead you want to stipple areas less than .01, you will need to >>> modify the contour levels in your res2 resource list to have a contour at >>> .01: >>> res2 at cnMinLevelValF = 0.00 ; set min contour level >>> res2 at cnMaxLevelValF = .10 ; set max contour level >>> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >>> >>> opt = True >>> opt at gsnShade*Low* = 17 ; stipple >>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >>> >>> Hope that answers your question. As always please respond to the >>> ncl-talk email list with further questions. >>> Best, >>> Adam >>> >>> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >>> ncl-talk at mailman.ucar.edu> wrote: >>> >>>> Dear all, >>>> >>>> I have tried to use the gsn_contour_shade function to show the >>>> significance level as a shaded pattern where p value is less than 0.05 and >>>> 0.01. However, the output map does not show the shaded areas (please see >>>> the attached map). >>>> >>>> Here are some codes: >>>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; >>>> ccr(lat,lon) >>>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>>> r_sst at long_name = "Correlation: r_sst" >>>> >>>> p = rtest(r_sst,42,0) >>>> copy_VarCoords(r_sst, p) >>>> >>>> opt = True >>>> opt at gsnShadeHigh = 17 ; stipple >>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>>> >>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>> >>>> plot2 = gsn_csm_contour(wks,p, res2) >>>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>>> >>>> overlay (plot, plot2) >>>> >>>> draw(plot) >>>> frame(wks) >>>> >>>> Please find the full code in the attached file. >>>> >>>> Any suggestions would be appreciated. >>>> >>>> -- >>>> >>>> *Regards,* >>>> >>>> Dr. Jalal >>>> >>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>> University, South Korea >>>> >>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>> >>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>> >>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>> >>>> *Founder and Director* of Research Society >>>> >>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>> >>>> *E-mails*: founder-and-director at researchsociety20.org >>>> >>>> jalal at jejunu.ac.kr >>>> >>>> 20205103002 at nuist.edu.cn >>>> _______________________________________________ >>>> 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/ >>> >>> >>> >> >> >> -- >> >> *Regards,* >> >> Dr. Jalal >> >> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >> University, South Korea >> >> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >> >> *M.Sc.* in Applied Meteorology, NUIST, China >> >> *B.Sc.* in Disaster Management, PSTU, Bangladesh >> >> *Founder and Director* of Research Society >> >> *Website*: https://researchsociety20.org/founder-and-director/ >> >> *E-mails*: founder-and-director at researchsociety20.org >> >> jalal at jejunu.ac.kr >> >> 20205103002 at nuist.edu.cn >> > > > -- > Adam Phillips > Associate Scientist IV, Climate Analysis Section > Climate and Global Dynamics Laboratory > National Center for Atmospheric Research > www.cgd.ucar.edu/staff/asphilli/ > > > -- *Regards,* Dr. Jalal *Postdoctoral Researcher,* Typhoon Research Center, Jeju National University, South Korea *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China *M.Sc.* in Applied Meteorology, NUIST, China *B.Sc.* in Disaster Management, PSTU, Bangladesh *Founder and Director* of Research Society *Website*: https://researchsociety20.org/founder-and-director/ *E-mails*: founder-and-director at researchsociety20.org jalal at jejunu.ac.kr 20205103002 at nuist.edu.cn -------------- next part -------------- An HTML attachment was scrubbed... URL: From dmjalal90 at gmail.com Fri Jan 19 03:42:31 2024 From: dmjalal90 at gmail.com (Md. Jalal Uddin) Date: Fri, 19 Jan 2024 19:42:31 +0900 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Thank you, sir, for your kind efforts. Your code works for hadisst.187001-202212.nc data. However, when I use my sst data from NOAA ( sst.mnmean.nc), the code does not work. I do not know the reason for this. Sorry to bother you. NOAA sst data: https://psl.noaa.gov/repository/entry/show?entryid=cac1c2a6-a864-4409-bb77-1fdead8eeb6e On Fri, 19 Jan 2024 at 15:32, Md. Jalal Uddin wrote: > This is great, thank you so much! > > On Wed, 17 Jan 2024 at 05:09, Adam Phillips wrote: > >> Hi Jalal, >> I was able to plug in a monthly SST dataset that I have here, so I was >> able to run your script. So I think what is happening is that the stippling >> is so small/faint/sparse that it is hard to see it. I made two changes to >> your script: >> 1 - I was getting stippling over land as the rtest function was seemingly >> not recognizing the _FillValue of the input array. (I find this strange, >> but have not used rtest before.) To fix this, I put in a where statement >> and added the _FillValue attribute to your p array: >> p = rtest(r_sst,42,0) >> * p at _FillValue = 1.e20* >> copy_VarCoords(r_sst, p) >> * p = where(ismissing(r_sst),p at _FillValue,p)* >> >> 2 - I modified/added the following gsn_contour_shade options to make the >> stippling more visible: >> >> * opt at gsnShadeFillScaleF = 0.25 ; add extra density >> opt at gsnShadeFillDotSizeF = 0.0015* >> >> I have attached the modified script, along with the plot that I am >> getting. >> Hope that helps! >> Best, >> Adam >> >> >> On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin >> wrote: >> >>> Thank you very much, sir. I used the following resources: >>> >>> opt = True >>> opt at gsnShadeFillType = "pattern" >>> opt at gsnShadeLow = 17 ; stipple >>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>> >>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>> plot2 = gsn_csm_contour(wks,p, res2) >>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ >>> >>> overlay (plot, plot2) >>> >>> draw(plot) >>> frame(wks) >>> >>> However, I did not get output as stipple areas where p value is less >>> than 0.05. I added the full code including Nino data in the attachment. >>> Should I send the SST data (sst.mon.mean.nc) via ftp? >>> >>> On Sat, 13 Jan 2024 at 06:55, Adam Phillips wrote: >>> >>>> Hi Jalal, >>>> As you are pattern filling, and the default is color filling, you need >>>> to set: >>>> opt at gsnShadeFillType = "pattern" >>>> I am also unclear what you mean by this statement: >>>> "...to show the significance level as a shaded pattern where p value is >>>> less than 0.05 and 0.01" >>>> >>>> Do you want to stipple areas less than .05? If so the following code >>>> should do that: >>>> opt = True >>>> opt at gsnShade*Low* = 17 ; stipple >>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >>>> >>>> If instead you want to stipple areas less than .01, you will need to >>>> modify the contour levels in your res2 resource list to have a contour at >>>> .01: >>>> res2 at cnMinLevelValF = 0.00 ; set min contour level >>>> res2 at cnMaxLevelValF = .10 ; set max contour level >>>> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >>>> >>>> opt = True >>>> opt at gsnShade*Low* = 17 ; stipple >>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >>>> >>>> Hope that answers your question. As always please respond to the >>>> ncl-talk email list with further questions. >>>> Best, >>>> Adam >>>> >>>> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >>>> ncl-talk at mailman.ucar.edu> wrote: >>>> >>>>> Dear all, >>>>> >>>>> I have tried to use the gsn_contour_shade function to show the >>>>> significance level as a shaded pattern where p value is less than 0.05 and >>>>> 0.01. However, the output map does not show the shaded areas (please see >>>>> the attached map). >>>>> >>>>> Here are some codes: >>>>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; >>>>> ccr(lat,lon) >>>>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>>>> r_sst at long_name = "Correlation: r_sst" >>>>> >>>>> p = rtest(r_sst,42,0) >>>>> copy_VarCoords(r_sst, p) >>>>> >>>>> opt = True >>>>> opt at gsnShadeHigh = 17 ; stipple >>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>>>> >>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>> >>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>>>> >>>>> overlay (plot, plot2) >>>>> >>>>> draw(plot) >>>>> frame(wks) >>>>> >>>>> Please find the full code in the attached file. >>>>> >>>>> Any suggestions would be appreciated. >>>>> >>>>> -- >>>>> >>>>> *Regards,* >>>>> >>>>> Dr. Jalal >>>>> >>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>> University, South Korea >>>>> >>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>>> >>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>> >>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>> >>>>> *Founder and Director* of Research Society >>>>> >>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>> >>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>> >>>>> jalal at jejunu.ac.kr >>>>> >>>>> 20205103002 at nuist.edu.cn >>>>> _______________________________________________ >>>>> 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/ >>>> >>>> >>>> >>> >>> >>> -- >>> >>> *Regards,* >>> >>> Dr. Jalal >>> >>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>> University, South Korea >>> >>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>> >>> *M.Sc.* in Applied Meteorology, NUIST, China >>> >>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>> >>> *Founder and Director* of Research Society >>> >>> *Website*: https://researchsociety20.org/founder-and-director/ >>> >>> *E-mails*: founder-and-director at researchsociety20.org >>> >>> jalal at jejunu.ac.kr >>> >>> 20205103002 at nuist.edu.cn >>> >> >> >> -- >> Adam Phillips >> Associate Scientist IV, Climate Analysis Section >> Climate and Global Dynamics Laboratory >> National Center for Atmospheric Research >> www.cgd.ucar.edu/staff/asphilli/ >> >> >> > > > -- > > *Regards,* > > Dr. Jalal > > *Postdoctoral Researcher,* Typhoon Research Center, Jeju National > University, South Korea > > *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China > > *M.Sc.* in Applied Meteorology, NUIST, China > > *B.Sc.* in Disaster Management, PSTU, Bangladesh > > *Founder and Director* of Research Society > > *Website*: https://researchsociety20.org/founder-and-director/ > > *E-mails*: founder-and-director at researchsociety20.org > > jalal at jejunu.ac.kr > > 20205103002 at nuist.edu.cn > -- *Regards,* Dr. Jalal *Postdoctoral Researcher,* Typhoon Research Center, Jeju National University, South Korea *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China *M.Sc.* in Applied Meteorology, NUIST, China *B.Sc.* in Disaster Management, PSTU, Bangladesh *Founder and Director* of Research Society *Website*: https://researchsociety20.org/founder-and-director/ *E-mails*: founder-and-director at researchsociety20.org jalal at jejunu.ac.kr 20205103002 at nuist.edu.cn -------------- next part -------------- An HTML attachment was scrubbed... URL: From setareh.rahimi at gmail.com Sat Jan 20 06:44:32 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Sat, 20 Jan 2024 17:14:32 +0330 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero Message-ID: Dear NCL users, I have a NetCDF file (attached) that contains some variables. I would like to change the values of all variables to zero for all grid points. So how can I do this using NCL, please? This file is an input file for the WRF-Chem model, and I do not want to change the format of the file, only want to change the values to zero. Please kindly advise me in this regard. Thanks in advance. Kind regards, -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: wrfbiochemi_d03 Type: application/octet-stream Size: 528760 bytes Desc: not available URL: From gthompsn at ucar.edu Sat Jan 20 09:01:26 2024 From: gthompsn at ucar.edu (Greg Thompson) Date: Sat, 20 Jan 2024 09:01:26 -0700 Subject: [ncl-talk] problem with PolyMarker opacity Message-ID: Dear NCL users, I am having a frustrating problem with using gsn_add_polymarker to add filled circle markers to a map object and getting it to acknowledge the gsMarkerOpacityF attribute. No matter what I try, it is not respecting any value of opacity, it is always 100% opaque. Yet there is no problem when It try running the NCL documented example code https://www.ncl.ucar.edu/Applications/Scripts/newcolor_3.ncl but that is not adding the polymarker to a map object. Can anyone suggest a workaround or potential solution to try? Thanks, Greg Thompson (UCAR/JCSDA, formerly NCAR-RAL) -------------- next part -------------- An HTML attachment was scrubbed... URL: From brownrig at ucar.edu Sat Jan 20 09:48:00 2024 From: brownrig at ucar.edu (Rick Brownrigg) Date: Sat, 20 Jan 2024 10:48:00 -0600 Subject: [ncl-talk] problem with PolyMarker opacity In-Reply-To: References: Message-ID: Hi Greg, Can you show us the code snippet you are using? Rick On Saturday, January 20, 2024, Greg Thompson via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > > Dear NCL users, > > I am having a frustrating problem with using > gsn_add_polymarker > to add filled circle markers to a map object and getting it to acknowledge > the > gsMarkerOpacityF > attribute. No matter what I try, it is not respecting any value of > opacity, it is always 100% opaque. > > Yet there is no problem when It try running the NCL documented example code > https://www.ncl.ucar.edu/Applications/Scripts/newcolor_3.ncl > but that is not adding the polymarker to a map object. > > Can anyone suggest a workaround or potential solution to try? > > Thanks, > > Greg Thompson (UCAR/JCSDA, formerly NCAR-RAL) > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From shea at ucar.edu Sun Jan 21 08:14:19 2024 From: shea at ucar.edu (Dennis Shea) Date: Sun, 21 Jan 2024 08:14:19 -0700 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: References: Message-ID: Attached is a script that does what you requested. On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > > Dear NCL users, > > I have a NetCDF file (attached) that contains some variables. I would like > to change the values of all variables to zero for all grid points. So how > can I do this using NCL, please? > > This file is an input file for the WRF-Chem model, and I do not want to > change the format of the file, only want to change the values to zero. > > Please kindly advise me in this regard. > > Thanks in advance. > Kind regards, > -- > 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: setareh2.ncl Type: application/octet-stream Size: 1042 bytes Desc: not available URL: From setareh.rahimi at gmail.com Mon Jan 22 01:19:15 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Mon, 22 Jan 2024 11:49:15 +0330 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: References: Message-ID: So many thanks for your help. I have another question please. If I want to do the same thing only for specific range of lat/lon, how can I do that? For example lat(30,40), and Lon (0,20)? Would you please kindly advise me in this regard? Kind regards. S.Rahimi On Sun, Jan 21, 2024 at 18:44 Dennis Shea wrote: > Attached is a script that does what you requested. > > > On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> >> Dear NCL users, >> >> I have a NetCDF file (attached) that contains some variables. I >> would like to change the values of all variables to zero for all grid >> points. So how can I do this using NCL, please? >> >> This file is an input file for the WRF-Chem model, and I do not want to >> change the format of the file, only want to change the values to zero. >> >> Please kindly advise me in this regard. >> >> Thanks in advance. >> Kind regards, >> -- >> 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 shea at ucar.edu Mon Jan 22 20:28:58 2024 From: shea at ucar.edu (Dennis Shea) Date: Mon, 22 Jan 2024 20:28:58 -0700 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: References: Message-ID: Hi Setareh, Use the seterah2.ncl script as a guide for what you wish to accomplish. Please read some NCL documentation. Then, try to write a test script. The following may be useful. *getind_latlon2d region_ind * *Good luck* *ncl-team* On Mon, Jan 22, 2024 at 1:19?AM Setareh Rahimi wrote: > So many thanks for your help. I have another question please. If I want to > do the same thing only for specific range of lat/lon, how can I do that? > For example lat(30,40), and Lon (0,20)? > > Would you please kindly advise me in this regard? > > Kind regards. > > S.Rahimi > > > > On Sun, Jan 21, 2024 at 18:44 Dennis Shea wrote: > >> Attached is a script that does what you requested. >> >> >> On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >>> >>> Dear NCL users, >>> >>> I have a NetCDF file (attached) that contains some variables. I >>> would like to change the values of all variables to zero for all grid >>> points. So how can I do this using NCL, please? >>> >>> This file is an input file for the WRF-Chem model, and I do not want to >>> change the format of the file, only want to change the values to zero. >>> >>> Please kindly advise me in this regard. >>> >>> Thanks in advance. >>> Kind regards, >>> -- >>> 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 Tue Jan 23 08:09:20 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Tue, 23 Jan 2024 18:39:20 +0330 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: References: Message-ID: Thank you very much for your feedback. On Tue, Jan 23, 2024 at 06:59 Dennis Shea wrote: > Hi Setareh, > > Use the seterah2.ncl script as a guide for what you wish to accomplish. > > Please read some NCL documentation. Then, try to write a test script. > > The following may be useful. > > > > *getind_latlon2d region_ind * > > > *Good luck* > > > *ncl-team* > > > > On Mon, Jan 22, 2024 at 1:19?AM Setareh Rahimi > wrote: > >> So many thanks for your help. I have another question please. If I want >> to do the same thing only for specific range of lat/lon, how can I do that? >> For example lat(30,40), and Lon (0,20)? >> >> Would you please kindly advise me in this regard? >> >> Kind regards. >> >> S.Rahimi >> >> >> >> On Sun, Jan 21, 2024 at 18:44 Dennis Shea wrote: >> >>> Attached is a script that does what you requested. >>> >>> >>> On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk < >>> ncl-talk at mailman.ucar.edu> wrote: >>> >>>> >>>> Dear NCL users, >>>> >>>> I have a NetCDF file (attached) that contains some variables. I >>>> would like to change the values of all variables to zero for all grid >>>> points. So how can I do this using NCL, please? >>>> >>>> This file is an input file for the WRF-Chem model, and I do not want to >>>> change the format of the file, only want to change the values to zero. >>>> >>>> Please kindly advise me in this regard. >>>> >>>> Thanks in advance. >>>> Kind regards, >>>> -- >>>> 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 dmjalal90 at gmail.com Wed Jan 24 01:53:28 2024 From: dmjalal90 at gmail.com (Md. Jalal Uddin) Date: Wed, 24 Jan 2024 17:53:28 +0900 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Dear Sir, Could I have this dataset hadisst.187001-202212.nc? On Fri, 19 Jan 2024 at 19:42, Md. Jalal Uddin wrote: > Thank you, sir, for your kind efforts. Your code works for > hadisst.187001-202212.nc data. However, when I use my sst data from NOAA ( > sst.mnmean.nc), the code does not work. I do not know the reason for > this. Sorry to bother you. > > NOAA sst data: > https://psl.noaa.gov/repository/entry/show?entryid=cac1c2a6-a864-4409-bb77-1fdead8eeb6e > > On Fri, 19 Jan 2024 at 15:32, Md. Jalal Uddin wrote: > >> This is great, thank you so much! >> >> On Wed, 17 Jan 2024 at 05:09, Adam Phillips wrote: >> >>> Hi Jalal, >>> I was able to plug in a monthly SST dataset that I have here, so I was >>> able to run your script. So I think what is happening is that the stippling >>> is so small/faint/sparse that it is hard to see it. I made two changes to >>> your script: >>> 1 - I was getting stippling over land as the rtest function was >>> seemingly not recognizing the _FillValue of the input array. (I find this >>> strange, but have not used rtest before.) To fix this, I put in a where >>> statement and added the _FillValue attribute to your p array: >>> p = rtest(r_sst,42,0) >>> * p at _FillValue = 1.e20* >>> copy_VarCoords(r_sst, p) >>> * p = where(ismissing(r_sst),p at _FillValue,p)* >>> >>> 2 - I modified/added the following gsn_contour_shade options to make the >>> stippling more visible: >>> >>> * opt at gsnShadeFillScaleF = 0.25 ; add extra density >>> opt at gsnShadeFillDotSizeF = 0.0015* >>> >>> I have attached the modified script, along with the plot that I am >>> getting. >>> Hope that helps! >>> Best, >>> Adam >>> >>> >>> On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin >>> wrote: >>> >>>> Thank you very much, sir. I used the following resources: >>>> >>>> opt = True >>>> opt at gsnShadeFillType = "pattern" >>>> opt at gsnShadeLow = 17 ; stipple >>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>> >>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>> plot2 = gsn_csm_contour(wks,p, res2) >>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ >>>> >>>> overlay (plot, plot2) >>>> >>>> draw(plot) >>>> frame(wks) >>>> >>>> However, I did not get output as stipple areas where p value is less >>>> than 0.05. I added the full code including Nino data in the attachment. >>>> Should I send the SST data (sst.mon.mean.nc) via ftp? >>>> >>>> On Sat, 13 Jan 2024 at 06:55, Adam Phillips wrote: >>>> >>>>> Hi Jalal, >>>>> As you are pattern filling, and the default is color filling, you need >>>>> to set: >>>>> opt at gsnShadeFillType = "pattern" >>>>> I am also unclear what you mean by this statement: >>>>> "...to show the significance level as a shaded pattern where p value >>>>> is less than 0.05 and 0.01" >>>>> >>>>> Do you want to stipple areas less than .05? If so the following code >>>>> should do that: >>>>> opt = True >>>>> opt at gsnShade*Low* = 17 ; stipple >>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >>>>> >>>>> If instead you want to stipple areas less than .01, you will need to >>>>> modify the contour levels in your res2 resource list to have a contour at >>>>> .01: >>>>> res2 at cnMinLevelValF = 0.00 ; set min contour level >>>>> res2 at cnMaxLevelValF = .10 ; set max contour level >>>>> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >>>>> >>>>> opt = True >>>>> opt at gsnShade*Low* = 17 ; stipple >>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >>>>> >>>>> Hope that answers your question. As always please respond to the >>>>> ncl-talk email list with further questions. >>>>> Best, >>>>> Adam >>>>> >>>>> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >>>>> ncl-talk at mailman.ucar.edu> wrote: >>>>> >>>>>> Dear all, >>>>>> >>>>>> I have tried to use the gsn_contour_shade function to show the >>>>>> significance level as a shaded pattern where p value is less than 0.05 and >>>>>> 0.01. However, the output map does not show the shaded areas (please see >>>>>> the attached map). >>>>>> >>>>>> Here are some codes: >>>>>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; >>>>>> ccr(lat,lon) >>>>>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>>>>> r_sst at long_name = "Correlation: r_sst" >>>>>> >>>>>> p = rtest(r_sst,42,0) >>>>>> copy_VarCoords(r_sst, p) >>>>>> >>>>>> opt = True >>>>>> opt at gsnShadeHigh = 17 ; stipple >>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>>>>> >>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>> >>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>>>>> >>>>>> overlay (plot, plot2) >>>>>> >>>>>> draw(plot) >>>>>> frame(wks) >>>>>> >>>>>> Please find the full code in the attached file. >>>>>> >>>>>> Any suggestions would be appreciated. >>>>>> >>>>>> -- >>>>>> >>>>>> *Regards,* >>>>>> >>>>>> Dr. Jalal >>>>>> >>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>> University, South Korea >>>>>> >>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>> China >>>>>> >>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>> >>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>> >>>>>> *Founder and Director* of Research Society >>>>>> >>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>> >>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>> >>>>>> jalal at jejunu.ac.kr >>>>>> >>>>>> 20205103002 at nuist.edu.cn >>>>>> _______________________________________________ >>>>>> 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/ >>>>> >>>>> >>>>> >>>> >>>> >>>> -- >>>> >>>> *Regards,* >>>> >>>> Dr. Jalal >>>> >>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>> University, South Korea >>>> >>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>> >>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>> >>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>> >>>> *Founder and Director* of Research Society >>>> >>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>> >>>> *E-mails*: founder-and-director at researchsociety20.org >>>> >>>> jalal at jejunu.ac.kr >>>> >>>> 20205103002 at nuist.edu.cn >>>> >>> >>> >>> -- >>> Adam Phillips >>> Associate Scientist IV, Climate Analysis Section >>> Climate and Global Dynamics Laboratory >>> National Center for Atmospheric Research >>> www.cgd.ucar.edu/staff/asphilli/ >>> >>> >>> >> >> >> -- >> >> *Regards,* >> >> Dr. Jalal >> >> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >> University, South Korea >> >> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >> >> *M.Sc.* in Applied Meteorology, NUIST, China >> >> *B.Sc.* in Disaster Management, PSTU, Bangladesh >> >> *Founder and Director* of Research Society >> >> *Website*: https://researchsociety20.org/founder-and-director/ >> >> *E-mails*: founder-and-director at researchsociety20.org >> >> jalal at jejunu.ac.kr >> >> 20205103002 at nuist.edu.cn >> > > > -- > > *Regards,* > > Dr. Jalal > > *Postdoctoral Researcher,* Typhoon Research Center, Jeju National > University, South Korea > > *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China > > *M.Sc.* in Applied Meteorology, NUIST, China > > *B.Sc.* in Disaster Management, PSTU, Bangladesh > > *Founder and Director* of Research Society > > *Website*: https://researchsociety20.org/founder-and-director/ > > *E-mails*: founder-and-director at researchsociety20.org > > jalal at jejunu.ac.kr > > 20205103002 at nuist.edu.cn > -- *Regards,* Dr. Jalal *Postdoctoral Researcher,* Typhoon Research Center, Jeju National University, South Korea *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China *M.Sc.* in Applied Meteorology, NUIST, China *B.Sc.* in Disaster Management, PSTU, Bangladesh *Founder and Director* of Research Society *Website*: https://researchsociety20.org/founder-and-director/ *E-mails*: founder-and-director at researchsociety20.org jalal at jejunu.ac.kr 20205103002 at nuist.edu.cn -------------- next part -------------- An HTML attachment was scrubbed... URL: From asphilli at ucar.edu Wed Jan 24 10:12:46 2024 From: asphilli at ucar.edu (Adam Phillips) Date: Wed, 24 Jan 2024 10:12:46 -0700 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Hi Jalal, I am guessing the reason why my code didn't work with the NOAA/PSL data is because that data's time coordinate variable did not match that of the HadISST dataset. (In the script I sent you I used coordinate subscripting of the time dimension to select the time needed.) Also, the NOAA/PSL SST dataset you pointed me to only runs from 1981-2012, and your script was set to analyze years 1980-2021. I will write to you offline with how to download the HadISST data file. Best, Adam On Wed, Jan 24, 2024 at 1:53?AM Md. Jalal Uddin wrote: > Dear Sir, > > Could I have this dataset hadisst.187001-202212.nc? > > On Fri, 19 Jan 2024 at 19:42, Md. Jalal Uddin wrote: > >> Thank you, sir, for your kind efforts. Your code works for >> hadisst.187001-202212.nc data. However, when I use my sst data from NOAA >> (sst.mnmean.nc), the code does not work. I do not know the reason for >> this. Sorry to bother you. >> >> NOAA sst data: >> https://psl.noaa.gov/repository/entry/show?entryid=cac1c2a6-a864-4409-bb77-1fdead8eeb6e >> >> On Fri, 19 Jan 2024 at 15:32, Md. Jalal Uddin >> wrote: >> >>> This is great, thank you so much! >>> >>> On Wed, 17 Jan 2024 at 05:09, Adam Phillips wrote: >>> >>>> Hi Jalal, >>>> I was able to plug in a monthly SST dataset that I have here, so I was >>>> able to run your script. So I think what is happening is that the stippling >>>> is so small/faint/sparse that it is hard to see it. I made two changes to >>>> your script: >>>> 1 - I was getting stippling over land as the rtest function was >>>> seemingly not recognizing the _FillValue of the input array. (I find this >>>> strange, but have not used rtest before.) To fix this, I put in a where >>>> statement and added the _FillValue attribute to your p array: >>>> p = rtest(r_sst,42,0) >>>> * p at _FillValue = 1.e20* >>>> copy_VarCoords(r_sst, p) >>>> * p = where(ismissing(r_sst),p at _FillValue,p)* >>>> >>>> 2 - I modified/added the following gsn_contour_shade options to make >>>> the stippling more visible: >>>> >>>> * opt at gsnShadeFillScaleF = 0.25 ; add extra density >>>> opt at gsnShadeFillDotSizeF = 0.0015* >>>> >>>> I have attached the modified script, along with the plot that I am >>>> getting. >>>> Hope that helps! >>>> Best, >>>> Adam >>>> >>>> >>>> On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin >>>> wrote: >>>> >>>>> Thank you very much, sir. I used the following resources: >>>>> >>>>> opt = True >>>>> opt at gsnShadeFillType = "pattern" >>>>> opt at gsnShadeLow = 17 ; stipple >>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>> >>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ >>>>> >>>>> overlay (plot, plot2) >>>>> >>>>> draw(plot) >>>>> frame(wks) >>>>> >>>>> However, I did not get output as stipple areas where p value is less >>>>> than 0.05. I added the full code including Nino data in the attachment. >>>>> Should I send the SST data (sst.mon.mean.nc) via ftp? >>>>> >>>>> On Sat, 13 Jan 2024 at 06:55, Adam Phillips wrote: >>>>> >>>>>> Hi Jalal, >>>>>> As you are pattern filling, and the default is color filling, you >>>>>> need to set: >>>>>> opt at gsnShadeFillType = "pattern" >>>>>> I am also unclear what you mean by this statement: >>>>>> "...to show the significance level as a shaded pattern where p value >>>>>> is less than 0.05 and 0.01" >>>>>> >>>>>> Do you want to stipple areas less than .05? If so the following code >>>>>> should do that: >>>>>> opt = True >>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >>>>>> >>>>>> If instead you want to stipple areas less than .01, you will need to >>>>>> modify the contour levels in your res2 resource list to have a contour at >>>>>> .01: >>>>>> res2 at cnMinLevelValF = 0.00 ; set min contour level >>>>>> res2 at cnMaxLevelValF = .10 ; set max contour level >>>>>> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >>>>>> >>>>>> opt = True >>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >>>>>> >>>>>> Hope that answers your question. As always please respond to the >>>>>> ncl-talk email list with further questions. >>>>>> Best, >>>>>> Adam >>>>>> >>>>>> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >>>>>> ncl-talk at mailman.ucar.edu> wrote: >>>>>> >>>>>>> Dear all, >>>>>>> >>>>>>> I have tried to use the gsn_contour_shade function to show the >>>>>>> significance level as a shaded pattern where p value is less than 0.05 and >>>>>>> 0.01. However, the output map does not show the shaded areas (please see >>>>>>> the attached map). >>>>>>> >>>>>>> Here are some codes: >>>>>>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; >>>>>>> ccr(lat,lon) >>>>>>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>>>>>> r_sst at long_name = "Correlation: r_sst" >>>>>>> >>>>>>> p = rtest(r_sst,42,0) >>>>>>> copy_VarCoords(r_sst, p) >>>>>>> >>>>>>> opt = True >>>>>>> opt at gsnShadeHigh = 17 ; stipple >>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>>>>>> >>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>> >>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>>>>>> >>>>>>> overlay (plot, plot2) >>>>>>> >>>>>>> draw(plot) >>>>>>> frame(wks) >>>>>>> >>>>>>> Please find the full code in the attached file. >>>>>>> >>>>>>> Any suggestions would be appreciated. >>>>>>> >>>>>>> -- >>>>>>> >>>>>>> *Regards,* >>>>>>> >>>>>>> Dr. Jalal >>>>>>> >>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>>> University, South Korea >>>>>>> >>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>> China >>>>>>> >>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>> >>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>> >>>>>>> *Founder and Director* of Research Society >>>>>>> >>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>> >>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>> >>>>>>> jalal at jejunu.ac.kr >>>>>>> >>>>>>> 20205103002 at nuist.edu.cn >>>>>>> _______________________________________________ >>>>>>> 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/ >>>>>> >>>>>> >>>>>> >>>>> >>>>> >>>>> -- >>>>> >>>>> *Regards,* >>>>> >>>>> Dr. Jalal >>>>> >>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>> University, South Korea >>>>> >>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>>> >>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>> >>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>> >>>>> *Founder and Director* of Research Society >>>>> >>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>> >>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>> >>>>> jalal at jejunu.ac.kr >>>>> >>>>> 20205103002 at nuist.edu.cn >>>>> >>>> >>>> >>>> -- >>>> Adam Phillips >>>> Associate Scientist IV, Climate Analysis Section >>>> Climate and Global Dynamics Laboratory >>>> National Center for Atmospheric Research >>>> www.cgd.ucar.edu/staff/asphilli/ >>>> >>>> >>>> >>> >>> >>> -- >>> >>> *Regards,* >>> >>> Dr. Jalal >>> >>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>> University, South Korea >>> >>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>> >>> *M.Sc.* in Applied Meteorology, NUIST, China >>> >>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>> >>> *Founder and Director* of Research Society >>> >>> *Website*: https://researchsociety20.org/founder-and-director/ >>> >>> *E-mails*: founder-and-director at researchsociety20.org >>> >>> jalal at jejunu.ac.kr >>> >>> 20205103002 at nuist.edu.cn >>> >> >> >> -- >> >> *Regards,* >> >> Dr. Jalal >> >> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >> University, South Korea >> >> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >> >> *M.Sc.* in Applied Meteorology, NUIST, China >> >> *B.Sc.* in Disaster Management, PSTU, Bangladesh >> >> *Founder and Director* of Research Society >> >> *Website*: https://researchsociety20.org/founder-and-director/ >> >> *E-mails*: founder-and-director at researchsociety20.org >> >> jalal at jejunu.ac.kr >> >> 20205103002 at nuist.edu.cn >> > > > -- > > *Regards,* > > Dr. Jalal > > *Postdoctoral Researcher,* Typhoon Research Center, Jeju National > University, South Korea > > *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China > > *M.Sc.* in Applied Meteorology, NUIST, China > > *B.Sc.* in Disaster Management, PSTU, Bangladesh > > *Founder and Director* of Research Society > > *Website*: https://researchsociety20.org/founder-and-director/ > > *E-mails*: founder-and-director at researchsociety20.org > > jalal at jejunu.ac.kr > > 20205103002 at nuist.edu.cn > -- 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 dmjalal90 at gmail.com Wed Jan 24 18:31:53 2024 From: dmjalal90 at gmail.com (Md. Jalal Uddin) Date: Thu, 25 Jan 2024 10:31:53 +0900 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Hi Sir, I have now used the same data and code that you provided. However, the output does not show the significant shaded area (see attached figure). And I get the following warnings: *warning*:dim_rmvmean_n: 21375 rightmost sections of the input array contained all missing values *warning*:escorc: Non-fatal conditions encountered in series or xstd equals zero. Possibly, all values of a series are constant. *warning*:escorc: Most likely, one or more series consisted of all constant values I am waiting for your kind advice. On Thu, 25 Jan 2024 at 02:13, Adam Phillips wrote: > Hi Jalal, > I am guessing the reason why my code didn't work with the NOAA/PSL data is > because that data's time coordinate variable did not match that of the > HadISST dataset. (In the script I sent you I used coordinate subscripting > of the time dimension to select the time needed.) Also, the NOAA/PSL SST > dataset you pointed me to only runs from 1981-2012, and your script was set > to analyze years 1980-2021. > > I will write to you offline with how to download the HadISST data file. > Best, > Adam > > > > On Wed, Jan 24, 2024 at 1:53?AM Md. Jalal Uddin > wrote: > >> Dear Sir, >> >> Could I have this dataset hadisst.187001-202212.nc? >> >> On Fri, 19 Jan 2024 at 19:42, Md. Jalal Uddin >> wrote: >> >>> Thank you, sir, for your kind efforts. Your code works for >>> hadisst.187001-202212.nc data. However, when I use my sst data from >>> NOAA (sst.mnmean.nc), the code does not work. I do not know the reason >>> for this. Sorry to bother you. >>> >>> NOAA sst data: >>> https://psl.noaa.gov/repository/entry/show?entryid=cac1c2a6-a864-4409-bb77-1fdead8eeb6e >>> >>> On Fri, 19 Jan 2024 at 15:32, Md. Jalal Uddin >>> wrote: >>> >>>> This is great, thank you so much! >>>> >>>> On Wed, 17 Jan 2024 at 05:09, Adam Phillips wrote: >>>> >>>>> Hi Jalal, >>>>> I was able to plug in a monthly SST dataset that I have here, so I was >>>>> able to run your script. So I think what is happening is that the stippling >>>>> is so small/faint/sparse that it is hard to see it. I made two changes to >>>>> your script: >>>>> 1 - I was getting stippling over land as the rtest function was >>>>> seemingly not recognizing the _FillValue of the input array. (I find this >>>>> strange, but have not used rtest before.) To fix this, I put in a where >>>>> statement and added the _FillValue attribute to your p array: >>>>> p = rtest(r_sst,42,0) >>>>> * p at _FillValue = 1.e20* >>>>> copy_VarCoords(r_sst, p) >>>>> * p = where(ismissing(r_sst),p at _FillValue,p)* >>>>> >>>>> 2 - I modified/added the following gsn_contour_shade options to make >>>>> the stippling more visible: >>>>> >>>>> * opt at gsnShadeFillScaleF = 0.25 ; add extra density >>>>> opt at gsnShadeFillDotSizeF = 0.0015* >>>>> >>>>> I have attached the modified script, along with the plot that I am >>>>> getting. >>>>> Hope that helps! >>>>> Best, >>>>> Adam >>>>> >>>>> >>>>> On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin >>>>> wrote: >>>>> >>>>>> Thank you very much, sir. I used the following resources: >>>>>> >>>>>> opt = True >>>>>> opt at gsnShadeFillType = "pattern" >>>>>> opt at gsnShadeLow = 17 ; stipple >>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>> >>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ >>>>>> >>>>>> overlay (plot, plot2) >>>>>> >>>>>> draw(plot) >>>>>> frame(wks) >>>>>> >>>>>> However, I did not get output as stipple areas where p value is less >>>>>> than 0.05. I added the full code including Nino data in the attachment. >>>>>> Should I send the SST data (sst.mon.mean.nc) via ftp? >>>>>> >>>>>> On Sat, 13 Jan 2024 at 06:55, Adam Phillips >>>>>> wrote: >>>>>> >>>>>>> Hi Jalal, >>>>>>> As you are pattern filling, and the default is color filling, you >>>>>>> need to set: >>>>>>> opt at gsnShadeFillType = "pattern" >>>>>>> I am also unclear what you mean by this statement: >>>>>>> "...to show the significance level as a shaded pattern where p value >>>>>>> is less than 0.05 and 0.01" >>>>>>> >>>>>>> Do you want to stipple areas less than .05? If so the following code >>>>>>> should do that: >>>>>>> opt = True >>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >>>>>>> >>>>>>> If instead you want to stipple areas less than .01, you will need to >>>>>>> modify the contour levels in your res2 resource list to have a contour at >>>>>>> .01: >>>>>>> res2 at cnMinLevelValF = 0.00 ; set min contour level >>>>>>> res2 at cnMaxLevelValF = .10 ; set max contour level >>>>>>> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >>>>>>> >>>>>>> opt = True >>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >>>>>>> >>>>>>> Hope that answers your question. As always please respond to the >>>>>>> ncl-talk email list with further questions. >>>>>>> Best, >>>>>>> Adam >>>>>>> >>>>>>> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >>>>>>> ncl-talk at mailman.ucar.edu> wrote: >>>>>>> >>>>>>>> Dear all, >>>>>>>> >>>>>>>> I have tried to use the gsn_contour_shade function to show the >>>>>>>> significance level as a shaded pattern where p value is less than 0.05 and >>>>>>>> 0.01. However, the output map does not show the shaded areas (please see >>>>>>>> the attached map). >>>>>>>> >>>>>>>> Here are some codes: >>>>>>>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; >>>>>>>> ccr(lat,lon) >>>>>>>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>>>>>>> r_sst at long_name = "Correlation: r_sst" >>>>>>>> >>>>>>>> p = rtest(r_sst,42,0) >>>>>>>> copy_VarCoords(r_sst, p) >>>>>>>> >>>>>>>> opt = True >>>>>>>> opt at gsnShadeHigh = 17 ; stipple >>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>>>>>>> >>>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>>> >>>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>>>>>>> >>>>>>>> overlay (plot, plot2) >>>>>>>> >>>>>>>> draw(plot) >>>>>>>> frame(wks) >>>>>>>> >>>>>>>> Please find the full code in the attached file. >>>>>>>> >>>>>>>> Any suggestions would be appreciated. >>>>>>>> >>>>>>>> -- >>>>>>>> >>>>>>>> *Regards,* >>>>>>>> >>>>>>>> Dr. Jalal >>>>>>>> >>>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>>>> University, South Korea >>>>>>>> >>>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>>> China >>>>>>>> >>>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>>> >>>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>>> >>>>>>>> *Founder and Director* of Research Society >>>>>>>> >>>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>>> >>>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>>> >>>>>>>> jalal at jejunu.ac.kr >>>>>>>> >>>>>>>> 20205103002 at nuist.edu.cn >>>>>>>> _______________________________________________ >>>>>>>> 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/ >>>>>>> >>>>>>> >>>>>>> >>>>>> >>>>>> >>>>>> -- >>>>>> >>>>>> *Regards,* >>>>>> >>>>>> Dr. Jalal >>>>>> >>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>> University, South Korea >>>>>> >>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>> China >>>>>> >>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>> >>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>> >>>>>> *Founder and Director* of Research Society >>>>>> >>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>> >>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>> >>>>>> jalal at jejunu.ac.kr >>>>>> >>>>>> 20205103002 at nuist.edu.cn >>>>>> >>>>> >>>>> >>>>> -- >>>>> Adam Phillips >>>>> Associate Scientist IV, Climate Analysis Section >>>>> Climate and Global Dynamics Laboratory >>>>> National Center for Atmospheric Research >>>>> www.cgd.ucar.edu/staff/asphilli/ >>>>> >>>>> >>>>> >>>> >>>> >>>> -- >>>> >>>> *Regards,* >>>> >>>> Dr. Jalal >>>> >>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>> University, South Korea >>>> >>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>> >>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>> >>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>> >>>> *Founder and Director* of Research Society >>>> >>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>> >>>> *E-mails*: founder-and-director at researchsociety20.org >>>> >>>> jalal at jejunu.ac.kr >>>> >>>> 20205103002 at nuist.edu.cn >>>> >>> >>> >>> -- >>> >>> *Regards,* >>> >>> Dr. Jalal >>> >>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>> University, South Korea >>> >>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>> >>> *M.Sc.* in Applied Meteorology, NUIST, China >>> >>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>> >>> *Founder and Director* of Research Society >>> >>> *Website*: https://researchsociety20.org/founder-and-director/ >>> >>> *E-mails*: founder-and-director at researchsociety20.org >>> >>> jalal at jejunu.ac.kr >>> >>> 20205103002 at nuist.edu.cn >>> >> >> >> -- >> >> *Regards,* >> >> Dr. Jalal >> >> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >> University, South Korea >> >> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >> >> *M.Sc.* in Applied Meteorology, NUIST, China >> >> *B.Sc.* in Disaster Management, PSTU, Bangladesh >> >> *Founder and Director* of Research Society >> >> *Website*: https://researchsociety20.org/founder-and-director/ >> >> *E-mails*: founder-and-director at researchsociety20.org >> >> jalal at jejunu.ac.kr >> >> 20205103002 at nuist.edu.cn >> > > > -- > Adam Phillips > Associate Scientist IV, Climate Analysis Section > Climate and Global Dynamics Laboratory > National Center for Atmospheric Research > www.cgd.ucar.edu/staff/asphilli/ > > > -- *Regards,* Dr. Jalal *Postdoctoral Researcher,* Typhoon Research Center, Jeju National University, South Korea *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China *M.Sc.* in Applied Meteorology, NUIST, China *B.Sc.* in Disaster Management, PSTU, Bangladesh *Founder and Director* of Research Society *Website*: https://researchsociety20.org/founder-and-director/ *E-mails*: founder-and-director at researchsociety20.org jalal at jejunu.ac.kr 20205103002 at nuist.edu.cn -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: sst_SepNov_cor.png Type: image/png Size: 525449 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: adam_sst_corr.ncl Type: application/octet-stream Size: 4613 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: AMM_DMI_NINO_PDO_ELI_SepNovAvg_1980_2021.csv Type: text/csv Size: 2586 bytes Desc: not available URL: From setareh.rahimi at gmail.com Thu Jan 25 07:46:55 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Thu, 25 Jan 2024 18:16:55 +0330 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: References: Message-ID: Dear Dennis, I noticed an issue. When I use the attached script it converts all values to zero, but the data is no georefrenced any more. So how can I fix this please? Thanks in advance. Best regards, S.Rahimi On Sun, Jan 21, 2024 at 18:44 Dennis Shea wrote: > Attached is a script that does what you requested. > > > On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> >> Dear NCL users, >> >> I have a NetCDF file (attached) that contains some variables. I >> would like to change the values of all variables to zero for all grid >> points. So how can I do this using NCL, please? >> >> This file is an input file for the WRF-Chem model, and I do not want to >> change the format of the file, only want to change the values to zero. >> >> Please kindly advise me in this regard. >> >> Thanks in advance. >> Kind regards, >> -- >> 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 asphilli at ucar.edu Thu Jan 25 09:37:48 2024 From: asphilli at ucar.edu (Adam Phillips) Date: Thu, 25 Jan 2024 09:37:48 -0700 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Hi Jalal, The error messages are happening as your SST array has missing data over land, so there are grid points that are completely missing. In this case the error messages can be ignored. With regards to the stippling: When I run the exact script you sent, I do get stippling in numerous areas. I altered the gsnShadeFillScaleF resource to .45 to make it more apparent, and outlined the areas that are significant in the attached plot. If when you run the script adam_sst_corr.ncl you are not getting any stippling in the areas I outlined, I am puzzled. I am running NCL v6.6.2 which is the latest version of NCL. What version are you running? Best, Adam On Wed, Jan 24, 2024 at 6:32?PM Md. Jalal Uddin wrote: > Hi Sir, > > I have now used the same data and code that you provided. However, the > output does not show the significant shaded area (see attached figure). And > I get the following warnings: > *warning*:dim_rmvmean_n: 21375 rightmost sections of the input array > contained all missing values > *warning*:escorc: Non-fatal conditions encountered in series or xstd > equals zero. > Possibly, all values of a series are constant. > *warning*:escorc: Most likely, one or more series consisted of all > constant values > > I am waiting for your kind advice. > > On Thu, 25 Jan 2024 at 02:13, Adam Phillips wrote: > >> Hi Jalal, >> I am guessing the reason why my code didn't work with the NOAA/PSL data >> is because that data's time coordinate variable did not match that of the >> HadISST dataset. (In the script I sent you I used coordinate subscripting >> of the time dimension to select the time needed.) Also, the NOAA/PSL SST >> dataset you pointed me to only runs from 1981-2012, and your script was set >> to analyze years 1980-2021. >> >> I will write to you offline with how to download the HadISST data file. >> Best, >> Adam >> >> >> >> On Wed, Jan 24, 2024 at 1:53?AM Md. Jalal Uddin >> wrote: >> >>> Dear Sir, >>> >>> Could I have this dataset hadisst.187001-202212.nc? >>> >>> On Fri, 19 Jan 2024 at 19:42, Md. Jalal Uddin >>> wrote: >>> >>>> Thank you, sir, for your kind efforts. Your code works for >>>> hadisst.187001-202212.nc data. However, when I use my sst data from >>>> NOAA (sst.mnmean.nc), the code does not work. I do not know the reason >>>> for this. Sorry to bother you. >>>> >>>> NOAA sst data: >>>> https://psl.noaa.gov/repository/entry/show?entryid=cac1c2a6-a864-4409-bb77-1fdead8eeb6e >>>> >>>> On Fri, 19 Jan 2024 at 15:32, Md. Jalal Uddin >>>> wrote: >>>> >>>>> This is great, thank you so much! >>>>> >>>>> On Wed, 17 Jan 2024 at 05:09, Adam Phillips wrote: >>>>> >>>>>> Hi Jalal, >>>>>> I was able to plug in a monthly SST dataset that I have here, so I >>>>>> was able to run your script. So I think what is happening is that the >>>>>> stippling is so small/faint/sparse that it is hard to see it. I made two >>>>>> changes to your script: >>>>>> 1 - I was getting stippling over land as the rtest function was >>>>>> seemingly not recognizing the _FillValue of the input array. (I find this >>>>>> strange, but have not used rtest before.) To fix this, I put in a where >>>>>> statement and added the _FillValue attribute to your p array: >>>>>> p = rtest(r_sst,42,0) >>>>>> * p at _FillValue = 1.e20* >>>>>> copy_VarCoords(r_sst, p) >>>>>> * p = where(ismissing(r_sst),p at _FillValue,p)* >>>>>> >>>>>> 2 - I modified/added the following gsn_contour_shade options to make >>>>>> the stippling more visible: >>>>>> >>>>>> * opt at gsnShadeFillScaleF = 0.25 ; add extra density >>>>>> opt at gsnShadeFillDotSizeF = 0.0015* >>>>>> >>>>>> I have attached the modified script, along with the plot that I am >>>>>> getting. >>>>>> Hope that helps! >>>>>> Best, >>>>>> Adam >>>>>> >>>>>> >>>>>> On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin >>>>>> wrote: >>>>>> >>>>>>> Thank you very much, sir. I used the following resources: >>>>>>> >>>>>>> opt = True >>>>>>> opt at gsnShadeFillType = "pattern" >>>>>>> opt at gsnShadeLow = 17 ; stipple >>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>> >>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ >>>>>>> >>>>>>> overlay (plot, plot2) >>>>>>> >>>>>>> draw(plot) >>>>>>> frame(wks) >>>>>>> >>>>>>> However, I did not get output as stipple areas where p value is less >>>>>>> than 0.05. I added the full code including Nino data in the attachment. >>>>>>> Should I send the SST data (sst.mon.mean.nc) via ftp? >>>>>>> >>>>>>> On Sat, 13 Jan 2024 at 06:55, Adam Phillips >>>>>>> wrote: >>>>>>> >>>>>>>> Hi Jalal, >>>>>>>> As you are pattern filling, and the default is color filling, you >>>>>>>> need to set: >>>>>>>> opt at gsnShadeFillType = "pattern" >>>>>>>> I am also unclear what you mean by this statement: >>>>>>>> "...to show the significance level as a shaded pattern where p >>>>>>>> value is less than 0.05 and 0.01" >>>>>>>> >>>>>>>> Do you want to stipple areas less than .05? If so the following >>>>>>>> code should do that: >>>>>>>> opt = True >>>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >>>>>>>> >>>>>>>> If instead you want to stipple areas less than .01, you will need >>>>>>>> to modify the contour levels in your res2 resource list to have a contour >>>>>>>> at .01: >>>>>>>> res2 at cnMinLevelValF = 0.00 ; set min contour level >>>>>>>> res2 at cnMaxLevelValF = .10 ; set max contour level >>>>>>>> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >>>>>>>> >>>>>>>> opt = True >>>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >>>>>>>> >>>>>>>> Hope that answers your question. As always please respond to the >>>>>>>> ncl-talk email list with further questions. >>>>>>>> Best, >>>>>>>> Adam >>>>>>>> >>>>>>>> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >>>>>>>> ncl-talk at mailman.ucar.edu> wrote: >>>>>>>> >>>>>>>>> Dear all, >>>>>>>>> >>>>>>>>> I have tried to use the gsn_contour_shade function to show the >>>>>>>>> significance level as a shaded pattern where p value is less than 0.05 and >>>>>>>>> 0.01. However, the output map does not show the shaded areas (please see >>>>>>>>> the attached map). >>>>>>>>> >>>>>>>>> Here are some codes: >>>>>>>>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; >>>>>>>>> ccr(lat,lon) >>>>>>>>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>>>>>>>> r_sst at long_name = "Correlation: r_sst" >>>>>>>>> >>>>>>>>> p = rtest(r_sst,42,0) >>>>>>>>> copy_VarCoords(r_sst, p) >>>>>>>>> >>>>>>>>> opt = True >>>>>>>>> opt at gsnShadeHigh = 17 ; stipple >>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>>>>>>>> >>>>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>>>> >>>>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>>>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>>>>>>>> >>>>>>>>> overlay (plot, plot2) >>>>>>>>> >>>>>>>>> draw(plot) >>>>>>>>> frame(wks) >>>>>>>>> >>>>>>>>> Please find the full code in the attached file. >>>>>>>>> >>>>>>>>> Any suggestions would be appreciated. >>>>>>>>> >>>>>>>>> -- >>>>>>>>> >>>>>>>>> *Regards,* >>>>>>>>> >>>>>>>>> Dr. Jalal >>>>>>>>> >>>>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>>>>> University, South Korea >>>>>>>>> >>>>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>>>> China >>>>>>>>> >>>>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>>>> >>>>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>>>> >>>>>>>>> *Founder and Director* of Research Society >>>>>>>>> >>>>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>>>> >>>>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>>>> >>>>>>>>> jalal at jejunu.ac.kr >>>>>>>>> >>>>>>>>> 20205103002 at nuist.edu.cn >>>>>>>>> _______________________________________________ >>>>>>>>> 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/ >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>> >>>>>>> >>>>>>> -- >>>>>>> >>>>>>> *Regards,* >>>>>>> >>>>>>> Dr. Jalal >>>>>>> >>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>>> University, South Korea >>>>>>> >>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>> China >>>>>>> >>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>> >>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>> >>>>>>> *Founder and Director* of Research Society >>>>>>> >>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>> >>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>> >>>>>>> jalal at jejunu.ac.kr >>>>>>> >>>>>>> 20205103002 at nuist.edu.cn >>>>>>> >>>>>> >>>>>> >>>>>> -- >>>>>> Adam Phillips >>>>>> Associate Scientist IV, Climate Analysis Section >>>>>> Climate and Global Dynamics Laboratory >>>>>> National Center for Atmospheric Research >>>>>> www.cgd.ucar.edu/staff/asphilli/ >>>>>> >>>>>> >>>>>> >>>>> >>>>> >>>>> -- >>>>> >>>>> *Regards,* >>>>> >>>>> Dr. Jalal >>>>> >>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>> University, South Korea >>>>> >>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>>> >>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>> >>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>> >>>>> *Founder and Director* of Research Society >>>>> >>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>> >>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>> >>>>> jalal at jejunu.ac.kr >>>>> >>>>> 20205103002 at nuist.edu.cn >>>>> >>>> >>>> >>>> -- >>>> >>>> *Regards,* >>>> >>>> Dr. Jalal >>>> >>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>> University, South Korea >>>> >>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>> >>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>> >>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>> >>>> *Founder and Director* of Research Society >>>> >>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>> >>>> *E-mails*: founder-and-director at researchsociety20.org >>>> >>>> jalal at jejunu.ac.kr >>>> >>>> 20205103002 at nuist.edu.cn >>>> >>> >>> >>> -- >>> >>> *Regards,* >>> >>> Dr. Jalal >>> >>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>> University, South Korea >>> >>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>> >>> *M.Sc.* in Applied Meteorology, NUIST, China >>> >>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>> >>> *Founder and Director* of Research Society >>> >>> *Website*: https://researchsociety20.org/founder-and-director/ >>> >>> *E-mails*: founder-and-director at researchsociety20.org >>> >>> jalal at jejunu.ac.kr >>> >>> 20205103002 at nuist.edu.cn >>> >> >> >> -- >> Adam Phillips >> Associate Scientist IV, Climate Analysis Section >> Climate and Global Dynamics Laboratory >> National Center for Atmospheric Research >> www.cgd.ucar.edu/staff/asphilli/ >> >> >> > > > -- > > *Regards,* > > Dr. Jalal > > *Postdoctoral Researcher,* Typhoon Research Center, Jeju National > University, South Korea > > *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China > > *M.Sc.* in Applied Meteorology, NUIST, China > > *B.Sc.* in Disaster Management, PSTU, Bangladesh > > *Founder and Director* of Research Society > > *Website*: https://researchsociety20.org/founder-and-director/ > > *E-mails*: founder-and-director at researchsociety20.org > > jalal at jejunu.ac.kr > > 20205103002 at nuist.edu.cn > -- 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: -------------- next part -------------- A non-text attachment was scrubbed... Name: sst_SepNov_cor.save.png Type: image/png Size: 647615 bytes Desc: not available URL: From ehsantaghizadeh at yahoo.com Thu Jan 25 11:06:45 2024 From: ehsantaghizadeh at yahoo.com (Ehsan Taghizadeh) Date: Thu, 25 Jan 2024 18:06:45 +0000 (UTC) Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: References: Message-ID: <633146882.489741.1706206005585@mail.yahoo.com> Dear Setareh,I suggest you read the variables from the fout file again. For example, if there is a variable named "V10" in fout, you just need to read that variable from fout, as follows:V10_zero=fout->V10 This V10_zero has the same coordinates as V10 in the original file, but its values are zero.I hope this is correct. ----------------------------------------------------------------------------- Sincerely,Ehsan Taghizade On Thursday, January 25, 2024 at 06:17:12 PM GMT+3:30, Setareh Rahimi via ncl-talk wrote: Dear?Dennis,I noticed an issue. When I use the attached script it converts all values to zero, but the data is no georefrenced any more. So how can I?fix this please?Thanks in advance.Best regards,? S.Rahimi On Sun, Jan 21, 2024 at 18:44 Dennis Shea wrote: Attached is a script that does what you requested. On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk wrote: Dear NCL users, I have a NetCDF file (attached) that?contains?some variables. I would?like to change the values of all variables?to zero for all grid points. So how can I do this using?NCL, please? This file is an input file for the WRF-Chem model, and I do not want to change?the format?of the file, only want to change the values to zero. Please kindly advise me in this regard. Thanks in advance.Kind regards,-- 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 _______________________________________________ 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 dave.allured at noaa.gov Thu Jan 25 12:48:27 2024 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Thu, 25 Jan 2024 12:48:27 -0700 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: <633146882.489741.1706206005585@mail.yahoo.com> References: <633146882.489741.1706206005585@mail.yahoo.com> Message-ID: The problem is that Dennis's program did not copy the file global attributes, which contain part of the georeferencing. Just add this line near the top of the program, after the line "print (vNames)": copy_VarAtts (f, fout) On Thu, Jan 25, 2024 at 11:07?AM Ehsan Taghizadeh via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear Setareh, > I suggest you read the variables from the fout file again. For example, if > there is a variable named "V10" in fout, you just need to read that > variable from fout, as follows: > V10_zero=fout->V10 > > This V10_zero has the same coordinates as V10 in the original file, but > its values are zero. > I hope this is correct. > > > On Thursday, January 25, 2024 at 06:17:12 PM GMT+3:30, Setareh Rahimi via > ncl-talk wrote: > > Dear Dennis, > I noticed an issue. When I use the attached script it converts all values > to zero, but the data is no georefrenced any more. So how can I fix this > please? > Thanks in advance. > Best regards, > > S.Rahimi > > > On Sun, Jan 21, 2024 at 18:44 Dennis Shea wrote: > > Attached is a script that does what you requested. > > > On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > > > Dear NCL users, > > I have a NetCDF file (attached) that contains some variables. I would like > to change the values of all variables to zero for all grid points. So how > can I do this using NCL, please? > > This file is an input file for the WRF-Chem model, and I do not want to > change the format of the file, only want to change the values to zero. > > Please kindly advise me in this regard. > > Thanks in advance. > Kind regards, > -- > S.Rahimi > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From setareh.rahimi at gmail.com Thu Jan 25 13:08:35 2024 From: setareh.rahimi at gmail.com (Setareh Rahimi) Date: Thu, 25 Jan 2024 23:38:35 +0330 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: References: <633146882.489741.1706206005585@mail.yahoo.com> Message-ID: Dear Dave, So many thanks for your help. I just checked it and it worked. May I ask you please to advise me on how to only convert the values of all variables to zero for a specific region and not whole the domain? (for example lat[0,10], and lon[10,30])? Thanks in advance. Kind regards, On Thu, Jan 25, 2024 at 11:18?PM Dave Allured - NOAA Affiliate < dave.allured at noaa.gov> wrote: > The problem is that Dennis's program did not copy the file global > attributes, which contain part of the georeferencing. Just add this line > near the top of the program, after the line "print (vNames)": > > copy_VarAtts (f, fout) > > > On Thu, Jan 25, 2024 at 11:07?AM Ehsan Taghizadeh via ncl-talk < > ncl-talk at mailman.ucar.edu> wrote: > >> Dear Setareh, >> I suggest you read the variables from the fout file again. For example, >> if there is a variable named "V10" in fout, you just need to read that >> variable from fout, as follows: >> V10_zero=fout->V10 >> >> This V10_zero has the same coordinates as V10 in the original file, but >> its values are zero. >> I hope this is correct. >> >> >> On Thursday, January 25, 2024 at 06:17:12 PM GMT+3:30, Setareh Rahimi via >> ncl-talk wrote: >> >> Dear Dennis, >> I noticed an issue. When I use the attached script it converts all values >> to zero, but the data is no georefrenced any more. So how can I fix this >> please? >> Thanks in advance. >> Best regards, >> >> S.Rahimi >> >> >> On Sun, Jan 21, 2024 at 18:44 Dennis Shea wrote: >> >> Attached is a script that does what you requested. >> >> >> On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >> >> Dear NCL users, >> >> I have a NetCDF file (attached) that contains some variables. I >> would like to change the values of all variables to zero for all grid >> points. So how can I do this using NCL, please? >> >> This file is an input file for the WRF-Chem model, and I do not want to >> change the format of the file, only want to change the values to zero. >> >> Please kindly advise me in this regard. >> >> Thanks in advance. >> Kind regards, >> -- >> S.Rahimi >> >> -- S.Rahimi -------------- next part -------------- An HTML attachment was scrubbed... URL: From dave.allured at noaa.gov Thu Jan 25 15:18:08 2024 From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate) Date: Thu, 25 Jan 2024 15:18:08 -0700 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: References: <633146882.489741.1706206005585@mail.yahoo.com> Message-ID: Setareh, with simple 1-D coordinate variables, you would easily do this with coordinate range subscripting. Read about that in the NCL ref manual. Here is an example that writes directly into an existing file, rather than making a new file. f->MSEBIO_ISOP( {0,10}, {10,30} ) = 0 However, the coordinates in your particular file are 3-D, so you have to use extra code. Here is the simplest way. This is also the slowest, but if you are only working with a small amount of data, it is better to keep the code simple. This example also writes directly into the existing file. Not tested: lat1 = 0 lat2 = 10 lon1 = 10 lon2 = 30 do t = 0, ntimes-1 do i = 0, nlats-1 do j = 0, nlons-1 if ( (XLAT(t,i,j) .ge. lat1) .and. (XLAT(t,i,j) .le. lat2) \ .and. (XLONG(t,i,j) .ge. lon1) .and. (XLONG(t,i,j) .le. lon2) ) then f->MSEBIO_ISOP(t,i,j) = 0 ! You can add more variables here end if end do end do end do On Thu, Jan 25, 2024 at 1:04?PM Setareh Rahimi wrote: > Dear Dave, > > So many thanks for your help. I just checked it and it worked. > > May I ask you please to advise me on how to only convert the values of all > variables to zero for a specific region and not whole the domain? (for > example lat[0,10], and lon[10,30])? > > Thanks in advance. > Kind regards, > > On Thu, Jan 25, 2024 at 11:18?PM Dave Allured - NOAA Affiliate < > dave.allured at noaa.gov> wrote: > >> The problem is that Dennis's program did not copy the file global >> attributes, which contain part of the georeferencing. Just add this line >> near the top of the program, after the line "print (vNames)": >> >> copy_VarAtts (f, fout) >> >> >> On Thu, Jan 25, 2024 at 11:07?AM Ehsan Taghizadeh via ncl-talk < >> ncl-talk at mailman.ucar.edu> wrote: >> >>> Dear Setareh, >>> I suggest you read the variables from the fout file again. For example, >>> if there is a variable named "V10" in fout, you just need to read that >>> variable from fout, as follows: >>> V10_zero=fout->V10 >>> >>> This V10_zero has the same coordinates as V10 in the original file, but >>> its values are zero. >>> I hope this is correct. >>> >>> >>> On Thursday, January 25, 2024 at 06:17:12 PM GMT+3:30, Setareh Rahimi >>> via ncl-talk wrote: >>> >>> Dear Dennis, >>> I noticed an issue. When I use the attached script it converts all >>> values to zero, but the data is no georefrenced any more. So how can I fix >>> this please? >>> Thanks in advance. >>> Best regards, >>> >>> S.Rahimi >>> >>> >>> On Sun, Jan 21, 2024 at 18:44 Dennis Shea wrote: >>> >>> Attached is a script that does what you requested. >>> >>> >>> On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk < >>> ncl-talk at mailman.ucar.edu> wrote: >>> >>> >>> Dear NCL users, >>> >>> I have a NetCDF file (attached) that contains some variables. I >>> would like to change the values of all variables to zero for all grid >>> points. So how can I do this using NCL, please? >>> >>> This file is an input file for the WRF-Chem model, and I do not want to >>> change the format of the file, only want to change the values to zero. >>> >>> Please kindly advise me in this regard. >>> >>> Thanks in advance. >>> Kind regards, >>> -- >>> S.Rahimi >>> >>> -------------- next part -------------- An HTML attachment was scrubbed... URL: From dmjalal90 at gmail.com Thu Jan 25 19:39:14 2024 From: dmjalal90 at gmail.com (Md. Jalal Uddin) Date: Fri, 26 Jan 2024 11:39:14 +0900 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Thank you, sir. I am using NCL v6.6.2 in Ubuntu (Windows Subsystem for Linux). On Fri, 26 Jan 2024 at 01:38, Adam Phillips wrote: > Hi Jalal, > The error messages are happening as your SST array has missing data over > land, so there are grid points that are completely missing. In this case > the error messages can be ignored. With regards to the stippling: When I > run the exact script you sent, I do get stippling in numerous areas. I > altered the gsnShadeFillScaleF resource to .45 to make it more apparent, > and outlined the areas that are significant in the attached plot. > > If when you run the script adam_sst_corr.ncl you are not getting any > stippling in the areas I outlined, I am puzzled. I am running NCL v6.6.2 > which is the latest version of NCL. What version are you running? > Best, > Adam > > > On Wed, Jan 24, 2024 at 6:32?PM Md. Jalal Uddin > wrote: > >> Hi Sir, >> >> I have now used the same data and code that you provided. However, the >> output does not show the significant shaded area (see attached figure). And >> I get the following warnings: >> *warning*:dim_rmvmean_n: 21375 rightmost sections of the input array >> contained all missing values >> *warning*:escorc: Non-fatal conditions encountered in series or xstd >> equals zero. >> Possibly, all values of a series are constant. >> *warning*:escorc: Most likely, one or more series consisted of all >> constant values >> >> I am waiting for your kind advice. >> >> On Thu, 25 Jan 2024 at 02:13, Adam Phillips wrote: >> >>> Hi Jalal, >>> I am guessing the reason why my code didn't work with the NOAA/PSL data >>> is because that data's time coordinate variable did not match that of the >>> HadISST dataset. (In the script I sent you I used coordinate subscripting >>> of the time dimension to select the time needed.) Also, the NOAA/PSL SST >>> dataset you pointed me to only runs from 1981-2012, and your script was set >>> to analyze years 1980-2021. >>> >>> I will write to you offline with how to download the HadISST data file. >>> Best, >>> Adam >>> >>> >>> >>> On Wed, Jan 24, 2024 at 1:53?AM Md. Jalal Uddin >>> wrote: >>> >>>> Dear Sir, >>>> >>>> Could I have this dataset hadisst.187001-202212.nc? >>>> >>>> On Fri, 19 Jan 2024 at 19:42, Md. Jalal Uddin >>>> wrote: >>>> >>>>> Thank you, sir, for your kind efforts. Your code works for >>>>> hadisst.187001-202212.nc data. However, when I use my sst data from >>>>> NOAA (sst.mnmean.nc), the code does not work. I do not know the >>>>> reason for this. Sorry to bother you. >>>>> >>>>> NOAA sst data: >>>>> https://psl.noaa.gov/repository/entry/show?entryid=cac1c2a6-a864-4409-bb77-1fdead8eeb6e >>>>> >>>>> On Fri, 19 Jan 2024 at 15:32, Md. Jalal Uddin >>>>> wrote: >>>>> >>>>>> This is great, thank you so much! >>>>>> >>>>>> On Wed, 17 Jan 2024 at 05:09, Adam Phillips >>>>>> wrote: >>>>>> >>>>>>> Hi Jalal, >>>>>>> I was able to plug in a monthly SST dataset that I have here, so I >>>>>>> was able to run your script. So I think what is happening is that the >>>>>>> stippling is so small/faint/sparse that it is hard to see it. I made two >>>>>>> changes to your script: >>>>>>> 1 - I was getting stippling over land as the rtest function was >>>>>>> seemingly not recognizing the _FillValue of the input array. (I find this >>>>>>> strange, but have not used rtest before.) To fix this, I put in a where >>>>>>> statement and added the _FillValue attribute to your p array: >>>>>>> p = rtest(r_sst,42,0) >>>>>>> * p at _FillValue = 1.e20* >>>>>>> copy_VarCoords(r_sst, p) >>>>>>> * p = where(ismissing(r_sst),p at _FillValue,p)* >>>>>>> >>>>>>> 2 - I modified/added the following gsn_contour_shade options to make >>>>>>> the stippling more visible: >>>>>>> >>>>>>> * opt at gsnShadeFillScaleF = 0.25 ; add extra density >>>>>>> opt at gsnShadeFillDotSizeF = 0.0015* >>>>>>> >>>>>>> I have attached the modified script, along with the plot that I am >>>>>>> getting. >>>>>>> Hope that helps! >>>>>>> Best, >>>>>>> Adam >>>>>>> >>>>>>> >>>>>>> On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin >>>>>>> wrote: >>>>>>> >>>>>>>> Thank you very much, sir. I used the following resources: >>>>>>>> >>>>>>>> opt = True >>>>>>>> opt at gsnShadeFillType = "pattern" >>>>>>>> opt at gsnShadeLow = 17 ; stipple >>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>> >>>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ >>>>>>>> >>>>>>>> overlay (plot, plot2) >>>>>>>> >>>>>>>> draw(plot) >>>>>>>> frame(wks) >>>>>>>> >>>>>>>> However, I did not get output as stipple areas where p value is >>>>>>>> less than 0.05. I added the full code including Nino data in the >>>>>>>> attachment. Should I send the SST data (sst.mon.mean.nc) via ftp? >>>>>>>> >>>>>>>> On Sat, 13 Jan 2024 at 06:55, Adam Phillips >>>>>>>> wrote: >>>>>>>> >>>>>>>>> Hi Jalal, >>>>>>>>> As you are pattern filling, and the default is color filling, you >>>>>>>>> need to set: >>>>>>>>> opt at gsnShadeFillType = "pattern" >>>>>>>>> I am also unclear what you mean by this statement: >>>>>>>>> "...to show the significance level as a shaded pattern where p >>>>>>>>> value is less than 0.05 and 0.01" >>>>>>>>> >>>>>>>>> Do you want to stipple areas less than .05? If so the following >>>>>>>>> code should do that: >>>>>>>>> opt = True >>>>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >>>>>>>>> >>>>>>>>> If instead you want to stipple areas less than .01, you will need >>>>>>>>> to modify the contour levels in your res2 resource list to have a contour >>>>>>>>> at .01: >>>>>>>>> res2 at cnMinLevelValF = 0.00 ; set min contour level >>>>>>>>> res2 at cnMaxLevelValF = .10 ; set max contour level >>>>>>>>> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >>>>>>>>> >>>>>>>>> opt = True >>>>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >>>>>>>>> >>>>>>>>> Hope that answers your question. As always please respond to the >>>>>>>>> ncl-talk email list with further questions. >>>>>>>>> Best, >>>>>>>>> Adam >>>>>>>>> >>>>>>>>> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >>>>>>>>> ncl-talk at mailman.ucar.edu> wrote: >>>>>>>>> >>>>>>>>>> Dear all, >>>>>>>>>> >>>>>>>>>> I have tried to use the gsn_contour_shade function to show the >>>>>>>>>> significance level as a shaded pattern where p value is less than 0.05 and >>>>>>>>>> 0.01. However, the output map does not show the shaded areas (please see >>>>>>>>>> the attached map). >>>>>>>>>> >>>>>>>>>> Here are some codes: >>>>>>>>>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) ; >>>>>>>>>> ccr(lat,lon) >>>>>>>>>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>>>>>>>>> r_sst at long_name = "Correlation: r_sst" >>>>>>>>>> >>>>>>>>>> p = rtest(r_sst,42,0) >>>>>>>>>> copy_VarCoords(r_sst, p) >>>>>>>>>> >>>>>>>>>> opt = True >>>>>>>>>> opt at gsnShadeHigh = 17 ; stipple >>>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>>>>>>>>> >>>>>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>>>>> >>>>>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>>>>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>>>>>>>>> >>>>>>>>>> overlay (plot, plot2) >>>>>>>>>> >>>>>>>>>> draw(plot) >>>>>>>>>> frame(wks) >>>>>>>>>> >>>>>>>>>> Please find the full code in the attached file. >>>>>>>>>> >>>>>>>>>> Any suggestions would be appreciated. >>>>>>>>>> >>>>>>>>>> -- >>>>>>>>>> >>>>>>>>>> *Regards,* >>>>>>>>>> >>>>>>>>>> Dr. Jalal >>>>>>>>>> >>>>>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju >>>>>>>>>> National University, South Korea >>>>>>>>>> >>>>>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>>>>> China >>>>>>>>>> >>>>>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>>>>> >>>>>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>>>>> >>>>>>>>>> *Founder and Director* of Research Society >>>>>>>>>> >>>>>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>>>>> >>>>>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>>>>> >>>>>>>>>> jalal at jejunu.ac.kr >>>>>>>>>> >>>>>>>>>> 20205103002 at nuist.edu.cn >>>>>>>>>> _______________________________________________ >>>>>>>>>> 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/ >>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> -- >>>>>>>> >>>>>>>> *Regards,* >>>>>>>> >>>>>>>> Dr. Jalal >>>>>>>> >>>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>>>> University, South Korea >>>>>>>> >>>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>>> China >>>>>>>> >>>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>>> >>>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>>> >>>>>>>> *Founder and Director* of Research Society >>>>>>>> >>>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>>> >>>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>>> >>>>>>>> jalal at jejunu.ac.kr >>>>>>>> >>>>>>>> 20205103002 at nuist.edu.cn >>>>>>>> >>>>>>> >>>>>>> >>>>>>> -- >>>>>>> Adam Phillips >>>>>>> Associate Scientist IV, Climate Analysis Section >>>>>>> Climate and Global Dynamics Laboratory >>>>>>> National Center for Atmospheric Research >>>>>>> www.cgd.ucar.edu/staff/asphilli/ >>>>>>> >>>>>>> >>>>>>> >>>>>> >>>>>> >>>>>> -- >>>>>> >>>>>> *Regards,* >>>>>> >>>>>> Dr. Jalal >>>>>> >>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>> University, South Korea >>>>>> >>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>> China >>>>>> >>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>> >>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>> >>>>>> *Founder and Director* of Research Society >>>>>> >>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>> >>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>> >>>>>> jalal at jejunu.ac.kr >>>>>> >>>>>> 20205103002 at nuist.edu.cn >>>>>> >>>>> >>>>> >>>>> -- >>>>> >>>>> *Regards,* >>>>> >>>>> Dr. Jalal >>>>> >>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>> University, South Korea >>>>> >>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>>> >>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>> >>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>> >>>>> *Founder and Director* of Research Society >>>>> >>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>> >>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>> >>>>> jalal at jejunu.ac.kr >>>>> >>>>> 20205103002 at nuist.edu.cn >>>>> >>>> >>>> >>>> -- >>>> >>>> *Regards,* >>>> >>>> Dr. Jalal >>>> >>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>> University, South Korea >>>> >>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>> >>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>> >>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>> >>>> *Founder and Director* of Research Society >>>> >>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>> >>>> *E-mails*: founder-and-director at researchsociety20.org >>>> >>>> jalal at jejunu.ac.kr >>>> >>>> 20205103002 at nuist.edu.cn >>>> >>> >>> >>> -- >>> Adam Phillips >>> Associate Scientist IV, Climate Analysis Section >>> Climate and Global Dynamics Laboratory >>> National Center for Atmospheric Research >>> www.cgd.ucar.edu/staff/asphilli/ >>> >>> >>> >> >> >> -- >> >> *Regards,* >> >> Dr. Jalal >> >> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >> University, South Korea >> >> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >> >> *M.Sc.* in Applied Meteorology, NUIST, China >> >> *B.Sc.* in Disaster Management, PSTU, Bangladesh >> >> *Founder and Director* of Research Society >> >> *Website*: https://researchsociety20.org/founder-and-director/ >> >> *E-mails*: founder-and-director at researchsociety20.org >> >> jalal at jejunu.ac.kr >> >> 20205103002 at nuist.edu.cn >> > > > -- > Adam Phillips > Associate Scientist IV, Climate Analysis Section > Climate and Global Dynamics Laboratory > National Center for Atmospheric Research > www.cgd.ucar.edu/staff/asphilli/ > > > -- *Regards,* Dr. Jalal *Postdoctoral Researcher,* Typhoon Research Center, Jeju National University, South Korea *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China *M.Sc.* in Applied Meteorology, NUIST, China *B.Sc.* in Disaster Management, PSTU, Bangladesh *Founder and Director* of Research Society *Website*: https://researchsociety20.org/founder-and-director/ *E-mails*: founder-and-director at researchsociety20.org jalal at jejunu.ac.kr 20205103002 at nuist.edu.cn -------------- next part -------------- An HTML attachment was scrubbed... URL: From shea at ucar.edu Thu Jan 25 20:05:58 2024 From: shea at ucar.edu (Dennis Shea) Date: Thu, 25 Jan 2024 20:05:58 -0700 Subject: [ncl-talk] How to convert the value of a NetCDF file to zero In-Reply-To: References: <633146882.489741.1706206005585@mail.yahoo.com> Message-ID: Dave ... As always .... THX ------------------------- Attached is an *untested* version of my original script with Dave's suggested modifications. Again, I do not have a local NCL version so I can not test it. Note: because the WRF coordinate dimensions are curvilinear, the subscripts may return a larger area than 1D coordinate variables. You may have to come up with your own version. Good Luck D On Thu, Jan 25, 2024 at 3:18?PM Dave Allured - NOAA Affiliate via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Setareh, with simple 1-D coordinate variables, you would easily do this > with coordinate range subscripting. Read about that in the NCL ref > manual. Here is an example that writes directly into an existing file, > rather than making a new file. > > f->MSEBIO_ISOP( {0,10}, {10,30} ) = 0 > > However, the coordinates in your particular file are 3-D, so you have to > use extra code. Here is the simplest way. This is also the slowest, but > if you are only working with a small amount of data, it is better to keep > the code simple. This example also writes directly into the existing > file. Not tested: > > lat1 = 0 > lat2 = 10 > lon1 = 10 > lon2 = 30 > > do t = 0, ntimes-1 > do i = 0, nlats-1 > do j = 0, nlons-1 > if ( (XLAT(t,i,j) .ge. lat1) .and. (XLAT(t,i,j) .le. lat2) \ > .and. (XLONG(t,i,j) .ge. lon1) .and. (XLONG(t,i,j) .le. > lon2) ) then > f->MSEBIO_ISOP(t,i,j) = 0 > ! You can add more variables here > end if > end do > end do > end do > > > On Thu, Jan 25, 2024 at 1:04?PM Setareh Rahimi > wrote: > >> Dear Dave, >> >> So many thanks for your help. I just checked it and it worked. >> >> May I ask you please to advise me on how to only convert the values of >> all variables to zero for a specific region and not whole the domain? (for >> example lat[0,10], and lon[10,30])? >> >> Thanks in advance. >> Kind regards, >> >> On Thu, Jan 25, 2024 at 11:18?PM Dave Allured - NOAA Affiliate < >> dave.allured at noaa.gov> wrote: >> >>> The problem is that Dennis's program did not copy the file global >>> attributes, which contain part of the georeferencing. Just add this line >>> near the top of the program, after the line "print (vNames)": >>> >>> copy_VarAtts (f, fout) >>> >>> >>> On Thu, Jan 25, 2024 at 11:07?AM Ehsan Taghizadeh via ncl-talk < >>> ncl-talk at mailman.ucar.edu> wrote: >>> >>>> Dear Setareh, >>>> I suggest you read the variables from the fout file again. For example, >>>> if there is a variable named "V10" in fout, you just need to read that >>>> variable from fout, as follows: >>>> V10_zero=fout->V10 >>>> >>>> This V10_zero has the same coordinates as V10 in the original file, but >>>> its values are zero. >>>> I hope this is correct. >>>> >>>> >>>> On Thursday, January 25, 2024 at 06:17:12 PM GMT+3:30, Setareh Rahimi >>>> via ncl-talk wrote: >>>> >>>> Dear Dennis, >>>> I noticed an issue. When I use the attached script it converts all >>>> values to zero, but the data is no georefrenced any more. So how can I fix >>>> this please? >>>> Thanks in advance. >>>> Best regards, >>>> >>>> S.Rahimi >>>> >>>> >>>> On Sun, Jan 21, 2024 at 18:44 Dennis Shea wrote: >>>> >>>> Attached is a script that does what you requested. >>>> >>>> >>>> On Sat, Jan 20, 2024 at 6:40?AM Setareh Rahimi via ncl-talk < >>>> ncl-talk at mailman.ucar.edu> wrote: >>>> >>>> >>>> Dear NCL users, >>>> >>>> I have a NetCDF file (attached) that contains some variables. I >>>> would like to change the values of all variables to zero for all grid >>>> points. So how can I do this using NCL, please? >>>> >>>> This file is an input file for the WRF-Chem model, and I do not want to >>>> change the format of the file, only want to change the values to zero. >>>> >>>> Please kindly advise me in this regard. >>>> >>>> Thanks in advance. >>>> Kind regards, >>>> -- >>>> 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: setareh3.ncl Type: application/octet-stream Size: 1989 bytes Desc: not available URL: From asphilli at ucar.edu Fri Jan 26 16:52:04 2024 From: asphilli at ucar.edu (Adam Phillips) Date: Fri, 26 Jan 2024 16:52:04 -0700 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Hi Jalal, As you are running the exact version of NCL as I am, I do not have an answer for you as to why the stippling is showing up for me but not for you. (That should not happen.) If you replace the draw(plot) line with draw(plot2) do you get a plot that shows the stippling? The only other idea I have is to bypass using gsn_contour_shade, which will allow you to color the stippling a color other than black. I have attached a new script along with the resulting graphic. Best, Adam On Thu, Jan 25, 2024 at 7:39?PM Md. Jalal Uddin wrote: > Thank you, sir. I am using NCL v6.6.2 in Ubuntu (Windows Subsystem for > Linux). > > On Fri, 26 Jan 2024 at 01:38, Adam Phillips wrote: > >> Hi Jalal, >> The error messages are happening as your SST array has missing data over >> land, so there are grid points that are completely missing. In this case >> the error messages can be ignored. With regards to the stippling: When I >> run the exact script you sent, I do get stippling in numerous areas. I >> altered the gsnShadeFillScaleF resource to .45 to make it more apparent, >> and outlined the areas that are significant in the attached plot. >> >> If when you run the script adam_sst_corr.ncl you are not getting any >> stippling in the areas I outlined, I am puzzled. I am running NCL v6.6.2 >> which is the latest version of NCL. What version are you running? >> Best, >> Adam >> >> >> On Wed, Jan 24, 2024 at 6:32?PM Md. Jalal Uddin >> wrote: >> >>> Hi Sir, >>> >>> I have now used the same data and code that you provided. However, the >>> output does not show the significant shaded area (see attached figure). And >>> I get the following warnings: >>> *warning*:dim_rmvmean_n: 21375 rightmost sections of the input array >>> contained all missing values >>> *warning*:escorc: Non-fatal conditions encountered in series or xstd >>> equals zero. >>> Possibly, all values of a series are constant. >>> *warning*:escorc: Most likely, one or more series consisted of all >>> constant values >>> >>> I am waiting for your kind advice. >>> >>> On Thu, 25 Jan 2024 at 02:13, Adam Phillips wrote: >>> >>>> Hi Jalal, >>>> I am guessing the reason why my code didn't work with the NOAA/PSL data >>>> is because that data's time coordinate variable did not match that of the >>>> HadISST dataset. (In the script I sent you I used coordinate subscripting >>>> of the time dimension to select the time needed.) Also, the NOAA/PSL SST >>>> dataset you pointed me to only runs from 1981-2012, and your script was set >>>> to analyze years 1980-2021. >>>> >>>> I will write to you offline with how to download the HadISST data file. >>>> Best, >>>> Adam >>>> >>>> >>>> >>>> On Wed, Jan 24, 2024 at 1:53?AM Md. Jalal Uddin >>>> wrote: >>>> >>>>> Dear Sir, >>>>> >>>>> Could I have this dataset hadisst.187001-202212.nc? >>>>> >>>>> On Fri, 19 Jan 2024 at 19:42, Md. Jalal Uddin >>>>> wrote: >>>>> >>>>>> Thank you, sir, for your kind efforts. Your code works for >>>>>> hadisst.187001-202212.nc data. However, when I use my sst data from >>>>>> NOAA (sst.mnmean.nc), the code does not work. I do not know the >>>>>> reason for this. Sorry to bother you. >>>>>> >>>>>> NOAA sst data: >>>>>> https://psl.noaa.gov/repository/entry/show?entryid=cac1c2a6-a864-4409-bb77-1fdead8eeb6e >>>>>> >>>>>> On Fri, 19 Jan 2024 at 15:32, Md. Jalal Uddin >>>>>> wrote: >>>>>> >>>>>>> This is great, thank you so much! >>>>>>> >>>>>>> On Wed, 17 Jan 2024 at 05:09, Adam Phillips >>>>>>> wrote: >>>>>>> >>>>>>>> Hi Jalal, >>>>>>>> I was able to plug in a monthly SST dataset that I have here, so I >>>>>>>> was able to run your script. So I think what is happening is that the >>>>>>>> stippling is so small/faint/sparse that it is hard to see it. I made two >>>>>>>> changes to your script: >>>>>>>> 1 - I was getting stippling over land as the rtest function was >>>>>>>> seemingly not recognizing the _FillValue of the input array. (I find this >>>>>>>> strange, but have not used rtest before.) To fix this, I put in a where >>>>>>>> statement and added the _FillValue attribute to your p array: >>>>>>>> p = rtest(r_sst,42,0) >>>>>>>> * p at _FillValue = 1.e20* >>>>>>>> copy_VarCoords(r_sst, p) >>>>>>>> * p = where(ismissing(r_sst),p at _FillValue,p)* >>>>>>>> >>>>>>>> 2 - I modified/added the following gsn_contour_shade options to >>>>>>>> make the stippling more visible: >>>>>>>> >>>>>>>> * opt at gsnShadeFillScaleF = 0.25 ; add extra density >>>>>>>> opt at gsnShadeFillDotSizeF = 0.0015* >>>>>>>> >>>>>>>> I have attached the modified script, along with the plot that I am >>>>>>>> getting. >>>>>>>> Hope that helps! >>>>>>>> Best, >>>>>>>> Adam >>>>>>>> >>>>>>>> >>>>>>>> On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin < >>>>>>>> dmjalal90 at gmail.com> wrote: >>>>>>>> >>>>>>>>> Thank you very much, sir. I used the following resources: >>>>>>>>> >>>>>>>>> opt = True >>>>>>>>> opt at gsnShadeFillType = "pattern" >>>>>>>>> opt at gsnShadeLow = 17 ; stipple >>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>> >>>>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ >>>>>>>>> >>>>>>>>> overlay (plot, plot2) >>>>>>>>> >>>>>>>>> draw(plot) >>>>>>>>> frame(wks) >>>>>>>>> >>>>>>>>> However, I did not get output as stipple areas where p value is >>>>>>>>> less than 0.05. I added the full code including Nino data in the >>>>>>>>> attachment. Should I send the SST data (sst.mon.mean.nc) via ftp? >>>>>>>>> >>>>>>>>> On Sat, 13 Jan 2024 at 06:55, Adam Phillips >>>>>>>>> wrote: >>>>>>>>> >>>>>>>>>> Hi Jalal, >>>>>>>>>> As you are pattern filling, and the default is color filling, you >>>>>>>>>> need to set: >>>>>>>>>> opt at gsnShadeFillType = "pattern" >>>>>>>>>> I am also unclear what you mean by this statement: >>>>>>>>>> "...to show the significance level as a shaded pattern where p >>>>>>>>>> value is less than 0.05 and 0.01" >>>>>>>>>> >>>>>>>>>> Do you want to stipple areas less than .05? If so the following >>>>>>>>>> code should do that: >>>>>>>>>> opt = True >>>>>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >>>>>>>>>> >>>>>>>>>> If instead you want to stipple areas less than .01, you will need >>>>>>>>>> to modify the contour levels in your res2 resource list to have a contour >>>>>>>>>> at .01: >>>>>>>>>> res2 at cnMinLevelValF = 0.00 ; set min contour level >>>>>>>>>> res2 at cnMaxLevelValF = .10 ; set max contour level >>>>>>>>>> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >>>>>>>>>> >>>>>>>>>> opt = True >>>>>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>>> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >>>>>>>>>> >>>>>>>>>> Hope that answers your question. As always please respond to the >>>>>>>>>> ncl-talk email list with further questions. >>>>>>>>>> Best, >>>>>>>>>> Adam >>>>>>>>>> >>>>>>>>>> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >>>>>>>>>> ncl-talk at mailman.ucar.edu> wrote: >>>>>>>>>> >>>>>>>>>>> Dear all, >>>>>>>>>>> >>>>>>>>>>> I have tried to use the gsn_contour_shade function to show the >>>>>>>>>>> significance level as a shaded pattern where p value is less than 0.05 and >>>>>>>>>>> 0.01. However, the output map does not show the shaded areas (please see >>>>>>>>>>> the attached map). >>>>>>>>>>> >>>>>>>>>>> Here are some codes: >>>>>>>>>>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) >>>>>>>>>>> ; ccr(lat,lon) >>>>>>>>>>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>>>>>>>>>> r_sst at long_name = "Correlation: r_sst" >>>>>>>>>>> >>>>>>>>>>> p = rtest(r_sst,42,0) >>>>>>>>>>> copy_VarCoords(r_sst, p) >>>>>>>>>>> >>>>>>>>>>> opt = True >>>>>>>>>>> opt at gsnShadeHigh = 17 ; stipple >>>>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>>>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>>>>>>>>>> >>>>>>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>>>>>> >>>>>>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>>>>>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>>>>>>>>>> >>>>>>>>>>> overlay (plot, plot2) >>>>>>>>>>> >>>>>>>>>>> draw(plot) >>>>>>>>>>> frame(wks) >>>>>>>>>>> >>>>>>>>>>> Please find the full code in the attached file. >>>>>>>>>>> >>>>>>>>>>> Any suggestions would be appreciated. >>>>>>>>>>> >>>>>>>>>>> -- >>>>>>>>>>> >>>>>>>>>>> *Regards,* >>>>>>>>>>> >>>>>>>>>>> Dr. Jalal >>>>>>>>>>> >>>>>>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju >>>>>>>>>>> National University, South Korea >>>>>>>>>>> >>>>>>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, >>>>>>>>>>> NUIST, China >>>>>>>>>>> >>>>>>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>>>>>> >>>>>>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>>>>>> >>>>>>>>>>> *Founder and Director* of Research Society >>>>>>>>>>> >>>>>>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>>>>>> >>>>>>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>>>>>> >>>>>>>>>>> jalal at jejunu.ac.kr >>>>>>>>>>> >>>>>>>>>>> 20205103002 at nuist.edu.cn >>>>>>>>>>> _______________________________________________ >>>>>>>>>>> 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/ >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>>> -- >>>>>>>>> >>>>>>>>> *Regards,* >>>>>>>>> >>>>>>>>> Dr. Jalal >>>>>>>>> >>>>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>>>>> University, South Korea >>>>>>>>> >>>>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>>>> China >>>>>>>>> >>>>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>>>> >>>>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>>>> >>>>>>>>> *Founder and Director* of Research Society >>>>>>>>> >>>>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>>>> >>>>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>>>> >>>>>>>>> jalal at jejunu.ac.kr >>>>>>>>> >>>>>>>>> 20205103002 at nuist.edu.cn >>>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> -- >>>>>>>> Adam Phillips >>>>>>>> Associate Scientist IV, Climate Analysis Section >>>>>>>> Climate and Global Dynamics Laboratory >>>>>>>> National Center for Atmospheric Research >>>>>>>> www.cgd.ucar.edu/staff/asphilli/ >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>> >>>>>>> >>>>>>> -- >>>>>>> >>>>>>> *Regards,* >>>>>>> >>>>>>> Dr. Jalal >>>>>>> >>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>>> University, South Korea >>>>>>> >>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>> China >>>>>>> >>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>> >>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>> >>>>>>> *Founder and Director* of Research Society >>>>>>> >>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>> >>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>> >>>>>>> jalal at jejunu.ac.kr >>>>>>> >>>>>>> 20205103002 at nuist.edu.cn >>>>>>> >>>>>> >>>>>> >>>>>> -- >>>>>> >>>>>> *Regards,* >>>>>> >>>>>> Dr. Jalal >>>>>> >>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>> University, South Korea >>>>>> >>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>> China >>>>>> >>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>> >>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>> >>>>>> *Founder and Director* of Research Society >>>>>> >>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>> >>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>> >>>>>> jalal at jejunu.ac.kr >>>>>> >>>>>> 20205103002 at nuist.edu.cn >>>>>> >>>>> >>>>> >>>>> -- >>>>> >>>>> *Regards,* >>>>> >>>>> Dr. Jalal >>>>> >>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>> University, South Korea >>>>> >>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>>> >>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>> >>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>> >>>>> *Founder and Director* of Research Society >>>>> >>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>> >>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>> >>>>> jalal at jejunu.ac.kr >>>>> >>>>> 20205103002 at nuist.edu.cn >>>>> >>>> >>>> >>>> -- >>>> Adam Phillips >>>> Associate Scientist IV, Climate Analysis Section >>>> Climate and Global Dynamics Laboratory >>>> National Center for Atmospheric Research >>>> www.cgd.ucar.edu/staff/asphilli/ >>>> >>>> >>>> >>> >>> >>> -- >>> >>> *Regards,* >>> >>> Dr. Jalal >>> >>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>> University, South Korea >>> >>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>> >>> *M.Sc.* in Applied Meteorology, NUIST, China >>> >>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>> >>> *Founder and Director* of Research Society >>> >>> *Website*: https://researchsociety20.org/founder-and-director/ >>> >>> *E-mails*: founder-and-director at researchsociety20.org >>> >>> jalal at jejunu.ac.kr >>> >>> 20205103002 at nuist.edu.cn >>> >> >> >> -- >> Adam Phillips >> Associate Scientist IV, Climate Analysis Section >> Climate and Global Dynamics Laboratory >> National Center for Atmospheric Research >> www.cgd.ucar.edu/staff/asphilli/ >> >> >> > > > -- > > *Regards,* > > Dr. Jalal > > *Postdoctoral Researcher,* Typhoon Research Center, Jeju National > University, South Korea > > *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China > > *M.Sc.* in Applied Meteorology, NUIST, China > > *B.Sc.* in Disaster Management, PSTU, Bangladesh > > *Founder and Director* of Research Society > > *Website*: https://researchsociety20.org/founder-and-director/ > > *E-mails*: founder-and-director at researchsociety20.org > > jalal at jejunu.ac.kr > > 20205103002 at nuist.edu.cn > -- 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: -------------- next part -------------- A non-text attachment was scrubbed... Name: adam_sst_corr.v2.ncl Type: application/octet-stream Size: 4717 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: sst_SepNov_cor.png Type: image/png Size: 533865 bytes Desc: not available URL: From dmjalal90 at gmail.com Sun Jan 28 19:05:20 2024 From: dmjalal90 at gmail.com (Md. Jalal Uddin) Date: Mon, 29 Jan 2024 11:05:20 +0900 Subject: [ncl-talk] Do not show significance level in a correlation map In-Reply-To: References: Message-ID: Thank you so much, sir, for your kind efforts. The code did not work on my computer. Alternatively, I have shaded the correlation results where p-value is less than 0.10 (90% confidence level). new_r_sst = where(p.lt.0.10,r_sst,r_sst at _FillValue) copy_VarCoords(r_sst,new_r_sst) On Sat, 27 Jan 2024 at 08:52, Adam Phillips wrote: > Hi Jalal, > As you are running the exact version of NCL as I am, I do not have an > answer for you as to why the stippling is showing up for me but not for > you. (That should not happen.) If you replace the draw(plot) line with > draw(plot2) do you get a plot that shows the stippling? The only other idea > I have is to bypass using gsn_contour_shade, which will allow you to color > the stippling a color other than black. I have attached a new script along > with the resulting graphic. > Best, > Adam > > On Thu, Jan 25, 2024 at 7:39?PM Md. Jalal Uddin > wrote: > >> Thank you, sir. I am using NCL v6.6.2 in Ubuntu (Windows Subsystem for >> Linux). >> >> On Fri, 26 Jan 2024 at 01:38, Adam Phillips wrote: >> >>> Hi Jalal, >>> The error messages are happening as your SST array has missing data over >>> land, so there are grid points that are completely missing. In this case >>> the error messages can be ignored. With regards to the stippling: When I >>> run the exact script you sent, I do get stippling in numerous areas. I >>> altered the gsnShadeFillScaleF resource to .45 to make it more apparent, >>> and outlined the areas that are significant in the attached plot. >>> >>> If when you run the script adam_sst_corr.ncl you are not getting any >>> stippling in the areas I outlined, I am puzzled. I am running NCL v6.6.2 >>> which is the latest version of NCL. What version are you running? >>> Best, >>> Adam >>> >>> >>> On Wed, Jan 24, 2024 at 6:32?PM Md. Jalal Uddin >>> wrote: >>> >>>> Hi Sir, >>>> >>>> I have now used the same data and code that you provided. However, the >>>> output does not show the significant shaded area (see attached figure). And >>>> I get the following warnings: >>>> *warning*:dim_rmvmean_n: 21375 rightmost sections of the input array >>>> contained all missing values >>>> *warning*:escorc: Non-fatal conditions encountered in series or xstd >>>> equals zero. >>>> Possibly, all values of a series are constant. >>>> *warning*:escorc: Most likely, one or more series consisted of all >>>> constant values >>>> >>>> I am waiting for your kind advice. >>>> >>>> On Thu, 25 Jan 2024 at 02:13, Adam Phillips wrote: >>>> >>>>> Hi Jalal, >>>>> I am guessing the reason why my code didn't work with the NOAA/PSL >>>>> data is because that data's time coordinate variable did not match that of >>>>> the HadISST dataset. (In the script I sent you I used coordinate >>>>> subscripting of the time dimension to select the time needed.) Also, the >>>>> NOAA/PSL SST dataset you pointed me to only runs from 1981-2012, and your >>>>> script was set to analyze years 1980-2021. >>>>> >>>>> I will write to you offline with how to download the HadISST data file. >>>>> Best, >>>>> Adam >>>>> >>>>> >>>>> >>>>> On Wed, Jan 24, 2024 at 1:53?AM Md. Jalal Uddin >>>>> wrote: >>>>> >>>>>> Dear Sir, >>>>>> >>>>>> Could I have this dataset hadisst.187001-202212.nc? >>>>>> >>>>>> On Fri, 19 Jan 2024 at 19:42, Md. Jalal Uddin >>>>>> wrote: >>>>>> >>>>>>> Thank you, sir, for your kind efforts. Your code works for >>>>>>> hadisst.187001-202212.nc data. However, when I use my sst data from >>>>>>> NOAA (sst.mnmean.nc), the code does not work. I do not know the >>>>>>> reason for this. Sorry to bother you. >>>>>>> >>>>>>> NOAA sst data: >>>>>>> https://psl.noaa.gov/repository/entry/show?entryid=cac1c2a6-a864-4409-bb77-1fdead8eeb6e >>>>>>> >>>>>>> On Fri, 19 Jan 2024 at 15:32, Md. Jalal Uddin >>>>>>> wrote: >>>>>>> >>>>>>>> This is great, thank you so much! >>>>>>>> >>>>>>>> On Wed, 17 Jan 2024 at 05:09, Adam Phillips >>>>>>>> wrote: >>>>>>>> >>>>>>>>> Hi Jalal, >>>>>>>>> I was able to plug in a monthly SST dataset that I have here, so I >>>>>>>>> was able to run your script. So I think what is happening is that the >>>>>>>>> stippling is so small/faint/sparse that it is hard to see it. I made two >>>>>>>>> changes to your script: >>>>>>>>> 1 - I was getting stippling over land as the rtest function was >>>>>>>>> seemingly not recognizing the _FillValue of the input array. (I find this >>>>>>>>> strange, but have not used rtest before.) To fix this, I put in a where >>>>>>>>> statement and added the _FillValue attribute to your p array: >>>>>>>>> p = rtest(r_sst,42,0) >>>>>>>>> * p at _FillValue = 1.e20* >>>>>>>>> copy_VarCoords(r_sst, p) >>>>>>>>> * p = where(ismissing(r_sst),p at _FillValue,p)* >>>>>>>>> >>>>>>>>> 2 - I modified/added the following gsn_contour_shade options to >>>>>>>>> make the stippling more visible: >>>>>>>>> >>>>>>>>> * opt at gsnShadeFillScaleF = 0.25 ; add extra density >>>>>>>>> opt at gsnShadeFillDotSizeF = 0.0015* >>>>>>>>> >>>>>>>>> I have attached the modified script, along with the plot that I am >>>>>>>>> getting. >>>>>>>>> Hope that helps! >>>>>>>>> Best, >>>>>>>>> Adam >>>>>>>>> >>>>>>>>> >>>>>>>>> On Mon, Jan 15, 2024 at 8:58?PM Md. Jalal Uddin < >>>>>>>>> dmjalal90 at gmail.com> wrote: >>>>>>>>> >>>>>>>>>> Thank you very much, sir. I used the following resources: >>>>>>>>>> >>>>>>>>>> opt = True >>>>>>>>>> opt at gsnShadeFillType = "pattern" >>>>>>>>>> opt at gsnShadeLow = 17 ; stipple >>>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>>> >>>>>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt)\ >>>>>>>>>> >>>>>>>>>> overlay (plot, plot2) >>>>>>>>>> >>>>>>>>>> draw(plot) >>>>>>>>>> frame(wks) >>>>>>>>>> >>>>>>>>>> However, I did not get output as stipple areas where p value is >>>>>>>>>> less than 0.05. I added the full code including Nino data in the >>>>>>>>>> attachment. Should I send the SST data (sst.mon.mean.nc) via ftp? >>>>>>>>>> >>>>>>>>>> On Sat, 13 Jan 2024 at 06:55, Adam Phillips >>>>>>>>>> wrote: >>>>>>>>>> >>>>>>>>>>> Hi Jalal, >>>>>>>>>>> As you are pattern filling, and the default is color filling, >>>>>>>>>>> you need to set: >>>>>>>>>>> opt at gsnShadeFillType = "pattern" >>>>>>>>>>> I am also unclear what you mean by this statement: >>>>>>>>>>> "...to show the significance level as a shaded pattern where p >>>>>>>>>>> value is less than 0.05 and 0.01" >>>>>>>>>>> >>>>>>>>>>> Do you want to stipple areas less than .05? If so the following >>>>>>>>>>> code should do that: >>>>>>>>>>> opt = True >>>>>>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>>>> plot2 = gsn_contour_shade(plot2,0.05,999.,opt) >>>>>>>>>>> >>>>>>>>>>> If instead you want to stipple areas less than .01, you will >>>>>>>>>>> need to modify the contour levels in your res2 resource list to have a >>>>>>>>>>> contour at .01: >>>>>>>>>>> res2 at cnMinLevelValF = 0.00 ; set min contour level >>>>>>>>>>> res2 at cnMaxLevelValF = .10 ; set max contour level >>>>>>>>>>> res2 at cnLevelSpacingF = 0.01 ; set contour spacing >>>>>>>>>>> >>>>>>>>>>> opt = True >>>>>>>>>>> opt at gsnShade*Low* = 17 ; stipple >>>>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>>>> plot2 = gsn_contour_shade(plot2,0.01,999.,opt)\ >>>>>>>>>>> >>>>>>>>>>> Hope that answers your question. As always please respond to the >>>>>>>>>>> ncl-talk email list with further questions. >>>>>>>>>>> Best, >>>>>>>>>>> Adam >>>>>>>>>>> >>>>>>>>>>> On Fri, Jan 5, 2024 at 6:15?AM Md. Jalal Uddin via ncl-talk < >>>>>>>>>>> ncl-talk at mailman.ucar.edu> wrote: >>>>>>>>>>> >>>>>>>>>>>> Dear all, >>>>>>>>>>>> >>>>>>>>>>>> I have tried to use the gsn_contour_shade function to show the >>>>>>>>>>>> significance level as a shaded pattern where p value is less than 0.05 and >>>>>>>>>>>> 0.01. However, the output map does not show the shaded areas (please see >>>>>>>>>>>> the attached map). >>>>>>>>>>>> >>>>>>>>>>>> Here are some codes: >>>>>>>>>>>> r_sst = escorc(ano_nino,ano_sstSON(lat|:,lon|:,time|:)) >>>>>>>>>>>> ; ccr(lat,lon) >>>>>>>>>>>> copy_VarCoords(ano_sstSON(0,:,:), r_sst) >>>>>>>>>>>> r_sst at long_name = "Correlation: r_sst" >>>>>>>>>>>> >>>>>>>>>>>> p = rtest(r_sst,42,0) >>>>>>>>>>>> copy_VarCoords(r_sst, p) >>>>>>>>>>>> >>>>>>>>>>>> opt = True >>>>>>>>>>>> opt at gsnShadeHigh = 17 ; stipple >>>>>>>>>>>> opt at gsnShadeFillScaleF = 0.5 ; add extra density >>>>>>>>>>>> opt at gsnShadeFillDotSizeF = 0.002 ; make dots larger >>>>>>>>>>>> >>>>>>>>>>>> plot = gsn_csm_contour_map(wks,r_sst, res) ; create plot >>>>>>>>>>>> >>>>>>>>>>>> plot2 = gsn_csm_contour(wks,p, res2) >>>>>>>>>>>> plot2 = gsn_contour_shade(plot2,0.01,0.05,opt) >>>>>>>>>>>> >>>>>>>>>>>> overlay (plot, plot2) >>>>>>>>>>>> >>>>>>>>>>>> draw(plot) >>>>>>>>>>>> frame(wks) >>>>>>>>>>>> >>>>>>>>>>>> Please find the full code in the attached file. >>>>>>>>>>>> >>>>>>>>>>>> Any suggestions would be appreciated. >>>>>>>>>>>> >>>>>>>>>>>> -- >>>>>>>>>>>> >>>>>>>>>>>> *Regards,* >>>>>>>>>>>> >>>>>>>>>>>> Dr. Jalal >>>>>>>>>>>> >>>>>>>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju >>>>>>>>>>>> National University, South Korea >>>>>>>>>>>> >>>>>>>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, >>>>>>>>>>>> NUIST, China >>>>>>>>>>>> >>>>>>>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>>>>>>> >>>>>>>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>>>>>>> >>>>>>>>>>>> *Founder and Director* of Research Society >>>>>>>>>>>> >>>>>>>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>>>>>>> >>>>>>>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>>>>>>> >>>>>>>>>>>> jalal at jejunu.ac.kr >>>>>>>>>>>> >>>>>>>>>>>> 20205103002 at nuist.edu.cn >>>>>>>>>>>> _______________________________________________ >>>>>>>>>>>> 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/ >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> -- >>>>>>>>>> >>>>>>>>>> *Regards,* >>>>>>>>>> >>>>>>>>>> Dr. Jalal >>>>>>>>>> >>>>>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju >>>>>>>>>> National University, South Korea >>>>>>>>>> >>>>>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>>>>> China >>>>>>>>>> >>>>>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>>>>> >>>>>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>>>>> >>>>>>>>>> *Founder and Director* of Research Society >>>>>>>>>> >>>>>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>>>>> >>>>>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>>>>> >>>>>>>>>> jalal at jejunu.ac.kr >>>>>>>>>> >>>>>>>>>> 20205103002 at nuist.edu.cn >>>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>>> -- >>>>>>>>> Adam Phillips >>>>>>>>> Associate Scientist IV, Climate Analysis Section >>>>>>>>> Climate and Global Dynamics Laboratory >>>>>>>>> National Center for Atmospheric Research >>>>>>>>> www.cgd.ucar.edu/staff/asphilli/ >>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> -- >>>>>>>> >>>>>>>> *Regards,* >>>>>>>> >>>>>>>> Dr. Jalal >>>>>>>> >>>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>>>> University, South Korea >>>>>>>> >>>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>>> China >>>>>>>> >>>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>>> >>>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>>> >>>>>>>> *Founder and Director* of Research Society >>>>>>>> >>>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>>> >>>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>>> >>>>>>>> jalal at jejunu.ac.kr >>>>>>>> >>>>>>>> 20205103002 at nuist.edu.cn >>>>>>>> >>>>>>> >>>>>>> >>>>>>> -- >>>>>>> >>>>>>> *Regards,* >>>>>>> >>>>>>> Dr. Jalal >>>>>>> >>>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>>> University, South Korea >>>>>>> >>>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>>> China >>>>>>> >>>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>>> >>>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>>> >>>>>>> *Founder and Director* of Research Society >>>>>>> >>>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>>> >>>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>>> >>>>>>> jalal at jejunu.ac.kr >>>>>>> >>>>>>> 20205103002 at nuist.edu.cn >>>>>>> >>>>>> >>>>>> >>>>>> -- >>>>>> >>>>>> *Regards,* >>>>>> >>>>>> Dr. Jalal >>>>>> >>>>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>>>> University, South Korea >>>>>> >>>>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, >>>>>> China >>>>>> >>>>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>>>> >>>>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>>>> >>>>>> *Founder and Director* of Research Society >>>>>> >>>>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>>>> >>>>>> *E-mails*: founder-and-director at researchsociety20.org >>>>>> >>>>>> jalal at jejunu.ac.kr >>>>>> >>>>>> 20205103002 at nuist.edu.cn >>>>>> >>>>> >>>>> >>>>> -- >>>>> Adam Phillips >>>>> Associate Scientist IV, Climate Analysis Section >>>>> Climate and Global Dynamics Laboratory >>>>> National Center for Atmospheric Research >>>>> www.cgd.ucar.edu/staff/asphilli/ >>>>> >>>>> >>>>> >>>> >>>> >>>> -- >>>> >>>> *Regards,* >>>> >>>> Dr. Jalal >>>> >>>> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >>>> University, South Korea >>>> >>>> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >>>> >>>> *M.Sc.* in Applied Meteorology, NUIST, China >>>> >>>> *B.Sc.* in Disaster Management, PSTU, Bangladesh >>>> >>>> *Founder and Director* of Research Society >>>> >>>> *Website*: https://researchsociety20.org/founder-and-director/ >>>> >>>> *E-mails*: founder-and-director at researchsociety20.org >>>> >>>> jalal at jejunu.ac.kr >>>> >>>> 20205103002 at nuist.edu.cn >>>> >>> >>> >>> -- >>> Adam Phillips >>> Associate Scientist IV, Climate Analysis Section >>> Climate and Global Dynamics Laboratory >>> National Center for Atmospheric Research >>> www.cgd.ucar.edu/staff/asphilli/ >>> >>> >>> >> >> >> -- >> >> *Regards,* >> >> Dr. Jalal >> >> *Postdoctoral Researcher,* Typhoon Research Center, Jeju National >> University, South Korea >> >> *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China >> >> *M.Sc.* in Applied Meteorology, NUIST, China >> >> *B.Sc.* in Disaster Management, PSTU, Bangladesh >> >> *Founder and Director* of Research Society >> >> *Website*: https://researchsociety20.org/founder-and-director/ >> >> *E-mails*: founder-and-director at researchsociety20.org >> >> jalal at jejunu.ac.kr >> >> 20205103002 at nuist.edu.cn >> > > > -- > Adam Phillips > Associate Scientist IV, Climate Analysis Section > Climate and Global Dynamics Laboratory > National Center for Atmospheric Research > www.cgd.ucar.edu/staff/asphilli/ > > > -- *Regards,* Dr. Jalal *Postdoctoral Researcher,* Typhoon Research Center, Jeju National University, South Korea *PhD* in Atmospheric Physics and Atmospheric Environment, NUIST, China *M.Sc.* in Applied Meteorology, NUIST, China *B.Sc.* in Disaster Management, PSTU, Bangladesh *Founder and Director* of Research Society *Website*: https://researchsociety20.org/founder-and-director/ *E-mails*: founder-and-director at researchsociety20.org jalal at jejunu.ac.kr 20205103002 at nuist.edu.cn -------------- next part -------------- An HTML attachment was scrubbed... URL: