Hi AnsleyMy observation data is spatial data, anom[time,lat,lon] and they are satellite observation data and my model data consists of 8 models (which I have aggregated, as per initial directions given here).I will ideally want to check if the models agree to the sign of change of the observation data.From the idea from Saurabh, that would work if I were to check that the models agree to the sign of change from model mean (which I understand), However, how will I be able to modify the commands from Saurabh, so that my script will compare the the models to the mean of the observation data?.Hope I’ve expressed myself with some clarity.Thanks once again.KwesiOn 22 March 2019 at 00:40:44, Ansley C. Manke (ansley.b.manke@xxxxxxxx) wrote:
I'd suggest making plots of the different expressions you are using to see if they're doing what you expect.
I can't tell what is "shape" of your observations? Are they timeseries at different stations?
A different idea would be to have a look at the sample* functions. Would it make sense to sample the model data at the times / locations of the observations?
yes? show function sample*
On 3/21/2019 7:36 AM, Kwesi Quagraine wrote:Hello Saurabh, Hello FerretersThanks for the script and guidelines.However, I want to know, how I can compare these models (8 in my case), say, precipitation to observation data (chirps), what would be the approach?Because here, my guess is I'll be comparing that to the observation mean. And how to do that is my frustration.RegardsKwesi
On Thu, Mar 21, 2019, 14:47 saurabh rathore <rohitsrb2020@xxxxxxxxx> wrote: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
--
REGARDS
Saurabh RathoreResearch Scholar (PhD.)Centre For Oceans, Rivers, Atmosphere & Land Science TechnologyIndian Institute Of Technology, Kharagpurcontact :- 91- 8345984434------------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