1*-----------------------------------------------------------------------
2*     MAP PROJECTION (CENTRAL CYLINDRICAL)          2007-10-14 E. TOYODA
3*-----------------------------------------------------------------------
4*     Copyright (C) 2000-2007 GFD Dennou Club. All rights reserved.
5*-----------------------------------------------------------------------
6      SUBROUTINE MPFCYC(XLON, YLAT, X, Y)
7
8      PARAMETER (EPSL = 1.0E-5)
9      EXTERNAL  RFPI, XMPLON
10
11      PI = RFPI()
12      X = XMPLON(XLON)
13      YY = MAX(-PI * 0.5 + EPSL, MIN(PI * 0.5 - EPSL, YLAT))
14      Y = TAN(YY)
15
16      RETURN
17*-----------------------------------------------------------------------
18      ENTRY MPICYC(X, Y, XLON, YLAT)
19
20      PI = RFPI()
21      XLON = X
22      YLAT = ATAN(Y)
23      IF (ABS(XLON) .LE. PI .AND. ABS(YLAT) .LE. PI/2) RETURN
24
25      CALL GLRGET('RUNDEF',RNA)
26      XLON = RNA
27      YLAT = RNA
28
29      END
30