|
Last change
on this file since 5186 was
5159,
checked in by abarral, 3 months ago
|
|
Put dimensions.h and paramet.h into modules
|
|
File size:
711 bytes
|
| Line | |
|---|
| 1 | ! This module replaces grid/fxy_sin.h |
|---|
| 2 | |
|---|
| 3 | MODULE lmdz_fxy_sin |
|---|
| 4 | USE comconst_mod, ONLY: pi |
|---|
| 5 | USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm |
|---|
| 6 | IMPLICIT NONE; PRIVATE |
|---|
| 7 | PUBLIC fx, fxprim, fy, fyprim |
|---|
| 8 | |
|---|
| 9 | CONTAINS |
|---|
| 10 | |
|---|
| 11 | REAL FUNCTION fx(ri) |
|---|
| 12 | REAL, INTENT(IN) :: ri |
|---|
| 13 | fx = 2. * pi / REAL(iim) * (ri - 0.5 * REAL(iim) - 1.) |
|---|
| 14 | END FUNCTION fx |
|---|
| 15 | |
|---|
| 16 | REAL FUNCTION fxprim(ri) |
|---|
| 17 | REAL, INTENT(IN) :: ri |
|---|
| 18 | fxprim = 2. * pi / REAL(iim) |
|---|
| 19 | END FUNCTION fxprim |
|---|
| 20 | |
|---|
| 21 | REAL FUNCTION fy(rj) |
|---|
| 22 | REAL, INTENT(IN) :: rj |
|---|
| 23 | fy = ASIN(1. + 2. * ((1. - rj) / REAL(jjm))) |
|---|
| 24 | END FUNCTION fy |
|---|
| 25 | |
|---|
| 26 | REAL FUNCTION fyprim(rj) |
|---|
| 27 | REAL, INTENT(IN) :: rj |
|---|
| 28 | fyprim = 1. / SQRT((rj - 1.) * (jjm + 1. - rj)) |
|---|
| 29 | END FUNCTION fyprim |
|---|
| 30 | END MODULE lmdz_fxy_sin |
|---|
Note: See
TracBrowser
for help on using the repository browser.
