da_set_merc.inc

1 subroutine da_set_merc(proj)
2   
3    !--------------------------------------------------------------------------
4    ! Purpose: Sets up the remaining basic elements for the mercator projection
5    !--------------------------------------------------------------------------
6 
7    implicit none
8 
9    type(proj_info), intent(inout)       :: proj
10    real                                 :: clain
11 
12    !  Preliminary variables
13 
14    clain = COS(rad_per_deg*proj%truelat1)
15    proj%dlon = proj%dx / (earth_radius_m * clain)
16 
17    ! Compute distance from equator to origin, and store in the 
18    ! proj%rsw tag.
19 
20    proj%rsw = 0.0
21    if (proj%lat1 .NE. 0.0) then
22       proj%rsw = (alog(tan(0.5*((proj%lat1+90.)*rad_per_deg))))/proj%dlon
23    end if
24 
25 end subroutine da_set_merc
26 
27