Last change
on this file since 3425 was
3415,
checked in by jgipsl, 7 years ago
|
Added script to create wood harvest forcing files for ORCHIDEE. Done by N Vuichard, LSCE
|
File size:
760 bytes
|
Line | |
---|
1 | use "/home/basc/dnarayan/LandUsemap_ORCHIDEE/LUH2v2h/staticData_quarterdeg.nc" |
---|
2 | use "/home/basc/dnarayan/LandUsemap_ORCHIDEE/LUH2v2h/states.nc" |
---|
3 | use "/home/basc/dnarayan/LandUsemap_ORCHIDEE/LUH2v2h/transitions.nc" |
---|
4 | |
---|
5 | ! Bioh |
---|
6 | ! Total harvest wood biomass |
---|
7 | ! Conversion from kg to g |
---|
8 | let bioh=(primf_bioh+primn_bioh+secmf_bioh+secyf_bioh+secnf_bioh)*1000 |
---|
9 | |
---|
10 | ! Land |
---|
11 | ! Total land fraction over which wood harvest acts |
---|
12 | let land=primf+primn+secdf+secdn |
---|
13 | |
---|
14 | |
---|
15 | ! Total Grid-cell Area |
---|
16 | ! Conversion from km2 to m2 |
---|
17 | let area=carea*1e6 |
---|
18 | |
---|
19 | ! Define harvest wood biomass per m2 of forest |
---|
20 | let woodharvest=if(land[d=2]) then bioh[d=3]/(land[d=2]*area[d=1]) else 0. |
---|
21 | |
---|
22 | |
---|
23 | |
---|
24 | |
---|
25 | repeat/l=1001:1165 save/file="/home/orchidee03/vuichard/WOODHARVEST_LUH2_v2h/woodharvest_`l+849`.nc" woodharvest[l=`l`] |
---|
26 | |
---|
Note: See
TracBrowser
for help on using the repository browser.