Re: [ferret_users] calcuting RMSE

[afwande juliet <afwandej965@xxxxxxxxx>]

now when I do like this i dont get error, but my main aim is to get one value for RSME for very location within my domain, My domain has about 4 boxes where I intend to calcuate RSME over . But when I save the file as .dat , i get so many values instead of expected one value. I have attached the .dat file for help

 set region/x=15E:54E/y=13n:18S
yes? let p=pre[d=1,l=1:3@ave]-pr[d=3,l=1:3@ave]
yes?  let q = p^2
yes? let s=q/3
yes? let m = s^0.5
yes? fill m;go land 1 "" 1
yes? list/clobber/file= test1.dat m
 LISTing to file test1.dat

On Mon, Oct 10, 2016 at 6:31 PM, afwande juliet <afwandej965@xxxxxxxxx> wrote:

My data has 3months I. E 3 time step for each file, so as you can see my l=1:3@ave, I regrided using cdo remapbil command, I used the cru data to regrid the wrf output

On Oct 10, 2016 6:04 PM, "Nicolas Freychet" <nicolas.freychet@xxxxxxxxx> wrote:

Hi Juliet,

It just seems that your 2 datasets do not have the same calendar or time record. If they only have 1 timestep, you can use the command "L=1" to precise that each variable should be associated to its first time:

let p=pre[d=2,L=1]-pr[d=1,L=1]

(if they have more than one timestep you may need to precise the time you consider differently, for example averaging the 6 first times L=1:6@ave.....)

2016-10-10 15:35 GMT+01:00 afwande juliet <afwandej965@xxxxxxxxx>:

Dear ALL, I am trying to calculate RMSE for two datasets over a box for MAM seasons for 1997 only. I have re-gridded the data so that have same grid information. But I get the error below. Could some help me. Also attached are two datasets am using

yes?  set region/x=32E:40E/y=8N:18S
yes? let p=pre[d=2]-pr[d=1]
yes? let q = p^2
yes? let r = q[l=1:3@sum]
yes? let s = r/3
yes? let m = s^0.5
yes? list m
           *** NOTE: Ambiguous coordinates on T axis: PRE[D=2]-PR[D=1]
 **ERROR: illegal limits: PR is not in the range T=35488:35579
          Axis extremes are T=19970301:19970532
yes?  set region/x=32E:40E/y=8N:18S
yes? let p=pre[d=2]-pr[d=1]
yes? let q = p^2
yes? let r = q[l=1:3@sum]
yes? let s = r/3
yes? let m = s^0.5
yes? list m
           *** NOTE: Ambiguous coordinates on T axis: PRE[D=2]-PR[D=1]
 **ERROR: illegal limits: PR does not exist at T=19970516
          Axis extremes are T=19970516

--

----------------------------------------------------

Nicolas Freychet

PDRA, School of Geosciences
University of Edinburgh
----------------------------------------------------

Attachment: test1.dat
Description: Binary data