vibeta and uv2dvF for moisture divergence

Hi, everyone
     I'd like to calculate moisture flux divergence from WRF output , in
which Vibeta function is used for integrating the Qflux firstly , and
uv2dvF function is used for the divergence calculating. My questions are
like:
1. If you have current function to make it ,or that's the only way I can do
it ;
2. Could I can use vibeta function directly with data at a sigma level
instead of getting a sigma to P conversion ?If so, how I can do it?
3. What the linlog =2 (nonlinear, log interploation) will be used in vibeta
function?
4. How people process the missing data in vibeta function, I get some
missing data in surface layers at and below 850mb?
5. I have a specific question about the dismatched dimension problem about
"x" and psfc, could you help me to fix it?(error message: fatal: Number of
dimensions on right hand side do not match number of dimension in left hand
side
fatal: Execute: Error occured at or near line 82 in file quv.test.ncl)

Thanks very much for help.

 YAN

PS: the script I made ,U, V, Q have interploated into p-level.

....

p=(/1000.,925.,....50,30, 20,10) ;17 ncep level
u= a->U(:,:,0:114,0:126); 17L
v= a->V(:,:,0:114,0:126);17L
q= a->Q(:,:,0:114,0:126);17L
psfc = a->PSFC(:,0:114,0:126)
linlog =1
qu =q*u*1000
qv =q*v*1000
ptop =300.
pbot =1000.

82 vint_qu=vibeta(p,qu, linlog,psfc,pbot,ptop)/9.8 !!!problem line
83 vint_qv=vibeta(p,qv, linlog,psfc,pbot,ptop)/9.8 !!!problem line

dv=uv2dvF(vint_qu, vint_qv)
uvd =dv2uvF(dv)

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Nov 24 14:47:15 2009