ParNCLDeveloperInfo/SimpleParNCLScripts: vort_homme.ncl

File vort_homme.ncl, 342 bytes (added by jayesh, 7 years ago)

Calculating vorticity from wind components on a HOMME grid

Line 
1begin
2fili = systemfunc("ls /media/wubi/home/jayesh/parvis/pargal-trunk/GridFiles/CAM/HOMME/datasets/homme26x3458.t[0-2].nc")  ; 3 files     
3
4setfileoption("nc","SuppressClose",False)   ; alter default mode
5
6f = addfiles(fili, "r") ; add the files
7
8u = f[:]->U ; read U
9v = f[:]->V ; read V
10
11vrt = uv2vrA(u,v) ; calc vort
12
13print(vrt)
14
15end
16