search

NorESMhub / CTSM

Community Terrestrial Systems Model (includes the Community Land Model of CESM)
http://www.cesm.ucar.edu/models/cesm2.0/land/
Other
2 stars 9 forks source link

Create fire emission parameters in FATES parameter file #47

Closed rosiealice closed 3 months ago

rosiealice commented 3 months ago

Need to add

https://github.com/NGEET/fates/blob/main/parameter_files/fates_params_default.cd

The default CTSM emissions factors live here, on glade. /glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/firedata/fire_emis_factors_c140116.nc

rosiealice commented 3 months ago

Noting for the sake of creating the emission factors, that it is necessary to map from the existing set of emission factors, indexed by PFT, into the FATES PFT space

I made the mapping in the following way, listing the FATES PFT (number, name), the HLM PFT (number, name) and then the 13 emission factors. This mapping did not require any subjective choices as where HLM PFTs needed to be combined (e.g. temperate and boreal EGT), their emission factors were already the same.

1 "broadleaf_evergreen_tropical_tree ", 4 broadleaf_evergreen_tropical_tree 1580, 103.2, 6.8, 8.1, 3.8, 1.85, 0.2, 9.1, 8.5, 6, 5.2, 0.63, 0.5\ 7, 2 "needleleaf_evergreen_extratrop_tree ", 1 needleleaf_evergreen_temperate_tree 1569, 106.7, 4.7, 5.7, 1.81, 3, 0.26, 13, 17.6, 8.3, 9.1, 0.56,\ 1, 3 "needleleaf_colddecid_extratrop_tree ", 3 needleleaf_deciduous_boreal_tree 1569, 106.7, 4.7, 5.7, 1.81, 3, 0.26, 13, 17.6, 8.3, 9.1, 0.56, 1,

4 "broadleaf_evergreen_extratrop_tree ", 5 broadleaf_evergreen_temperate_tree 1569, 106.7, 4.7, 5.7, 1.81, 3, 0.26, 13, 17.6, 8.3, 9.1, 0.56, 1\ , 5 "broadleaf_hydrodecid_tropical_tree ", 6 broadleaf_deciduous_tropical_tree 1580, 103.2, 6.8, 8.1, 3.8, 1.85, 0.2, 9.1, 8.5, 6, 5.2, 0.63, 0.5\ 7, 6 "broadleaf_colddecid_extratrop_tree ", 7 broadleaf_deciduous_temperate_tree 1569, 106.7, 4.7, 5.7, 1.81, 3, 0.26, 13, 17.6, 8.3, 9.1, 0.56, 1\ , 7 "broadleaf_evergreen_extratrop_shrub ", 9 broadleaf_evergreen_shrub 1663, 61.6, 2.2, 3.4, 0.99, 2.35, 0.21, 4.9, 8.5, 3.7, 3.2, 0.46, 0.71,

8 "broadleaf_hydrodecid_extratrop_shrub ",10 broadleaf_deciduous_temperate_shrub 1663, 61.6, 2.2, 3.4, 0.99, 2.35, 0.21, 4.9, 8.5, 3.7, 3.2, 0.4\ 6, 0.71, 9 "broadleaf_colddecid_extratrop_shrub ",10 broadleaf_deciduous_temperate_shrub 1663, 61.6, 2.2, 3.4, 0.99, 2.35, 0.21, 4.9, 8.5, 3.7, 3.2, 0.4\ 6, 0.71, 10 "arctic_c3_grass ", 12 c3_arctic_grass 1663, 61.6, 2.2, 3.4, 0.99, 2.35, 0.21, 4.9, 8.5, 3.7, 3.2, 0.46, 0.71,

11 "cool_c3_grass ", 13 c3_non-arctic_grass 1663, 61.6, 2.2, 3.4, 0.99, 2.35, 0.21, 4.9, 8.5, 3.7, 3.2, 0.46, 0.71,

12 "c4_grass " ; 14 c4_grass 1663, 61.6, 2.2, 3.4, 0.99, 2.35, 0.21, 4.9, 8.5, 3.7, 3.2, 0.46, 0.71,

rosiealice commented 3 months ago

For the injection heights, we needed to map the hardwired code from CTSM into FATES.

The code looks like this: ! temperate and boreal evergreen needleleaf trees else if ( veg_type == ndllf_evr_tmp_tree .or. veg_type == ndllf_evr_brl_tree ) then ztop = 4.e3_r8 ! m ! temperate and boreal trees else if ( veg_type == ndllf_dcd_brl_tree .or. veg_type == nbrdlf_evr_tmp_tree .or. & veg_type == nbrdlf_dcd_tmp_tree .or. veg_type == nbrdlf_dcd_brl_tree ) then ztop = 3.e3_r8 ! m ! tropical broadleaf trees (evergreen and decidious) else if ( veg_type == nbrdlf_evr_trp_tree .or. veg_type == nbrdlf_dcd_trp_tree ) then ztop = 2.5e3_r8 ! m ! shrubs else if ( veg_type >= nbrdlf_evr_shrub .and. veg_type <= nbrdlf_dcd_brl_shrub ) then ztop = 2.e3_r8 ! m ! grasses else if ( veg_type >= nc3_arctic_grass .and. veg_type <= nc3_nonarctic_grass ) then ztop = 1.e3_r8 ! m ! generic unmanaged crops else if ( veg_type == nc3crop .or. veg_type <= nc3irrig ) then ztop = 1.e3_r8 ! m ! Prognostic crops else if ( veg_type >= npcropmin .and. veg_type <= npcropmax ) then ztop = 1.e3_r8 ! m else

and we end up with this:

"broadleaf_evergreen_tropical_tree ", 2500 "needleleaf_evergreen_extratrop_tree ", 4000 "needleleaf_colddecid_extratrop_tree " , 3000 "broadleaf_evergreen_extratrop_tree ", 3000 "broadleaf_hydrodecid_tropical_tree ", 2500 "broadleaf_colddecid_extratrop_tree " , 3000 "broadleaf_evergreen_extratrop_shrub ", 2000 "broadleaf_hydrodecid_extratrop_shrub ", 2000 "broadleaf_colddecid_extratrop_shrub " , 2000 "arctic_c3_grass ", 1000 "cool_c3_grass ", 1000 "c4_grass ", 1000 ;

rosiealice commented 3 months ago

These are now introduced here:

https://github.com/NorESMhub/fates/blob/e3d231af66bb5c51f05b5ffaa27fd2b3e04c30c2/parameter_files/fates_params_default.cdl#L269