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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 7585|回复: 6

[作图] [求助]ncl画海冰趋势,海岸线附近不连续

[复制链接]

新浪微博达人勋

发表于 2017-6-17 10:04:19 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 dongfang16 于 2017-6-17 10:06 编辑

求助:
画1971-2000年北极海冰在9,10,11个月的平均趋势。
用ncl画图,在海岸线附近图不连续,可以看到很明显的拐角,是不是画错了?

lat=192,lon=256


;==========================================================================
;==========================================================================
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"

begin

MODEL = 18
start = 1850
model = "MIROC-ESM-CHEM"

AD = 8
season = "SON"

  f1 = addfile("sic_OImon_MIROC-ESM-CHEM_historical_r1i1p1_185001-200512.nc","r")

  xa = f1->sic(0,:,:)
  xb = f1->sic(0,:,:)
  xc = f1->sic(0,:,:)

;==========================================================================
    aa    = 1971
    dd    = 1985

    x1    = f1->sic((aa-start)*12+AD:(dd-start)*12+AD:12,:,:)
    x2    = f1->sic((aa-start)*12+1+AD:(dd-start)*12+1+AD:12,:,:)
    x3    = f1->sic((aa-start)*12+2+AD:(dd-start)*12+2+AD:12,:,:)

    aa    = 1986
    dd    = 2000

    x4    = f1->sic((aa-start)*12+AD:(dd-start)*12+AD:12,:,:)
    x5    = f1->sic((aa-start)*12+1+AD:(dd-start)*12+1+AD:12,:,:)
    x6    = f1->sic((aa-start)*12+2+AD:(dd-start)*12+2+AD:12,:,:)

;==========================================================================

        xa = (dim_avg_n_Wrap(x1,0)+dim_avg_n_Wrap(x2,0)+dim_avg_n_Wrap(x3,0))/45
        xb = (dim_avg_n_Wrap(x4,0)+dim_avg_n_Wrap(x5,0)+dim_avg_n_Wrap(x6,0))/45
      
xc = xb - xa

;==========================================================================

    res            = True

    res@gsnPolar        = "NH"
    res@mpMinLatF        = 60.
    res@mpCenterLonF    = 90.
    res@mpGridLonSpacingF    = 30.
    res@mpGridLatSpacingF    = 10.
    res@mpOutlineBoundarySets    ="National"
    res@mpDataBaseVersion    = "MediumRes"

    res@cnFillOn        = True   
    res@cnLinesOn        = False  

    res@gsnLeftString    = season
    res@gsnCenterString    = ""
    res@gsnRightString    = "Trend"
    res@tiMainString    = model

    res@pmLabelBarHeightF    = 0.1  
    res@pmLabelBarWidthF    = 0.5
    res@cnLevelSelectionMode= "ManualLevels"
    res@cnMinLevelValF      = -1.5
    res@cnMaxLevelValF    = 1.5
    res@cnLevelSpacingF     = 0.15
    res@lbLabelFontHeightF = 0.015

    res@gsnMaximize        = True
  
;==========================================================================
wks    = gsn_open_wks("png",MODEL+"-S1T1-trd-sic-SON-test")

  gsn_define_colormap(wks,"BlueWhiteOrangeRed")
  gsn_reverse_colormap(wks)  
    plot    = new(1,graphic)   
    plot    = gsn_csm_contour_map_polar(wks,xc,res)
        
end
18-S1T1-trd-sic-SON-test.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-17 11:33:03 | 显示全部楼层
没做过,全当我瞎说哈

看这个样子,回不回和本身数据有关系?
我看到图整体样子还是蛮好看的,海陆也基本有对应,是不是数据本身分辨率有问题还是怎么?
可能可以直接画某一个时次的数据出来看看,先确认下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-17 12:53:22 | 显示全部楼层
MIROC 的格点大小大概2deg,这样的不平滑很正常
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-23 12:04:08 | 显示全部楼层
请问1986-2000   减去1970-1985  的9,10,11三个月的平均应该叫什么的趋势?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-26 12:56:31 | 显示全部楼层
藕茭白空心菜 发表于 2017-6-23 12:04
请问1986-2000   减去1970-1985  的9,10,11三个月的平均应该叫什么的趋势?

距平
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-26 12:56:35 | 显示全部楼层
藕茭白空心菜 发表于 2017-6-23 12:04
请问1986-2000   减去1970-1985  的9,10,11三个月的平均应该叫什么的趋势?

距平
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-7-10 14:42:25 | 显示全部楼层

距平不对吧,如果前后15年都是9,10,11月的季节平均的话,这个应该是合成分析之类的前后时间段的对比?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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