Ticket #400: genere_woodharvest.jnl

File genere_woodharvest.jnl, 760 bytes (added by nvuilsce, 4 years ago)

script for generating the wood harvest data

Line 
1use "/home/basc/dnarayan/LandUsemap_ORCHIDEE/LUH2v2h/staticData_quarterdeg.nc"
2use "/home/basc/dnarayan/LandUsemap_ORCHIDEE/LUH2v2h/states.nc"
3use "/home/basc/dnarayan/LandUsemap_ORCHIDEE/LUH2v2h/transitions.nc"
4
5! Bioh
6! Total harvest wood biomass
7! Conversion from kg to g
8let bioh=(primf_bioh+primn_bioh+secmf_bioh+secyf_bioh+secnf_bioh)*1000
9
10! Land
11! Total land fraction over which wood harvest acts
12let land=primf+primn+secdf+secdn
13
14
15! Total Grid-cell Area
16! Conversion from km2 to m2
17let area=carea*1e6
18
19! Define harvest wood biomass per m2 of forest
20let woodharvest=if(land[d=2]) then bioh[d=3]/(land[d=2]*area[d=1]) else 0.
21
22
23
24
25repeat/l=1001:1165 save/file="/home/orchidee03/vuichard/WOODHARVEST_LUH2_v2h/woodharvest_`l+849`.nc" woodharvest[l=`l`]
26