17 |
zthe_2d, zpic_2d, zval_2d) |
zthe_2d, zpic_2d, zval_2d) |
18 |
|
|
19 |
use conf_dat2d_m, only: conf_dat2d |
use conf_dat2d_m, only: conf_dat2d |
20 |
use comgeom, only: rlatu, rlonv |
use dynetat0_m, only: rlatu, rlonv |
21 |
use grid_noro_m, only: grid_noro |
use grid_noro_m, only: grid_noro |
22 |
use indicesol, only: epsfra |
use indicesol, only: epsfra |
23 |
use netcdf, only: nf90_nowrite |
use netcdf, only: nf90_nowrite |
53 |
INTEGER iml_rel |
INTEGER iml_rel |
54 |
INTEGER jml_rel |
INTEGER jml_rel |
55 |
INTEGER ncid, varid |
INTEGER ncid, varid |
56 |
REAL, pointer:: relief(:, :) |
REAL, ALLOCATABLE:: relief(:, :) ! in m |
57 |
REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:) |
REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:) |
58 |
REAL, pointer:: lon_ini(:), lat_ini(:) |
REAL, ALLOCATABLE:: lon_ini(:), lat_ini(:) |
59 |
|
|
60 |
!----------------------------------- |
!----------------------------------- |
61 |
|
|
92 |
|
|
93 |
CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, relief , & |
CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, relief , & |
94 |
interbar=.FALSE.) |
interbar=.FALSE.) |
|
deallocate(lon_ini, lat_ini) ! pointers |
|
95 |
|
|
96 |
print *, 'Compute all the parameters needed for the gravity wave drag code' |
print *, 'Compute all the parameters needed for the gravity wave drag code' |
|
|
|
97 |
CALL grid_noro(lon_rad, lat_rad, relief, rlonv, rlatu, phis, zmea_2d, & |
CALL grid_noro(lon_rad, lat_rad, relief, rlonv, rlatu, phis, zmea_2d, & |
98 |
zstd_2d, zsig_2d, zgam_2d, zthe_2d, zpic_2d, zval_2d, mask) |
zstd_2d, zsig_2d, zgam_2d, zthe_2d, zpic_2d, zval_2d, mask) |
|
deallocate(relief) ! pointer |
|
99 |
phis(:, :) = phis(:, :) * 9.81 |
phis(:, :) = phis(:, :) * 9.81 |
100 |
|
|
101 |
mask(2:, 1) = mask(1, 1) ! north pole |
mask(2:, 1) = mask(1, 1) ! north pole |