Dear KwesiG'dayI tried to use previous suggestion given by Russ, Ryo and other on my data. SO my data consists of 11 CMIP5 models and I want to show atleast 7 or 8 models agreement in sign. So I did it in this way. Slope is my linear trend from 11 cmip5 models. But before applying this method I combined all my slope from all the different models into a single netcdf file (attached jnl script). So that I can do the computation in my desired dimension which is 'm' here. Here 'm' represents the total number of models or ensembles.I hope this will help.let tsign = if slope ge 0 then 1 else (-1)let agree80 = if abs(tsign[m=1:11@sum]) ge 8 then 1fill agree80*slope[m=@ave]On Wed, Mar 20, 2019 at 8:13 PM Kwesi A. Quagraine <starskykwesi@xxxxxxxxx> wrote:Hello FerretersI am trying to make a plot with model agreement.I want to plot model agreement over some observation. Hence, how the models agree to the sign of change.Attached is a snippet of my script. I would kindly want to know if I am doing it right and if not, the best way to approach it.Thanks for your help in advance.Kwesi------------Try not to become a man of success but rather a man of value- Albert EinsteinKwesi A. QuagraineDepartment of PhysicsSchool of Physical SciencesCollege of Agriculture and Natural SciencesUniversity of Cape CoastCape Coast, GhanaAlt. Email: kwesi@xxxxxxxxxxxxxxOffice: +27 21 650 3164Skype: quagraine_cwasi--REGARDSSaurabh RathoreResearch Scholar (PhD.)Centre For Oceans, Rivers, Atmosphere & Land Science TechnologyIndian Institute Of Technology, Kharagpurcontact :- 91- 8345984434