Moisture flux convergence (fwd)

[Arpita Mandal]

hi
i am calculating moisture flux vector using the foll script which
reinit
sdfopen oct2001uwindecm.nc
sdfopen oct2001sphumecm.nc
sdfopen oct2001slpecm.nc
set fwrite zfl.1998060100.288
set gxout fwrite
set x 1 145
set y 1 73
set t 1
define aa=vint(msl.3(z=1)/100.,q.2*u,300)
d aa
disable fwrite
reinit
sdfopen oct2001vwindecm.nc
sdfopen oct2001sphumecm.nc
sdfopen oct2001slpecm.nc
set fwrite mfl.1998100100.288
set gxout fwrite
set x 1 145
set y 1 73
set t 1
define bb=vint(msl.3(z=1)/100.,q.2*v,300)
d bb
disable fwrite
it works with ncep/ncar data, but when i use it with ecmwf data, it forve
s errors, as
world grid coordinates beyond data , error with q.2.
i have downloaded the files from ecmwf with levels 1000-10 for u,vwind and
1000-300 for sphum.
kindly help me

also when i open the files separately in grads. for u and v wind it shows
lev set from 10 10
but for sphum it shows lev set from 300 300.
i have downloaded the data in sequence but the levels are from high to low
both should show levels 1000 1000 as with necp


     WITH WARM REGARDS

     Dr Arpita Mandal
     Post Doctoral Fellow
     Planetary and Geosciences Division
     Physical Research Laboratory
     Navrangpura, Ahmedabad-380009
     Gujarat, India
     Tel No:07926314260
     Mob No:09427321259
     Email:arpita at prl.ernet.in
           arpita_mandal2000 at yahoo.com