[gradsusr] ave function error

[Eric Altshuler]

Hello Toshiya, 


Your increment value of 4.0 (degrees) for xdef and ydef is much too large. The resulting domain will cover 4000x1000 degrees, which is absurd (the entire globe is 360x180). Try using a much smaller value such as 0.02 (domain will be roughly 20x5 degrees). Since you have a Cartesian grid that's not associated with any geographic coordinate system, the domain's location and map background won't really matter. 


Best regards, 


Eric L. Altshuler 
Research Associate 
Center for Ocean-Land-Atmosphere Studies 
253 Research Hall, Mail Stop 6C5 
George Mason University 
4400 University Drive 
Fairfax, VA 22030 USA 

E-mail: ela at cola.iges.org 
Phone: (703) 993-5725 
Fax: (703) 993-5770 

----- Original Message -----
From: "吉田 敏哉" <yoshida_t at storm.dpri.kyoto-u.ac.jp> 
To: gradsusr at gradsusr.org 
Sent: Thursday, June 11, 2015 3:04:55 AM 
Subject: [gradsusr] ave function error 


Dear users, 


I am using CFD output computed on Cartesian coordinates. 
My ctl.file is as follows: 


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

DSET ^data10 
TITLE LES 
UNDEF -9.99E33 
OPTIONS byteswapped 
XDEF 1000 LINEAR 0.000 4.000 
YDEF 250 LINEAR 2.000 4.000 
ZDEF 105 LEVELS 
0.500 1.550 2.705 3.976 ... 
TDEF 1 LINEAR 00:01z01jan2000 1mn 
VARS 3 
u 105 99 u winds 
w 105 99 w winds 
v 105 99 v winds 
ENDVARS 
-------------------------------------------------------- 


When I use ave function for Y-axis direction, I have the following error: 



ga-> d ave(u,y=24,y=250) 
Averaging. dim = 1, start = 24, end = 250 
Error from AVE: Internal logic check 100 
Operation Error: Error from ave function 
Error ocurred at column 1 
DISPLAY error: Invalid expression 
Expression = ave(u,y=24,y=250) 


If I change "ave(u,y=24,y=250)" to "ave(u,y=1,y=250)" , it work normally. 
In addition, "ave(u,y=1,y=250)" and "ave(u,y=1,y=23)" are exactly the same value so 
I think ave function isn't working normally more than y=23. 


I confirmed the CFD output was correct with GrADS and 

this error only occur when I use ave function for Y-axis direction. 


I tried using 2.0.21,2.0.2, 2.0.a9 and 2.1.a2 versions but all of them don't work normally. 
I didn't know what to do, do you have any ideas? 






Thank you for being patient with my English. 
Toshiya 
_______________________________________________ 
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/20150611/0867a7f7/attachment.html