1 |
program test_inter_barxy |
program test_inter_barxy |
2 |
|
|
3 |
use comconst, only: dtvr, daysec, iniconst |
use comconst, only: iniconst |
4 |
use comgeom, only: rlonu, rlatv, inigeom |
use comdissnew, only: read_comdissnew |
5 |
use conf_gcm_m, only: conf_gcm, day_step |
use comgeom, only: inigeom |
6 |
use dimens_m, only: iim, jjm |
use conf_gcm_m, only: conf_gcm |
7 |
use disvert_m, only: pa |
use dimensions, only: iim, jjm |
8 |
|
USE dynetat0_m, only: rlonu, rlatv, rlatu, rlatu1, rlatu2, rlonv, xprimm025, & |
9 |
|
xprimp025, xprimu, xprimv, yprimu1, yprimu2, read_serre |
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 |
|
|
34 |
allocate(var_ana3d(iml_dyn, jml_dyn)) |
allocate(var_ana3d(iml_dyn, jml_dyn)) |
35 |
|
|
36 |
CALL conf_gcm |
CALL conf_gcm |
37 |
dtvr = daysec / real(day_step) |
call read_comdissnew |
|
print *, 'dtvr = ', dtvr |
|
|
pa = 5e4 |
|
38 |
CALL iniconst |
CALL iniconst |
39 |
|
call read_serre |
40 |
|
CALL fyhyp(rlatu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1) |
41 |
|
CALL fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025) |
42 |
|
|
43 |
|
rlatu(1) = pi / 2. |
44 |
|
rlatu(jjm + 1) = -rlatu(1) |
45 |
|
|
46 |
CALL inigeom |
CALL inigeom |
47 |
|
|
48 |
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)/) |