rich-iannone / puffr Goto Github PK

Dear Richard,
I'm trying to use your excellent system PuffR for an odour exposure simulation near to a waste treatment plant in Italy (43.126324, 12.409735).
Following your instructions, that I found on slideshare, I see that your example for 'vancouver' works fine, whereas when I try to set my parameters I'm not able to run the function (I stopped at the calmet_define_geophys function, see https://dl.dropboxusercontent.com/u/74598336/puffR_snap1.PNG).
Can you help me, please, and tell me where I'm wrong?
Moreover, I would use if possible my own files for DEM, land cover and meteorological parameters (you can find them on https://dl.dropboxusercontent.com/u/74598336/input_data.zip, if you want).
Could you explain me how to correctly set the package functions to do that?
I thank you in advance,
best regards,
Paolo Rossini
Rimini, Italy

Hi @rich-iannone Thanks for your package. I need some help, I have a question.

I have compiled Calpuff and mmif, hence Calmet it is not necessary because I create the Calmet output from WRF using mmif. Now I only need to create the emission inputs files. According to calpuff.inp, the emission inputs files are:

Emission Files

PTEMARB.DAT input * PTDAT = *
VOLEMARB.DAT input * VOLDAT = *
BAEMARB.DAT input * ARDAT = *
LNEMARB.DAT input * LNDAT = *

How do I generate this files? I've navigated into the calpuff.for file but I still cant find the proper Fortran format.

Does PuffR generate this emission files?

I've used your function to generate emissions files, calpuff_add_line_sources, but it generated only one line. I tried also calpuff_add_area_sources, but I don't know if I use it correctly. Anyways, is the output these functions, the information required for the emission files *.DAT in Calpuff?

Thanks

Hi, rich-iannone, thank you for your great PuffR. I'm intrested in calpuff recently and I really want to use it in my Rstudio, However, there are troubles all the time. I want to send this question to your for help.
The "srtm_UTM_resampled.df” in the calmet_define_geophys founction is defined to c("x", "y", "z"), but from my results is c("z", "x", "y"). I'm not sure of it, is it right?
And the ggsave for "srtm_UTM_resampled.df" terrain--location_name-...... made a mistake in my R. The error is "Error in dev(filename = filename, width = dim[1], height = dim[2], ...) :
unused argument (filename = filename)". Thus, I deleted the "device = pdf, width = 8, height = 8, units = "in"", it became ok. But I still don't know the reason of it.
Because I haven't finish the first step of calmet_define_geophys, I'm worried about the next step for the project. I really want to get your help. Thank you very much.