Commit 59eccd23 authored by Alan Wai Hou Lio's avatar Alan Wai Hou Lio
Browse files

licoreim wp3 bug fix

parent 41ce0133
Pipeline #22468 passed with stage
in 4 minutes and 26 seconds
...@@ -804,16 +804,15 @@ subroutine update_regulation(array1, array2) bind(c,name='update_regulation') ...@@ -804,16 +804,15 @@ subroutine update_regulation(array1, array2) bind(c,name='update_regulation')
TTfa_vel = array1(16) TTfa_vel = array1(16)
! !
! LICOREIM ! LICOREIM
Licovar%gain = 1.00_mk
if ((time .gt. Licovar%time_detected).and.(time .le.(Licovar%time_detected+Licovar%x0+Licovar%sigma*4))) then if ((time .gt. Licovar%time_detected).and.(time .le.(Licovar%time_detected+Licovar%x0+Licovar%sigma*4))) then
! if ((time .gt. 100.00_mk) .and. (time .le. 110.00_mk)) then ! if ((time .gt. 100.00_mk) .and. (time .le. 110.00_mk)) then
write(6,*) "*** Re-writing the RPM gain at time = ", time, " s ***"
Licovar%gain = 1.0000_mk+ Licovar%amp*exp( - ((time - (Licovar%time_detected+Licovar%x0))**2/(2*Licovar%sigma**2))**Licovar%p) Licovar%gain = 1.0000_mk+ Licovar%amp*exp( - ((time - (Licovar%time_detected+Licovar%x0))**2/(2*Licovar%sigma**2))**Licovar%p)
else
Licovar%gain = 1.000000000_mk
endif endif
GenSpeed_unmod = GenSpeed GenSpeed_unmod = GenSpeed
if (time .gt. 5.0_mk) then if ((time .gt. 5.0_mk) .and. (Licovar%gain.gt.1.00_mk)) then
GenSpeed = GenSpeed_unmod*Licovar%gain GenSpeed = GenSpeed_unmod*Licovar%gain
write(6,*) "*** Re-writing the RPM gain at time = ", time, " s ***"
endif endif
!*********************************************************************************************** !***********************************************************************************************
! Safety system ! Safety system
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment