请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3629|回复: 5

[作图] ncl出不了图

[复制链接]

新浪微博达人勋

发表于 2023-3-9 20:09:29 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
用ncl处理nc文件,数据都可以正常输出,也不报错,就是出不了图,有没有大佬帮我看一下哪里出错了
loadscript("$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl")
loadscript("$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl")
loadscript("$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl")

begin
        f   = addfile ("D:/dataPRE/PRE_MON_1961-2020.nc","r")
        pre = f->PRE_MON(:,:,:)             ;读取1961-2020的降水数据

        pre_JJA = month_to_season(pre,"JJA")   ;青藏高原地区夏季平均降水
;        pre_JJA = (pre,"JJA")   ;青藏高原地区夏季平均降水
        pre_sum = pre_JJA*3      ;青藏高原地区夏季累计降水,单位化为mm

  printVarSummary(pre_sum)
        copy_VarCoords(pre_JJA,pre_sum);复制参数
        time = ispan(1961,2020,1);读取时间变量数据
       
        pre_avg = new((/60/),"float");新建维度存放区域夏季降水量时间序列
        pre_jp  = new((/60/),"float");新建维度存放区域夏季降水距平时间序列

    latS     = 25         
    latN     = 42
    lonW     = 73
    lonE     = 106
    lat     = f->lat
    lon     = f->lon
   

   pre_avg = dim_avg_n_Wrap(pre_sum(:,:,:),(/1,2/))

        jp = dim_rmvmean_Wrap(pre_avg)

     ;   opt = True
      ; opt@spi_type = 3 ; calculate using Pearson III distribution

       jp    = dim_spi_n(pre_avg, 12, False, 0)

   
      ; dmean =  avg(pre_avg(:49))
       ; jp = (pre_avg - dmean)/dmean*100
        jp@units="mm"
        print(jp)
        print(pre_avg)
     ;  print("dmean=="+dmean)

        wks=gsn_open_wks("pdf","TP-pre_time");设置绘图空间
        res=True
        res@gsnXYBarChart=True;绘制柱状图
        res@gsnMaximize = True
        res@tmXBMode="Explicit";坐标轴标注
        res@tmXBValues=time
        res@tmXBLabels=time
        res@tmXBLabelStride=5;每5个刻度添加标注
        res@tmXBLabelFontHeightF=0.1;标注字号
        res@tmYLLabelFontHeightF=0.1
        res@tmXBLabelDirection="Down"
        res@tmXBLabelFontHeightF=0.1
        res@tmXBLabels=time
        res@gsnScale = True   ;xy轴统一尺寸            
        res@vpHeightF = 1.5;坐标轴高度
        res@vpWidthF  = 6;坐标轴宽度
        ;res@trYMinF   = -12;y轴范围
        ;res@trYMaxF   = 12
        res@tiYAxisString = "Precipitation Anomaly/mm" ; 标题
        res@tiXAxisString = "year"      
        res@tiMainString  = "Summer Precipitation Anomaly in TP of China"
        res@tiMainFontHeightF=0.19;标题字号
        res@tiXAxisFontHeightF=0.1
        res@tiYAxisFontHeightF=0.1
        res@gsnYRefLine = 0 ; 参考线数值  
        res@gsnAboveYRefLineColor = "red"           ; 参考线上下颜色
        res@gsnBelowYRefLineColor = "blue"         
        ;res@xyLineColors = (/"red","blue"/)         ; 折线颜色
         ;res@pmLegendDisplayMode="Always";图例显示方式
         ;res@xyExplicitLegendLabels=(/"Data","Trend"/);设置图例标注文本
         ;res@pmLegendHeightF=0.05 ;图例高度
         ;res@pmLegendWidthF=0.1 ;图例宽度
         ;res@pmLegendOrthogonalPosF=-1.4;图例垂直移动
         ;res@pmLegendParallelPosF=0.1;图例水平移动
         plot=gsn_csm_xy(wks,time,jp,res);绘制距平直方图
         
      end
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2023-3-10 09:21:18 | 显示全部楼层
检查一下ncl打开路径的文件夹里有没有pdf格式的图
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2023-3-10 16:26:03 | 显示全部楼层
guoguohh 发表于 2023-3-10 09:21
检查一下ncl打开路径的文件夹里有没有pdf格式的图

一开始是画出来了一张图,然后发现了点错误,改了一下标题名字,结果就出不了新图,也不知道出了什么问题
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2023-3-10 20:37:19 | 显示全部楼层
synt 发表于 2023-3-10 16:26
一开始是画出来了一张图,然后发现了点错误,改了一下标题名字,结果就出不了新图,也不知道出了什么问题 ...

新的图是不是把原来的图覆盖了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2023-3-11 13:53:01 | 显示全部楼层
guoguohh 发表于 2023-3-10 20:37
新的图是不是把原来的图覆盖了

没有覆盖,还是原来的图
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2023-3-11 14:04:45 | 显示全部楼层
guoguohh 发表于 2023-3-10 20:37
新的图是不是把原来的图覆盖了

找到新图了,在ncl的其中一个文件夹里,之前都没发现,总之谢谢大佬了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表