! Program Name: ! Author(s)/Contact(s): ! Abstract: ! History Log: ! ! Usage: ! Parameters: ! Input Files: ! ! Output Files: ! ! ! Condition codes: ! ! If appropriate, descriptive troubleshooting instructions or ! likely causes for failures could be mentioned here with the ! appropriate error code ! ! User controllable options: subroutine lccone (fsplat,ssplat,sign,confac) real, parameter :: conv=0.01745329251994 integer :: sign real :: fsplat,ssplat,confac if (abs(fsplat-ssplat).lt.1.E-2) then confac = sin(fsplat*conv) else confac = log10(cos(fsplat*conv))-log10(cos(ssplat*conv)) confac = confac/(log10(tan((45.-float(sign)*fsplat/2.)*conv))-& log10(tan((45.-float(sign)*ssplat/2.)*conv))) endif end subroutine lccone