| 1 | MODULE YOMLEG |
|---|
| 2 | |
|---|
| 3 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
|---|
| 4 | |
|---|
| 5 | IMPLICIT NONE |
|---|
| 6 | |
|---|
| 7 | SAVE |
|---|
| 8 | |
|---|
| 9 | ! ------------------------------------------------------------------ |
|---|
| 10 | !* Description of Legendre polynomials |
|---|
| 11 | |
|---|
| 12 | ! RW : weights of the gaussian quadrature |
|---|
| 13 | ! RMU : mu sin(theta) |
|---|
| 14 | ! R1MU2 : 1.-MU*MU cos(theta)**2 |
|---|
| 15 | ! R1MUI : 1./R1MU2 1/cos(theta)**2 |
|---|
| 16 | ! R1MUA : 1./R1MU2/RA 1/(a*cos(theta)**2) |
|---|
| 17 | ! RSQM2 : SQRT(R1MU2) cos(theta) |
|---|
| 18 | ! R1QM2 : 1./SQRT(R1MU2) 1/cos(theta) |
|---|
| 19 | ! RACTHE : 1./SQRT(R1MU2) 1/(a*cos(theta)) |
|---|
| 20 | ! RLATIG : arcsin(mu) theta GLOBAL VIEW |
|---|
| 21 | ! RLATI : arcsin(mu) theta |
|---|
| 22 | ! RIPI0 : bi-cubic interpolation coefficients |
|---|
| 23 | ! RIPI1 : bi-cubic interpolation coefficients |
|---|
| 24 | ! RIPI2 : bi-cubic interpolation coefficients |
|---|
| 25 | |
|---|
| 26 | REAL(KIND=JPRB),ALLOCATABLE:: RW(:) |
|---|
| 27 | REAL(KIND=JPRB),ALLOCATABLE:: RMU(:) |
|---|
| 28 | REAL(KIND=JPRB),ALLOCATABLE:: R1MU2(:) |
|---|
| 29 | REAL(KIND=JPRB),ALLOCATABLE:: R1MUI(:) |
|---|
| 30 | REAL(KIND=JPRB),ALLOCATABLE:: R1MUA(:) |
|---|
| 31 | REAL(KIND=JPRB),ALLOCATABLE:: RSQM2(:) |
|---|
| 32 | REAL(KIND=JPRB),ALLOCATABLE:: R1QM2(:) |
|---|
| 33 | REAL(KIND=JPRB),ALLOCATABLE:: RACTHE(:) |
|---|
| 34 | REAL(KIND=JPRB),ALLOCATABLE:: RLATIG(:) |
|---|
| 35 | REAL(KIND=JPRB),ALLOCATABLE:: RLATI(:) |
|---|
| 36 | REAL(KIND=JPRB),ALLOCATABLE:: RIPI0(:) |
|---|
| 37 | REAL(KIND=JPRB),ALLOCATABLE:: RIPI1(:) |
|---|
| 38 | REAL(KIND=JPRB),ALLOCATABLE:: RIPI2(:) |
|---|
| 39 | |
|---|
| 40 | !$OMP THREADPRIVATE(r1mu2,r1mua,r1mui,r1qm2,racthe,ripi0,ripi1,ripi2,rlati,rlatig,rmu,rsqm2,rw) |
|---|
| 41 | END MODULE YOMLEG |
|---|
