|
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.
