explicit declaration in modules use
This commit is contained in:
parent
88d0bffa22
commit
f8c7aaf924
@ -1,6 +1,8 @@
|
|||||||
module calcei_mod
|
module calcei_mod
|
||||||
|
|
||||||
contains
|
implicit none
|
||||||
|
|
||||||
|
contains
|
||||||
|
|
||||||
! ======================================================================
|
! ======================================================================
|
||||||
! nist guide to available math software.
|
! nist guide to available math software.
|
||||||
|
@ -34,6 +34,7 @@
|
|||||||
REAL(wp_), PARAMETER :: zero = 0.0_wp_
|
REAL(wp_), PARAMETER :: zero = 0.0_wp_
|
||||||
REAL(wp_), PARAMETER :: one = 1.0_wp_
|
REAL(wp_), PARAMETER :: one = 1.0_wp_
|
||||||
real(wp_), parameter :: pi = 3.141592653589793_wp_ ! 3.141592653589793238462643383280
|
real(wp_), parameter :: pi = 3.141592653589793_wp_ ! 3.141592653589793238462643383280
|
||||||
|
real(wp_), parameter :: pihalf = 1.57079632679489661923_wp_
|
||||||
REAL(wp_), PARAMETER :: sqrt_pi = 1.772453850905516_wp_
|
REAL(wp_), PARAMETER :: sqrt_pi = 1.772453850905516_wp_
|
||||||
! REAL(wp_), PARAMETER :: sqrt_2 = 1.414213562373095_wp_
|
! REAL(wp_), PARAMETER :: sqrt_2 = 1.414213562373095_wp_
|
||||||
REAL(wp_), PARAMETER :: degree = pi/180.0_wp_
|
REAL(wp_), PARAMETER :: degree = pi/180.0_wp_
|
||||||
|
@ -1,7 +1,7 @@
|
|||||||
module dispersion
|
module dispersion
|
||||||
!
|
!
|
||||||
use const_and_precisions, only : wp_,zero,one,im,czero,cunit,pi,sqrt_pi
|
use const_and_precisions, only : wp_,zero,one,im,czero,cunit,pi,sqrt_pi
|
||||||
use calcei_mod
|
use calcei_mod, only : calcei3
|
||||||
implicit none
|
implicit none
|
||||||
! local constants
|
! local constants
|
||||||
integer, parameter :: npts=500
|
integer, parameter :: npts=500
|
||||||
|
@ -1,6 +1,7 @@
|
|||||||
c
|
c
|
||||||
c
|
c
|
||||||
c Integration routine dqags.f from quadpack and dependencies: BEGIN
|
c Integration routine dqags.f from quadpack and dependencies: BEGIN
|
||||||
|
c Modified version for functions f(x,yi) with more than one variable
|
||||||
c
|
c
|
||||||
c
|
c
|
||||||
subroutine dqagsmv(f,a,b,apar,np,epsabs,epsrel,result,abserr,
|
subroutine dqagsmv(f,a,b,apar,np,epsabs,epsrel,result,abserr,
|
||||||
|
@ -881,7 +881,7 @@ c
|
|||||||
|
|
||||||
if (iprof.eq.1) then
|
if (iprof.eq.1) then
|
||||||
nprof=98
|
nprof=98
|
||||||
lprf=length(filenmprf)
|
lprf=len_trim(filenmprf)
|
||||||
open(file=filenmprf(1:lprf)//'.prf',
|
open(file=filenmprf(1:lprf)//'.prf',
|
||||||
. status= 'unknown',unit=nprof)
|
. status= 'unknown',unit=nprof)
|
||||||
call profdata
|
call profdata
|
||||||
@ -892,7 +892,7 @@ c read equilibrium data from file if iequil=2
|
|||||||
c
|
c
|
||||||
if (iequil.eq.2) then
|
if (iequil.eq.2) then
|
||||||
neqdsk=99
|
neqdsk=99
|
||||||
leqq=length(filenmeqq)
|
leqq=len_trim(filenmeqq)
|
||||||
open(file=filenmeqq(1:leqq)//'.eqdsk',
|
open(file=filenmeqq(1:leqq)//'.eqdsk',
|
||||||
. status= 'unknown', unit=neqdsk)
|
. status= 'unknown', unit=neqdsk)
|
||||||
call equidata
|
call equidata
|
||||||
@ -1081,7 +1081,7 @@ c
|
|||||||
c initial beam data are measured in mm -> transformed to cm
|
c initial beam data are measured in mm -> transformed to cm
|
||||||
c
|
c
|
||||||
nfbeam=97
|
nfbeam=97
|
||||||
lbm=length(filenmbm)
|
lbm=len_trim(filenmbm)
|
||||||
open(file=filenmbm(1:lbm)//'.txt',status= 'unknown',unit=nfbeam)
|
open(file=filenmbm(1:lbm)//'.txt',status= 'unknown',unit=nfbeam)
|
||||||
c
|
c
|
||||||
read(nfbeam,*) nisteer
|
read(nfbeam,*) nisteer
|
||||||
|
@ -8,7 +8,7 @@
|
|||||||
! the current drive value by adjoint approach
|
! the current drive value by adjoint approach
|
||||||
!
|
!
|
||||||
!########################################################################
|
!########################################################################
|
||||||
USE const_and_precisions
|
USE const_and_precisions, only : wp_
|
||||||
!-------
|
!-------
|
||||||
IMPLICIT NONE
|
IMPLICIT NONE
|
||||||
CHARACTER(Len=1), PRIVATE :: adj_appr(6) ! adjoint approach switcher
|
CHARACTER(Len=1), PRIVATE :: adj_appr(6) ! adjoint approach switcher
|
||||||
@ -170,6 +170,7 @@
|
|||||||
! K - Spitzer's function
|
! K - Spitzer's function
|
||||||
! dKdu = dK/du, i.e. its derivative over normalized momentum
|
! dKdu = dK/du, i.e. its derivative over normalized momentum
|
||||||
!=======================================================================
|
!=======================================================================
|
||||||
|
use const_and_precisions, only : comp_eps
|
||||||
IMPLICIT NONE
|
IMPLICIT NONE
|
||||||
REAL(wp_), INTENT(in) :: mu,Zeff,fc,u,q,gam
|
REAL(wp_), INTENT(in) :: mu,Zeff,fc,u,q,gam
|
||||||
REAL(wp_), INTENT(out) :: K,dKdu
|
REAL(wp_), INTENT(out) :: K,dKdu
|
||||||
@ -234,6 +235,7 @@
|
|||||||
! sfd(1),...,sfd(4) - coefficients of the polynomial expansion of the
|
! sfd(1),...,sfd(4) - coefficients of the polynomial expansion of the
|
||||||
! "Spitzer"-function (the same as in the Hirshman paper)
|
! "Spitzer"-function (the same as in the Hirshman paper)
|
||||||
!=======================================================================
|
!=======================================================================
|
||||||
|
use const_and_precisions, only : mc2_
|
||||||
IMPLICIT NONE
|
IMPLICIT NONE
|
||||||
REAL(wp_), INTENT(in) :: mu,Zeff,fc
|
REAL(wp_), INTENT(in) :: mu,Zeff,fc
|
||||||
INTEGER :: n,i,j
|
INTEGER :: n,i,j
|
||||||
@ -385,6 +387,7 @@
|
|||||||
! K - Spitzer function
|
! K - Spitzer function
|
||||||
! dKdu - its derivative
|
! dKdu - its derivative
|
||||||
!=======================================================================
|
!=======================================================================
|
||||||
|
use const_and_precisions, only : zero,one
|
||||||
IMPLICIT NONE
|
IMPLICIT NONE
|
||||||
REAL(wp_), INTENT(in) :: Zeff,fc,u,q,gam
|
REAL(wp_), INTENT(in) :: Zeff,fc,u,q,gam
|
||||||
REAL(wp_), INTENT(out) :: K,dKdu
|
REAL(wp_), INTENT(out) :: K,dKdu
|
||||||
|
Loading…
Reference in New Issue
Block a user