! modified for use of Seattle filenames<br>
<br>
use acom3-cell-thickness.50s-20n.cdf ! from topog_acom3.1.dta.nc<br>
<br>
use subdomain_notides-mean.50s-20n.cdf<br>
use subdomain_notides-mean-w.50s-20n.cdf<br>
<br>
! note that e435 (dhu/100 from original) is already on the (x,y) u grid<br>
let delh = e435[d=acom3-cell-thickness.50s-20n,gz=v_30d[d=subdomain_notides-mean.50s-20n]@asn]<br>
<br>
let u1 = missing(u_30d[d=subdomain_notides-mean.50s-20n]/100 , 0) ! convert to m/s<br>
let v1 = missing(v_30d[d=subdomain_notides-mean.50s-20n]/100 , 0)<br>
<br>
let ux = u1[x=@ddb]<br>
let vy = v1[y=@ddb]<br>
<br>
let uxbar = 0.5*(ux + ux[j=@shf:-1])<br>
let vybar = 0.5*(vy + vy[i=@shf:-1])<br>
<br>
let div = delh * (uxbar + vybar)<br>
let vzint = div[k=1:36@rsum]<br>
let vonw = vzint[gz=w_30d[d=subdomain_notides-mean-w.50s-20n]@asn]<br>