//WWBBAFSD JOB (WD23008M1H1020,WWB2-B6),BALLISH,TIME=1,REGION=8200K,
//    MSGCLASS=R,NOTIFY=$WD40BB
//*MAIN CLASS=XACHK,LINES=(12,C)
//*FORMAT PR,DDNAME=,DEST=NASJS3.WWR27
//*FORMAT PR,DDNAME=FT06F001,DEST=WWR27,COPIES=1
//*FORMAT PR,DDNAME=SYSPRINT,DEST=WWR27,COPIES=1
//*FORMAT PR,DDNAME=PLOTLOG,DEST=WWR27,COPIES=0
//*FORMAT PR,DDNAME=SYSMSG,DEST=WWR27,COPIES=1
//*FORMAT PR,DDNAME=SYSTERM,DEST=WWR27,COPIES=0
//*FORMAT PR,DDNAME=JESJCL,DEST=WWR27,COPIES=1
//*FORMAT PR,DEST=NASJS3.VTPLOT,DDNAME=RJERASTR,FORMS=PLOT,COPIES=1
//STEP        EXEC JS7F2CLG,SYSOUT='*'
//FORT.SYSIN  DD DISP=SHR,DSN=NWS.WD23.KAC.RTNEPH.SORC(PCHAR)
//        DD *
C.....   FOR BRAD BALLISH......20 MAR 92
C.....   CODE READS AF SNODEP DATA SET CREATED IN NESDIS MASTER MAP
C         FORMAT BY JOHN PRITCHARD....FEB 91
C         AND CHARACTER PLOTS THE SNOW,ICE ETC POINTS...
C.....   GLOBAL DATA FROM NH AND SH ARRAYS ... APRIL 91
C.....   FILLING HOLES IN THE ANALYSIS IF PT SURROUNDED BY 6 OR MORE
C          SNOW COVERED POINTS...,,,
C          SINCE OPEN SEA VS SEAICE IS TRUSTWORTHY (K.E.M.),TRIED TO
C          USE LAND/SEA FROM SSMI SNOW FILE AND ONLY FILL LAND PTS..
C          HOWEVER,SEEMS NOT SAME AS USED AT OFFUTT(SEAICE IN AF SNOW
C          NOT NECESSARILY OVER OPEN WATER OF SSMI DATA)..
C          SOOO....SINCE FILLING ONLY EFFECTS THE SMALL POLAR
C          INTERPOLATION REGION OF NMC GRID ANYHOW,
C          DECIDED TO CONTINUE TO FILL ALL PTS AS BEFORE...26APR91..
C
      PARAMETER (IIAF=512,JJAF=512,IIHAF=IIAF/2,JJHAF=JJAF/2)
      PARAMETER (JHAF1=JJHAF+1)
      INTEGER*2 IDATA(IIAF),CC(2),DD(2),MM(2),YY(2)
      DIMENSION SNODEP(IIAF,JJAF),SNOFIL(IIAF,JJAF)
      DIMENSION GRDLAT(IIAF,JJAF)
      CHARACTER*1 AFSNOW(IIAF,JJAF),SNOW,NOSNO,IPRT,SEAICE,IGLAC
      CHARACTER*1 ILAND,ISEA,ICOAST,QUEST
      CHARACTER*1 NUGTS,NUGTL,NULTS,NULTL,FILLD
      CHARACTER*1 KDAT1(8),INUM(10)
      CHARACTER*2 LLMO,NMO(12)
      CHARACTER*6  TITLE
      CHARACTER*30 LABL
      DIMENSION IDAYT(2),IMONT(2),IYRT(2)
      DIMENSION JDAYT(2),JMONT(2),JYRT(2)
      EQUIVALENCE (IDATA(5),KDAT1(1))
      EQUIVALENCE (LLMO,KDAT1(4))
C....... ABBREVIATIONS OF MONTHS..FROM J.PRITCHARD,NESDIS...JAN92
      DATA NMO/'JA','FE','MR','AP','MY','JN','JL','AU','SE','OC','NO',
     1         'DE'/
      DATA INUM/'1','2','3','4','5','6','7','8','9','0'/
      DATA SNOW/'S'/,NOSNO/' '/,SEAICE/'I'/,QUEST/'?'/,IGLAC/'G'/
      DATA ISEA/' '/,ILAND/'X'/,FILLD/'F'/,ICOAST/':'/
      DATA NUGTS/'+'/,NULTS/'-'/,NUGTL/'*'/,NULTL/':'/
      DATA LABL/'- A.F. - SNOW - ...........   '/
C....    BEGIN HERE ....
      IAF = IIAF
      JAF = IIAF
      IAF1 = IAF - 1
      JAF1 = JAF - 1
C.....    AF GRID ORIENTED WITH I=1 ON LEFT AND J=1 AT TOP
C.....      SO FOR NMC (WHERE) J=1 IS AT BOTTOM J AT POLE IS AFPOLE-1
      IPOLE = IIHAF + 1
      JPOLE = JJHAF
C     DO 1000 LHEM=1,2
C....      MHEM=1,-1 FOR NH,SH...
       LHEM = 1
       NHEM = 1
       IIN = 20
       IF (LHEM.EQ.2) THEN
        NHEM = -1
        IIN = 21
       END IF
       PRINT 200,NHEM
       REWIND IIN
       READ (IIN) IDATA
       PRINT 100,(IDATA(I),I=1,8)
C...   SAVE DATE (4 I*2 WORDS), FOR 17DE91, CC= 1,DD=7D,MM=E9,YY=1 ,,,
       CC(LHEM) = IDATA(5)
       DD(LHEM) = IDATA(6)
       MM(LHEM) = IDATA(7)
       YY(LHEM) = IDATA(8)
       TITLE(1:1)=KDAT1(2)
       TITLE(2:2)=KDAT1(3)
       TITLE(3:3)=KDAT1(4)
       TITLE(4:4)=KDAT1(5)
       TITLE(5:5)=KDAT1(6)
       TITLE(6:6)=KDAT1(7)
       WRITE(73,309) TITLE
  309  FORMAT(A6)
       DO 310 J=1,10
        JTENS = J
        IF (KDAT1(2).EQ.INUM(J)) GO TO 312
  310  CONTINUE
  312  IF (JTENS.EQ.10) JTENS=0
       DO 320 J=1,10
        JONES = J
        IF (KDAT1(3).EQ.INUM(J)) GO TO 314
  320  CONTINUE
  314  IF (JONES.EQ.10) JONES=0
       IDAYT(LHEM) = JTENS * 10 + JONES
       DO 410 J=1,10
        JTENS = J
        IF (KDAT1(6).EQ.INUM(J)) GO TO 412
  410  CONTINUE
  412  IF (JTENS.EQ.10) JTENS=0
       DO 420 J=1,10
        JONES = J
        IF (KDAT1(7).EQ.INUM(J)) GO TO 414
  420  CONTINUE
  414  IF (JONES.EQ.10) JONES=0
       IYRT(LHEM) = JTENS * 10 + JONES
       DO 422 J=1,12
        JMONTH = J
        IF (LLMO.EQ.NMO(J)) GO TO 424
  422  CONTINUE
  424  IMONT(LHEM) = JMONTH
       READ (IIN) IDATA
       PRINT 100,(IDATA(I),I=1,8)
C....   SNODEP IS IN TENTHS OF INCHES..(THOUGH HOLE INCHES IF GT 1 IN.)
C          J INDEX FLIPPING DONE TO CONFORM TO NMC CONVENTION.......
       DO 50 J=1,JAF
        READ (IIN) IDATA
        DO 45 I=1,IAF
          SNODEP(I,JAF+1-J) = IDATA(I)
          SNOFIL(I,JAF+1-J) = IDATA(I)
   45   CONTINUE
   50  CONTINUE
C....   BECAUSE OF SPECKLING HOLES IN AF DATA (CAN BE LAND OR WATER)
C             FILL THE NOSNO POINTS, WHICH ARE
C          SURROUNDED BY GE 7 SNO POINTS, WITH .1 INCH SNOW.....
C          (NOTE WE LOSE ICE DESIGNATION OF THE ERRONEOUS OPEN SEA PTS)
       DO 160 J=2,JAF1
        DO 160 I=2,IAF1
          ISUM = 0
          DO 159 L=1,3
           IF (SNOFIL(I+L-2,J-1).GT.0.) ISUM=ISUM+1
           IF (SNOFIL(I+L-2,J).GT.0.) ISUM=ISUM+1
           IF (SNOFIL(I+L-2,J+1).GT.0.) ISUM=ISUM+1
  159     CONTINUE
          IF (SNODEP(I,J).LE.0..AND.ISUM.GE.6) SNODEP(I,J)=5000.
  160  CONTINUE
       SMAX = SNODEP(1,1)
       SMIN = SNODEP(1,1)
C....     WE FLIP THE CHARACTER ARRAY(AFSNOW) BACK TO AF CONVENTION
       DO 60 J=1,JAF
        DO 60 I=1,IAF
          AFSNOW(I,J) = NOSNO
          IF (SNODEP(I,JAF+1-J).GT.0.) AFSNOW(I,J) = SNOW
          IF (SNODEP(I,JAF+1-J).GE.4090..AND.SNODEP(I,JAF+1-J).LT.4999.)
     1      AFSNOW(I,J) = NUGTS
          IF (SNODEP(I,JAF+1-J).GT.0..AND.SNODEP(I,JAF+1-J).LT.10.)
     1      AFSNOW(I,J) = NUGTL
          IF (SNODEP(I,JAF+1-J).GE.5000.) THEN
            AFSNOW(I,J) = FILLD
            SNODEP(I,JAF+1-J) = 1.
          END IF
          IF (SNODEP(I,J).GT.SMAX) SMAX=SNODEP(I,J)
          IF (SNODEP(I,J).LT.SMIN) SMIN=SNODEP(I,J)
   60  CONTINUE
       PRINT  110,SMAX,SMIN
       KPRT = (IAF-1)/128 + 1
C      CALL PCHAR(AFSNOW,LABL,IAF,JAF,KPRT)
C....    1/8 BEDIENT GRID POLAR STEREOGRAPHIC GRID....
C          NEEDED TO GET LAT/LON OF EACH NEPH POINT
       XMESHL = 47.625
C...    ORIENTATION FOR N.H. / S.H. GRID IS 80. / 260.
       ORIENT = 80.0
       IF (NHEM.LT.0) ORIENT = 260.0
C.................  SPECIAL IN CASE S.HEMISPHERE    ................
C....    FOR USE WITH W3FB05,XMSH GT 0 = NH.../ LT 0 = SH ...
       XMSH = NHEM * XMESHL
C.................  SPECIAL IN CASE S.HEMISPHERE    ................
       DO 150 J=1,JAF
        DO 150 I=1,IAF
         GRDLAT(I,J) = -99.
C----    DETERMINE THE LAT OF THE A.F. POINTS AND PLACE DATA
C          IN PROPER DIAGNOSTIC GRID ROW.....
C----    GET APPROPRIATE I,J IN THE NMC SENSE, POLE AT (0,0)
         XI = I - IPOLE
         XJ = J - JPOLE
         CALL W3FB05 (XI,XJ,XMSH,ORIENT,ALAT,ALON)
         IF(SNODEP(I,J) .LE. 9.9)  GO TO 148
         IF(SNODEP(I,J) .GT. 4900.0)  GO TO 148
  147    CALL SNOUT(XI,XJ,SNODEP(I,J))
  148    CONTINUE
         IF (NHEM.GT.0) THEN
           IF (ALAT.LT.0.) GO TO 150
         ELSE
           IF (ALAT.GT.0.) GO TO 150
         END IF
         GRDLAT(I,J) = ABS(ALAT)
  150  CONTINUE
 1000 CONTINUE
C     PRINT 120,CC(1),DD(1),MM(1),YY(1),CC(2),DD(2),MM(2),YY(2)
C     PRINT 140,IDAYT(1),IMONT(1),IYRT(1),IDAYT(2),IMONT(2),IYRT(2)
      STOP
  100 FORMAT(1H ,' DATA REFERENCE= ',8A2)
  110 FORMAT(1H ,' SNOW MAX=',E15.5,'  SNOW MIN=',E15.5)
  120 FORMAT(1H ,' NH-D,M,Y-TODA =',4A2,'..SH D,M,Y-TODA =',4A2)
  140 FORMAT(1H ,' NH-D,M,Y-TODA =',3I3,'..SH D,M,Y-TODA =',3I3)
  200 FORMAT(1H ,'  NHEM =',I5)
      END
       SUBROUTINE SNOUT(XI,XJ,PP)
       CHARACTER *1 XN,XW
       CHARACTER *80 CARD
       DATA     CARD / 80*' ' /
       DATA XN/'N'/
       DATA XW/'W'/
       DATA IC/1/
C      DO CORNERS PLUS CENTER OF BOXES
       I1=2
       XII=XI+.5*FLOAT(I1-1)
       J1=2
       IC=IC+1
       XJJ=XJ+.5*FLOAT(J1-1)
       CALL W3FB05(XII,XJJ,47.625,80.0,ALAT,ALON)
C      WRITE(72,105) CARD
  105  FORMAT(A)
       IPPF=INT(PP/10.0)
       IF(IPPF .GT. 9) IPPF=9
       WRITE(72,115) IC,ALAT,ALON,IPPF
  115  FORMAT(I5,1X,F5.2,'N ',F6.2,'W ',I3)
C      PRINT 100,ALAT,XN,ALON,XW,ID
   10  CONTINUE
C 100  FORMAT('       ',F5.2,A1,F7.2,A1,'     ',I5)
       RETURN
       END
/*
//*
//GO.FT06F001 DD SYSOUT=*
//GO.FT20F001 DD DISP=SHR,DSN=COM.SPPROD.SNODEPH.NHMAMAP
//GO.FT21F001 DD DISP=SHR,DSN=COM.SPPROD.SNODEPH.SHMAMAP
//GO.FT72F001 DD DSN=&&SNOGRID,DISP=(NEW,PASS),UNIT=SYSDA,
//            DCB=(RECFM=FB,LRECL=80,BLKSIZE=24000),SPACE=(TRK,150)
//GO.FT73F001 DD DSN=&&AFTITLE,DISP=(NEW,PASS),UNIT=SYSDA,
//            DCB=(RECFM=FB,LRECL=80,BLKSIZE=80),SPACE=(TRK,01)
//*
//PRFSPLOT EXEC NMCIDAS
//DRIVE.FT06F001 DD SYSOUT=*
//DRIVE.PLOTPARM DD *
 &PLOT MODEL=8224,XMAX=4800.,YMAX=4608.,RESERV=10.,UNITS=170.,MODE=5,
  ID=1,REGION=200.,LBLK=8000,MSGLVL=1,IOMASK=1,SCALE=1.0,
  IOPT=2,LYNES=2500 &END
/*
//DRIVE.SYSIN DD *
IN
SP SAMNXT 1 1
DM -50 90080 -140  1 3
PM NEXT
SL UNIT 66 4 CENTER 2 0 (STA=Y TXTCOL=1 TXTSIZ=7 PRT=N
PT 228 * AIRFORCE SNOW ANALYSIS    00Z  ^91,1
END FRAME
/*
//FT66F001 DD DSN=&&SNOGRID,DISP=(OLD,DELETE)
//FT91F001 DD DSN=&&AFTITLE,DISP=(OLD,DELETE)
//
//*  SECOND JOB NO LONGER RUNS VIA K. MITCHELL (1/3/96)
//*
//**BBSNOD JOB (WD23008M1H1020,WWB2-B6),BALLISH,TIME=1,REGION=8200K,
//    MSGCLASS=R,NOTIFY=$WD40BB
//*MAIN CLASS=XACHK,LINES=(12,C)
//*FORMAT PR,DDNAME=,DEST=NASJS3.WWR27
//*FORMAT PR,DDNAME=FT06F001,DEST=WWR27,COPIES=1
//*FORMAT PR,DDNAME=SYSPRINT,DEST=WWR27,COPIES=1
//*FORMAT PR,DDNAME=PLOTLOG,DEST=WWR27,COPIES=0
//*FORMAT PR,DDNAME=SYSMSG,DEST=WWR27,COPIES=1
//*FORMAT PR,DDNAME=SYSTERM,DEST=WWR27,COPIES=0
//*FORMAT PR,DDNAME=JESJCL,DEST=WWR27,COPIES=1
//*FORMAT PR,DEST=NASJS3.VTPLOT,DDNAME=RJERASTR,FORMS=PLOT,COPIES=1
//STEP        EXEC JS7F2CLG,SYSOUT='*'
//FORT.SYSIN  DD DISP=SHR,DSN=NWS.WD23.KAC.RTNEPH.SORC(PCHAR)
//        DD *
C.....   FOR BRAD BALLISH......20 MAR 92
C.....   CODE READS AF SNODEP DATA SET CREATED IN NESDIS MASTER MAP
C         FORMAT BY JOHN PRITCHARD....FEB 91
C         AND CHARACTER PLOTS THE SNOW,ICE ETC POINTS...
C.....   GLOBAL DATA FROM NH AND SH ARRAYS ... APRIL 91
C.....   FILLING HOLES IN THE ANALYSIS IF PT SURROUNDED BY 6 OR MORE
C          SNOW COVERED POINTS...,,,
C          SINCE OPEN SEA VS SEAICE IS TRUSTWORTHY (K.E.M.),TRIED TO
C          USE LAND/SEA FROM SSMI SNOW FILE AND ONLY FILL LAND PTS..
C          HOWEVER,SEEMS NOT SAME AS USED AT OFFUTT(SEAICE IN AF SNOW
C          NOT NECESSARILY OVER OPEN WATER OF SSMI DATA)..
C          SOOO....SINCE FILLING ONLY EFFECTS THE SMALL POLAR
C          INTERPOLATION REGION OF NMC GRID ANYHOW,
C          DECIDED TO CONTINUE TO FILL ALL PTS AS BEFORE...26APR91..
C
      PARAMETER (IIAF=512,JJAF=512,IIHAF=IIAF/2,JJHAF=JJAF/2)
      PARAMETER (JHAF1=JJHAF+1)
      INTEGER*2 IDATA(IIAF),CC(2),DD(2),MM(2),YY(2)
      DIMENSION SNODEP(IIAF,JJAF),SNOFIL(IIAF,JJAF)
      DIMENSION GRDLAT(IIAF,JJAF)
      CHARACTER*1 AFSNOW(IIAF,JJAF),SNOW,NOSNO,IPRT,SEAICE,IGLAC
      CHARACTER*1 ILAND,ISEA,ICOAST,QUEST
      CHARACTER*1 NUGTS,NUGTL,NULTS,NULTL,FILLD
      CHARACTER*1 KDAT1(8),INUM(10)
      CHARACTER*2 LLMO,NMO(12)
      CHARACTER*6  TITLE
      CHARACTER*30 LABL
      DIMENSION IDAYT(2),IMONT(2),IYRT(2)
      DIMENSION JDAYT(2),JMONT(2),JYRT(2)
      EQUIVALENCE (IDATA(5),KDAT1(1))
      EQUIVALENCE (LLMO,KDAT1(4))
C....... ABBREVIATIONS OF MONTHS..FROM J.PRITCHARD,NESDIS...JAN92
      DATA NMO/'JA','FE','MR','AP','MY','JN','JL','AU','SE','OC','NO',
     1         'DE'/
      DATA INUM/'1','2','3','4','5','6','7','8','9','0'/
      DATA SNOW/'S'/,NOSNO/' '/,SEAICE/'I'/,QUEST/'?'/,IGLAC/'G'/
      DATA ISEA/' '/,ILAND/'X'/,FILLD/'F'/,ICOAST/':'/
      DATA NUGTS/'+'/,NULTS/'-'/,NUGTL/'*'/,NULTL/':'/
      DATA LABL/'- A.F. - SNOW - ...........   '/
C....    BEGIN HERE ....
      IAF = IIAF
      JAF = IIAF
      IAF1 = IAF - 1
      JAF1 = JAF - 1
C.....    AF GRID ORIENTED WITH I=1 ON LEFT AND J=1 AT TOP
C.....      SO FOR NMC (WHERE) J=1 IS AT BOTTOM J AT POLE IS AFPOLE-1
      IPOLE = IIHAF + 1
      JPOLE = JJHAF
C     DO 1000 LHEM=1,2
C....      MHEM=1,-1 FOR NH,SH...
       LHEM = 1
       NHEM = 1
       IIN = 20
       IF (LHEM.EQ.2) THEN
        NHEM = -1
        IIN = 21
       END IF
       PRINT 200,NHEM
       REWIND IIN
       READ (IIN) IDATA
       PRINT 100,(IDATA(I),I=1,8)
C...   SAVE DATE (4 I*2 WORDS), FOR 17DE91, CC= 1,DD=7D,MM=E9,YY=1 ,,,
       CC(LHEM) = IDATA(5)
       DD(LHEM) = IDATA(6)
       MM(LHEM) = IDATA(7)
       YY(LHEM) = IDATA(8)
       TITLE(1:1)=KDAT1(2)
       TITLE(2:2)=KDAT1(3)
       TITLE(3:3)=KDAT1(4)
       TITLE(4:4)=KDAT1(5)
       TITLE(5:5)=KDAT1(6)
       TITLE(6:6)=KDAT1(7)
       WRITE(73,309) TITLE
  309  FORMAT(A6)
       DO 310 J=1,10
        JTENS = J
        IF (KDAT1(2).EQ.INUM(J)) GO TO 312
  310  CONTINUE
  312  IF (JTENS.EQ.10) JTENS=0
       DO 320 J=1,10
        JONES = J
        IF (KDAT1(3).EQ.INUM(J)) GO TO 314
  320  CONTINUE
  314  IF (JONES.EQ.10) JONES=0
       IDAYT(LHEM) = JTENS * 10 + JONES
       DO 410 J=1,10
        JTENS = J
        IF (KDAT1(6).EQ.INUM(J)) GO TO 412
  410  CONTINUE
  412  IF (JTENS.EQ.10) JTENS=0
       DO 420 J=1,10
        JONES = J
        IF (KDAT1(7).EQ.INUM(J)) GO TO 414
  420  CONTINUE
  414  IF (JONES.EQ.10) JONES=0
       IYRT(LHEM) = JTENS * 10 + JONES
       DO 422 J=1,12
        JMONTH = J
        IF (LLMO.EQ.NMO(J)) GO TO 424
  422  CONTINUE
  424  IMONT(LHEM) = JMONTH
       READ (IIN) IDATA
       PRINT 100,(IDATA(I),I=1,8)
C....   SNODEP IS IN TENTHS OF INCHES..(THOUGH HOLE INCHES IF GT 1 IN.)
C          J INDEX FLIPPING DONE TO CONFORM TO NMC CONVENTION.......
       DO 50 J=1,JAF
        READ (IIN) IDATA
        DO 45 I=1,IAF
          SNODEP(I,JAF+1-J) = IDATA(I)
          SNOFIL(I,JAF+1-J) = IDATA(I)
   45   CONTINUE
   50  CONTINUE
C....   BECAUSE OF SPECKLING HOLES IN AF DATA (CAN BE LAND OR WATER)
C             FILL THE NOSNO POINTS, WHICH ARE
C          SURROUNDED BY GE 7 SNO POINTS, WITH .1 INCH SNOW.....
C          (NOTE WE LOSE ICE DESIGNATION OF THE ERRONEOUS OPEN SEA PTS)
       DO 160 J=2,JAF1
        DO 160 I=2,IAF1
          ISUM = 0
          DO 159 L=1,3
           IF (SNOFIL(I+L-2,J-1).GT.0.) ISUM=ISUM+1
           IF (SNOFIL(I+L-2,J).GT.0.) ISUM=ISUM+1
           IF (SNOFIL(I+L-2,J+1).GT.0.) ISUM=ISUM+1
  159     CONTINUE
          IF (SNODEP(I,J).LE.0..AND.ISUM.GE.6) SNODEP(I,J)=5000.
  160  CONTINUE
       SMAX = SNODEP(1,1)
       SMIN = SNODEP(1,1)
C....     WE FLIP THE CHARACTER ARRAY(AFSNOW) BACK TO AF CONVENTION
       DO 60 J=1,JAF
        DO 60 I=1,IAF
          AFSNOW(I,J) = NOSNO
          IF (SNODEP(I,JAF+1-J).GT.0.) AFSNOW(I,J) = SNOW
          IF (SNODEP(I,JAF+1-J).GE.4090..AND.SNODEP(I,JAF+1-J).LT.4999.)
     1      AFSNOW(I,J) = NUGTS
          IF (SNODEP(I,JAF+1-J).GT.0..AND.SNODEP(I,JAF+1-J).LT.10.)
     1      AFSNOW(I,J) = NUGTL
          IF (SNODEP(I,JAF+1-J).GE.5000.) THEN
            AFSNOW(I,J) = FILLD
            SNODEP(I,JAF+1-J) = 1.
          END IF
          IF (SNODEP(I,J).GT.SMAX) SMAX=SNODEP(I,J)
          IF (SNODEP(I,J).LT.SMIN) SMIN=SNODEP(I,J)
   60  CONTINUE
       PRINT  110,SMAX,SMIN
       KPRT = (IAF-1)/128 + 1
C      CALL PCHAR(AFSNOW,LABL,IAF,JAF,KPRT)
C....    1/8 BEDIENT GRID POLAR STEREOGRAPHIC GRID....
C          NEEDED TO GET LAT/LON OF EACH NEPH POINT
       XMESHL = 47.625
C...    ORIENTATION FOR N.H. / S.H. GRID IS 80. / 260.
       ORIENT = 80.0
       IF (NHEM.LT.0) ORIENT = 260.0
C.................  SPECIAL IN CASE S.HEMISPHERE    ................
C....    FOR USE WITH W3FB05,XMSH GT 0 = NH.../ LT 0 = SH ...
       XMSH = NHEM * XMESHL
C.................  SPECIAL IN CASE S.HEMISPHERE    ................
       DO 150 J=1,JAF
        DO 150 I=1,IAF
         GRDLAT(I,J) = -99.
C----    DETERMINE THE LAT OF THE A.F. POINTS AND PLACE DATA
C          IN PROPER DIAGNOSTIC GRID ROW.....
C----    GET APPROPRIATE I,J IN THE NMC SENSE, POLE AT (0,0)
         XI = I - IPOLE
         XJ = J - JPOLE
         CALL W3FB05 (XI,XJ,XMSH,ORIENT,ALAT,ALON)
         IF(SNODEP(I,J) .LE. 9.9)  GO TO 148
         IF(SNODEP(I,J) .GT. 4900.0)  GO TO 148
  147    CALL SNOUT(XI,XJ,SNODEP(I,J))
  148    CONTINUE
         IF (NHEM.GT.0) THEN
           IF (ALAT.LT.0.) GO TO 150
         ELSE
           IF (ALAT.GT.0.) GO TO 150
         END IF
         GRDLAT(I,J) = ABS(ALAT)
  150  CONTINUE
 1000 CONTINUE
C     PRINT 120,CC(1),DD(1),MM(1),YY(1),CC(2),DD(2),MM(2),YY(2)
C     PRINT 140,IDAYT(1),IMONT(1),IYRT(1),IDAYT(2),IMONT(2),IYRT(2)
      STOP
  100 FORMAT(1H ,' DATA REFERENCE= ',8A2)
  110 FORMAT(1H ,' SNOW MAX=',E15.5,'  SNOW MIN=',E15.5)
  120 FORMAT(1H ,' NH-D,M,Y-TODA =',4A2,'..SH D,M,Y-TODA =',4A2)
  140 FORMAT(1H ,' NH-D,M,Y-TODA =',3I3,'..SH D,M,Y-TODA =',3I3)
  200 FORMAT(1H ,'  NHEM =',I5)
      END
       SUBROUTINE SNOUT(XI,XJ,PP)
       CHARACTER *1 XN,XW
       CHARACTER *80 CARD
       DATA     CARD / 80*' ' /
       DATA XN/'N'/
       DATA XW/'W'/
       DATA IC/1/
C      DO CORNERS PLUS CENTER OF BOXES
       I1=2
       XII=XI+.5*FLOAT(I1-1)
       J1=2
       IC=IC+1
       XJJ=XJ+.5*FLOAT(J1-1)
       CALL W3FB05(XII,XJJ,47.625,80.0,ALAT,ALON)
C      WRITE(72,105) CARD
  105  FORMAT(A)
       ISNO=INT(PP/10.0)
       IPPF=1
       IF(ISNO .GE. 6) IPPF=6
       WRITE(72,115) IC,ALAT,ALON,IPPF
  115  FORMAT(I5,1X,F5.2,'N ',F6.2,'W ',I3)
C      PRINT 100,ALAT,XN,ALON,XW,ID
   10  CONTINUE
C 100  FORMAT('       ',F5.2,A1,F7.2,A1,'     ',I5)
       RETURN
       END
/*
//*
//GO.FT06F001 DD SYSOUT=*
//GO.FT20F001 DD DISP=SHR,DSN=COM.SPPROD.SNODEPH.NHMAMAP
//GO.FT21F001 DD DISP=SHR,DSN=COM.SPPROD.SNODEPH.SHMAMAP
//GO.FT72F001 DD DSN=&&SNOGRID,DISP=(NEW,PASS),UNIT=SYSDA,
//            DCB=(RECFM=FB,LRECL=80,BLKSIZE=24000),SPACE=(TRK,150)
//GO.FT73F001 DD DSN=&&AFTITLE,DISP=(NEW,PASS),UNIT=SYSDA,
//            DCB=(RECFM=FB,LRECL=80,BLKSIZE=80),SPACE=(TRK,01)
//*
//PRFSPLOT EXEC NMCIDAS
//DRIVE.FT06F001 DD SYSOUT=*
//DRIVE.PLOTPARM DD *
 &PLOT MODEL=8224,XMAX=4800.,YMAX=4608.,RESERV=10.,UNITS=170.,MODE=5,
  ID=1,REGION=200.,LBLK=8000,MSGLVL=1,IOMASK=1,SCALE=1.0,
  IOPT=2,LYNES=2500 &END
/*
//DRIVE.SYSIN DD *
IN
SP SAMNXT 2 2
DM -50 45095 -40  1 3
PM NEXT
SL UNIT 66 4 CENTER 2 0 (STA=Y TXTCOL=1 TXTSIZ=7 PRT=N
PT 228 * AIRFORCE SNOW DEPTH       00Z  ^91,1
END FRAME
/*
//FT66F001 DD DSN=&&SNOGRID,DISP=(OLD,DELETE)
//FT91F001 DD DSN=&&AFTITLE,DISP=(OLD,DELETE)
//
//
//
//
//