Initialize Profiles at GP!!!
SUBROUTINE InitProfilesGP() ! MODULES USE MODgvec_MPI , ONLY: par_Bcast USE MODgvec_MHD3D_Vars, ONLY: pres_profile, chi_profile, Phi_profile IMPLICIT NONE !----------------------------------------------------------------------------------------------------------------------------------- ! INPUT VARIABLES !----------------------------------------------------------------------------------------------------------------------------------- ! OUTPUT VARIABLES !----------------------------------------------------------------------------------------------------------------------------------- ! LOCAL VARIABLES INTEGER :: iGP !=================================================================================================================================== IF(MPIroot)THEN !$OMP PARALLEL DO & !$OMP SCHEDULE(STATIC) DEFAULT(SHARED) PRIVATE(iGP) DO iGP=1,nGP chiPrime_GP( iGP) = chi_profile%eval_at_rho(s_GP(iGP),deriv=1) pres_GP( iGP) = pres_profile%eval_at_rho(s_GP(iGP)) PhiPrime_GP( iGP) = Phi_profile%eval_at_rho(s_GP(iGP),deriv=1) PhiPrime2_GP(iGP) = PhiPrime_GP( iGP)**2 END DO !iGP !$OMP END PARALLEL DO END IF !MPIroot CALL par_Bcast(chiPrime_GP,0) CALL par_Bcast(pres_GP,0) CALL par_Bcast(PhiPrime_GP,0) CALL par_Bcast(PhiPrime2_GP,0) END SUBROUTINE InitProfilesGP