fc2k_power_modulation.F90

Go to the documentation of this file.

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

  !-------------------------------------------------------!
  !     This routine provides power array and time array  !
  !     for power modulation                              !
  !-------------------------------------------------------!
  !     Source:       ---                                 !
  !     Developers:   D.Kalupin                           !
  !     Kontacts:     Denis.Kalupin@efda.org              !
  !                                                       !
  !     Comments:     created for ETS                     !
  !                                                       !
  !-------------------------------------------------------!

    SUBROUTINE fc2k_power_modulation(POWER_IN,  NLEVEL, TSTART, TEND, FREQUENCY, TAU, POWER_OUT, TIMES_OUT)                      

    USE itm_constants


    IMPLICIT NONE


! +++ i/o parameters:
    REAL (R8), INTENT (IN)     :: power_in(1000)      
    REAL (R8), INTENT (OUT)    :: power_out(1000)      
    REAL (R8), INTENT (OUT)    :: times_out(1000)      
    REAL (R8), INTENT (IN)     :: tstart, tend, frequency, tau     
    REAL (R8)                  :: tstep     

    INTEGER                    :: i, k
    INTEGER                    :: nstep,  nlevel


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

! +++ Number of steps:
    nstep               = int((tend-tstart)*frequency)
    k                   = 0
    tstep               = (tend-tstart)/(nstep)

    times_out           = 1.d10

    DO i = 1, min(500,nstep)
       k                = k+1

       times_out(i*2-1) = tstart + tstep*(i-1)
       times_out(i*2)   = tstart + tstep*i - min(tau,tstep)*0.01_r8

       power_out(i*2-1) = power_in(k)
       power_out(i*2)   = power_in(k)

       IF (k.GE.nlevel) k = 0
    END DO


 10 RETURN


    END SUBROUTINE fc2k_power_modulation