fc2k_workaround_eqstart.f90

Go to the documentation of this file.

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

     SUBROUTINE eqstart(PROF_FLAG,    J0_FLAG,         Q0_FLAG,        &
                        coreprof_in,  equilibrium_in,  toroidfield_in, &
                        coreprof_out, equilibrium_out) 

     USE euitm_schemas

     USE equilibrium_start

     IMPLICIT NONE

! +++ CPO derived types:
     TYPE (type_equilibrium), POINTER  :: equilibrium_in(:)   !input CPO with geometry quantities from previous time
     TYPE (type_equilibrium), POINTER  :: equilibrium_out(:)  !output CPO with geometry quantities from previous iteration
     TYPE (type_toroidfield), POINTER  :: toroidfield_in(:)   !toroidal field, major radius, total current
     TYPE (type_coreprof),    POINTER  :: coreprof_in(:)      !input CPO with plasma profiles   
     TYPE (type_coreprof),    POINTER  :: coreprof_iter(:)    !iterration CPO with plasma profiles   
     TYPE (type_coreprof),    POINTER  :: coreprof_out(:)     !output CPO with plasma profiles   
     TYPE (type_coreprof),    POINTER  :: coreprof_new(:)     !output CPO with plasma profiles   


! +++ Dimensions:
     INTEGER                           :: prof_flag, flag     !determines primary profile
     INTEGER                           :: j0_flag             !flag for negative current, if J0_FLAG>0 negative current will be cut off
     INTEGER                           :: q0_flag             !Flag for positive dq/drho: 0-allowed, >0-cut off
     INTEGER,               PARAMETER  :: ext_equil = 2


     CALL start_profiles_consistency(prof_flag,    j0_flag,        q0_flag,        ext_equil,  &
                                       coreprof_in,  equilibrium_in, toroidfield_in,             &
                                       coreprof_out, equilibrium_out) 

     RETURN


     END SUBROUTINE eqstart