Menu
▾
▴
#3 Advection of fluff and initialisation, getm_bio.f90
I am currently using ERGOM (9 compartments + fluff), thus the total
number of compartments numc=10. In getm_bio.f90 the do_advection routine
runs from i=1,numc. Therefore, fluff is also advected, which is wrong. I
had to change the code to i=1,numc-1. This only works, because I know I
am using a system with fluff.
Currently I dont know if there is a nice work around, instead of hard
coding.
Further, during initialisation of the compartments, the fluff array is
3d and fully read in. Because, only the first layer is used, the rest of
the array can be set to 0, which is currently not the case. This does
not change the biology, but might lead to some confusion during
postprocessing.
