fc2k_power_from_array.F90

Go to the documentation of this file.

! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  

  !-------------------------------------------------------!
  !     This routine selects the power from the time      !
  !     sequence                                          !
  !-------------------------------------------------------!
  !     Source:       ---                                 !
  !     Developers:   D.Kalupin                           !
  !     Kontacts:     Denis.Kalupin@efda.org              !
  !                                                       !
  !     Comments:     created for ETS                     !
  !                                                       !
  !-------------------------------------------------------!

    SUBROUTINE fc2k_power_from_array(POWER_AR,  TIME_AR, TIME, POWER)                      

    USE itm_constants


    IMPLICIT NONE


! +++ i/o parameters:
    REAL (R8), INTENT (IN)     :: power_ar(1000)      
    REAL (R8), INTENT (IN)     :: time_ar(1000)     
    REAL (R8), INTENT (OUT)    :: power      
    REAL (R8), INTENT (IN)     :: time

    INTEGER                    :: i

!+++++++++++++++++++++++++++++++++++++++++

    DO i = 1, 999
       IF (time_ar(i).LE.time .AND. time.LE.time_ar(i+1)) THEN
          power = power_ar(i) + (power_ar(i+1)- power_ar(i))/(time_ar(i+1)-time_ar(i))*(time-time_ar(i))
          !write(*,*) 'power from NBI', POWER, ' for time', time 
          goto 10
       END IF

       IF (time_ar(i).GT.time_ar(i+1)) THEN
          WRITE (*,*)'NONMONOTHONIC TIME HAS BEEN SPECIFIED FOR POWER SEQUENCE'
          stop
       END IF
    END DO


 10 RETURN


    END SUBROUTINE fc2k_power_from_array