#      PROGRAM SCPKT3

nx=73; ny=37
xmin=0; xmax=360; ymin=-90; ymax=90

p1 = []; nx.times{p1.push []}
p2 = []; ny.times{p2.push []}

#### EXTERNAL  ISGTRC

iws = 1

d1 = open('t810630.dat').read
d2 = open('t811231.dat').read
d1.split(" ").filter{|x| x.to_f}

      DO 20 r=190,245,5
        amin=R
        amax=r+5
*       idx=(r-180)*1.4*1000+999
        idx=int((r-170)*1.25)*1000+999
        CALL uestlv(AMIN,AMAX,IDX)
   20 CONTINUE 

      CALL sgopn( IWS )

      CALL sglset( 'LSOFTF', .FALSE. )
      CALL umlset( 'LGRIDMJ', .FALSE. )
      CALL umrset( 'DGRIDMN', 30.0 )
      CALL umiset( 'ITYPEMN', 1 )
      CALL umiset( 'ITYPEOUT', 1 )
      CALL SGFRM

      CALL sglset( 'L2TO3', .TRUE. )
*     CALL scseye(-1.0, -1.5,  2.0 )
      CALL scsobj( 0.5,  0.5,  0.0 )
      CALL scspln( 1, 2, -0.3)
      CALL SCSPRJ

      CALL sgssim( 0.12, 0.0, 0.0 )
      CALL sgsmpl( 180.0, 90.0, 0.0 )
      CALL sgsvpt( 0.1, 0.9, 0.1, 0.9 )
      CALL sgswnd( XMIN, XMAX, YMIN, YMAX )
      CALL sgstxy( -180.0, 180.0, -90.0, 90.0 )
      CALL sgstrn( isgtrc('MER') )
      CALL SGSTRF

      CALL sglset( 'LCLIP', .TRUE. )
      CALL slpwwr( 1 )

      CALL uetone( P1, NX, NX, NY )
      CALL umpmap( 'coast_world' )
      CALL UMPGLB
      CALL slpvpr( 1 )

      CALL sglset( 'LCLIP', .FALSE. )
      CALL sgtxzv( 0.5, 0.95, 'Temp. : 10 hPa : Jun 30, 1981',
     +             0.025, 0, 0, 3 )

      CALL scspln( 1, 2, 0.4)
*     CALL sgrset( 'YOFF3', -0.25 )
      CALL SCSPRJ

      CALL sglset( 'LCLIP', .TRUE. )
      CALL slpwwr( 1 )

      CALL uetone( P2, NX, NX, NY )
      CALL umpmap( 'coast_world' )
      CALL UMPGLB
      CALL slpvpr( 1 )

      CALL sglset( 'LCLIP', .FALSE. )
      CALL sgtxzv( 0.5, 0.95, 'Temp. : 10 hPa : Dec 31, 1981',
     +             0.025, 0, 0, 3 )

      CALL SGCLS

      END