Last change
on this file since 5293 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
|
Rev | Line | |
---|
[4774] | 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.