Open robertladwig opened 5 years ago
Tried the different settling options using GLM 3.0.5 and AED 1.3.2. Still, using Stokes Law (settling =3 ) for particulate organic matter relies on sedimentation velocity and doesn't recognize changes in the density and diameter.
POP with settling = 1 (constant) and velocity set to 0
POP with settling = 3 (Stokes) and velocity set to 0
POP with settling = 1 (constant) and velocity set to 1
POP with settling = 3 (Stokes) and velocity set to 1
POP with settling = 3, velocity = 1 and density to 2000, so it should sink like a stone
POP with settling = 3, velocity = 1 and density to 0.1, so it should float
Hey @robertladwig Can you plot OGM_Psed_pop? In fact OGM_Psed_pop/OGM_pop should be vvel in m/s
@robertladwig
Can you also make suremobility_off = .false.
in glm.nml &wq_setup
though I doubt this is wrong as otherwise changing velocity wouldn't work
@matthipsey mobility_off is set to .false. And OGM_Psed_pop seems to be null in my simulation. Does the Psed-Variable only include advective transport?
Currently running GLM 3.0.0 with AED2 1.3.0 and tested this on a Mac and a Windows machine. It seems that the settling configuration in the aed2_mobility subroutines is not functioning as intended. Changing the settling configuration won't affect the model outcome, as it seems that GLM is always using the DEFAULT CASE with vvel = data%w_pom and vvel_cpom = data%w_cpom (instead of for instance Stoke's Law).
My guess is that the keywords from aed2.h aren't read in properly, as AED2 doesn't recognize the settling/migration keyword.