MSFPS.f

Go to the documentation of this file.

 
      SUBROUTINE msfps(LAT,TRUELAT1,MSF)
 
 
! Computes the map scale factor for a Polar Stereographic grid at a given
! latitude.
 
      IMPLICIT NONE
 
! Define some private constants
!
      REAL, PARAMETER   :: pi = 3.1415927
      REAL, PARAMETER   :: rad_per_deg = pi / 180.
 
      REAL, INTENT(IN)           :: lat  ! latitude where msf is requested
      REAL, INTENT(IN)           :: truelat1
      REAL, INTENT(OUT)          :: msf
 
      REAL                       :: psi1, psix, pole
 
      IF (truelat1 >= 0.) THEN
        psi1 = (90. - truelat1) * rad_per_deg
        pole =90.
      ELSE
        psi1 = (90. + truelat1) * rad_per_deg
        pole = -90.
      ENDIF
 
      psix = (pole - lat)*rad_per_deg
      msf = ((1.+cos(psi1))/(1.0 + cos(psix)))
      RETURN
 
      END SUBROUTINE msfps