|
◆ 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 864 del file geo_proj_class.F90.
867 ELEMENTAL SUBROUTINE proj_stretched_ll(lon,lat,x,y, &
868 longitude_stretch_pole, latitude_stretch_pole, stretch_factor)
869 DOUBLE PRECISION, INTENT(in) :: lon,lat
870 DOUBLE PRECISION, INTENT(out) :: x,y
871 DOUBLE PRECISION, INTENT(in) :: longitude_stretch_pole, latitude_stretch_pole, &
874 DOUBLE PRECISION :: csq
876 csq = stretch_factor**2
878 y = raddeg*asin((1.0d0 - csq + (1.0d0 + csq)*sin(degrad*lat)) / &
879 (1.0d0 + csq + (1.0d0 - csq)*sin(degrad*lat)))
881 END SUBROUTINE proj_stretched_ll
883 ELEMENTAL SUBROUTINE unproj_stretched_ll(x,y,lon,lat,&
884 longitude_stretch_pole, latitude_stretch_pole, stretch_factor)
885 DOUBLE PRECISION, INTENT(in) :: x,y
886 DOUBLE PRECISION, INTENT(out) :: lon,lat
887 DOUBLE PRECISION, INTENT(in) :: longitude_stretch_pole, latitude_stretch_pole, &
890 DOUBLE PRECISION :: csq
892 csq = stretch_factor**2
895 lat = raddeg*asin((csq - 1.0d0 + (csq + 1.0d0)*sin(degrad*y)) / &
896 (csq + 1.0d0 + (csq - 1.0d0)*sin(degrad*y)))
898 END SUBROUTINE unproj_stretched_ll
|