!! $HeadURL$
!! $Date$
!! $Revision$

use "/home/basc/dnarayan/LandUsemap_ORCHIDEE/LUH2v2h/staticData_quarterdeg.nc"
use "/home/basc/dnarayan/LandUsemap_ORCHIDEE/LUH2v2h/states.nc"
use "/home/basc/dnarayan/LandUsemap_ORCHIDEE/LUH2v2h/transitions.nc"

! Bioh
! Total harvest wood biomass
! Conversion from kg to g
let bioh=(primf_bioh+primn_bioh+secmf_bioh+secyf_bioh+secnf_bioh)*1000

! Land
! Total land fraction over which wood harvest acts
let land=primf+primn+secdf+secdn


! Total Grid-cell Area 
! Conversion from km2 to m2
let area=carea*1e6

! Define harvest wood biomass per m2 of forest 
let woodharvest=if(land[d=2]) then bioh[d=3]/(land[d=2]*area[d=1]) else 0.




repeat/l=1001:1165 save/file="/home/orchidee03/vuichard/WOODHARVEST_LUH2_v2h/woodharvest_`l+849`.nc" woodharvest[l=`l`]