sign change in mean flow control

This commit is contained in:
ignis 2014-07-15 14:23:52 +09:00
parent fc92088b1a
commit 949fbe102d

View file

@ -289,11 +289,12 @@ module m_fdm_calc
write(St_data,'(f10.5,10e25.15)')fdmtime,dummyu_,sum_wrate,fl_location/(REAL(nx-1)*hx)*100.
! Control the inflow mean velocity, dummyu_
write(*,'(a20,2f12.8)') ' ** oldsumc, sumc = ',oldsumc,sumc
if(swadtv.ne.0.and.mod((fdmcyc+1),swadtv).eq.0) then
time_int=fdmtime-time_int
delu=(oldsumc-sumc)/time_int
delu=(sumc-oldsumc)/time_int
if ((fdmtime).gt.startad.and.oldsumc.ne.0.) then
dummyu_=dummyu_-delu
dummyu_=dummyu_+delu
endif
time_int=fdmtime
endif