[gradsusr] use "set gxout fwrite" to show 2D (t, z) average wind vectors

[Yi-Chih Huang]

Hello,

    I used "set gxout fwrite" to show 2D (t,z) average wind vectors as
below.  But the vertical wind field changes abruptly that is not a right
vertical wind field.  I have a feeling that the wind vectors were arranged
horizontal, then vertically -- similar to the sequence of an array in
Fortran.  But even I switched z and t, the pattern is the same style.
Could someone tell me how to show show 2D (t,z) average wind vectors?

    Thanks much,

                    Yi-Chih

###############
"reinit"
"open Morakot_d01.ctl"
"set fwrite -sq aveuv.dat"

tt=17
while(tt<=57)
"set t "tt""
zz=8
while(zz<=18)

"Set z "zz""

"Set x 65 91"
"Set y 37 65"
*"set x 1"
*"Set y 1"

"set gxout print"
"d aave(u,x=65,x=91,y=37,y=65)"
line = sublin(result,2)
u_value=subwrd(line,1)

say result
"d aave(v,x=65,x=91,y=37,y=65)"
line = sublin(result,2)
v_value=subwrd(line,1)


"set gxout fwrite"
"set x 1"
"set y 1"
"define uu="u_value""
"define vv="v_value""
"d uu"
"d vv"


"undefine uu"
"undefine vv"
zz=zz+1
endwhile
tt=tt+2

endwhile

"disable fwrite"
"close 1"

"reinit"
"open aveuv.ctl"
"set x 1"
"Set y 1"
"set z 1 11"
"set t 1 21"
"d uu;vv"
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111206/267cfde5/attachment-0003.html