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