5 |
use comgeom, only: inigeom |
use comgeom, only: inigeom |
6 |
use conf_gcm_m, only: conf_gcm |
use conf_gcm_m, only: conf_gcm |
7 |
use dimensions, only: iim, jjm |
use dimensions, only: iim, jjm |
8 |
USE dynetat0_m, only: rlonu, rlatv, rlatu, rlatu1, rlatu2, rlonv, xprimm025, & |
USE dynetat0_m, only: rlonu, rlatv, fyhyp, fxhyp |
9 |
xprimp025, xprimu, xprimv, yprimu1, yprimu2, read_serre |
use dynetat0_chosen_m, only: read_serre |
|
use fxhyp_m, only: fxhyp |
|
|
use fyhyp_m, only: fyhyp |
|
10 |
use inter_barxy_m, only: inter_barxy |
use inter_barxy_m, only: inter_barxy |
11 |
USE nr_util, ONLY : pi |
USE nr_util, ONLY : pi |
12 |
|
|
35 |
call read_comdissnew |
call read_comdissnew |
36 |
CALL iniconst |
CALL iniconst |
37 |
call read_serre |
call read_serre |
38 |
CALL fyhyp(rlatu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1) |
CALL fyhyp |
39 |
CALL fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025) |
CALL fxhyp |
|
|
|
|
rlatu(1) = pi / 2. |
|
|
rlatu(jjm + 1) = -rlatu(1) |
|
|
|
|
40 |
CALL inigeom |
CALL inigeom |
41 |
|
|
42 |
lon_ini = - pi + 2 * pi / iml_dyn * (/(i, i = 0, iml_dyn - 1)/) |
lon_ini = - pi + 2 * pi / iml_dyn * (/(i, i = 0, iml_dyn - 1)/) |