fc2k_power_from_array2.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                     !
  !                dy   add arbitrary array size          !
  !                                                       !
  !-------------------------------------------------------!

    SUBROUTINE fc2k_power_from_array2(POWER_AR,  TIME_AR, TIME, POWER)                      

    USE itm_constants


    IMPLICIT NONE


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

    INTEGER                    :: i
    integer                    :: np,nt
    real (r8)                  :: teps

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

   teps=1.0e-7 !0.1 mks
   np=size(power_ar)
   nt=size(time_ar)

   if (np.ne.nt) then
      write(*,*) 'PowerFromArray: size of power and time arrays are notequal'
      stop
   end if

   if (np.eq.0) then
      write(*,*) 'PowerFromArray: no elements in power array'
      stop
    end if
    power=0.0
    DO i = 1, np
       IF (time_ar(i).LE.time .AND. time.LE.time_ar(i+1)) THEN
          if (abs(time_ar(i)-time_ar(i+1)).lt.teps) then
             power = power_ar(i)
             write (*,*) 'time values are equal taking power at that point'
          else
             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 
          end if
          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_array2