[gradsusr] RMSE for accumulated Rainfall

[Lyndon Mark Olaguera]

ow...I forgot to remove the parenthesis..I have corrected it in my original
script..
Thanks


On Fri, May 17, 2013 at 7:22 AM, Eric Altshuler <ela at cola.iges.org> wrote:

> There is still a sign error in your calculation of rn1 and rn2. It should
> be:
>
> 'define rn1=(rainc+rainnc)-(rainc(t='Time_prev')+rainnc(t='Time_prev'))'
> 'define rn2=(rainc.2+rainnc.2)-(rainc.2(t='Time_prev')+
> rainnc.2(t='Time_prev'))'
>
> ------------------------------
> *From: *"Lyndon Mark Olaguera" <olagueralyndonmark429 at gmail.com>
> *To: *"GrADS Users Forum" <gradsusr at gradsusr.org>
> *Sent: *Wednesday, May 15, 2013 9:47:48 PM
> *Subject: *Re: [gradsusr] RMSE for accumulated Rainfall
>
>
> thanks for the kindest reply..This is a big help..
> Here's my modified script:
>
> 'open 1D_nudging/surfnudge/surface_grid_nudging.ctl'
> 'open 1D_nudging/model/fnl_nonudging.ctl'
>
> prompt 'Max Time: '
> pull Max_time
> prompt 'Title Header: '
> pull Title_header
> Time_val=2
> Time_prev=1
> while(Time_val<=Max_time)
>
> 'set t 'Time_val
> 'define rn1=(rainc+rainnc)-(rainc(t='Time_prev')-rainnc(t='Time_prev'))'
> 'define
> rn2=(rainc.2+rainnc.2)-(rainc.2(t='Time_prev')-rainnc.2(t='Time_prev'))'
> 'define dif=rn1-rn2'
> 'define sq=dif*dif'
> 'define am=amean(sq,lon=115,lon=135,lat=4,lat=25)'
> 'define rmse=sqrt(am)'
> 'd rmse'
>
> rmse=subwrd(result,4)
> output=rmse
> res=write(rmse_rain.csv,output,append)
>
> Time_val=Time_val+1
> Time_prev=Time_prev+1
> endwhile
> quit
>
>
>  This is working now..Thanks for your suggestions.
>
>
>
>
> On Thu, May 16, 2013 at 12:20 AM, Jeff Duda <jeffduda319 at gmail.com> wrote:
>
>> Yes, you can choose to leave the expression in that form.  What I have
>> done is simply to distribute the minus sign.  The expressions are
>> mathematically identical.
>>
>> Jeff
>>
>>
>> On Wed, May 15, 2013 at 2:53 AM, Lyndon Mark Olaguera <
>> olagueralyndonmark429 at gmail.com> wrote:
>>
>>>
>>> On Wed, May 15, 2013 at 1:06 PM, Jeff Duda <jeffduda319 at gmail.com>wrote:
>>>
>>>> 'define fcstrain = rainc.1 + rainnc.1 - rainc.1(t-1) - rainnc.1(t-1)'
>>>> 'define obsrain = rainc.2 + rainnc.2 - rainc.2(t-1) - rainnc.2(t-1)'
>>>> 'define rainerror = fcstrain - obsrain'
>>>>
>>>
>>> Hello,
>>> Thanks for the reply..
>>> I have another question on this line..
>>>
>>> 'define fcstrain = rainc.1 + rainnc.1 - rainc.1(t-1) - rainnc.1(t-1)'
>>> 'define obsrain = rainc.2 + rainnc.2 - rainc.2(t-1) - rainnc.2(t-1)'
>>> 'define rainerror = fcstrain - obsrain'
>>>
>>> should this be
>>>
>>> 'define fcstrain = (rainc.1 + rainnc.1) - (rainc.1(t-1) + rainnc.1(t-1))'
>>> 'define obsrain =( rainc.2 + rainnc.2) - (rainc.2(t-1) + rainnc.2(t-1))'
>>>
>>> plus not minus?
>>>
>>>
>>>
>>> I will try this later
>>> Lyndon
>>>
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>>>
>>
>>
>> --
>> Jeff Duda
>> Graduate research assistant
>> University of Oklahoma School of Meteorology
>> Center for Analysis and Prediction of Storms
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130517/ca06600a/attachment-0003.html