## 13.7 Examples

                PROGRAM EX13_3
                DIMENSION NXA(4),NYA(4),XPOL(4),YPOL(4)
                CHARACTER*60 CTIT
                DATA NXA/200,1150,200,1150/NYA/1600,1600,2700,2700/
                DATA XPOL/0.,0.,0.,0./YPOL/0.,45.,90.,-45./

                CTIT='Azimuthal Lambert Projections'

                CALL SETPAG('DA4P')
                CALL DISINI
                CALL PAGERA
                CALL COMPLX

                CALL HEIGHT(50)
                NL=NLMESS(CTIT)
                NX=(2250-NL)/2.
                CALL MESSAG(CTIT,NX,300)

                CALL AXSLEN(900,900)
                CALL PROJCT('LAMBERT')

                DO I=1,4
                  CALL AXSPOS(NXA(I),NYA(I))
                  CALL MAPPOL(XPOL(I),YPOL(I))
                  CALL GRAFMP(-180.,180.,-180.,30.,-90.,90.,-90.,30.)

                  CALL COLOR('GREEN')
                  CALL WORLD
                  CALL COLOR('FORE')
                  CALL GRIDMP(1,1)
                  CALL ENDGRF
                END DO

                CALL DISFIN
                END

Figure 13.3: Azimuthal Lambert Projections