Re: [ferret_users] Calculting model agreement

[saurabh rathore <rohitsrb2020@xxxxxxxxx>]

Dear Kwesi

G'day

I 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. 

use slope_histruss.nc

let tsign = if slope ge 0 then 1 else (-1)

let agree80 = if abs(tsign[m=1:11@sum]) ge 8 then 1

fill agree80*slope[m=@ave]

On Wed, Mar 20, 2019 at 8:13 PM Kwesi A. Quagraine <starskykwesi@xxxxxxxxx> wrote:

Hello Ferreters

I 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 Einstein

Kwesi A. Quagraine

Department of Physics

School of Physical Sciences

College of Agriculture and Natural Sciences

University of Cape Coast

Cape Coast, Ghana

Alt. Email: kwesi@xxxxxxxxxxxxxx 

Web: http://www.recycleupghana.org/ 

Office: +27 21 650 3164

Skype: quagraine_cwasi

--

REGARDS

Saurabh Rathore

Research Scholar (PhD.)

Centre For Oceans, Rivers, Atmosphere & Land Science Technology

Indian Institute Of Technology, Kharagpur

contact :- 91- 8345984434

Attachment: trend_histensemble.jnl
Description: Binary data