From xi.chang01 at gmail.com Tue May 9 18:22:29 2023 From: xi.chang01 at gmail.com (Xi Chang) Date: Tue, 9 May 2023 17:22:29 -0700 Subject: [ncl-talk] sub drveof: ier,jopt= 10 1 returned from vcmssm/crmssm Message-ID: Hi everyone, I cannot figure this out, I am trying todo EOF analysis by stacking the data in latitude axis. Everything looks fine to me but I still get this error (see below). There is no missing value in the data and I am using NCL 6.6.2. This is the error: === Variable: x Type: float Total Size: 253281600 bytes 63320400 values Number of Dimensions: 3 Dimensions and sizes: [time | 5863] x [lat | 75] x [lon | 144] Coordinates: Number Of Attributes: 1 _FillValue : -9.96921e+36 (0) Calculating EOF.. sub drveof: ier,jopt= 10 1 returned from vcmssm/crmssm warning:eofunc_n: 10 eigenvectors failed to converge Variable: eof Type: float Total Size: 86400 bytes 21600 values Number of Dimensions: 3 Dimensions and sizes: [evn | 2] x [lat | 75] x [lon | 144] Coordinates: evn: [1..2] Number Of Attributes: 5 eval : ( -9.96921e+36, -9.96921e+36 ) pcvar : ( -9.96921e+36, -9.96921e+36 ) matrix : correlation method : no transpose _FillValue : -9.96921e+36 (0) (0) min=-9.96921e+36 max=-9.96921e+36 -------------- next part -------------- An HTML attachment was scrubbed... URL: From M.Lazenby at sussex.ac.uk Wed May 10 09:52:07 2023 From: M.Lazenby at sussex.ac.uk (Melissa Lazenby) Date: Wed, 10 May 2023 15:52:07 +0000 Subject: [ncl-talk] Bar Chart Plotting Issues Message-ID: Dear NCL Users I am trying to plot a bar chart with year and months on the X-axis and am struggling to do so correctly. I would also like to choose some years and months to be grey instead of black to show certainty levels. I am trying to make a similar plot to the one below: [cid:0ba8b8f1-2012-4b28-8d87-b13e4017b864] This is what I currently have: [cid:292cb285-6cd1-464e-b5d0-933539d55edb] Any help would be greatly appreciated! Kindest Regards Melissa Code: ;************************************************ load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" ;************************************************ begin ;************************************************ ; plotting ;************************************************ wks = gsn_open_wks("X11","Mpwapwa_Index") res = True res at gsnMaximize = True res at tiXAxisString = "Year & Month" ; x-axis label res at tiYAxisString = "Precipitation Index" ; y-axis label res at gsnXYBarChart = True res at gsnXYBarChartBarWidth = 0.1 res at gsnYRefLine = 0 ; reference line res at tiMainString = "Mpwapwa" res at trYMinF = -3.5 ; min level res at trYMaxF = 3.5 ; max level res at vpWidthF = 1.5 ; Make the plot wider than it res at vpHeightF = 0.75 ; is high. res at gsnDraw = False ; do not draw the plot res at gsnFrame = False ;Years and Months x = (/185501,?????185502,?????185503,?????185504,?????185505,?????185506,?????185507,?????185508,?????185509,?????185510,?????185511,?????185512,?????185601,?????185602,?????185603,?????185604,?????185605,?????185606,?????185607,?????185608,?????185609,?????185610,?????185611,?????185612,?????185701,?????185702,?????185703,?????185704,?????185705,?????185706,?????185707,?????185708,?????185709,?????185710,?????185711,?????185712,?????185801,?????185802,?????185803,?????185804,?????185805,?????185806,?????185807,?????185808,?????185809,?????185810,?????185811,?????185812,?????185901,?????185902,?????185903,?????185904,?????185905,?????185906,?????185907,?????185908,?????185909,?????185910,?????185911,?????185912,?????186001,?????186002,?????186003,?????186004,?????186005,?????186006,?????186007,?????186008,?????186009,?????186010,?????186011,?????186012,?????186101,?????186102,?????186103,?????186104,?????186105,?????186106,?????186107,?????186108,?????186109,?????186110,?????186111,?????186112,?????186201,?????186202,?????186203,?????186204,?????186205,?????186206,?????186207,?????186208,?????186209,?????186210,?????186211,?????186212,?????186301,?????186302,?????186303,?????186304,?????186305,?????186306,?????186307,?????186308,?????186309,?????186310,?????186311,?????186312,?????186401,?????186402,?????186403,?????186404,?????186405,?????186406,?????186407,?????186408,?????186409,?????186410,?????186411,?????186412,?????186501,?????186502,?????186503,?????186504,?????186505,?????186506,?????186507,?????186508,?????186509,?????186510,?????186511,?????186512,?????186601,?????186602,?????186603,?????186604,?????186605,?????186606,?????186607,?????186608,?????186609,?????186610,?????186611,?????186612,?????186701,?????186702,?????186703,?????186704,?????186705,?????186706,?????186707,?????186708,?????186709,?????186710,?????186711,?????186712,?????186801,?????186802,?????186803,?????186804,?????186805,?????186806,?????186807,?????186808,?????186809,?????186810,?????186811,?????186812,?????186901,?????186902,?????186903,?????186904,?????186905,?????186906,?????186907,?????186908,?????186909,?????186910,?????186911,?????186912,?????187001,?????187002,?????187003,?????187004,?????187005,?????187006,?????187007,?????187008,?????187009,?????187010,?????187011,?????187012,?????187101,?????187102,?????187103,?????187104,?????187105,?????187106,?????187107,?????187108,?????187109,?????187110,?????187111,?????187112,?????187201,?????187202,?????187203,?????187204,?????187205,?????187206,?????187207,?????187208,?????187209,?????187210,?????187211,?????187212,?????187301,?????187302,?????187303,?????187304,?????187305,?????187306,?????187307,?????187308,?????187309,?????187310,?????187311,?????187312,?????187401,?????187402,?????187403,?????187404,?????187405,?????187406,?????187407,?????187408,?????187409,?????187410,?????187411,?????187412,?????187501,?????187502,?????187503,?????187504,?????187505,?????187506,?????187507,?????187508,?????187509,?????187510,?????187511,?????187512,?????187601,?????187602,?????187603,?????187604,?????187605,?????187606,?????187607,?????187608,?????187609,?????187610,?????187611,?????187612,?????187701,?????187702,?????187703,?????187704,?????187705,?????187706,?????187707,?????187708,?????187709,?????187710,?????187711,?????187712,?????187801,?????187802,?????187803,?????187804,?????187805,?????187806,?????187807,?????187808,?????187809,?????187810,?????187811,?????187812,?????187901,?????187902,?????187903,?????187904,?????187905,?????187906,?????187907,?????187908,?????187909,?????187910,?????187911,?????187912,?????188001,?????188002,?????188003,?????188004,?????188005,?????188006,?????188007,?????188008,?????188009,?????188010,?????188011,?????188012,?????188101,?????188102,?????188103,?????188104,?????188105,?????188106,?????188107,?????188108,?????188109,?????188110,?????188111,?????188112,?????188201,?????188202,?????188203,?????188204,?????188205,?????188206,?????188207,?????188208,?????188209,?????188210,?????188211,?????188212,?????188301,?????188302,?????188303,?????188304,?????188305,?????188306,?????188307,?????188308,?????188309,?????188310,?????188311,?????188312,?????188401,?????188402,?????188403,?????188404,?????188405,?????188406,?????188407,?????188408,?????188409,?????188410,?????188411,?????188412,?????188501,?????188502,?????188503,?????188504,?????188505,?????188506,?????188507,?????188508,?????188509,?????188510,?????188511,?????188512,?????188601,?????188602,?????188603,?????188604,?????188605,?????188606,?????188607,?????188608,?????188609,?????188610,?????188611,?????188612,?????188701,?????188702,?????188703,?????188704,?????188705,?????188706,?????188707,?????188708,?????188709,?????188710,?????188711,?????188712,?????188801,?????188802,?????188803,?????188804,?????188805,?????188806,?????188807,?????188808,?????188809,?????188810,?????188811,?????188812,?????188901,?????188902,?????188903,?????188904,?????188905,?????188906,?????188907,?????188908,?????188909,?????188910,?????188911,?????188912,?????189001,?????189002,?????189003,?????189004,?????189005,?????189006,?????189007,?????189008,?????189009,?????189010,?????189011,?????189012/) ;Index y = (/ 0,???0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????3,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????-1,???0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????-3,???-3,???-3,???-3,???0,????0,????0,????0,????0,????0,????0,????0,????2,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0.05,?0,????0,????0,????0,????0,????0,????2,????1,????1,????1,????1,????1,????2,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????-1,???1,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????-1,???-1,???-1,???-1,???-1,???0,????0,????0,????0,????0,????0,????0,????-2,???-2,???-2,???-2,???-2,???0,????0,????0,????0,????0,????0,????1,????1,????1,????1,????1,????1,????0.05,?0,????0,????0,????0,????1,????0.05,?0.05,?0.05,?0.05,?0.05,?2,????0,????0,????0,????0,????0,????0,????0,????-1,???-2,???-3,???-2,???-2,???0,????0,????0,????0,????0,????0,????0,????-2,???0,????-1,???0,????0,????0,????0,????0,????0,????0,????0,????2,????1,????0.05,?1,????0.05,?0.05,?0.05,?0,????0,????0,????0,????0,????0,????-1,???-1,???0.05,?0.05,?0,????0,????0,????0,????0,????0,????0,????0,????-3,???0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0.05,?0.05,?0.05,?0.05,?0.05,?0,????0,????0,????0,????0,????0,????0,????0.05,?0.05,?0.05,?0.05,?0.05,?0,????0,????0,????0,????0,????0,????0.05,?0.05,?0.05,?0.05,?0.05,?0.05,?0,????0,????0,????0,????0,????0,????0,????-2,???-3,???-3,???-3,???-2,???0,????0,????0,????0,????0,????0.05,?0,????0.05,?0.05,?0.05,?0.05,?0.05,?0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0,????0/) res at tmXBMode = "Manual"????? res at tmXBTickStartF = 185501 res at tmXBTickEndF = 189012 ;res at tmXBTickSpacingF= 120 res at tmXBLabelStride = 1 bc = new(432,string) bc= "black" time = (/185501,185502,185503,??????185504,?????185505,?????185506,?????185507,?????185508,?????185509,?????185510,?????185511,?????185512,?????185601,?????185602,?????185603,?????185604,?????185605,?????185606,?????185607,?????185608,?????185609,?????185610,?????185611,?????185612,?????185701,?????185702,?????185703,?????185704,?????185705,?????185706,?????185707,?????185708,?????185709,?????185710,?????185711,?????185712,?????185801,?????185802,?????185803,?????185804,?????185805,?????185806,?????185807,?????185808,?????185809,?????185810,?????185811,?????185812,?????185901,?????185902,?????185903,?????185904,?????185905,?????185906,?????185907,?????185908,?????185909,?????185910,?????185911,?????185912,?????186001,?????186002,?????186003,?????186004,?????186005,?????186006,?????186007,?????186008,?????186009,?????186010,?????186011,?????186012,?????186101,?????186102,?????186103,?????186104,?????186105,?????186106,?????186107,?????186108,?????186109,?????186110,?????186111,?????186112,?????186201,?????186202,?????186203,?????186204,?????186205,?????186206,?????186207,?????186208,?????186209,?????186210,?????186211,?????186212,?????186301,?????186302,?????186303,?????186304,?????186305,?????186306,?????186307,?????186308,?????186309,?????186310,?????186311,?????186312,?????186401,?????186402,?????186403,?????186404,?????186405,?????186406,?????186407,?????186408,?????186409,?????186410,?????186411,?????186412,?????186501,?????186502,?????186503,?????186504,?????186505,?????186506,?????186507,?????186508,?????186509,?????186510,?????186511,?????186512,?????186601,?????186602,?????186603,?????186604,?????186605,?????186606,?????186607,?????186608,?????186609,?????186610,?????186611,?????186612,?????186701,?????186702,?????186703,?????186704,?????186705,?????186706,?????186707,?????186708,?????186709,?????186710,?????186711,?????186712,?????186801,?????186802,?????186803,?????186804,?????186805,?????186806,?????186807,?????186808,?????186809,?????186810,?????186811,?????186812,?????186901,?????186902,?????186903,?????186904,?????186905,?????186906,?????186907,?????186908,?????186909,?????186910,?????186911,?????186912,?????187001,?????187002,?????187003,?????187004,?????187005,?????187006,?????187007,?????187008,?????187009,?????187010,?????187011,?????187012,?????187101,?????187102,?????187103,?????187104,?????187105,?????187106,?????187107,?????187108,?????187109,?????187110,?????187111,?????187112,?????187201,?????187202,?????187203,?????187204,?????187205,?????187206,?????187207,?????187208,?????187209,?????187210,?????187211,?????187212,?????187301,?????187302,?????187303,?????187304,?????187305,?????187306,?????187307,?????187308,?????187309,?????187310,?????187311,?????187312,?????187401,?????187402,?????187403,?????187404,?????187405,?????187406,?????187407,?????187408,?????187409,?????187410,?????187411,?????187412,?????187501,?????187502,?????187503,?????187504,?????187505,?????187506,?????187507,?????187508,?????187509,?????187510,?????187511,?????187512,?????187601,?????187602,?????187603,?????187604,?????187605,?????187606,?????187607,?????187608,?????187609,?????187610,?????187611,?????187612,?????187701,?????187702,?????187703,?????187704,?????187705,?????187706,?????187707,?????187708,?????187709,?????187710,?????187711,?????187712,?????187801,?????187802,?????187803,?????187804,?????187805,?????187806,?????187807,?????187808,?????187809,?????187810,?????187811,?????187812,?????187901,?????187902,?????187903,?????187904,?????187905,?????187906,?????187907,?????187908,?????187909,?????187910,?????187911,?????187912,?????188001,?????188002,?????188003,?????188004,?????188005,?????188006,?????188007,?????188008,?????188009,?????188010,?????188011,?????188012,?????188101,?????188102,?????188103,?????188104,?????188105,?????188106,?????188107,?????188108,?????188109,?????188110,?????188111,?????188112,?????188201,?????188202,?????188203,?????188204,?????188205,?????188206,?????188207,?????188208,?????188209,?????188210,?????188211,?????188212,?????188301,?????188302,?????188303,?????188304,?????188305,?????188306,?????188307,?????188308,?????188309,?????188310,?????188311,?????188312,?????188401,?????188402,?????188403,?????188404,?????188405,?????188406,?????188407,?????188408,?????188409,?????188410,?????188411,?????188412,?????188501,?????188502,?????188503,?????188504,?????188505,?????188506,?????188507,?????188508,?????188509,?????188510,?????188511,?????188512,?????188601,?????188602,?????188603,?????188604,?????188605,?????188606,?????188607,?????188608,?????188609,?????188610,?????188611,?????188612,?????188701,?????188702,?????188703,?????188704,?????188705,?????188706,?????188707,?????188708,?????188709,?????188710,?????188711,?????188712,?????188801,?????188802,?????188803,?????188804,?????188805,?????188806,?????188807,?????188808,?????188809,?????188810,?????188811,?????188812,?????188901,?????188902,?????188903,?????188904,?????188905,?????188906,?????188907,?????188908,?????188909,?????188910,?????188911,?????188912,?????189001,?????189002,?????189003,?????189004,?????189005,?????189006,?????189007,?????189008,?????189009,?????189010,?????189011,?????189012/) selyr = (/185709,185812,185912,186001, 186002,18603,186012, 187105,187211,187412,187506,187601,188001,188102,188204/) do gg = 0,dimsizes(selyr)-1 bc = where(time.eq.selyr(gg),"grey",bc) end do res at gsnXYBarChartColors = "bc" plot=gsn_csm_xy(wks,x,y,res) draw(plot) frame(wks) end Dr. Melissa Lazenby Lecturer in Climate Change Department of Geography Chichester 1 C150 University of Sussex "Education is the most powerful weapon which you can use to change the world" Nelson Mandela [cid:72ad72ef-975a-4887-a128-63e5bd225f2b] -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image.png Type: image/png Size: 64237 bytes Desc: image.png URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image.png Type: image/png Size: 24966 bytes Desc: image.png URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Outlook-gkrj01ui.png Type: image/png Size: 24937 bytes Desc: Outlook-gkrj01ui.png URL: From francesco.trotta4 at unibo.it Thu May 11 10:22:55 2023 From: francesco.trotta4 at unibo.it (Francesco Trotta) Date: Thu, 11 May 2023 16:22:55 +0000 Subject: [ncl-talk] install ncl with conda Message-ID: Dear Ncl-User I'm trying to install ncl in my Linux VM inside my new macbook M2 but I got these message below How can I solve it? best Francesco conda create -n ncl_stable -c conda-forge ncl Collecting package metadata (current_repodata.json): done Solving environment: failed with repodata from current_repodata.json, will retry with next repodata source. Collecting package metadata (repodata.json): done Solving environment: failed PackagesNotFoundError: The following packages are not available from current channels: - ncl Current channels: - https://conda.anaconda.org/conda-forge/linux-aarch64 - https://conda.anaconda.org/conda-forge/noarch - https://repo.anaconda.com/pkgs/main/linux-aarch64 - https://repo.anaconda.com/pkgs/main/noarch - https://repo.anaconda.com/pkgs/r/linux-aarch64 - https://repo.anaconda.com/pkgs/r/noarch To search for alternate channels that may provide the conda package you're looking for, navigate to https://anaconda.org and use the search bar at the top of the page. -------------- next part -------------- An HTML attachment was scrubbed... URL: From asphilli at ucar.edu Thu May 11 17:38:32 2023 From: asphilli at ucar.edu (Adam Phillips) Date: Thu, 11 May 2023 17:38:32 -0600 Subject: [ncl-talk] Bar Chart Plotting Issues In-Reply-To: References: Message-ID: Hi Melissa, I think there are two things going on. 1- Your input x-axis values (= your x array) go up by 1 for every month of the year, and then jump by 89. This will result in large blank areas in the plot as NCL will allocate space for (say) 185513-185600. The solution is to create a new x-axis array after you've defined x: x_alt = ispan(0,dimsizes(x)-1,1) manually set the x-axis tick marks and labels (to show say every 5th year): res at tmXBLabelMode = "ExplicitLabels" res at tmXBValues = x_alt(0::60). ; create x-axis tick mark every 5 years res at tmXBLabels = x(0::60) ; and label them appropriately res at tmXBMinorValues = x_alt(0::6) create minor tick marks every 6 months; optional and then pass in your x_alt array as the x-axis values to plot against your y array: plot=gsn_csm_xy(wks,x_alt,y,res) 2 - I would reset your bar width to be set to 1, or, if you don't want successive month bars touching then set it to say 0.8: res at gsnXYBarChartBarWidth = 0.8 If you have any further questions on this please respond to ncl-talk. Best, Adam On Wed, May 10, 2023 at 9:52?AM Melissa Lazenby via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear NCL Users > > I am trying to plot a bar chart with year and months on the X-axis and am > struggling to do so correctly. > I would also like to choose some years and months to be grey instead of > black to show certainty levels. I am trying to make a similar plot to the > one below: > > > This is what I currently have: > > > Any help would be greatly appreciated! > > Kindest Regards > Melissa > > > Code: > > ;************************************************ > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" > ;************************************************ > begin > > ;************************************************ > ; plotting > ;************************************************ > wks = gsn_open_wks("X11","Mpwapwa_Index") > > res = True > res at gsnMaximize = True > res at tiXAxisString = "Year & Month" ; x-axis label > res at tiYAxisString = "Precipitation Index" ; y-axis label > res at gsnXYBarChart = True > res at gsnXYBarChartBarWidth = 0.1 > res at gsnYRefLine = 0 ; reference line > res at tiMainString = "Mpwapwa" > > res at trYMinF = -3.5 ; min level > res at trYMaxF = 3.5 ; max level > > res at vpWidthF = 1.5 ; Make the plot wider than it > res at vpHeightF = 0.75 ; is high. > > res at gsnDraw = False ; do not draw the plot > res at gsnFrame = False > > ;Years and Months > x = (/185501, 185502, 185503, 185504, 185505, 185506, 185507, 185508, > 185509, 185510, 185511, 185512, 185601, 185602, 185603, 185604, 185605, > 185606, 185607, 185608, 185609, 185610, 185611, 185612, 185701, 185702, > 185703, 185704, 185705, 185706, 185707, 185708, 185709, 185710, 185711, > 185712, 185801, 185802, 185803, 185804, 185805, 185806, 185807, 185808, > 185809, 185810, 185811, 185812, 185901, 185902, 185903, 185904, 185905, > 185906, 185907, 185908, 185909, 185910, 185911, 185912, 186001, 186002, > 186003, 186004, 186005, 186006, 186007, 186008, 186009, 186010, 186011, > 186012, 186101, 186102, 186103, 186104, 186105, 186106, 186107, 186108, > 186109, 186110, 186111, 186112, 186201, 186202, 186203, 186204, 186205, > 186206, 186207, 186208, 186209, 186210, 186211, 186212, 186301, 186302, > 186303, 186304, 186305, 186306, 186307, 186308, 186309, 186310, 186311, > 186312, 186401, 186402, 186403, 186404, 186405, 186406, 186407, 186408, > 186409, 186410, 186411, 186412, 186501, 186502, 186503, 186504, 186505, > 186506, 186507, 186508, 186509, 186510, 186511, 186512, 186601, 186602, > 186603, 186604, 186605, 186606, 186607, 186608, 186609, 186610, 186611, > 186612, 186701, 186702, 186703, 186704, 186705, 186706, 186707, 186708, > 186709, 186710, 186711, 186712, 186801, 186802, 186803, 186804, 186805, > 186806, 186807, 186808, 186809, 186810, 186811, 186812, 186901, 186902, > 186903, 186904, 186905, 186906, 186907, 186908, 186909, 186910, 186911, > 186912, 187001, 187002, 187003, 187004, 187005, 187006, 187007, 187008, > 187009, 187010, 187011, 187012, 187101, 187102, 187103, 187104, 187105, > 187106, 187107, 187108, 187109, 187110, 187111, 187112, 187201, 187202, > 187203, 187204, 187205, 187206, 187207, 187208, 187209, 187210, 187211, > 187212, 187301, 187302, 187303, 187304, 187305, 187306, 187307, 187308, > 187309, 187310, 187311, 187312, 187401, 187402, 187403, 187404, 187405, > 187406, 187407, 187408, 187409, 187410, 187411, 187412, 187501, 187502, > 187503, 187504, 187505, 187506, 187507, 187508, 187509, 187510, 187511, > 187512, 187601, 187602, 187603, 187604, 187605, 187606, 187607, 187608, > 187609, 187610, 187611, 187612, 187701, 187702, 187703, 187704, 187705, > 187706, 187707, 187708, 187709, 187710, 187711, 187712, 187801, 187802, > 187803, 187804, 187805, 187806, 187807, 187808, 187809, 187810, 187811, > 187812, 187901, 187902, 187903, 187904, 187905, 187906, 187907, 187908, > 187909, 187910, 187911, 187912, 188001, 188002, 188003, 188004, 188005, > 188006, 188007, 188008, 188009, 188010, 188011, 188012, 188101, 188102, > 188103, 188104, 188105, 188106, 188107, 188108, 188109, 188110, 188111, > 188112, 188201, 188202, 188203, 188204, 188205, 188206, 188207, 188208, > 188209, 188210, 188211, 188212, 188301, 188302, 188303, 188304, 188305, > 188306, 188307, 188308, 188309, 188310, 188311, 188312, 188401, 188402, > 188403, 188404, 188405, 188406, 188407, 188408, 188409, 188410, 188411, > 188412, 188501, 188502, 188503, 188504, 188505, 188506, 188507, 188508, > 188509, 188510, 188511, 188512, 188601, 188602, 188603, 188604, 188605, > 188606, 188607, 188608, 188609, 188610, 188611, 188612, 188701, 188702, > 188703, 188704, 188705, 188706, 188707, 188708, 188709, 188710, 188711, > 188712, 188801, 188802, 188803, 188804, 188805, 188806, 188807, 188808, > 188809, 188810, 188811, 188812, 188901, 188902, 188903, 188904, 188905, > 188906, 188907, 188908, 188909, 188910, 188911, 188912, 189001, 189002, > 189003, 189004, 189005, 189006, 189007, 189008, 189009, 189010, 189011, > 189012/) > > ;Index > y = (/ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, > 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, > -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -3, -3, -3, -3, 0, 0, 0, 0, 0, 0, 0, > 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, > 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, > 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, > 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, > 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, > 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.05, 0, 0, 0, 0, 0, 0, 2, 1, 1, 1, 1, > 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, > 0, 0, 0, 0, 0, 0, 0, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0, 0, -2, -2, -2, > -2, -2, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0.05, 0, 0, 0, 0, 1, 0.05, > 0.05, 0.05, 0.05, 0.05, 2, 0, 0, 0, 0, 0, 0, 0, -1, -2, -3, -2, -2, 0, 0, > 0, 0, 0, 0, 0, -2, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 2, 1, 0.05, 1, 0.05, > 0.05, 0.05, 0, 0, 0, 0, 0, 0, -1, -1, 0.05, 0.05, 0, 0, 0, 0, 0, 0, 0, 0, > -3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.05, 0.05, 0.05, 0.05, 0.05, 0, 0, 0, > 0, 0, 0, 0, 0.05, 0.05, 0.05, 0.05, 0.05, 0, 0, 0, 0, 0, 0, 0.05, 0.05, > 0.05, 0.05, 0.05, 0.05, 0, 0, 0, 0, 0, 0, 0, -2, -3, -3, -3, -2, 0, 0, 0, > 0, 0, 0.05, 0, 0.05, 0.05, 0.05, 0.05, 0.05, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, > 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/) > > res at tmXBMode = "Manual" > res at tmXBTickStartF = 185501 > res at tmXBTickEndF = 189012 > ;res at tmXBTickSpacingF= 120 > res at tmXBLabelStride = 1 > > > bc = new(432,string) > bc= "black" > time = (/185501,185502,185503, 185504, 185505, 185506, 185507, 185508, > 185509, 185510, 185511, 185512, 185601, 185602, 185603, 185604, 185605, > 185606, 185607, 185608, 185609, 185610, 185611, 185612, 185701, 185702, > 185703, 185704, 185705, 185706, 185707, 185708, 185709, 185710, 185711, > 185712, 185801, 185802, 185803, 185804, 185805, 185806, 185807, 185808, > 185809, 185810, 185811, 185812, 185901, 185902, 185903, 185904, 185905, > 185906, 185907, 185908, 185909, 185910, 185911, 185912, 186001, 186002, > 186003, 186004, 186005, 186006, 186007, 186008, 186009, 186010, 186011, > 186012, 186101, 186102, 186103, 186104, 186105, 186106, 186107, 186108, > 186109, 186110, 186111, 186112, 186201, 186202, 186203, 186204, 186205, > 186206, 186207, 186208, 186209, 186210, 186211, 186212, 186301, 186302, > 186303, 186304, 186305, 186306, 186307, 186308, 186309, 186310, 186311, > 186312, 186401, 186402, 186403, 186404, 186405, 186406, 186407, 186408, > 186409, 186410, 186411, 186412, 186501, 186502, 186503, 186504, 186505, > 186506, 186507, 186508, 186509, 186510, 186511, 186512, 186601, 186602, > 186603, 186604, 186605, 186606, 186607, 186608, 186609, 186610, 186611, > 186612, 186701, 186702, 186703, 186704, 186705, 186706, 186707, 186708, > 186709, 186710, 186711, 186712, 186801, 186802, 186803, 186804, 186805, > 186806, 186807, 186808, 186809, 186810, 186811, 186812, 186901, 186902, > 186903, 186904, 186905, 186906, 186907, 186908, 186909, 186910, 186911, > 186912, 187001, 187002, 187003, 187004, 187005, 187006, 187007, 187008, > 187009, 187010, 187011, 187012, 187101, 187102, 187103, 187104, 187105, > 187106, 187107, 187108, 187109, 187110, 187111, 187112, 187201, 187202, > 187203, 187204, 187205, 187206, 187207, 187208, 187209, 187210, 187211, > 187212, 187301, 187302, 187303, 187304, 187305, 187306, 187307, 187308, > 187309, 187310, 187311, 187312, 187401, 187402, 187403, 187404, 187405, > 187406, 187407, 187408, 187409, 187410, 187411, 187412, 187501, 187502, > 187503, 187504, 187505, 187506, 187507, 187508, 187509, 187510, 187511, > 187512, 187601, 187602, 187603, 187604, 187605, 187606, 187607, 187608, > 187609, 187610, 187611, 187612, 187701, 187702, 187703, 187704, 187705, > 187706, 187707, 187708, 187709, 187710, 187711, 187712, 187801, 187802, > 187803, 187804, 187805, 187806, 187807, 187808, 187809, 187810, 187811, > 187812, 187901, 187902, 187903, 187904, 187905, 187906, 187907, 187908, > 187909, 187910, 187911, 187912, 188001, 188002, 188003, 188004, 188005, > 188006, 188007, 188008, 188009, 188010, 188011, 188012, 188101, 188102, > 188103, 188104, 188105, 188106, 188107, 188108, 188109, 188110, 188111, > 188112, 188201, 188202, 188203, 188204, 188205, 188206, 188207, 188208, > 188209, 188210, 188211, 188212, 188301, 188302, 188303, 188304, 188305, > 188306, 188307, 188308, 188309, 188310, 188311, 188312, 188401, 188402, > 188403, 188404, 188405, 188406, 188407, 188408, 188409, 188410, 188411, > 188412, 188501, 188502, 188503, 188504, 188505, 188506, 188507, 188508, > 188509, 188510, 188511, 188512, 188601, 188602, 188603, 188604, 188605, > 188606, 188607, 188608, 188609, 188610, 188611, 188612, 188701, 188702, > 188703, 188704, 188705, 188706, 188707, 188708, 188709, 188710, 188711, > 188712, 188801, 188802, 188803, 188804, 188805, 188806, 188807, 188808, > 188809, 188810, 188811, 188812, 188901, 188902, 188903, 188904, 188905, > 188906, 188907, 188908, 188909, 188910, 188911, 188912, 189001, 189002, > 189003, 189004, 189005, 189006, 189007, 189008, 189009, 189010, 189011, > 189012/) > > selyr = (/185709,185812,185912,186001, 186002,18603,186012, > 187105,187211,187412,187506,187601,188001,188102,188204/) > do gg = 0,dimsizes(selyr)-1 > bc = where(time.eq.selyr(gg),"grey",bc) > end do > > res at gsnXYBarChartColors = "bc" > > > plot=gsn_csm_xy(wks,x,y,res) > > draw(plot) > frame(wks) > > end > > > > > > > > Dr. Melissa Lazenby > Lecturer in Climate Change > Department of Geography > Chichester 1 C150 > University of Sussex > > *"Education is the most powerful weapon which you can use to change the > world" Nelson Mandela* > > > _______________________________________________ > 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: -------------- next part -------------- A non-text attachment was scrubbed... Name: image.png Type: image/png Size: 64237 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image.png Type: image/png Size: 24966 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Outlook-gkrj01ui.png Type: image/png Size: 24937 bytes Desc: not available URL: From czarzycki at psu.edu Fri May 12 09:52:23 2023 From: czarzycki at psu.edu (Colin Zarzycki) Date: Fri, 12 May 2023 11:52:23 -0400 Subject: [ncl-talk] install ncl with conda In-Reply-To: References: Message-ID: Hi Francesco, See this discussion on the NCL github repository: https://github.com/NCAR/ncl/issues/187 Essentially, there is no "bottled" NCL with conda and the arm64 architecture, so it currently must be built from source using something like Macports. -Colin On Thu, May 11, 2023 at 12:23?PM Francesco Trotta via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear Ncl-User > > I'm trying to install ncl in my Linux VM inside my new macbook M2 > but I got these message below > > How can I solve it? > > best > Francesco > > conda create -n ncl_stable -c conda-forge ncl > Collecting package metadata (current_repodata.json): done > Solving environment: failed with repodata from current_repodata.json, will > retry with next repodata source. > Collecting package metadata (repodata.json): done > Solving environment: failed > > PackagesNotFoundError: The following packages are not available from > current channels: > > - ncl > > Current channels: > > - https://conda.anaconda.org/conda-forge/linux-aarch64 > - https://conda.anaconda.org/conda-forge/noarch > - https://repo.anaconda.com/pkgs/main/linux-aarch64 > - https://repo.anaconda.com/pkgs/main/noarch > - https://repo.anaconda.com/pkgs/r/linux-aarch64 > - https://repo.anaconda.com/pkgs/r/noarch > > To search for alternate channels that may provide the conda package you're > looking for, navigate to > > https://anaconda.org > > and use the search bar at the top of the page. > > _______________________________________________ > ncl-talk mailing list > ncl-talk at mailman.ucar.edu > List instructions, subscriber options, unsubscribe: > https://mailman.ucar.edu/mailman/listinfo/ncl-talk > -- Colin M. Zarzycki, Ph.D. Assistant Professor of Meteorology and Climate Dynamics Department of Meteorology and Atmospheric Science Pennsylvania State University Web: http://www.colinzarzycki.com My working day may not be your working day. Please do not feel obliged to reply to this email outside of your normal working hours. -------------- next part -------------- An HTML attachment was scrubbed... URL: From xi.chang01 at gmail.com Sat May 13 08:52:24 2023 From: xi.chang01 at gmail.com (Xi Chang) Date: Sat, 13 May 2023 07:52:24 -0700 Subject: [ncl-talk] scatter plot with a connected line Message-ID: Hi NCL, Is there any example how to create the attached plot using NCL? It is a scatter plot but with one connected line. Best, Chang -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Screen Shot 2023-05-13 at 7.50.13 AM.png Type: image/png Size: 83874 bytes Desc: not available URL: From olagueralyndonmark429 at gmail.com Fri May 19 21:48:17 2023 From: olagueralyndonmark429 at gmail.com (Lyndz) Date: Sat, 20 May 2023 11:48:17 +0800 Subject: [ncl-talk] How to plot significance of fractional variance (MJO) Message-ID: Dear NCL experts, I would like to plot the significance of the MJO (30-60 ISO) band similar to Figure 8 of Hartmann and Michelsen 1989 ( https://journals.ametsoc.org/view/journals/atsc/46/18/1520-0469_1989_046_2838_ipiir_2_0_co_2.xml ) [image: image.png] I found an example in the NCL webpage that plots the fractional variance: https://www.ncl.ucar.edu/Applications/Scripts/mjoclivar_3.ncl How can I plot the significance from here (above the red noise?) I'll appreciate any help on this. Sincerely, Lyndz -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image.png Type: image/png Size: 80286 bytes Desc: not available URL: From olagueralyndonmark429 at gmail.com Fri May 19 23:27:04 2023 From: olagueralyndonmark429 at gmail.com (Lyndz) Date: Sat, 20 May 2023 13:27:04 +0800 Subject: [ncl-talk] smoothed daily anomalies Message-ID: Dear NCL experts, Is there a function that I can use to remove the smoothed (2 harmonics) annual cycle from a one dimensional time series? I saw this example but for 3 dimensional time series only: https://www.ncl.ucar.edu/Document/Functions/Contributed/calcDayAnomTLL.shtml Any suggestions on how to do this in NCL? Attached is the sample data that I have. -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: test_rainfall.xlsx Type: application/vnd.openxmlformats-officedocument.spreadsheetml.sheet Size: 315079 bytes Desc: not available URL: From olagueralyndonmark429 at gmail.com Wed May 24 01:49:28 2023 From: olagueralyndonmark429 at gmail.com (Lyndz) Date: Wed, 24 May 2023 15:49:28 +0800 Subject: [ncl-talk] Plotting wavelet Message-ID: Dear NCL-Experts, I'm following example 1 of plotting the wavelet of daily rainfall from here: https://www.ncl.ucar.edu/Applications/wavelet.shtml How do I cut the y-axis to 128? Attached is the sample plot and script. Thanks in advance for the help! -Lyndz -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Screenshot from 2023-05-24 15-43-50.png Type: image/png Size: 200899 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: wavelet.ncl Type: application/octet-stream Size: 3980 bytes Desc: not available URL: From asphilli at ucar.edu Wed May 24 18:24:33 2023 From: asphilli at ucar.edu (Adam Phillips) Date: Wed, 24 May 2023 18:24:33 -0600 Subject: [ncl-talk] Plotting wavelet In-Reply-To: References: Message-ID: Hi Lyndz, Try setting res at trYMaxF = 128. (See XY example #5. ) If you have any further questions please respond to the ncl-talk email list. Adam On Wed, May 24, 2023 at 2:02?AM Lyndz via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dear NCL-Experts, > > I'm following example 1 of plotting the wavelet of daily rainfall from > here: > https://www.ncl.ucar.edu/Applications/wavelet.shtml > > How do I cut the y-axis to 128? > Attached is the sample plot and script. > > > Thanks in advance for the help! > > -Lyndz > _______________________________________________ > 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 G.Monte at isac.cnr.it Tue May 30 07:41:35 2023 From: G.Monte at isac.cnr.it (G.Monte at isac.cnr.it) Date: Tue, 30 May 2023 15:41:35 +0200 Subject: [ncl-talk] How to read CloudSat-GPM data Message-ID: Dears all, I have a question for you: can NCL read CloudSat-GPM data? I work with version 6.5. Thanks in advance, Giulio Monte From shea at ucar.edu Tue May 30 08:38:06 2023 From: shea at ucar.edu (Dennis Shea) Date: Tue, 30 May 2023 08:38:06 -0600 Subject: [ncl-talk] How to read CloudSat-GPM data In-Reply-To: References: Message-ID: %> ncl_filedump CloudSatFile | less With the appropriate file extension On Tue, May 30, 2023 at 7:41 AM G.Monte--- via ncl-talk < ncl-talk at mailman.ucar.edu> wrote: > Dears all, > > I have a question for you: can NCL read CloudSat-GPM data? I work with > version 6.5. > > Thanks in advance, > Giulio Monte > > _______________________________________________ > 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: