Planetary Boundary Layer Height

[Berry]

I need PBL height at a specific area in Eastern Europe with hourly resolution for a duration of a month.Is it possible to be done with Hysplit?
So far I managed to estimate it only for a certain hour of the day its time.
I used reanalysis meteorological data - "Trajectory" - "Total run time (hrs)" = 0 - "Simple Listing".
Thank you in advance.

[Berry]

Even though I waited for a hysplit advice, that worked nice for me.
If anyone has anything else to add as far as Hysplit is concered I would appreciate it. I wanted to something like that at the last slide but calculating PBL height for many more hours and days and with std (http://www.google.gr/url?sa=t&rct=j&q=& ... qE4kV0NpEg )
By the way, meteoinfo is pretty useful tool. Thank you mentioning it.

[gracy.hm]

thanks for such a nice forum. I also need to calculate the Boundary layer height(BLH) for my analysis I need to check the daily evolution of BLH also I want the average picture of BLH for month as well as season.
Can I get the MeteoInfo script for that.

[gracy.hm]

Thanks for such a nice forum, I also need boundary layer calculation, can I get the sample MeteoInfo script .

[hdelattre]

I’m interested to get Planetary Boundary Layer Height to a specific station located South of France (Latitude: 43.52N, Longitude: 4.64E) in July and August 2011. Following the previous discussions on the forum, I downloaded GDAS1 data and MeteoInfo software in order to extract PBLH data.

I can successfully perform One Dimension Plot but don’t manage to run the script “demo_pbl_station.py” to extract the data. It seems that the loop never ends and stays on July data. The control panel indicates in a continuous loop:

C:\Users\hdelattre\Documents\Doctorat\GlobalDataAssimilationSystem\gdas1.jul11.w1
0
C:\Users\hdelattre\Documents\Doctorat\GlobalDataAssimilationSystem\gdas1.jul11.w2
0
C:\Users\hdelattre\Documents\Doctorat\GlobalDataAssimilationSystem\gdas1.jul11.w3
0
C:\Users\hdelattre\Documents\Doctorat\GlobalDataAssimilationSystem\gdas1.jul11.w4
0
C:\Users\hdelattre\Documents\Doctorat\GlobalDataAssimilationSystem\gdas1.jul11.w5
0
C:\Users\hdelattre\Documents\Doctorat\GlobalDataAssimilationSystem\gdas1.jul11.w1
0
etc.

When I attempt to access to the software, it does not respond.

As it’s the first time I work with python, my problem is probably a coding issue. Below is the modified script I tried to run and here are the modifications I performed:
- Line 18, dataDir value
- Line 21, sTime value
- Line 22, eTime value
- Line 28, lon value
- Line 29, lat value
- Lines 40, 41, 55 and 57, I changed the indentation because the software returned the message “Syntax Error: unindent does not match any outer indentation level”

Can somebody help me ?

The modified script:
#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2012-11-29
# Purpose: Derive station PBL data from ARL meteorological data
# Note: Sample
#-----------------------------------------------------------
import clr
from System import *
from System.Globalization import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *
from MeteoInfoC.Data.MeteoData import *

import os.path

#---- Set directories
dataDir = "C:\\Users\\hdelattre\\Documents\\Doctorat\\GlobalDataAssimilationSystem\\"

#---- Set times - this case show seasonal variation
sTime = DateTime.Parse("2011-07-01")
eTime = DateTime.Parse("2011-08-31")
aTime = sTime

#---- Create a MeteoDataInfo object
mydata = MeteoDataInfo()
#---- Set lon/lat of the location
lon = 44.0
lat = 4.0
#---- Set output text file
outpath = dataDir + 'PBL_data.txt'
outf = open(outpath, 'w')
outf.write('Time,PBLH')
#---- Loop
i = 0
sum = 0.0
while aTime <= eTime:
for w in range(5):
inFile = dataDir + "gdas1." + aTime.ToString("MMM", CultureInfo("en-us")).ToLower() + aTime.ToString("yy") + ".w" + (w + 1).ToString()
print inFile
if os.path.isfile(inFile):
mydata.OpenARLData(inFile) #Open ARL data file
tNum = mydata.GetTimeNumber()
print tNum
for t in range(tNum):
dTime = mydata.GetTime(t)
print dTime.ToString("yyyy-MM-dd HH:00")
mydata.TimeIndex = t #Set time index
gData = mydata.GetGridData("PBLH") #Get PBLH grid data
pblh = gData.ToStation(lon, lat) #Intepolate PBLH data to the location
print 'PBLH = %10.2f' %pblh
outf.write(os.linesep)
outf.write(dTime.ToString("yyyyMMddHH") + ',%-7.2f ' %pblh)
sum += pblh
i += 1

aTime = aTime.AddMonths(1)

outf.close()
#---- Calculate seasonal average PBLH value
ave = sum / i
print 'Average PBLH: %10.2f' %ave

print 'Finished!'

[ariel.stein]

There is also a program called vmixing that is part of the linux distribution that uses a CONTROL and a SETUP.CFG file to read the meteorology and extract several diagnostic stability parameters.