使用NCL沿着任意两点做垂直剖面图

kisscate

Mark一下，学习了

wawatou83

占位…………

wawatou83

占位…………

meng333

{:eb502:}

圆梓君

zhangxiao 发表于 2016-8-12 16:16
你看一下这个帖子 http://bbs.06climate.com/forum.php?mod=viewthread&tid=21983&extra=&page=1

请问有画过垂直剖面风场吗？不知道怎么叠加感觉总是不对哇

圆梓君

吴阿好大人 发表于 2016-8-12 16:46
师兄用WRF后处理的数据直接就有NCL函数进行斜剖面绘图，但是如果用别的数据就没办法了。。。

请问你师兄有没有画过风场剖面哇。。。感觉这个坎迈不过去了

342901wzx

;计算绘制剖面图的数据,；(注：坐标超过地图范围，则无法绘图)
        leftlat  = 30.75
        rightlat = 33.75
        leftlon  = 115.25
        rightlon = 118.25

        npts     = 100 ; number of points in resulting transect
        dist     = gc_latlon(leftlat,leftlon,rightlat,rightlon,npts,2) ;计算剖面直线的坐标点(经度、纬度)
        points   = ispan(0,npts-1,1)*1.0 ;计算剖面直线点的序列
        ;printVarSummary(dist)
        ;print("("+dist@gclat + " , " + dist@gclon + ")")
      
        ;计算显示剖面的数据
       ;trans   = linint2_points(t&lon_t,t&lat_t,t,True,dist@gclon,dist@gclat,2)  ;计算的数据不含坐标
        Tt   = linint2_points_Wrap(lon({leftlon:rightlon}),lat({leftlat:rightlat}),t(:,:,{leftlat:rightlat},{leftlon:rightlon}),True,dist@gclon,dist@gclat,2)  ;计算的数据含坐标
        Tu   = linint2_points_Wrap(lon({leftlon:rightlon}),lat({leftlat:rightlat}),u(:,:,{leftlat:rightlat},{leftlon:rightlon}),True,dist@gclon,dist@gclat,2)
        Tv   = linint2_points_Wrap(lon({leftlon:rightlon}),lat({leftlat:rightlat}),v(:,:,{leftlat:rightlat},{leftlon:rightlon}),True,dist@gclon,dist@gclat,2)

。。。。。。
          uvplot4      =  gsn_csm_vector(wks,   Tv(Lst,:,:), Tw1(Lst,:,:), wvres)   ;垂直风场
          dplot4       =  gsn_csm_contour(wks,  Td(Lst,:,:), dres)                  ;散度
          dplot4       =  ColorNegDashZeroPosContour(dplot4,"blue","black","brown")

  
          plots(4)    =  dplot4
          overlay(plots(4) ,uvplot4)
。。。。
gsn_panel(wks, plots, (/2,2/), panres)

342901wzx

;计算绘制剖面图的数据,；(注：坐标超过地图范围，则无法绘图)
         leftlat  = 30.75
         rightlat = 33.75
         leftlon  = 115.25
         rightlon = 118.25

         npts     = 100 ; number of points in resulting transect
         dist     = gc_latlon(leftlat,leftlon,rightlat,rightlon,npts,2) ;计算剖面直线的坐标点(经度、纬度)
         points   = ispan(0,npts-1,1)*1.0 ;计算剖面直线点的序列
        ;printVarSummary(dist)
         ;print("("+dist@gclat + " , " + dist@gclon + ")")
        
         ;计算显示剖面的数据
       ;trans   = linint2_points(t&lon_t,t&lat_t,t,True,dist@gclon,dist@gclat,2)  ;计算的数据不含坐标
        Tt   = linint2_points_Wrap(lon({leftlon:rightlon}),lat({leftlat:rightlat}),t(:,:,{leftlat:rightlat},{leftlon:rightlon}),True,dist@gclon,dist@gclat,2)  ;计算的数据含坐标
        Tu   = linint2_points_Wrap(lon({leftlon:rightlon}),lat({leftlat:rightlat}),u(:,:,{leftlat:rightlat},{leftlon:rightlon}),True,dist@gclon,dist@gclat,2)
         Tv   = linint2_points_Wrap(lon({leftlon:rightlon}),lat({leftlat:rightlat}),v(:,:,{leftlat:rightlat},{leftlon:rightlon}),True,dist@gclon,dist@gclat,2)

。。。。。。
          uvplot4      =  gsn_csm_vector(wks,   Tv(Lst,:,:), Tw1(Lst,:,:), wvres)   ;垂直风场
          dplot4       =  gsn_csm_contour(wks,  Td(Lst,:,:), dres)                  ;散度
          dplot4       =  ColorNegDashZeroPosContour(dplot4,"blue","black","brown")

   
           plots(4)    =  dplot4
           overlay(plots(4) ,uvplot4)
。。。。
gsn_panel(wks, plots, (/2,2/), panres)

以上程序，基本涵盖所有任意剖面绘图，结合官方函数，搞定没问题

lulunuist

342901wzx 发表于 2020-12-12 09:45
;计算绘制剖面图的数据,；(注：坐标超过地图范围，则无法绘图)
         leftlat  = 30.75
         righ ...

楼主你好：
uvplot4      =  gsn_csm_vector(wks,   Tv(Lst,:,:), Tw1(Lst,:,:), wvres)   ;垂直风场
这个命令只画了v和w，不能反映真实的垂直环流吧?
如果再加上考虑u风不知道应该怎么处理

Eudora99

圆梓君 发表于 2018-6-30 07:36
请问有画过垂直剖面风场吗？不知道怎么叠加感觉总是不对哇

请问你解决这个问题了么？我的垂直风场画上去也不对，求告知正确的程序段。