Skip to content

Implement crown fire flame length in kitral #119

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
matildeRivas wants to merge 3 commits into
base: main
Choose a base branch
from
Draft

Implement crown fire flame length in kitral #119

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
ffdcc94
Read CBD and CBH from raster layers in Kitral
matildeRivas Jan 27, 2025
882a18e
Calculate crown intensity and flamen length in kitral
matildeRivas Jan 28, 2025
a5251fc
[bug] Use real tree height
matildeRivas Jan 30, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 6 additions & 6 deletions Cell2Fire/Cell2Fire.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -682,7 +682,7 @@ Cell2Fire::reset(int rnumber, double rnumber2, int simExt = 1)
this->surfaceIntensityFolder = this->args.OutFolder + "SurfaceIntensity" + separator();
}
// Crown Byram Intensity Folder
if ((this->args.OutIntensity) && (this->args.AllowCROS) && (this->args.Simulator == "S"))
if ((this->args.OutIntensity) && (this->args.AllowCROS) && (this->args.Simulator != "C"))
{
CSVWriter CSVFolder("", "");
this->crownIntensityFolder = this->args.OutFolder + "CrownIntensity";
Expand All @@ -698,15 +698,15 @@ Cell2Fire::reset(int rnumber, double rnumber2, int simExt = 1)
this->surfaceFlameLengthFolder = this->args.OutFolder + "SurfaceFlameLength" + separator();
}
// Crown Flame Length Folder
if ((this->args.OutFl) && (this->args.AllowCROS) && (this->args.Simulator == "S"))
if ((this->args.OutFl) && (this->args.AllowCROS) && (this->args.Simulator != "C"))
{
CSVWriter CSVFolder("", "");
this->crownFlameLengthFolder = this->args.OutFolder + "CrownFlameLength";
CSVFolder.MakeDir(this->crownFlameLengthFolder);
this->crownFlameLengthFolder = this->args.OutFolder + "CrownFlameLength" + separator();
}
// max Flame Length Folder
if ((this->args.OutFl) && (this->args.AllowCROS) && (this->args.Simulator == "S"))
if ((this->args.OutFl) && (this->args.AllowCROS) && (this->args.Simulator != "C"))
{
CSVWriter CSVFolder("", "");
this->maxFlameLengthFolder = this->args.OutFolder + "MaxFlameLength";
Expand Down Expand Up @@ -1963,7 +1963,7 @@ Cell2Fire::Results()
}

// Crown Intensity
if ((this->args.OutIntensity) && (this->args.AllowCROS) && (this->args.Simulator == "S"))
if ((this->args.OutIntensity) && (this->args.AllowCROS) && (this->args.Simulator != "C"))
{
this->crownIntensityFolder = this->args.OutFolder + "CrownIntensity" + separator();
std::string intensityName;
Expand Down Expand Up @@ -1999,7 +1999,7 @@ Cell2Fire::Results()
}

// Crown Flame length
if ((this->args.OutFl) && (this->args.AllowCROS) && (this->args.Simulator == "S"))
if ((this->args.OutFl) && (this->args.AllowCROS) && (this->args.Simulator != "C"))
{
this->crownFlameLengthFolder = this->args.OutFolder + "CrownFlameLength" + separator();
std::string fileName;
Expand All @@ -2017,7 +2017,7 @@ Cell2Fire::Results()
}

// Max Flame length
if ((this->args.OutFl) && (this->args.AllowCROS) && (this->args.Simulator == "S"))
if ((this->args.OutFl) && (this->args.AllowCROS) && (this->args.Simulator != "C"))
{
this->maxFlameLengthFolder = this->args.OutFolder + "MaxFlameLength" + separator();
std::string fileName;
Expand Down
2 changes: 1 addition & 1 deletion Cell2Fire/Cells.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -736,7 +736,7 @@ Cells::manageFire(int period,
surfFraction[nb] = metrics.sfc;
SurfaceFlameLengths[this->realId - 1] = mainstruct.fl;
SurfaceFlameLengths[nb - 1] = metrics.fl;
if ((args->AllowCROS) && (args->Simulator == "S"))
if ((args->AllowCROS) && (args->Simulator != "C"))
{
float comp_zero = 0;
MaxFlameLengths[this->realId - 1]
Expand Down
125 changes: 115 additions & 10 deletions Cell2Fire/FuelModelKitral.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -680,13 +680,79 @@ byram_intensity(inputs* data, main_outs* at)
return ib; // unidad de medida
}

/**
* Calculates byram fire intensity when there is active crown fire.
* In order for this to be calculated, the input folder must contain
* files with CBD, CBH and tree height data for each cell.
* @param at Structure containing the cell's output data.
* @param data Structure containing the cell's input data.
* @return Fire intensity.
*/

float
crown_byram_intensity_k(main_outs* at, inputs* data)
{
float cbd, cbh;
if (std::isnan(data->cbh))
{
cbh = cbhs[data->nftype][0];
}
else
{
cbh = data->cbh;
}
// CBH is 0,1012 * Height^(1,4822)
float canopy_height = std::pow(cbh / 0.1012, 1 / 1.4822) - cbh;
if (canopy_height < 0)
{
throw std::runtime_error("Tree height is lower than canopy base height, "
"please provide valid files.");
}
if (std::isnan(data->cbd))
{
cbd = cbds[data->nftype][0];
}
else
{
cbd = data->cbd;
}
return std::ceil((hs[data->nftype][0] / 60) * cbd * canopy_height * at->ros_active * 100.0) / 100.0;
}

/**
* @brief Calculates the flame length of a cell when there is crown fire.
* @param intensity Byram intensity for crown fires
* @return the flame length
*/

float
crown_flame_length_k(float intensity)
{
float fl = 0.1 * pow(intensity, 0.5);
if (fl < 0.01)
{
return 0;
}
else
{
return std::ceil(fl * 100.0) / 100.0;
}
}

bool
fire_type(inputs* data, main_outs* at, int FMC)
{
float intensity, critical_intensity, cbh;
bool crownFire = false;
intensity = at->sfi;
cbh = cbhs[data->nftype][0];
if (std::isnan(data->cbh))
{
cbh = cbhs[data->nftype][0];
}
else
{
cbh = data->cbh;
}
// cbh = data->cbh;
critical_intensity = pow((0.01 * cbh * (460 + 25.9 * FMC)), 1.5);
if ((intensity > critical_intensity) && cbh != 0)
Expand All @@ -698,11 +764,26 @@ float
crownfractionburn(inputs* data, main_outs* at, int FMC)
{ // generar output de cfb
float a, cbd, ros, ros0, H, wa, i0, cbh, cfb;
cbh = cbhs[data->nftype][0];
if (std::isnan(data->cbh))
{
cbh = cbhs[data->nftype][0];
}
else
{
cbh = data->cbh;
}
i0 = pow((0.01 * cbh * (460 + 25.9 * FMC)), 1.5);
H = hs[data->nftype][0];
wa = fls_david[data->nftype][0];
cbd = cbds[data->nftype][0];
if (std::isnan(data->cbd))
{
cbd = cbds[data->nftype][0];
}
else
{
cbd = data->cbd;
}

ros0 = 60 * i0 / (H * wa);
ros = at->rss;
if (cbd != 0)
Expand Down Expand Up @@ -767,15 +848,26 @@ checkActive(inputs* data, main_outs* at, int FMC) // En KITRAL SE USA PL04
{
float ros_critical, cbd, H, wa, i0, cbh;
bool active;
cbh = cbhs[data->nftype][0];
if (std::isnan(data->cbh))
{
cbh = cbhs[data->nftype][0];
}
else
{
cbh = data->cbh;
}
i0 = pow((0.01 * cbh * (460 + 25.9 * FMC)), 1.5);
H = hs[data->nftype][0];
wa = fls_david[data->nftype][0];
cbd = cbds[data->nftype][0];
ros_critical = 60 * i0 / (H * wa);

cbd = cbds[data->nftype][0];

if (std::isnan(data->cbd))
{
cbd = cbds[data->nftype][0];
}
else
{
cbd = data->cbd;
}
active = cbd * ros_critical >= 3;
return active;
}
Expand Down Expand Up @@ -828,7 +920,14 @@ calculate_k(inputs* data,
FMC = args->FMC;
ptr->nftype = data->nftype;
ptr->fmc = fmcs[data->nftype][0];
ptr->cbh = cbhs[data->nftype][0];
if (isnan(data->cbh))
{
ptr->cbh = cbhs[data->nftype][0];
}
else
{
ptr->cbh = data->cbh;
}
// cout << " cbh " << ptr->cbh << "\n";

ptr->fl = fls_david[data->nftype][0];
Expand Down Expand Up @@ -870,7 +969,7 @@ calculate_k(inputs* data,

// Step 10: Criterion for Crown Fire Initiation (no init if user does not
// want to include it)
if (args->AllowCROS && cbhs[data->nftype][0] != 0)
if (args->AllowCROS && (data->cbh > 0 || cbhs[data->nftype][0] != 0))
{
if (activeCrown)
{ // si el fuego esta activo en copas chequeamos condiciones
Expand Down Expand Up @@ -904,6 +1003,8 @@ calculate_k(inputs* data,
at->rss = hptr->ros;
bptr->ros = backfire_ros10_k(hptr, sec);
fptr->ros = flankfire_ros_k(hptr->ros, bptr->ros, sec->lb);
at->crown_intensity = crown_byram_intensity_k(at, data);
at->crown_flame_length = crown_flame_length_k(at->crown_intensity);

if (args->verbose)
{
Expand All @@ -926,6 +1027,8 @@ calculate_k(inputs* data,
at->rss = hptr->ros;
bptr->ros = backfire_ros10_k(hptr, sec);
fptr->ros = flankfire_ros_k(hptr->ros, bptr->ros, sec->lb);
at->crown_intensity = crown_byram_intensity_k(at, data);
at->crown_flame_length = crown_flame_length_k(at->crown_intensity);

if (args->verbose)
{
Expand Down Expand Up @@ -1012,6 +1115,8 @@ determine_destiny_metrics_k(inputs* data, fuel_coefs* ptr, arguments* args, main
if (crownFire)
{
metrics->cfb = crownfractionburn(data, metrics, FMC);
metrics->crown_intensity = crown_byram_intensity_k(metrics, data);
metrics->crown_flame_length = crown_flame_length_k(metrics->crown_intensity);
}
if (args->verbose)
{
Expand Down
6 changes: 6 additions & 0 deletions Cell2Fire/FuelModelKitral.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -32,6 +32,9 @@ float backfire_ros_k(main_outs* at, snd_outs* sec);
// Flame length [m])
float flame_length(inputs* data, main_outs* at);

// Crown Flame length [m])
float crown_flame_length_k(float intensity);

// Angle of the flame w.r.t. horizontal surface (Putnam's)
float angleFL(inputs* data, main_outs* at);

Expand All @@ -41,6 +44,9 @@ float flame_height(inputs* data, main_outs* at);
// byram intensity
float byram_intensity(inputs* data, main_outs* at);

// crown byram intensity
float crown_byram_intensity_k(main_outs* at, inputs* data);

// Type of fire (Crown = 1, SROS = 0)
bool fire_type(inputs* data, main_outs* atr, int FMC);
// Active crown
Expand Down
5 changes: 5 additions & 0 deletions Cell2Fire/FuelModelSpain.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,6 +35,8 @@ float backfire_ros_s(main_outs* at, snd_outs* sec);

// Flame length [m])
float flame_length(inputs* data, fuel_coefs* ptr);
// Crown Flame length [m])
float crown_flame_length(float intensity);

// Angle of the flame w.r.t. horizontal surface (Putnam's)
float angleFL(inputs* data, fuel_coefs* ptr);
Expand All @@ -45,6 +47,9 @@ float flame_height(inputs* data, fuel_coefs* ptr);
// byram intensity
float byram_intensity(inputs* data, fuel_coefs* ptr);

// crown byram intensity
float crown_byram_intensity(main_outs* at, inputs* data);

// Type of fire (Crown = 1, SROS = 0)
bool fire_type(inputs* data, fuel_coefs* ptr);

Expand Down