NCL计算涡度平流

shoddy

请教各位大神，NCL有没有自带计算涡度平流的函数？
如没有，跪求计算涡度平流的代码……

vip545

这个不太清楚

TIANYU

dcx

请问楼主解决了吗？涡度怎么计算的

shoddy

dcx 发表于 2015-7-11 13:16
请问楼主解决了吗？涡度怎么计算的

涡度是用NCL自带的函数算的

当年明月

shoddy 发表于 2017-7-28 09:40
涡度是用NCL自带的函数算的

请问是哪个函数呀？

shoddy

当年明月 发表于 2017-11-13 23:59
请问是哪个函数呀？

"uv2vrF"。。。

liuxiaoyue

这个是ncl计算平流的函数
需要注意的是：变量坐标的纬度必须是单调递增的，如果不是需要调整一下坐标
像这样u = u(:,::-1,:)楼主先计算涡度，再使用这个函数计算涡度平流就行了~

shoddy

liuxiaoyue 发表于 2017-11-18 12:16
这个是ncl计算平流的函数
需要注意的是：变量坐标的纬度必须是单调递增的，如果不是需要调整一下坐标
像 ...

赞~~~~~~~~~

盈星

u         = short2flt(f1->u(a(x),{700:850},{22:38},{99:116}))
    v         = short2flt(f1->v(a(x),{700:850},{22:38},{99:116}))
        lat       = f1->latitude({22:38})
       
        f         = (2.*2.*3.1415926/86400)*sin(lat*3.1415926/180.)
        dx        = 6371000.*cos(lat*3.1415926/180.)*3.1415926/180.*0.5
        dy        = 6371000.*3.1415926/180.*0.5
       
        advvor    = new((/6,65,69/),float)
        advvor1   = new((/6,65,69/),float)
        advvor2   = new((/6,65,69/),float)
        advvor    = 0.0
        advvor1   = 0.0
        advvor2   = 0.0
        do z = 0,5,1
          do i = 1,63,1
            do j = 1,67,1
             advvor1(z,i,j) = -u(z,i,j)*((vor(z,i,j+1)-vor(z,i,j-1))/dx(i))-v(z,i,j)*((vor(z,i+1,j)-vor(z,i-1,j))/dy)
             advvor2(z,i,j) = -v(z,i,j)*(f(i+1)-f(i-1))/dy
            end do
          end do
        end do
        advvor = (advvor1 + advvor2)*1.e9
        copy_VarMeta(vor,advvor)
        printMinMax(advvor,True)