ivosuter
ivosuter
1) If you have a switch 'lneutral', it's easy to set the wall temperature equal to the fluid temperature 2) there is no real solution: maybe set the Ri or...
If you want it to be 'neutral' skip the whole wfuno routine and set Cd directly based on wall roughness length? From subroutine unom it seems that in this case...
hmm. It seems it's still quite a bit of coding. So, I would recommend to make a wfmneutral: copy wf_uno.f90, rename delete cases 'x2' and function 'unoh' & 'unom' delete...
yes usually its 'ctm=unom(...)', followed by 'dummy=u^2*ctm'; sometimes its directly 'dummy=u^2*unom()'. it's pretty much the same thing. wherever there is a unom() you can just replace it with 'fkar2/(logdz**2)' _"We...
I said I would have a last look at it, but that we believe it's done. So gonna check wfmneutral now
Hmm it appears to me there are indeed some mistakes in both, wfm_neutral and wf_uno in the edge cases. On the staggered grid the tangential velocities on a surface are...
ui-1 ui ui+1 ui+2 o-----o-----o-----o |...................................| |...................................| |...................................| ui-1 and ui+2 only experience half a cell
wfm_neutral: do we need w bottom cases at all? since blocks are never floating the bottom is always the floor? If we do, shouldn't it be (on lines 142,251,352 and...
Hello, it is currently not possible to read in large scale forcing at given time-steps and the code would need to be adapted. Currently there is also only soil-moisture in...