画两坐标连线这一垂直剖面的TN通量

DZJ0019

Tevent1 = f1->Tevent1   ;(day|41,level|27,lat|181,lon|451)
    Hevent1 = f2->Hevent1
    WAFy = f3->wafy1
    WAFz = f4->wafz1

    leftlat  =  38
    rightlat =  81.85

    leftlon  =  -60
    rightlon =  30

    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

    trans1   = linint2_points(Tevent1&lon,Tevent1&lat,Tevent1,True,dist@gclon,dist@gclat,2)

    copy_VarAtts(Tevent1,trans1)

    trans1!0      = "day"         
    trans1&day    = Tevent1&day
    trans1!1      = "level"         
    trans1&level    = Tevent1&level

    trans2   = linint2_points(Hevent1&lon,Hevent1&lat,Hevent1,True,dist@gclon,dist@gclat,2)

    copy_VarAtts(Hevent1,trans2)

    trans2!0      = "day"         
    trans2&day    = Hevent1&day
    trans2!1      = "level"         
    trans2&level    = Hevent1&level

    trans3   = linint2_points(WAFy&lon,WAFy&lat,WAFy,True,dist@gclon,dist@gclat,2)

    copy_VarAtts(WAFy,trans3)

    trans3!0      = "day"         
    trans3&day    = WAFy&day
    trans3!1      = "level"         
    trans3&level    = WAFy&level

    trans4   = linint2_points(WAFz&lon,WAFz&lat,WAFz,True,dist@gclon,dist@gclat,2)

    copy_VarAtts(WAFz,trans4)

    trans4!0      = "day"         
    trans4&day    = WAFz&day
    trans4!1      = "level"         
    trans4&level    = WAFz&level

    wks = gsn_open_wks("pdf","North_East")
    res1                     = True         
    res1@tmXBMode            = "Explicit"   
    res1@tmXBValues          = (/points(0),points(npts-1)/)

    res1@tmXBLabels          = (/"lat:"+leftlat +", "+"lon:"+leftlon,"lat:"+rightlat+", "+"lon:"+rightlon/)

    res1@cnFillOn            = True         
    res1@cnFillPalette       = "NCV_blue_red"
    res1@gsnDraw = False
    res1@gsnFrame = False
    res1@cnLineLabelsOn = False
    res1@gsnAddCyclic = False
    res1@cnLinesOn           = False
    res1@lbLabelBarOn = True
    res1@cnInfoLabelOn = False      
    res1@pmLabelBarOrthogonalPosF = -0.05        
    res1@trYReverse          = True      
    ;res1@tmYRMode             = "Automatic"  
    res1@tiYAxisString = "Pressure(hpa)"
    res1@gsnLeftString = "TN01"
    res1@gsnRightString = "Temperature/Hgt"
    ;res1@cnLevelSpacingF     = 1.0
    res1@cnLevelSelectionMode = "ExplicitLevels"
    res1@cnLevels = ispan(-5,5,1)         

    res2 = True
    res2@cnFillOn = False
    res2@gsnDraw = False
    res2@gsnFrame = False
    res2@cnLinesOn = True
    res2@gsnContourNegLineDashPattern = 2
    res2@gsnContourPosLineDashPattern = 0
    res2@cnLineLabelsOn = True
    res2@cnLineLabelAngleF = 0.0
    res2@cnLevelSelectionMode = "ExplicitLevels"
    res2@cnLevels = ispan(-120,120,10)

    Tes2 = True     
    Tes2@gsnDraw = False
    Tes2@gsnFrame = False
    Tes2@vcRefMagnitudeF = 0.01           
    Tes2@vcRefLengthF = 0.04            
    Tes2@vcGlyphStyle = "CurlyVector"   
    Tes2@vcMinDistanceF = 0.01
    Tes2@gsnAddCyclic = False

    do i = 0,40,1
        res1@tiMainString = "day " + (i-30)
        plot1 = gsn_csm_pres_hgt(wks,trans1(i,:,:),res1)
        plot2 = gsn_csm_pres_hgt(wks,trans2(i,:,:),res2)
        plot3 = gsn_csm_vector(wks,trans3(i,:,:),trans4(i,:,:),Tes2)
        overlay(plot1,plot2)
        overlay(plot1,plot3)
        draw(plot1)
        frame(wks)
    end do不知道哪出问题了，垂直剖面上只有温度场和高度场的值，没有TN通量的值，求大佬指点