1 |
program test_inter_barxy |
program test_inter_barxy |
2 |
|
|
3 |
use comconst, only: dtvr, daysec, iniconst |
use comconst, only: iniconst |
4 |
use comdissnew, only: read_comdissnew |
use comdissnew, only: read_comdissnew |
5 |
use comgeom, only: inigeom |
use comgeom, only: inigeom |
6 |
use conf_gcm_m, only: conf_gcm, day_step |
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 |
USE dynetat0_m, only: rlonu, rlatv, rlatu, rlatu1, rlatu2, rlonv, xprimm025, & |
9 |
use disvert_m, only: pa |
xprimp025, xprimu, xprimv, yprimu1, yprimu2 |
10 |
|
use fxhyp_m, only: fxhyp |
11 |
|
use fyhyp_m, only: fyhyp |
12 |
use inter_barxy_m, only: inter_barxy |
use inter_barxy_m, only: inter_barxy |
13 |
USE nr_util, ONLY : pi |
USE nr_util, ONLY : pi |
14 |
use read_serre_m, only: read_serre |
use read_serre_m, only: read_serre |
36 |
|
|
37 |
CALL conf_gcm |
CALL conf_gcm |
38 |
call read_comdissnew |
call read_comdissnew |
|
dtvr = daysec / real(day_step) |
|
|
print *, 'dtvr = ', dtvr |
|
|
pa = 5e4 |
|
39 |
CALL iniconst |
CALL iniconst |
40 |
call read_serre |
call read_serre |
41 |
|
CALL fyhyp(rlatu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1) |
42 |
|
CALL fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025) |
43 |
|
|
44 |
|
rlatu(1) = pi / 2. |
45 |
|
rlatu(jjm + 1) = -rlatu(1) |
46 |
|
|
47 |
CALL inigeom |
CALL inigeom |
48 |
|
|
49 |
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)/) |