1*-----------------------------------------------------------------------
2*     MAP PROJECTION (BRAUN PERSPECTIVE CYLINDRICAL) 2007-10-14 E. TOYODA
3*-----------------------------------------------------------------------
4*     Copyright (C) 2000-2007 GFD Dennou Club. All rights reserved.
5*-----------------------------------------------------------------------
6      SUBROUTINE MPFCYB(XLON, YLAT, X, Y)
7
8      PARAMETER (EPSL = 1.0E-5)
9      EXTERNAL  RFPI, XMPLON
10      SAVE C
11      DATA C /1.0/
12
13        PI = RFPI()
14        X = XMPLON(XLON)
15        Y = 2.0 * C * TAN(YLAT * 0.5)
16
17      RETURN
18*-----------------------------------------------------------------------
19      ENTRY MPICYB(X, Y, XLON, YLAT)
20
21        PI = RFPI()
22        XLON = X
23        YLAT = ATAN(Y * 0.5 / C) * 2.0
24        IF (ABS(XLON) .LE. PI) RETURN
25
26        CALL GLRGET('RUNDEF',RNA)
27        XLON = RNA
28        YLAT = RNA
29      RETURN
30*-----------------------------------------------------------------------
31      ENTRY MPSCYB(STDLAT)
32
33        C = COS(STDLAT)
34
35      END
36