28 |
! the US Navy gridpoint area and the surface of the model gridpoint |
! the US Navy gridpoint area and the surface of the model gridpoint |
29 |
! area. See "grid_noto.txt". |
! area. See "grid_noto.txt". |
30 |
|
|
31 |
use dimens_m, only: iim, jjm |
use dimensions, only: iim, jjm |
32 |
use mva9_m, only: mva9 |
use mva9_m, only: mva9 |
33 |
use nr_util, only: assert, pi |
use nr_util, only: assert, pi |
34 |
|
|
36 |
REAL, intent(in):: xdata(:) ! (iusn) |
REAL, intent(in):: xdata(:) ! (iusn) |
37 |
REAL, intent(in):: ydata(:) ! (jusn) |
REAL, intent(in):: ydata(:) ! (jusn) |
38 |
|
|
39 |
REAL, intent(in):: zdata(:, :) ! (iusn, jusn) input field |
REAL, intent(in):: zdata(:, :) ! (iusn, jusn) input field, in m |
40 |
REAL, intent(in):: x(:), y(:) ! coordinates of output field |
REAL, intent(in):: x(:), y(:) ! coordinates of output field |
41 |
|
|
42 |
! Correlations of US Navy orography gradients: |
! Correlations of US Navy orography gradients: |
43 |
REAL, intent(out):: zphi(:, :) ! (iim + 1, jjm + 1) orography not smoothed |
|
44 |
|
REAL, intent(out):: zphi(:, :) ! (iim + 1, jjm + 1) |
45 |
|
! geoptential height of orography, not smoothed, in m |
46 |
|
|
47 |
real, intent(out):: zmea(:, :) ! (iim + 1, jjm + 1) smoothed orography |
real, intent(out):: zmea(:, :) ! (iim + 1, jjm + 1) smoothed orography |
48 |
real, intent(out):: zstd(:, :) ! (iim + 1, jjm + 1) Standard deviation |
real, intent(out):: zstd(:, :) ! (iim + 1, jjm + 1) Standard deviation |
49 |
REAL, intent(out):: zsig(:, :) ! (iim + 1, jjm + 1) Slope |
REAL, intent(out):: zsig(:, :) ! (iim + 1, jjm + 1) Slope |
57 |
|
|
58 |
real, intent(out):: mask(:, :) ! (iim + 1, jjm + 1) fraction of land |
real, intent(out):: mask(:, :) ! (iim + 1, jjm + 1) fraction of land |
59 |
|
|
60 |
! Variables local to the procedure: |
! Local: |
61 |
|
|
62 |
! In this version it is assumed that the input data come from |
! In this version it is assumed that the input data come from |
63 |
! the US Navy dataset: |
! the US Navy dataset: |
64 |
integer, parameter:: iusn = 2160, jusn = 1080 |
integer, parameter:: iusn = 2160, jusn = 1080 |
65 |
integer, parameter:: iext = 216 |
integer, parameter:: iext = 216 |
66 |
REAL xusn(iusn + 2 * iext), yusn(jusn + 2) |
REAL xusn(iusn + 2 * iext), yusn(jusn + 2) |
67 |
REAL zusn(iusn + 2 * iext, jusn + 2) |
REAL zusn(iusn + 2 * iext, jusn + 2) ! in m |
68 |
|
|
69 |
! Intermediate fields (correlations of orography gradient) |
! Intermediate fields (correlations of orography gradient) |
70 |
REAL, dimension(iim + 1, jjm + 1):: ztz, zxtzx, zytzy, zxtzy, weight |
REAL, dimension(iim + 1, jjm + 1):: ztz, zxtzx, zytzy, zxtzy, weight |