【已解决】NCL FNL 850hPa温度场和风场

channam

用2020年1月份FNL数据做850hPa温度场和10m风速平均，发现温度场北边温度比南边高，感觉有问题，而且平均风速偏小（2m/s的参考矢量），不知道原因在哪里，另外想问一下FNL数据有缺测值吗，我在绘图时候没考虑过缺测的问题。  求助大家吖~

部分相关脚本在这里：

;;read;;
        dir   = "./fnl/"
        files = systemfunc("ls " + dir + "fnl_202001*.grib2")     
        f   = addfiles(files,"r")
        ListSetType (f, "join")   
        u = f[:]->UGRD_P0_L103_GLL0(:,{10},:,:)
        v = f[:]->VGRD_P0_L103_GLL0(:,{10},:,:)
        hgt =f[:]->HGT_P0_L100_GLL0(:,{85000},:,:)
        tmp =f[:]->TMP_P0_L100_GLL0(:,{85000},:,:)
        tmp = tmp - 273.15
        ..........

;;plot;;
        map     = gsn_csm_map(wks,mpres)
        contour1 = gsn_csm_contour(wks,dim_avg_n(hgt,0),hgtres)   
        contour2 = gsn_csm_contour(wks,dim_avg_n(tmp,0),tmpres)
        vector  = gsn_csm_vector(wks,dim_avg_n(u,0),dim_avg_n(v,0),res_vc)   

ncl之后出现这种warning，感觉不影响，不知道和这个有关吗。
warning:ContourPlotPreDraw: out of range coordinates encountered; standard  rendering method may be unreliable;
consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values
warning:ContourPlotDraw: out of range coordinates encountered; standard AreaFill rendering method may be unreliable;
consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values
warning:ContourPlotDraw: out of range coordinates encountered; standard  rendering method may be unreliable;
consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values

后来自己找到的解决方法：dim_avg_n变为dim_avg_n_Wrap，保留原数据信息包括经纬度。

;;plot;;
        map      = gsn_csm_map(wks,mpres)
        contour1 = gsn_csm_contour(wks,dim_avg_n_Wrap(tmp,0),tmpres)  ;使用dim_avg_n_Wrap做平均可保留原数据信息包括经纬度等不易出错
        contour2 = gsn_csm_contour(wks,dim_avg_n_Wrap(hgt,0),hgtres)  ;对变量第一维做平均
        vector   = gsn_csm_vector(wks,dim_avg_n_Wrap(u,0),dim_avg_n_Wrap(v,0),res_vc)

珮瑶瑶

LZ你好，完整的脚本分享一下吧，谢谢

channam

哈哈哈哈或或或 发表于 2020-5-11 23:11
我想知道你的色板用的是什么呀

用的 cmp_b2r

哈哈哈哈或或或

我想知道你的色板用的是什么呀

@me

想问一下2和箭头怎么放到图形的右上方？

淡定冷静

可以看一下你完整的脚本吗

临水照花

请问方便分享完整代码吗？在网络上都找不到画风场的

1297712690

ncar的fnl数据好像要转置一下？我拿matlab读的时候是这样，不知道楼主是什么情况

独角兽

请问怎么解决的呢