var/da/da_tools/da_lc_cone.inc

   1 subroutine da_lc_cone(truelat1, truelat2, cone)
   2
   3    !------------------------------------------------------------------------
   4    ! Purpose: compute the cone factor of a Lambert Conformal projection
   5    !------------------------------------------------------------------------
   6
   7    implicit none
   8
   9    real, intent(in)             :: truelat1  ! (-90 -> 90 degrees N)
  10    real, intent(in)             :: truelat2  !   "   "  "   "     "
  11    real, intent(out)            :: cone
  12
  13    if (trace_use_dull) call da_trace_entry("da_lc_cone")
  14
  15    ! First, see if this is a secant or tangent projection.  For tangent
  16    ! projections, truelat1 = truelat2 and the cone is tangent to the
  17    ! Earth's surface at this latitude.  For secant projections, the cone
  18    ! intersects the Earth's surface at each of the distinctly different
  19    ! latitudes
  20    if (abs(truelat1-truelat2) > 0.1) then
  21       cone = alog10(cos(truelat1*rad_per_deg)) - &
  22              alog10(cos(truelat2*rad_per_deg))
  23       cone = cone /(alog10(tan((45.0 - abs(truelat1)/2.0) * rad_per_deg)) - &
  24              alog10(tan((45.0 - abs(truelat2)/2.0) * rad_per_deg)))
  25    else
  26       cone = sin(abs(truelat1)*rad_per_deg)
  27    end if
  28
  29    if (trace_use_dull) call da_trace_exit("da_lc_cone")
  30
  31 end subroutine da_lc_cone
  32
  33