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(:, :) |
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 |
|
|
98 |
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, & |
99 |
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 |
|
100 |
phis(:, :) = phis(:, :) * 9.81 |
phis(:, :) = phis(:, :) * 9.81 |
101 |
|
|
102 |
mask(2:, 1) = mask(1, 1) ! north pole |
mask(2:, 1) = mask(1, 1) ! north pole |