|
◆ geo_proj_proj()
elemental subroutine geo_proj_proj |
( |
type(geo_proj), intent(in) |
this, |
|
|
double precision, intent(in) |
lon, |
|
|
double precision, intent(in) |
lat, |
|
|
double precision, intent(out) |
x, |
|
|
double precision, intent(out) |
y |
|
) |
| |
Computes and returns coordinates in the projected system given the geographical coordinates.
- Parametri
-
[in] | this | object to project |
[in] | lon | geographical coordinates |
[in] | lat | geographical coordinates |
[out] | x | projected coordinates |
[out] | y | projected coordinates |
Definizione alla linea 874 del file geo_proj_class.F90.
874 DOUBLE PRECISION, INTENT(out) :: lon,lat 875 DOUBLE PRECISION, INTENT(in) :: longitude_stretch_pole, latitude_stretch_pole, & 878 DOUBLE PRECISION :: csq 880 csq = stretch_factor**2 883 lat = raddeg*asin((csq - 1.0d0 + (csq + 1.0d0)*sin(degrad*y)) / & 884 (csq + 1.0d0 + (csq - 1.0d0)*sin(degrad*y))) 886 END SUBROUTINE unproj_stretched_ll 897 ELEMENTAL SUBROUTINE proj_lambert(lon,lat,x,y, & 898 latin1, latin2, lov, lad, projection_center_flag) 899 DOUBLE PRECISION, INTENT(in) :: lon,lat 900 DOUBLE PRECISION, INTENT(out) :: x,y 901 DOUBLE PRECISION, INTENT(in) :: latin1, latin2, lov, lad 902 INTEGER, INTENT(in) :: projection_center_flag 904 DOUBLE PRECISION :: n, f, ro0, ro, cs1, cs2, cs3, pollat, angle, cot 905 DOUBLE PRECISION, PARAMETER :: epsy = 1.0d-100
|