TNC-NMFO / NWLAND

carbon accounting model
0 stars 0 forks source link

write_caland_inputs() | confirm vestigial nature of these forest NPP values for regions #131

Open sbassett opened 2 years ago

sbassett commented 2 years ago

this code is only referenced in a section marked if(FALSE) {} so it seems like it can be commented out or removed

line 500 of write_caland_inputs

### forest regional npp values, all biomass (hudiburg et al 2009, appendix Table A2-A)
# units converted to MgC/ha/yr
# these are used to get relative regional breakdown to apply to the statewide input values
# north coast is from coast range
# eastside is from east cascades
# klamath is from klamath mountains
# sierra cascades is from sierra nevada
# deserts is from central basin
# south coast is from oak-chapparal
# central coast is from oak-chapparal
# Delta and central valley use the oak-chapparal because it is very close to the unweighted average of 4.76
reg_names = c("C08001", "C08003", "C08005",   "C08007", "C08009", "C08011",   "C08013", "C08014", "C08015",   "C08017", "C08019", "C08021", "C08023",   "C08025", "C08027", "C08029",   "C08031", "C08033", "C08035",   "C08037", "C08039", "C08041",   "C08043", "C08045", "C08047",   "C08049", "C08051", "C08053",   "C08055", "C08057", "C08059",   "C08061", "C08063", "C08065",   "C08067", "C08069", "C08071", "C08073",   "C08075", "C08077", "C08079",   "C08081", "C08083", "C08085",   "C08087", "C08089", "C08091",   "C08093", "C08095", "C08097", "C08099",   "C08101", "C08103", "C08105",   "C08107", "C08109", "C08111",   "C08113", "C08115", "C08117",   "C08119", "C08121", "C08123",   "C08125", "C35001", "C35003",   "C35005", "C35006", "C35007",   "C35009", "C35011", "C35013",   "C35015", "C35017", "C35019", "C35021",   "C35023", "C35025", "C35027",   "C35028", "C35029", "C35031",   "C35033", "C35035", "C35037",   "C35039", "C35041", "C35043", "C35045",   "C35047", "C35049", "C35051",   "C35053", "C35055", "C35057",   "C35059", "C35061")
reg_vals = c(2.7665805940976,   2.19413160973837, 2.76903527728246, 2.63244596906997, 2.73156886762252,   1.96065837175619, 2.48106007005976, 2.83053659684482, 2.50674633597612,   1.47423783008082, 2.17096150901139, 2.22751021638909, 2.35203682882593,   1.85780717383626, 2.37986719643171, 2.76533621517, 2.85744668150817,   2.68309393038726, 2.59211034661613, 2.47142310893281, 2.35526522105354,   2.58431796022273, 2.68467417266612, 2.60565594532341, 2.22044797695949,   2.16685554590379, 2.0840697968623, 2.10787609620547, 2.49165205065151,   1.9908450075498, 2.43398197009325, 1.65489526347314, 2.20635498322028,   2.01240348242764, 2.66346811441964, 2.42558636352397, 2.58752231967444,   1.39397801950559, 2.56120779924163, 2.73712517620926, 2.2270407002044,   2.55307919087443, 2.79548332072972, 2.81695358991931, 2.58784567162072,   1.98898701165137, 2.56885700422033, 2.21435714550648, 2.66865093751068,   2.23145460700705, 2.26391729179792, 2.44760661358592, 2.387488171171,   2.25899628552496, 2.20372980676083, 2.1553245025352, 2.10926918670168,   2.56945181569684, 1.9235928955523, 2.166836311221, 2.07863550401721,   2.57364838844398, 2.51801211172068, 2.3560541659734, 2.78747919811024,   2.97898795685335, 1.80974260382829, 3.022560474804, 2.16823688171128,   1.69593465713667, 1.4927280655417, 2.14845195977066, 2.44623622668831,   2.83133071142975, 1.68112506172512, 2.36289995081388, 2.21303591253923,   1.62845566428421, 2.88239113368413, 2.91395751850666, 2.38291015846412,   2.94560135363299, 2.47079253138599, 2.48706571203913, 1.34950548411691,   2.72613407232389, 1.42701414649999, 2.77707486008997, 2.68901161939305,   2.74181181089144, 2.78211152711749, 2.67879484012917, 2.57130741240519,   2.52834981330927, 2.6229441239555, 2.66147093432118, 2.61573418211572)
forest_npp = data.frame(Region=reg_names, npp=reg_vals)
forest_npp$Land_Type = "Forest"