QUIC-Fire-TT / ttrs_quicfire

ttrs_quicfire is a Python library to easily configure burn models for plots of land defined using shape files for the quicfire model.
MIT License
0 stars 1 forks source link

Fuel Topography and Parabolic dz Arrays #7

Closed cbonesteel closed 2 years ago

cbonesteel commented 2 years ago

Issues #2 and #3

Description

This branch implements the use of topography when gathering fast fuels data for input into the quicfire model. It makes the default build use topography, with topography not being used only when specified. It generates a topo.dat and topo.inp file with the appropriate data for the model.

This branch also implements a variable parabolic dz array that is built using either the fuel_height + relief * 3 where relief is the difference between the tallest and lowest point of the topography or fuel_height + 100. The later is the default, as when topography is not used, this is the better option since relief is 0.

Motivation and Context

These changes were made to improve the input files accuracy for any given plot of land since before the model was only running for flatlands. While this is a simple way to represent the burn plot, topography will greatly improve the accuracy of the runs that are generated.

Test Cases

For the topography generation, the outputs for the topo.dat were compared to that of the flatlands. The data was being generated and stored in the data file.

For the parabolic dz array, the following test case was used. It was compared to the original matlab script and was similar in results.

build_parabolic_dz_array(nz=22, Lz=350, n_surf=4, dz_surf=2.5)

Checklist