MeteoInfoLab脚本示例：计算气团轨迹

MeteoInfo

dtlyr 发表于 2016-6-23 10:46
脚本计算的轨迹不是像软件版的生成tgs文件吗？请问老师有没有关于meteinfolab脚本的教程？

参考这里：http://www.meteothinker.com/docs/meteoinfolab/index.html

453271826

MeteoInfo 发表于 2016-6-23 10:55
参考这里：http://www.meteothinker.com/docs/meteoinfolab/index.html

老师 有没有办法把脚本中轨迹计算的文件转化为tgs文件呢 或者 脚本计算出来的数据怎么在MeteoInfoLab中做聚类分析呢 TrajStat插件中的连续运算只能一个月一个月来 这边需要两年的后向轨迹数据 并进行聚类分析 王老师您有什么建议么

MeteoInfo

453271826 发表于 2016-11-30 17:11
老师 有没有办法把脚本中轨迹计算的文件转化为tgs文件呢 或者 脚本计算出来的数据怎么在MeteoInfoLab中做 ...

TrajStat里有转tgs文件的功能。本帖1楼第2个脚本可以一次计算多个月的轨迹。

453271826

MeteoInfo 发表于 2016-11-30 18:23
TrajStat里有转tgs文件的功能。本帖1楼第2个脚本可以一次计算多个月的轨迹。

好的 谢谢王老师 我再摸索一下

peteryu90

王老师，我按照您的代码在Lab里编辑完成后，如何才能运行呢，我是这方面的小白，一点都不懂，请您赐教

MeteoInfo

peteryu90 发表于 2017-4-10 16:52
王老师，我按照您的代码在Lab里编辑完成后，如何才能运行呢，我是这方面的小白，一点都不懂，请您赐教

工具栏里最后一个按钮（三角形）

peteryu90

MeteoInfo 发表于 2017-4-10 17:39
工具栏里最后一个按钮（三角形）

脚本里的数据有需要更改的地方吗？

littlelock

# Set working directory
metDir = 'I:\Tianhe\wrf_exp18'
outDir = 'I:\Tianhe\wrf_exp18'
workingDir = 'I:\Tianhe\wrf_exp18'
os.chdir(workingDir)
print 'Current directory: ' + os.getcwd()
# Set parameters
lon = '140.0'
lat = '30.0'
shour = '06'
heights = ['100.0','500.0','1000.0']
hnum = len(heights)
hours = '-12'
vertical = '0'
top = '10000.0'
# Set meteorological data files
fns = []
fn = 'test_wrfout1.arl'
fns.append(fn)
# Set start/end time
stime = datetime.datetime(2001,6,29)
# Write CONTROL file
ctFile = './CONTROL'
print stime.strftime('%Y-%m-%d ') + shour + ':00'
ctf = open(ctFile, 'w')
ctf.write(stime.strftime('%y %m %d ') + shour + "\n")
ctf.write(str(hnum) + '\n')
for j in range(0,hnum):
    ctf.write(lat + ' ' + lon + ' ' + heights[j] + '\n')
ctf.write(hours + '\n')
ctf.write(vertical + '\n')
ctf.write(top + '\n')
fnnum = len(fns)
ctf.write(str(fnnum) + '\n')
for j in range(0,fnnum):
    ctf.write(metDir + '/' + '\n')
    ctf.write(fns[j] + '\n')
ctf.write(outDir + '/' + '\n')
outfn = stime.strftime('traj_%Y%m%d')
ctf.write(outfn)
ctf.close()
# Calculate trajectories
os.system('c:/hysplit4/exec/hyts_std.exe')
print 'Finish...'

王老师，根据您提供的脚本，我把wrfout数据转化成arl，但对比了原来的wrfout和新的数据发现出现了数据的移位（附图），上面是我的脚本，请问这可能是什么原因呢？

MeteoInfo

littlelock 发表于 2017-4-18 21:48
# Set working directory
metDir = 'I:\Tianhe\wrf_exp18'
outDir = 'I:\Tianhe\wrf_exp18'

参考这里：http://www.meteothinker.com/exam ... ory/wrfout2arl.html

yydyoucan

王老师，这个脚本好像只能算gds1开头的同化数据，不能算RP开头同化数据