uataq / stilt

Stochastic Time-Inverted Lagrangian Transport model

Stochastic Time-Inverted Lagrangian Transport model

An open source Lagrangian particle dispersion model.

Docs

STILT documentation
Methods details

About

STILT would not be possible without the strong community of developers behind it. This distribution contains a completely redesigned STILT wrapper and proposes a centralized, collaborative platform for documentation and future development. Model development in the form of feature enhancements, documentation updates, bug fixes, or simple suggestions from the community are welcome. Contribution guidelines can be found here.

Relevant manuscripts

Loughner, C. P., Fasoli, B., Stein, A. F., Lin, J. C.: Incorporating features from the Stochastic Time-Inverted Lagrangian Transport (STILT) model into the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model: a unified dispersion model for time-forward and time-reversed applications, J. Appl. Meteorol. Climatol., 10.1175/JAMC-D-20-0158.1, 2021.

Fasoli, B., Lin, J. C., Bowling, D. R., Mitchell, L., and Mendoza, D.: Simulating atmospheric tracer concentrations for spatially distributed receptors: updates to the Stochastic Time-Inverted Lagrangian Transport model's R interface (STILT-R version 2), Geosci. Model Dev., 10.5194/gmd-11-2813-2018, 2018.

Stein, A. R., Draxler, R. R., Rolph, G. D., Stunder, B. J. B., and Cohen M. D.: NOAA’s HYSPLIT atmospheric transport and dispersion modeling system. Bull. Amer. Meteor. Soc., 10.1175/BAMS-D-14-00110.1, 2015.

Lin, J. C., Gerbig, C., Wofsy, S. C., Andrews, A. E., Daube, B. C., Davis, K. J. and Grainger, C. A.: A near-field tool for simulating the upstream influence of atmospheric observations: The Stochastic Time-Inverted Lagrangian Transport (STILT) model, J. Geophys. Res., 10.1029/2002JD003161, 2003.

stilt's Issues

problems with wrf-stilt

Hi,

I am running the STILT tutorial with wrf data (the wrfout in the tutorial folder) and I am getting the following error:
Using an existing decoding configuration: WRFDATA.CFG
NetCDF: Variable not found
ZNT
Variable not found: ZNT
Some defined variables not found in file.

Thank you for your help!
Best regards, Cristian Arango

Alternative Projections for STILT footprints

It would be nice to be able to save footprints in a variety of projections. My suggestion would be to use proj4 strings to define the projection. In the footprint configuration file, we can set the domain by defining the projection, the lower-left and upper-right corners, and the x and y resolution.

from foot_config.r:

#Footprint projection.
#Use proj4 convention ('http://proj4.org/faq.html') use '+proj=latlon' for latlon grid.
#projection <- '+proj=latlon'
projection <- '+proj=aea  +lat_0=48.128 +lon_0=11.575 +lat_1=40 +lat_2=60 +x_0=0 +y_0=0'
#Lower left and upper right corner lat and lon
ll_lon <- -73 #Degrees East
ur_lon <- -69.5
ll_lat <- 41.5   #Degrees North
ur_lat <- 43.5

#Resolution...match units with CRS! latlon: degrees, most other projections: meters
xres <- 1000
yres <- 1000

Then, in the code that calls the footprint maker, you convert all of the partcle lat/lon pairs to x/y in the projection before calling create_footprint. Note you will have to change the gaussian smoother to remove the latitude correction.

The trickiest part is getting the netCDF files CF compliant. It won't recognize the proj4 string, so in order to get it so that Panoply will correctly project it I made a clunky look up table (projection_dictionary.r) It also requires a bunch of other attributes.

I've attached my version of the footprint related files. Let me know how I can most effectively contribute. I don't want to make you sift through a bunch of peoples' poorly commented code!

footprint_maker.zip

stilt project

Hi, sorry to bother you again, I would like to ask you, can I change the projection in stilt, my wrfout file is Lambert projection and the projection in stilt is long-lat, can I change the projection in stilt?
And I tried this problem, I used proj4 and changed the projection method to projection <- '+proj=laea +lat_0=30 +lon_0=106', unfortunately, I got Error : cannot allocate vector of size 611239.8 Gb error I'd like to ask if there is any way to solve it
I seem to have found the cause of the error, may I ask if it is set to lambert, is the unit of xres in meters?

Thank you for your time and have a nice life!

STILT Model Applications on different source sizes and heights

In the application of the model for a period of time, it is found that the STILT model is suitable for point source pollution simulation. On the model introduction page, I can see the STILT model is a variant of the HYSPLIT model. However, the HYSPLIT model supports point, line, and area sources, as well as setting the scale (size) of pollution sources. If the STILT model is to support three types of pollution source settings, please tell me how to set it? At the same time, I would like to know if the model can set and simulate pollution diffusion results under different source sizes and vertical height layers, and how? Thank you!

Could not find 'bindrcpp'

Hi, complete newbie when it comes to R. Ran:

if (!require('devtools')) install.packages('devtools')
devtools::install_github('benfasoli/uataq')

and after installing a slew of packages, I got this:

• installing source package ‘uataq’ ...
  ** R
  ** byte-compile and prepare package for lazy loading
  Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) :
  there is no package called ‘bindrcpp’

Could be user error on my end. Also don't see a package named 'parallel' from a few of the CRAN source sites I looked at in the US.

Any help is appreciated with getting this up and running. Trying to run some models to determine PM dispersion to help argue against a factory they are trying to build near by.

-Thanks

Metrological Grid Nesting does not work

Hi the hymodelc binary built for Centos does not seem to allow for meteorological grid nesting. Instead of transitioning to using the meteorology from a lower resolution file at the boundary of a high resolution file the trajectories stop at the high resolution boundary. This feature is available in the hymodelc binary in the AER's STILT distribution and seems to be in Hysplit. In those cases the trajectories transition to using the lower resolution file at the boundary of the high resolution one. Is this a bug with UATAQ's stilt distribution?

STILT for lower PBL??

Dear Sir, I am trying to visualize the effect of local transport to the receptor site.
The observation station at our receptor site is located at very high altitude (hilly area) such a way that boundary layer lies very close to surface or sometimes lies below the observation station.

My doubt is in such a case, can I use STILT to understand the influence of local transport to the receptor site, if yes what parameters i need to be taken care of i run_stilt.r file.
What value of PBL i should use?? Will STILT be able to provide some significant results?

OR i should run WRF and then use the output for STILT??

Thank you in advance.

about emission inventory convolution

Hello!
I would like to ask you two questions
1、The emission inventory of wbb tutorial is a .rds file, is the file you are downloading an rds file, I am downloading an nc file, do I need to convert it to an rds file? How do I need to convert it, and how do I convert it please?
2. Is the concentration change value obtained after convolving the footprint file with the emission inventory the value of the latitude and longitude of the center of each grid?

Thank you for your reply and your time!
Have a nice life!

initialise tracer concentration

I have a query regarding replacing footprint units (ppm/umol/m2/s) with particle concentration (ug/m3).
I want to use the concentrations of a specie in (ug/m3) as a tracer concentration for the simulations. So that in the final result i would be able to visualize the particle concentration w.r.t to its trajectory. just like in following figure. such that I can insert the concentration (ug/m3) [measured from instruments] at the receptor into the STILT.

Is it feasible? What are the changes I need to make in the rstilt.r or which variable needs to be modified for this?

Bugs Report

Hi developers,

When I learn your WBB tutorial, all is OK, except for the html out step.
It hints the following information.

Loading required namespace: ncdf4
Warning message:
data_frame() is deprecated, use tibble().
This warning is displayed once per session.
Error in htmlwidgets::saveWidget(p, "timeseries.html") :
Saving a widget with selfcontained = TRUE requires pandoc. For details see:
https://github.com/rstudio/rmarkdown/blob/master/PANDOC.md
Execution halted

Regards,

ddlddl58

plotting trajectories using .netcdf outputs

Hii sir
Thanks for the wonderful tutorials you have uploaded. I was able to install and run the model successfully. I have got the outputs. I am interested to plot the trajectories (i.e. only path trajectory) using these outputs like as follows:

But since it has footprints i am getting an output like this

So, can you please help me how i can be able to plot trajectories only using STILT outputs.

stilt

HI!
1、 Can I use stilt mode to simulate the concentration between altitude layers, for example, simulating the concentration at 40~50m in a site area?
2、Does the stilt model provide a calculation of the background CO2 concentration?

Other meteorological data to ARL

Hi,
I'm struggling to convert ERA5 met data to ARL format. I have tried grib2arl (https://www.ready.noaa.gov/HYSPLIT_data2arl.php) and era52arl (https://github.com/amcz/hysplit_metdata), but neither of them have worked for me. Are you planning to offer support on converting ERA5 files to ARL, or are there any converters you recommend?

Discrepancies Between Self-Generated and Tutorial-Provided Time Series Graphs in Stationary Simulation Tutorials

The graph output did not appear as expected, as shown in the tutorial. Could you please help me identify potential causes for this discrepancy or how to solve this problem? I have checked every step but am still unable to find the reason. Thanks for your help.

Wbb tutorial.r Warning Messages

I am trying to run the first tutorial: and I received the following warning messages when I run the tutorial.r in the directory: 01-wbb:

Warning messages:
1: In showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj = prefer_proj) :
Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs +type=crs
2: In showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj = prefer_proj) :
Discarded datum World Geodetic System 1984 in Proj4 definition
Warning message:
In sum(new("RasterLayer", file = new(".RasterFile", name = "", datanotation = "FLT4S", :
Nothing to summarize if you provide a single RasterLayer; see cellStats

I was wondering if anyone was familiar with these warning messages and how I might fix them. I am still able to produce the timeseries.html, average_footprint.html, and average_contribution.html, but when I try to open them with the links or lynx command, I get a black page.

Using the STILT model for the monitoring of volcanic eruptions

Hello,

I am trying to adapt the available tutorials to model the volcanic eruptions and its pollutants dispersion in the atmosphere. I’ve studied the STILT model and worked the associated tutorials and I’m now trying to adapt them to simulate an emission of sulfur dioxide from the Etna and its dispersion.
I would like to use STILT with other inputs to investigate the variation of SO2 in an area centered on the volcano which will be the only point of the grid with a non-zero emission value. As STILT is used to show the backward trajectories, I planned to replace in the script run_stilt.r the file receptors.rds (the file available for the light-rail tutorial) with another raster file according to the same type (data.frame) that will be a (latitude, longitude) grid centered on the volcano. For the emission.rds file, I planned to recreate a RasterBrick file like the one in the tutorials, also with a grid (4 columns : lat, lon, SO2 emission value, time) centered on the volcano, but, as I previously said, with the whole "SO2 emission value" column equal to zero, except the line corresponding to the point (lat, lon) of the volcano. By doing that, I hope that the backward trajectories from each receptor point will illustrate the dispersion of the SO2 after being emitted by the volcano.

So, I have a couple questions concerning these two input files, in order to recreate them and be able to restart STILT for my simulation.
I tried to recreate a data.frame identical to the one of the receptors.rds file, but when I start the run_stilt.r script, I only obtain output files in the "out/by-id" directory. The "out/footprints" remains empty, and I don't understand why since when I run the script with the file of the light-rail tutorial, everything looks fine. I noticed that the footprints directory is completed by the simulation_step.r script, but I don't understand the difference between the files that disrupt STILT.
Concerning the emissions.rds file, which is a pretty much bigger file with a lot of information, I also tried to recreate its structure but I can't test it for now, because I need to run run_stilt.r correctly before I run the tutorial.r script (that need a emissions.rds file), but in cases where it doesn't work, I would like to have, if this is possible, more information about its creation, and the origin of the data used to complete it.

To finish with, I wonder about the ability of STILT to work on the « forward » mode, that can appear to be a choice maybe more intuitive in the context of my study around the volcanic eruptions.

If you have any idea concerning this adaptation of the STILT tutorials, I would be grateful for your help. I can of course precise anything that wasn't comprehensible enough in my explanations, or share my adapted scripts and inputs.

Thank you for your time

Maximum time period for Backward-forward trajectory

Dear sir,
I was wondering how long (maximum number of hours) can we run STILT model for backward - forward trajectory simulation? i.e. if i give value of n_hours for one month (24*30), will the model be able to provide good results.
And please give your valuable comments on the issue i created earlier https://github.com/uataq/stilt-tutorials/issues/4#issue-853197631. I am frequently getting this error for different simulations even after re-installing the all repositories.
Thanks.

Predicting receptor concentrations from convolution of footprint with emissions inventories

Hi,

I'm interested in comparing point sensor methane measurements to local emissions inventories by scaling the inventories with STILT-R footprints. This is motivated by the WBB Carbon Dioxide Github example where observed enhancements were estimated as:

dCO2 = sum(values(foot * emissions_subset), na.rm = T))

A series of questions have arisen while trying to apply this approach to my work.

First, the tutorial mentions that this approach "makes some overly-basic assumptions" - could you elaborate on what these assumptions are and when they do not hold?

Second, what is the proper way to change the simulation's spatial resolution when using the summation approach? For instance, the WBB example is run at 0.002 deg^2 (~222 m^2) resolution. Footprint and emissions units are ppm m^2 s umol^-1 and umol m^-2 s^-1 respectively. A reasonable dCO2 enhancement is calculated with this combination of spatial resolutions/units. However, the summation approach also means that changing the spatial resolution dramatically changes the estimated enhancement (more/fewer points increases/decreases the estimated enhancement). Do the footprint and emissions units need to be rescaled relative to the 0.002 deg^2 values? In other words, is there something "special" about this resolution which makes the footprint and emissions inventories give reasonable results?

Third, is there a reason to perform a summation instead of an integration? In principle these provide the same area under the curve, but in practice I've found integration to yield erroneous results.

Thank you for your help and your time,
Griffin

Incomplete Trajectory for higher altitudes

Hi sir
I have run STILT at zagl <- 4000 meter altitude and 120 hour backward mode. And in the output i cant see a complete trajectory of footprints to the receptor (RED DOT). Is there any explanation of this incomplete trajectory for higher altitudes? or is it something like model incapability for higher altitudes? Whereas we can see a complete trajectory of footprints for lower altitudes i.e. 100 meter in my run. Please clear my doubt.

4000 meter altitude output

100 Meter altitude output

Questions in using STILT model

Hello, professors：
I used the WRF output to driven the STILT model, I have 4 questions about it.

1. First question about the projection, the projection of the WRF output is lambert_conformal_conic, and I found the run_stilt.r projection <- '+proj=longlat', so, should the projection of the run_stilt.r must be consistent with the WRF output?
2. I found the spatial resolution in the run_stilt.r is degree, but in WRF output is m, the spatial resolution must consistent with the WRF output？

3. Receptor location(s)

lati <- 40.766189
long <- -111.847672
zagl <- 25
the receptor location only can set one in once simulation？
How it can set multiple locations in once simulation? Just like
lati <- 40.766189 38.0 …..
long <- -111.847672 -112.0 …
zagl <- 25 24 …
i got an error about it when set 2 locations
[1] "2021-04-01 12:00:00 UTC" "2021-04-01 13:00:00 UTC"
Error: unexpected numeric constant in "lati <- 43.88 44.0"
Execution halted
4. how to generate the emissions.rds file？

Running new STILT script

I am trying to run my own R script with STILT that I have created by modifying the tutorial R scripts. However, when I attempt the R script, I get the following error:
[ckbekker@login1 STILT_runs]$ Rscript run_stilt_mariposa_jul28.R Initializing STILT Number of receptors: 24 Number of parallel threads: 8 Error in load_libs("dplyr", "ncdf4", "parallel", "raster", "rslurm", lib.loc = lib.loc) : could not find function "load_libs" Calls: source -> withVisible -> eval -> eval Execution halted
Do you know what this error might be related to or how I can figure out?
Thank you.

run_stilt_mariposa_jul28.txt

./hycs_std >> stilt.log 2>&1' timed out after 3600s

What is the meaning of following errors?

emission.rds file

Hello, sorry to bother you, is the script for creating the in your tutorial.r there? I am using this program to read my emission nc file and I have the following problem, can you please help me， thank you for your precious time, thank you very much.
Uploading codefile.rar…

1. non-numeric argument to binary operator
2. 'subset' must be logical
   I can't upload my code file, can you contact me?
   thank you very much.

Error linking Homebrew gfortran installation in hycs_std (?)

When I try to run a tutorial example it returns the following shell abort message and R warning for each receptor. No output files (footprints etc) are produced.

sh: line 1:  3156 Abort trap: 6           ./hycs_std >> stilt.log 2>&1
Warning message:
In calc_trajectory(namelist, rundir, emisshrs, hnf_plume, met_files,  :
  Failed to output  /Users/ethan/Desktop/research/phd/em27/wbb-tutorial/out/by-id/201512100000_-111.847672_40.766189_25/PARTICLE_STILT.DAT

I have not been able to follow the source of the error.

I am running on Mac OS 10.14.16, R version 4.1.2, gfortran from Homebrew gcc 11.2.0_1.

I hope this is an appropriate venue to ask about this.

Thanks!

UPDATE----

I found the stilt.log file, and it said:

dyld: Library not loaded: /usr/local/gfortran/lib/libgfortran.5.dylib
  Referenced from: /Users/ethan/Desktop/research/phd/em27/wbb-tutorial/out/by-id/201512100000_-111.847672_40.766189_25/./hycs_std
  Reason: no suitable image found.  Did find:
	/usr/local/lib/libgfortran.5.dylib: stat() failed with errno=62
Failed to output  /Users/ethan/Desktop/research/phd/em27/wbb-tutorial/out/by-id/201512100000_-111.847672_40.766189_25/PARTICLE_STILT.DAT 

This seems to have to do with an earlier issue I had while running stilt_init(), which would result in this error when compiling:

Compiling footprint kernel aggregation subroutine...
gfortran -mmacosx-version-min=10.13 -fno-optimize-sibling-calls  -fPIC  -Wall -g -O2  -c  r/src/permute.f90 -o r/src/permute.o
clang -mmacosx-version-min=10.13 -dynamiclib -Wl,-headerpad_max_install_names -undefined dynamic_lookup -single_module -multiply_defined suppress -L/Library/Frameworks/R.framework/Resources/lib -L/usr/local/lib -o r/src/permute.so r/src/permute.o -L/usr/local/gfortran/lib/gcc/x86_64-apple-darwin18/8.2.0 -L/usr/local/gfortran/lib -lgfortran -lquadmath -lm -F/Library/Frameworks/R.framework/.. -framework R -Wl,-framework -Wl,CoreFoundation
ld: warning: directory not found for option '-L/usr/local/gfortran/lib/gcc/x86_64-apple-darwin18/8.2.0'
ld: warning: directory not found for option '-L/usr/local/gfortran/lib'
ld: library not found for -lgfortran
clang: error: linker command failed with exit code 1 (use -v to see invocation)
make: *** [r/src/permute.so] Error 1

I fixed this issue with the solution at this stackexchange answer about another R package (Rccp) having issues linking Homebrew installations of gfortran:
https://stackoverflow.com/questions/29992066/rcpp-warning-directory-not-found-for-option-l-usr-local-cellar-gfortran-4-8/29993906#29993906

I fixed the original error from hycs_std by linking libgfortran.X.dylib to where it says it wants it in stilt.log

I also had to link the library libquadmath.X.dylib

Everything seems to be working now!

It seems there is a widespread hard time to be had finding Homebrew installations of gfortran.

dmass correction

Hi Ben,

Excuse my ignorance here, but I've been comparing the wrappers from the original version of STILT (i.e. http://stilt-model.org) to yours, and I wondered if you could explain one of the differences for me. In the original wrappers there is an option to correct for violation of mass conservation in the met fields (dmassTF in Trajecfoot.r), which defaults to TRUE (at least in my copy). I couldn't find an equivalent option in your wrappers - have I missed it or did you decide not to implement this correction?

Cheers,
Joe

would like ability to have separate installation and output directories

We would like to keep our STILT installation on /home2, but write output to our /wrk directory. Right now, the /out directory is assumed to be under the STILT project directory, so it is hard to separate the two.

In order to get around this, we are currently using the following "hacks". (Please note that Stilt_Utah installed in /home2/smg4.)

r/dependencies.r
line 7, added: stilt_wd2 <- '/home2/smg4/stilt_Utah'
line 10, modified: rsc <- dir(file.path(stilt_wd2, 'r', 'src'), pattern = '.*\.r$', full.names = T)
line 20, modified: permute_exe <- file.path(stilt_wd2, 'r/src/permute.so')

r/src/simulation_step.r
line 44, added: stilt_wd2 <- '/home2/smg4/stilt_Utah'
line 45, modified: source(file.path(stilt_wd2, 'r/dependencies.r'), local = T)
line 81, modified: link_files <- dir(file.path(stilt_wd2, 'exe'))
line 82, modified: file.symlink(file.path(file.path(stilt_wd2, 'exe'), link_files), rundir)

r/run_stilt.r
line 10, set stilt_wd <- file.path('/wrk/smg4', project)
lines 153-154: setwd(stilt_wd)
source('/home2/smg4/stilt_Utah/r/dependencies.r')

proj_as_proj_string: Cannot find proj.db

Why this error is coming though i am getting output netcdf file?
proj_as_proj_string: Cannot find proj.db
proj_as_wkt: Cannot find proj.db
proj_create: Cannot find proj.db
proj_as_wkt: Cannot find proj.db

& Also Please tell the role of number of receptors. How should i decide that it should be 1 or 10 or more?
Thanks.

Error while using multiple nodes

Dear Benfasoli

Can you please have a look at the following error? I am currently trying to run the script on the HPC server, so I am putting multiple nodes for parallel tasks.
i.e. n_cores <- 6
n_nodes <- 2
I have just started working on HPC, so if you also provide any additional comments from your side on handling STILT codes on High Computing systems would be very beneficial for me.

Forward mode results in "PARTICLE_STILT.DAT does not contain any trajectory data" while backward mode works

Hello,

When I try to run STILT in forward mode with HRRR ARL's from ftp://ftp.arl.noaa.gov/archives/hrrr/ I have issues with forward in time simulations but not backward in time simulations. I have double checked that the necessary met files are in the met directory. The issue seems to arise when trying to cross HRRR file boundaries. If run_time is at 10:00 UTC, it exhibits the failing behavior. But if run_time is at 12:00 UTC, it works. This behavior is not entirely consistent however, and I cannot figure out what controls it. gfs0p25 ARL's seem to work fine.

Thanks!

Backward in time:

n_hours   <- -6 

Results in successful trajectory calculations:

R
-----------------------------------------
Initializing STILT
Commit ID: fc4daceee2f608650040688213fdb53554d37b62
Number of receptors: 1
Number of parallel workers: 1
Parallelization disabled. Executing simulations sequentially...
Running simulation ID:   202208061000_-103.5_31.5_10
>

Trying to go forward with the same meteorology, receptors, only changing n_hours:

n_hours   <- 6 

results in:

R
-----------------------------------------
Initializing STILT
Commit ID: fc4daceee2f608650040688213fdb53554d37b62
Number of receptors: 1
Number of parallel workers: 1
Parallelization disabled. Executing simulations sequentially...
Running simulation ID:   202208061000_-103.5_31.5_10
Warning message:
In calc_trajectory(namelist, rundir, emisshrs, hnf_plume, met_files,  
/n/home00/emanninen/apps/stilt/myproject/out/by-id/202208061000_-103.5_31.5_10
PARTICLE_STILT.DAT does not contain any trajectory data.

Examining stilt.log reveals:

stilt.log
--------------------------------------------
STOP 900
 HYSPLIT - Initialization
 HYSPLIT version: hysplit.v5.1.0
 Last Changed Date: 2021-05-13
  NOTICE   main: using namelist file - SETUP.CFG
  NOTICE   main: STILT emulation (mixing ratio, vary layer, mass sum)
 Calculation Started ... please be patient
 *ERROR* metpos: start point not within (x,y,t) any data file
  - start time before start of meteorology data
/n/home00/emanninen/apps/stilt/myproject/out/by-id/202208061000_-103.5_31.5_10/PARTICLE_STILT.DAT does not contain any trajectory data.

about emission inventory convolution

Hello!
I would like to ask you two questions
1、The emission inventory of wbb tutorial is a .rds file, is the file you are downloading an rds file, I am downloading an nc file, do I need to convert it to an rds file? How do I need to convert it, and how do I convert it please?
2. Is the concentration change value obtained after convolving the footprint file with the emission inventory the value of the latitude and longitude of the center of each grid?

Thank you for your reply and your time!
Have a nice life!

f951: error: unrecognized command-line option '-quiet'

When I type, Rscript -e "uataq::stilt_init('myproject')"
I got this,

fatal: destination path 'myproject' already exists and is not an empty directory.
Compiling footprint kernel aggregation subroutine...
/opt/R/arm64/bin/gfortran -mtune=native -fno-optimize-sibling-calls -fPIC -Wall -g -O2 -c r/src/permute.f90 -o r/src/permute.o
f951: error: unrecognized command-line option '-quiet'
make: *** [r/src/permute.o] Error 1

How I can solve it?

following problem when running multiple receptors

Hello, I have the following problem when running multiple receptors: Error in r_run_time + layers * interval : non-numeric argument to binary operator.
Here is my code，
receptors <- data.frame(
run_time=c("2015-12-10 00:00:00","2015-12-10 03:00:00"),
lati = c(40.766189,40.81235),
long = c(-111.847672,-111.754683),
zagl = c(25,25),
stringsAsFactors=FALSE
)

stilt

Can I use stilt mode to simulate the concentration between altitude layers, for example, simulating the concentration at 40~50m in a site area?

STILT Multiple receptors

Hello!
I need to simulate two receptor sites in stilt, is this feasible? How should I do it?
Thank you!

missing permute.so

Hi, when I run stilt on mac system, it shows me this error:

source('r/dependencies.r')
Error in eval(ei, envir) :
calc_footprint(): failed to find permute.so in r/src/

Then I got a permute.so file from others and put it under r/src/ folder, it shows me this error:

source('r/dependencies.r')
Error in dyn.load(permute_exe) :
unable to load shared object '/Users/mac/Documents/Research/stilt/r/src/permute.so':
dlopen(/Users/mac/Documents/Research/stilt/r/src/permute.so, 0x0006): tried: '/Users/mac/Documents/Research/stilt/r/src/permute.so' (not a mach-o file), '/System/Volumes/Preboot/Cryptexes/OS/Users/mac/Documents/Research/stilt/r/src/permute.so' (no such file), '/Users/mac/Documents/Research/stilt/r/src/permute.so' (not a mach-o file)

Do you have any suggestions?