额外变量输出

Meteorological Field Descriptions

https://ready.arl.noaa.gov/hypub/flddesc.html

FIELD ID CONVERSION FACTOR UNITS CONTOUR INTERVAL DESCRIPTION
3 NULL 1.00 0.00 MISSING DATA
1 TMPS 1.00 DEGC 0.00 SURFACE TEMPERATURE
1 VSBY .001 KM 10.00 VISIBILITY
1 MSLP 1.00 HPA 4.00 MEAN SEA-LEVEL PRESSURE
1 PRSS 1.00 HPA 0.00 SURFACE PRESSURE
1 ICWT 1.00 0--1 .00 ICE-COVERED WATER
1 SNOW 1.00 0--1 .00 SNOW COVERAGE
1 SNWD 1.00 CM 1.00 SNOW DEPTH
1 SHGT 1.00 M 0.00 SURFACE HEIGHT
1 TOPO 1.00 M 0.00 TOPOGRAPHY
1 CPPT 1.00 MM 5.00 CONVECTIVE PRECIPITATION
1 TPP1 1.00 M 1.00 1 HOUR ACCUMULATED PRECIPITATION
1 TPP3 1.00 MM 1.00 3-HOUR ACCUMULATED PRECIPITATION
1 CPP3 1.00 MM 1.00 3-HOUR ACC. CONVECTIVE PRECIPITATION
1 TPP6 1.00 MM 5.00 6-HOUR ACCUMULATED PRECIPITATION
1 CPP6 1.00 MM 5.00 6-HOUR ACC. CONVECTIVE PRECIPITATION
1 TPPT 1.00 MM 5.00 12-HOUR ACCUMULATED PRECIPITATION
1 TPPA 1.00 MM 5.00 TOTAL ACCUMULATED PRECIPITATION
1 PRAT 1.00 KGM2 0.00 PRECIPITATION RATE
1 HFLX 1.00 W/M2 25.00 LATENT HEAT FLUX
1 SHTF 1.00 W/M2 0.00 SENSIBLE HEAT NET FLUX
1 LHTF 1.00 W/M2 0.00 LATENT HEAT NET FLUX
1 ULWF 1.00 W/M2 0.00 UPWARD LONG WAVE RADIATION FLUX
1 DLWF 1.00 W/M2 0.00 DOWNWARD LONG WAVE RADIATION FLUX
1 DSWF 1.00 W/M2 0.00 DOWNWARD SHORT WAVE RADIATION FLUX
1 QSTR 1000.00 G/KG 1.00 FLUX MIXING RATIO
1 WFLX 1000.00 GM2S 1.00 WATER VAPOR FLUX
1 EXCO 1000.00 GM2S 0.00 EXCHANGE COEFFICIENT
3 MXLR 1.00 0--N 1.00 NUMBER OF MIXED SIGMA LAYERS
1 MXHT 1.00 M 0.00 MIXED LAYER DEPTH
1 PBLH 1.00 M 0.00 PLANETARY BOUNDARY LAYER HEIGHT
1 THTS 1.00 DEFC 5.00 SURFACE POTENTIAL TEMPERATURE
1 U10M 1.00 M/S .00 10 M U-WIND COMPONENT
1 V10M 1.00 M/S .00 10 M V-WIND COMPONENT
1 UMOF 1.00 N/M2 0.00 MOMENTUM FLUX, U-WIND COMPONENT
1 VMOF 1.00 N/M2 0.00 MOMENTUM FLUX, V-WIND COMPONENT
1 MOMF 1.00 N/M2 0.00 TOTAL MOMENTUM FLUX
1 DEPV 1.00 M/S 0.01 DEPOSITION VELOCITY
1 USTR 100.00 M/S 10.00 FRICTION VELOCITY
1 TSTR 1.00 DEGC 0.00 FLUX TEMPERATURE
1 VGTP 1.00 0--1 1.00 VEGETATION TYPE
1 T02M 1.00 DEGC 5.00 2 M TEMPERATURE
2 TEMP 1.00 DEGC 0.00 TEMPERATURE
2 UWND 1.00 M/S 5.00 U-WIND COMPONENT
2 VWND 1.00 M/S 5.00 V-WIND COMPONENT
2 FLAG 1.00 KNTS 0.00 WIND FLAGS
2 WVCT 1.00 KNTS 0.00 WIND VECTORS
2 STRM 1.00 KNTS 0.00 STREAMLINES
2 WSPD 1.00 KNTS 5.00 WIND SPEED
2 WWND 3600.00 MB/H 5.00 W-WIND COMPONENT
2 SPHU 1000.00 G/KG 1.00 SPECIFIC HUMIDITY
2 RELH 1.00 PCT 20.00 RELATIVE HUMIDITY
1 DP2M 1.00 DEGC 0.00 2 M DEW POINT
2 HGTS 0.10 DM 6.00 HEIGHT
2 TKEN 1.00 JOUL 0.00 TOTAL KINETIC ENERGY
2 PRES 1.00 HPA 4.00 PRESSURE
1 THKN 0.10 DM 6.00 THICKNESS
3 HGT1 1.00 M 60.00 1000 MB HEIGHT
3 HGT5 1.00 M 60.00 500 MB HEIGHT
3 ABSV 100000.00 /S 2.00 ABSOLUTE VORTICITY
3 CONC 1.00 /M3 .00 CONCENTRATION
1 SFCC 1.00 /M3 .00 SURFACE CONCENTRATION
1 DEPS 1.00 /M2 .00 SURFACE DEPOSITION
1 DEPO 1.00 /M2 .00 SURFACE DEPOSITION
1 OPPT 1.00 MM .25 OBSERVED PRECIPITATION
1 SOLW 1.00 KG/M2 .00 0 TO 200 CM SOIL MOISTURE CONTENT
1 SOLT 1.00 DEGC .00 SOIL TEMPERATURE
1 SOLM 1.00 KG/M2 .00 SOIL MOISTURE
1 RH2M 1.00 PCT 20.00 2 METER RELATIVE HUMIDITY
3 P10M 1.00 K .00 POTENTIAL TEMPERATURE
3 WESD 1.00 KG/M2 .00 WATER EQUIV. OF ACCUM. SNOW DEPTH
3 SNOC 1.00 PCT 10.00 SNOW COVERAGE
3 CSNO 1.00 0--1 .50 CATEGORIAL SNOW (YES=1/NO=0)
3 CICE 1.00 0--1 .50 CATEGORIAL ICE PELLETS (YES=1/NO=0)
3 CFZR 1.00 0--1 .50 CATEGORIAL FREEZING RAIN (YES=1/NO=0)
3 CRAI 1.00 0--1 .50 CATEGORIAL RAIN (YES=1/NO=0)
3 PTYP 1.00 1--5 1.00 PRECIPITATION TYPE (RA=1,TRW=2,ZR=3,ICE=4,SN=5)
1 RGHS 1.00 M .00 SURFACE ROUGHNESS
3 LCLD 1.00 PCT 10.00 LOW CLOUD COVER
3 MCLD 1.00 PCT 10.00 MEDIUM CLOUD COVER
3 HCLD 1.00 PCT 10.00 HIGH CLOUD COVER
3 TCLD 1.00 PCT 20.00 AVERAGE TOTAL CLOUD COVER
3 CLDB 1.00 M 0.00 CLOUD BOTTOM HEIGHT
3 CLDT 1.00 M 0.00 CLOUD TOP HEIGHT
3 CWMR 1.00 G/KG 0.00 CLOUD WATER MIXING RATIO
3 WTMP 1.00 DEGC .00 WATER TEMPERATURE
3 LISD 1.00 DEGC 3.00 STANDARD LIFTED INDEX
3 LIB4 1.00 K 3.00 BEST 4-LAYER LIFTED INDEX
3 CAPE 1.00 J/KG .00 CONVECTIVE AVAILABLE POTENTIAL ENERGY
3 CINH 1.00 J/KG .00 CONVECTIVE INHIBITION
3 REFC 1.00 dBZ 5.00 COMPOSITE REFLECTIVITY
3 LTNG 1.00 .00 LIGHTNING STRIKES

关于空气质量来源和HYSPLIT人口

显示其他沿Traj变量

我知道如何访问/修改SETUP文件。我怎样才能使HYSPLIT将潜热/显热通量变量显示为沿轨迹变量?这些变量已经在我正在使用的NARR ARL文件中(LHTF,SHTF)。
On the other hand, you could also use the GUI under Meteorology/Display Data/Text Profile to look at the value of the fluxes at fixed latitude and longitudes.

沿轨迹降雨率

For precip, at hours 00, 06, 12, and 18z, the GDAS has the 6-four forecast from the previous cycle.
At hours 03, 09, 15, and 21z, the GDAS has the 3-hour forecast from the current cycle.

如何沿着hysplit轨迹输出涡度。

您可以先计算涡度(http://www.meteothink.org/examples/me ... hdivg.html),然后沿轨迹提取涡度数据(http://www.meteothinker.com/examples/me。 .. _data.html

# Set working directory
trajDir = 'D:/Temp/HYSPLIT'
meteoDir = r'U:\data\ARL\2015'

# Open trjactory data file
print 'Open trajectory data file ...'
trajfn = os.path.join(trajDir, 'traj_20150331')
print 'Trajectory file: ' + trajfn
trajf = addfile_hytraj(trajfn)

# Create trajectory layer
trajLayer = trajf.trajlayer()

# Open meteorological data file
print 'Open meteorological data file...'
meteofn = os.path.join(meteoDir, 'gdas1.mar15.w5')
print 'Meteorological file: ' + meteofn
meteof = addfile(meteofn)

# Get meteorological data along trajectory
print 'Get meteorological data along trajectory...'
outfn = os.path.join(trajDir, 'traj.txt')
outf = open(outfn, 'w')
outf.write('Lon,Lat,Time,Heigh,pres,VORT\n')
#uwnd = meteof['UWND'][:,:,:,:]
#vwnd = meteof['VWND'][:,:,:,:]
uwnd = meteof['UWND'][:,:,[10, 80],[50,140]]
vwnd = meteof['VWND'][:,:,[10, 80],[50,140]]
vort = hcurl(uwnd, vwnd)
idx = 0
for tline in trajLayer.shapes():    
    t = trajLayer.cellvalue('Date', idx)
    h = trajLayer.cellvalue('Hour', idx)
    t.replace(hour=h)
    for ps in tline.getPoints():  
        lon = ps.X
        lat = ps.Y 
        z = ps.M
        pres = ps.Z
        v = vort.interpn([t,pres,lat,lon])
        print 'lon: %.3f; lat: %.3f; time: %s; height: %.1f; pres: %.1f; VORT: %f' % (lon,lat,t.strftime('%Y%m%d_%H:%M'),z,pres,v)
        line = '%.3f,%.3f,%s,%.1f,%.1f,%f' % (lon,lat,t.strftime('%Y%m%d_%H:%M'),z,pres,v)
        outf.write(line + '\n')
        t = t + datetime.timedelta(hours=-1)
    idx += 1
outf.close()
print 'Finish...'

Wind speed along a trajectory

Generate trajectory on other map

METEOINFO

Get meteorological data along trajectory


. After all the changes have been made to the setup menu, press the Save button whichcopies all the information to a special file called default_traj and then closes the menu. Nowthe model simulation can be started from the Trajectory / Run model menu tab. This buttoncauses two actions. First, the trajectory setup information is copied from default_traj to a filecalled CONTROL, which is the input name required by the HYSPLIT trajectory model executable (hyts_std.exe). Second, the GUI then launches the model.

image.png

find it SETUP.CFG

How to get meteorology output along the trajectory on LINUX

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 199,711评论 5 468
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 83,932评论 2 376
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 146,770评论 0 330
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 53,799评论 1 271
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,697评论 5 359
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,069评论 1 276
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,535评论 3 390
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,200评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,353评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,290评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,331评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,020评论 3 315
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,610评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,694评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,927评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,330评论 2 346
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 41,904评论 2 341

推荐阅读更多精彩内容