unit of emission rate

Topics about the HYSPLIT dispersion model.
shanky
Posts: 14
Joined: February 14th, 2015, 12:59 am
Registered HYSPLIT User: No

Re: unit of emission rate

Post by shanky »

Hi Ariel,
Sorry for the delayed reply. I am still not able to get the desired results.
I had used the following steps to run my trajectory and dispersion model.
Actually, I have concatenate the WRF-ARW files using NCO utility into a single wrf output file. The concatenated single wrf output file contains around 13 days data with timestep (00,03,06,09,12,15,18,21) per day. I had used ncview to check the single wrf output file ..the newly created file seems fine.
Then, I had used the arw2arl utility to convert this single wrf output file into ARL format and the newly created filename is wrf_new_d01.
While converting the file it shows the following message but the conversion status message says finished and successful.
"Using an existing decoding configuration: WRFDATA.CFG
Using an existing encoding configuration: ARLDATA.CFG
NOTICE pakset:
Number of index records = 1
Number of records /time = 265
NetCDF: Start+count exceeds dimension bound

NOTICE pakini: start initialization rec = 1
NOTICE pakini: end initialization rec = 265
NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

Initialized accumulated precip from file: WRFRAIN.BIN
Accumulated precip file not from the same simulation
Completed: 13 6 19 0 0

"
I have checked the converted output file using the Meteorology and check file utility of HYSPLIT... the output window shows the ARL file contents.

Then I'd set up the model for the trajectory backtracing and checked the trajectory output ....but the result was not as expected :( :( :( .......

I've used the HYSPLIT online facility and I set up the trajectory model with the same settings I'd used HYSPLIT on Yellowstone(UCAR/NCAR supercomputer). I'd run the online HYSPLIT utility using the archived datasets. The online HYSPLIT run showed the desired results but theHYSPLIT model I'd installed on Yellowstone does not give me the desired results.

Could you please help me in this. I am not sure what I am missing here on yellowstone.

Thank You
shanky
Posts: 14
Joined: February 14th, 2015, 12:59 am
Registered HYSPLIT User: No

Re: unit of emission rate

Post by shanky »

Hi Ariel,
I am attaching herewith my arw2arl conversion shell script along with the shell script output file (MESSAGE).

Thank You
ariel.stein
Posts: 660
Joined: November 7th, 2012, 3:14 pm
Registered HYSPLIT User: Yes

Re: unit of emission rate

Post by ariel.stein »

I did not get the attachment. Can you please try again?
Could you also send me a sample of your WRF output?
shanky
Posts: 14
Joined: February 14th, 2015, 12:59 am
Registered HYSPLIT User: No

Re: unit of emission rate

Post by shanky »

HI Ariel,
There is some problem to upload the files... So that I am pasting my WRF to ARL conversion shell script
#!/bin/bash

# Where the program is located

PGM="/glade/p/work/kurkute/HYS/trunk/exec"

### Location of the WRF data you would like to convert ###
WRFDIR="/glade/p/nhap0001/chenl/wrf_chenl/WRFV3/run_hgt_ctl/wrfout_ctl_10-23"

### Location that you would like to output HYSPLIT data to ###
HYDIR="/glade/p/work/kurkute/HYS/trunk/working.782"

### Name you would like to give the HYSPLIT data ###
HYDATA="wrf_hysplit_d01"



### Convert all WRF data in WRFDIR to Hysplit format:

#ln -sf ${WRFDIR}/wrfout_d01* .

ln -sf ${WRFDIR}/wrfout_d01_2013-06-{15..20}* .

FILES=$(ls -1 wrfout_d01*)

for FILE in ${FILES}
do

rm -f ARLDATA.BIN
${PGM}/arw2arl -i${FILE}
cat ARLDATA.BIN >> ${HYDIR}/${HYDATA}

done

rm -f ARLDATA.BIN

rm -f wrfout_d01*



Currently, I am uploading the original wrf output file which I need to convert into ARL format on the dropbox ..I am providing herewith the dropbox link of original WRF-ARW output file and the output of my arw2arl shell script link respt.

https://www.dropbox.com/s/wfvy1i2wxjjc8 ... toarl?dl=0
https://www.dropbox.com/s/c1wy0g79mlpcb ... %3A00?dl=0

Thank You
ariel.stein
Posts: 660
Joined: November 7th, 2012, 3:14 pm
Registered HYSPLIT User: Yes

Re: unit of emission rate

Post by ariel.stein »

Well, I was able to convert the file to the ARL format. Please try to do the conversion of the files one by one and then concatenate them.
shanky
Posts: 14
Joined: February 14th, 2015, 12:59 am
Registered HYSPLIT User: No

Re: unit of emission rate

Post by shanky »

Hi Ariel,
thank you for your reply. I have tried to convert the single WRF-ARW file into the ARL but still it throws the following message
Using an existing decoding configuration: WRFDATA.CFG
Using an existing encoding configuration: ARLDATA.CFG
NOTICE pakset:
Number of index records = 1
Number of records /time = 265
NetCDF: Start+count exceeds dimension bound

NOTICE pakini: start initialization rec = 1
NOTICE pakini: end initialization rec = 265
NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

NetCDF: Start+count exceeds dimension bound

Forecast time greater than zero, but no previous precip file found!
Completed: 13 6 19 3 0

I am not sure possibly what is going wrong with it ?

Thanking Yours,
ariel.stein
Posts: 660
Joined: November 7th, 2012, 3:14 pm
Registered HYSPLIT User: Yes

Re: unit of emission rate

Post by ariel.stein »

I am not sure what is wrong but everything points to your NetCDF libraries.
Do you have any other machine your can work with?
shanky
Posts: 14
Joined: February 14th, 2015, 12:59 am
Registered HYSPLIT User: No

Re: unit of emission rate

Post by shanky »

Hi Ariel,
I have complied HYSPLIT using Netcdf latest version and it works now. It does not give an error and the new results are satisfactory.

Thanking Yours,
Post Reply

Return to “Dispersion Model”