! ! ! ! !  CONVERSION FACTORS  ! ! ! ! ! ! ! 
! convert from per second to per day
let mult = 60*60*24

! convert from per day to per year
let mult1 = mult*365

! convert nitrogen to carbon
let ntc = 6.625

! convert mol C to grams C
let mtg = 12.011
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

use tavg.02015.01.06.nc 

!Model Primary Production
let NPP1 = (o_diatnpp+o_phytnpp+O_diaznpp+o_coccnpp)*ntc*mtg*mult1 !gC per y
let cocc = if o_cocc ge 0 then o_cocc else 0
let phyt = if o_phyt ge 0 then o_phyt else 0
let diat = if o_diat ge 0 then o_diat else 0
let diaz = if o_diaz ge 0 then o_diaz else 0


go prepareFigure hovmoeller_seasonal
go "$HOME/opt_postprocessing_scripts/go/ilandscape5x1.jnl"
ppl AXLSZE 0.15 0.15
set view v1
sha/set/pal=purple_red/nolab/lev=(0,10,1)(inf) NPP1[k=1,x=@ave,l=1:73]
ppl shakey 1,0,0.15,5, , , , ,`($ppl$ylen)*(0.08)`,`($ppl$ylen)*(0.13)`
ppl txlsze 0.15
ppl shade
!con/ov/nolab/lev=(0,13,1)(inf) NPP1[k=1,x=@ave,l=1:73]
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "NPP (gC m^-^2 y^-^1)"

set view v2
sha/set/pal=purple_red/nolab/lev=(0,80,10)(inf)/axes=1,1,1,1 (diat[k=1,x=@ave,l=1:73])*ntc*mtg*1000
ppl shakey 1,0,0.15,5, , , , ,`($ppl$ylen)*(0.08)`,`($ppl$ylen)*(0.13)`
ppl shade
!con/ov/nolab/lev=10/axes=1,1,1,1 (diat[k=1,x=@ave,l=1:73])*ntc*mtg*1000
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "DT (mgC m^-^2)"

set view v3
sha/set/pal=purple_red/nolab/lev=(0,70,10)(inf)/axes=1,1,1,1 (cocc[k=1,x=@ave,l=1:73])*ntc*mtg*1000
ppl shakey 1,0,0.15,5, , , , ,`($ppl$ylen)*(0.08)`,`($ppl$ylen)*(0.13)`
ppl txlsze 0.15
ppl shade
!con/ov/nolab/lev=(0,60,10)(inf)/axes=1,1,1,1 (cocc[k=1,x=@ave,l=1:73])*ntc*mtg*1000
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "CP (mgC m^-^2)"

set view v4
sha/set/pal=purple_red/nolab/lev=(0,50,10)(inf)/axes=1,1,1,1 (diaz[k=1,x=@ave,l=1:73])*ntc*mtg*1000
ppl shakey 1,0,0.15,5, , , , ,`($ppl$ylen)*(0.08)`,`($ppl$ylen)*(0.13)`
ppl txlsze 0.15
ppl shade
!con/ov/nolab/lev=10/axes=1,1,1,1 (diaz[k=1,x=@ave,l=1:73])*ntc*mtg*1000
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "DZ (mgC m^-^2)"

set view v5
sha/set/pal=purple_red/nolab/lev=(0,5,0.5)(inf)/axes=1,1,1,1 (phyt[k=1,x=@ave,l=1:73])*ntc*mtg*1000
ppl shakey 1,0,0.15,5, , , , ,`($ppl$ylen)*(0.08)`,`($ppl$ylen)*(0.13)`
ppl txlsze 0.15
ppl shade
!con/ov/nolab/lev=10/axes=1,1,1,1 (phyt[k=1,x=@ave,l=1:73])*ntc*mtg*1000
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "LP (mgC m^-^2)"

go finalizeFigure_psthicken