WRF UTC时间转化为BEIJING时间NCL脚本

870744170

第一次发帖！通过调用这两个函数，可以将WRF输出的UTC时间转化为beijing时间

hzwjy

  fi = addfile("./wrfout_d01.nc", "r")

  Times = wrf_user_getvar(fi,"Times",-1)
  TimeUTC = wrf_times_c(Times,1)
  TimeBJT = TimeUTC + 8
  copy_VarAtts(TimeUTC, TimeBJT)
  BJT = cd_calendar(TimeBJT,-3)

unice

好东西，赞一个

longlivehj

时间格式化可以用cd_string。

870744170

longlivehj 发表于 2014-11-26 01:39
时间格式化可以用cd_string。

这样就不需要substring函数了哈哈！

hhumuer

        seconds            = ut_inv_calendar(year,month,day,00,00,00,units,0)+(8+hour)*3600
敢问下面几行什么意思呢，看不太懂~~~

      p2 = dimsizes (main) - 2                ; go to last avail char in main string


        seconds@units   = units
        utc                = ut_calendar(seconds,0)
;---------------------------------------------------------------------------------------
                if((utc(0,1).lt.10).and.(utc(0,2).lt.10).and.(utc(0,3).lt.10))then
                        beijing = utc(0,0)+"-0"+utc(0,1)+"-0"+utc(0,2)+" 0"+utc(0,3)+":00:00"
                end if

870744170

hzwjy 发表于 2014-11-26 19:30
fi = addfile("./wrfout_d01.nc", "r")

  Times = wrf_user_getvar(fi,"Times",-1)

{:eb502:}{:eb502:}

870744170

hhumuer 发表于 2014-11-26 17:22
seconds            = ut_inv_calendar(year,month,day,00,00,00,units,0)+(8+hour)*3600
敢问下面几行什么意思 ...

ut_calendar函数根据输入的seconds，以数组的形式转换为日期,这里返回给utc这个变量。utc(0,0)代表的是年，utc(0,1)代表的是月。。。。为了保证日期以2001-09-02 00:00:00的形式输出，要判断utc(0,1)是否大于10....方法比较傻。。。还是参照下面的回复吧！

空谷幽竹

hzwjy 发表于 2014-11-26 19:30
fi = addfile("./wrfout_d01.nc", "r")

  Times = wrf_user_getvar(fi,"Times",-1)

wrf_times_c 函数无法识别啊？

空谷幽竹

空谷幽竹 发表于 2014-11-29 16:11
wrf_times_c 函数无法识别啊？

好吧，查了下知道了