cat > p.f <<\EOR
C=============================================================================C
C International Comprehensive Ocean-Atmosphere Data Set (ICOADS)  22 Aug 2012 C
C Filename:level: gts:01J                            Fortran 77 program+shell C
C Purpose: Convert NCDC GTS data to IMMA                     Author: S.Lubker C
C=============================================================================C
C Software Revision Information (previous version: 29 May 2012, level 01I):
C IMMA1 with attm5 including relative humidity (RH) and setting the relative
C humidity indicator to 1 (Relative humidity in whole Percentage, measured
C and originally reported).
C-----------------------------------------------------------------------3456789
C Software documentation for the (modifiable) example program {gts} and routine
C {locate}, and for the (invariant) user-interface routines {ship,buoy,cman}.
C
C As provided {gts}: (a) reads marine Global Telecommunication System (GTS)
C data from standard input, (b) outputs the data on standard output in IMMA
C format, and (c) outputs a summary of any errors encountered on fort.8.
C The recognized GTS data types are SHIP (FM 13), BUOY (FM 18) and Coastal-
C Marine Automated Network (C-MAN) data (modified FM 13).  These data types
C are distinguished by code identifiers BBXX, ZZYY, and CMAN, respectively.
C
C Each input report consist of an 18-character header followed by a space
C and the GTS string, e.g.:
C    BBXX01040220040212 HRRF 01234 99197 31492 46/// /1319 10280 40168
C where the header provides the code identifier and receipt date/time
C (assembled from bulletin header information) according to this format:
C    ZZYY01060320010923 Header example
C    ^^^^               Code identifier (BBXX, ZZYY, or CMAN)
C        ^^             Sequence number for split messages (obsolete; unused)
C          ^^           Month
C            ^^         Day
C              ^^^^     Year
C                  ^^   Hour
C                    ^^ Minute
C Or an 37-character header extended with the WMO "abbreviated heading"
C (assembled from bulletin header information):
C    ZZYY01060320010923SSVX13KWBC060923RRA
C                      TTAAii              Bulletin ID
C                            CCCC          Collecting centre
C                                YYGGgg    Day Hour Minute
C                                      BBB Delayed or Correction or Amendment
C
C The output IMMA data have a minimal IMMA structure (core+attm1+attm5+attm99)
C with deck (992-995) and SID (114) assigned appropriately for the NCDC GTS
C data, and are formatted according to ICOADS, instead of NCDC, conventions.
C The obsolete fixed NCDC structure (core+attm1+attm5+attm6+attm99) and NCDC
C conventions (signed longitude representation and leading zeros in some
C numeric fields) could if necessary be added via {qcimma}, as can settings
C of NID and C1 (e.g., for VOSClim data) and QC flags.
C
C External library required: {lmrlib}.
C Machine dependencies: none known.
C-----------------------------------------------------------------------3456789
      program gts
      implicit integer(a-e,g-z)
      character str*1019
  100 read(*,'(a)',end=900)str
      lth=index(str,' ')
           if (str(1:4).eq.'ZZYY') then
        call buoy(str(:lth),str(lth+1:))
      else if (str(1:4).eq.'BBXX') then
        call ship(str(:lth),str(lth+1:))
      else if (str(1:4).eq.'CMAN') then
cc        call ship(str(:lth),str(lth+1:))
        call cman(str(:lth+6),str(lth+7:))
      else
cc        print *,str
cc        stop 'STOP'
      endif
      goto 100
  900 call prnstr
      end
c-----------------------------------------------------------------------3456789
      subroutine locate(id,latit,longit)
c identify C-MAN locations using an external list
      character*(*) id
      real latit,longit
      parameter(nlist=9999)
      character list(nlist)*64
      character xml*999
      data n/0/
      save n
c
      if (n.eq.0) then
        open(13,file='activestations.xml')
        open(14,file='activestations.txt')
        write(14,'(a)')'stations missing in activestations.xml'
   10   read(13,'(a)',end=100)xml
        if (index(xml,'<station ').eq.0) goto 10
        if (n+1.gt.nlist) stop 'nlist'
        ndx=index(xml,' id="')+5
        if (ndx.eq.5) goto 10
        list(n+1)(:43)=xml(ndx:ndx+index(xml(ndx:),'"')-2)
        ndx=index(xml,' lat="')+6
        if (ndx.eq.6) goto 10
        list(n+1)(44:)=xml(ndx:ndx+index(xml(ndx:),'"')-2)
        ndx=index(xml,' lon="')+6
        if (ndx.eq.6) goto 10
        list(n+1)(51:)=xml(ndx:ndx+index(xml(ndx:),'"')-2)
        do j=1,16
          if (ichar(list(n+1)(j:j)).ge.ichar('a')
     +   .and.ichar(list(n+1)(j:j)).le.ichar('z')
     +    ) list(n+1)(j:j)=char(ichar(list(n+1)(j:j))-ichar(' '))
        enddo
        n=n+1
        goto 10
      endif
c
  100 do j=1,n
        if (id.eq.list(j)(:16)) then
          read(list(j),'(t44,f7.0,f8.0)')latit,longit
          latit=(nint(latit*1000.)+nint(sign(5.,latit)))/10/100.
          longit=(nint(longit*1000.)+nint(sign(5.,longit)))/10/100.
          longit=mod(longit+360.,360.)
          return
        endif
      enddo
      if (n+1.gt.nlist) stop 'nlist'
      write(list(n+1),'(a,t44,f7.0,f8.0)')id,latit,longit
      write(14,'(a)')list(n+1)
      n=n+1
      end
C=============================================================================C
C WARNING:  Code beyond this point should not require any modification.       C
C=============================================================================C
C-----------------------------------------------------------------------3456789
      subroutine ship(head,FM13)
      implicit integer(a-e,g-z)
      character*(*) head,FM13
     +,abbr*4,cvalu*8,imma*267,core*108,attm1*65,attm5*94,attm6*68
     +,tens*8,tenths*8
      parameter(YR=1,MO=2,DY=3,HR=4,LAT=5,LON=6,IM=7,ATTC=8,TI=9,LI=10
     +,DS=11,VS=12,NID=13,II=14,ID=15,C1=24,DI=26,D=27,WI=28,W=29,VI=30
     +,VV=31,WW=32,W1=33,SLP=34,A=35,PPP=36,IT=37,AT=38,WBTI=39,WBT=40
     +,DPTI=41,DPT=42,SI=43,SST=44,N=45,NH=46,CL=47,HI=48,H=49,CM=50
     +,CH=51,WD=52,WP=53,WH=54,SD=55,SP=56,SH=57,OS=58,OP=59,FM=60
     +,IMMV=61,IX=62,W2=63,WMI=64,SD2=65,SP2=66,SH2=67,IS=68,ES=69,RS=70
     +,IC1=71,IC2=72,IC3=73,IC4=74,IC5=75,IR=76,RRR=77,TR=78,NU=79
     +,QCI=80,QI1=81,QI2=82,QI3=83,QI4=84,QI5=85,QI6=86,QI7=87,QI8=88
     +,QI9=89,QI10=90,QI11=91,QI12=92,QI13=93,QI14=94,QI15=95,QI16=96
     +,QI17=97,QI18=98,QI19=99,QI20=100,QI21=101,HDG=102,COG=103,SOG=104
     +,SLL=105,SLHH=106,RWD=107,RWS=108,QI22=109,QI23=110,QI24=111
     +,QI25=112,QI26=113,QI27=114,QI28=115,QI29=116,RH=117,RHI=118
     +,AWSI=119,IMONO=120)
      dimension abbr(imono),cvalu(imono),fvalu(imono)
      data abbr/'yr','mo','dy','hr','lat','lon','im','attc','ti','li'
     +,'ds','vs','nid','ii','id','id+1','id+2','id+3','id+4','id+5'
     +,'id+6','id+7','id+8','c1','c1+1','di','d','wi','w','vi','vv','ww'
     +,'w1','slp','a','ppp','it','at','wbti','wbt','dpti','dpt','si'
     +,'sst','n','nh','cl','hi','h','cm','ch','wd','wp','wh','sd','sp'
     +,'sh','os','op','fm','immv','ix','w2','wmi','sd2','sp2','sh2','is'
     +,'es','rs','ic1','ic2','ic3','ic4','ic5','ir','rrr','tr','nu'
     +,'qci','qi1','qi2','qi3','qi4','qi5','qi6','qi7','qi8','qi9'
     +,'qi10','qi11','qi12','qi13','qi14','qi15','qi16','qi17','qi18'
     +,'qi19','qi20','qi21','hdg','cog','sog','sll','slhh','rwd','rws'
     +,'qi22','qi23','qi24','qi25','qi26','qi27','qi28','qi29','rh'
     +,'rhi','awsi','imono'/
     +,attm1/' 165'/,attm6/' 668'/,fmiss/-99999./
      external float
      real dwpt
C
      do 9 i=1,imono
        cvalu(i)=' '
        fvalu(i)=fmiss
    9 continue
      cvalu(yr)=head(11:14)
      fvalu(yr)=float(cvalu(yr),fmiss)
      if (fvalu(yr).ne.fmiss) then
        if (nint(fvalu(yr)).lt.1770
     +  .or.nint(fvalu(yr)).gt.2024) fvalu(yr)=fmiss
      endif
      cvalu(mo)=head(7:8)
      fvalu(mo)=float(cvalu(mo),fmiss)
      if (fvalu(mo).ne.fmiss) then
        if (nint(fvalu(mo)).lt.1
     +  .or.nint(fvalu(mo)).gt.12) fvalu(mo)=fmiss
      endif
      cvalu(dy)=head(9:10)
      fvalu(dy)=float(cvalu(dy),fmiss)
      if (fvalu(dy).ne.fmiss) then
        if (nint(fvalu(dy)).lt.1
     +  .or.nint(fvalu(dy)).gt.31) fvalu(dy)=fmiss
      endif
      if (head(:4).eq.'CMAN') then
        cvalu(hr)=head(15:16)
        fvalu(hr)=float(cvalu(hr),fmiss)
        if (fvalu(hr).ne.fmiss) then
          if (nint(fvalu(hr)).gt.23) then
            fvalu(hr)=fmiss
          else
            fvalu(ti)=0.
          endif
        endif
        cvalu(wi)='4'
        fvalu(wi)=4.
      endif
      fvalu(im)=1.
      fvalu(attc)=3.
C
C SECTION 0
C M(i)M(i)M(j)M(j)
C D....D
C A(1)b(w)n(b)n(b)n(b)
      j=index(FM13,' ')
      do 90 i=1,min(j-1,9)
        cvalu(id+i-1)=FM13(i:i)
        fvalu(id+i-1)=ichar(cvalu(id+i-1)(1:1))
   90 continue
      if (cvalu(wi).ne.' ') then
        attm1(11:16)='995114'
        call locate(FM13(:j),fvalu(lat),fvalu(lon))
        if (fvalu(lat).ne.fmiss) fvalu(li)=5.
      else
          if (j.eq.6
     +  .and. cvalu(id  ).ge.'0' .and. cvalu(id  ).le.'9'
     +  .and. cvalu(id+1).ge.'0' .and. cvalu(id+1).le.'9'
     +  .and. cvalu(id+2).ge.'0' .and. cvalu(id+2).le.'9'
     +  .and. cvalu(id+3).ge.'0' .and. cvalu(id+3).le.'9'
     +  .and. cvalu(id+4).ge.'0' .and. cvalu(id+4).le.'9') then
          attm1(11:16)='993114'
        else
          attm1(11:16)='992114'
        endif
C YYGGi(w)
      if (index(FM13(j+1:),' ').ne.6) goto 700
      day=fvalu(dy)
      cvalu(dy)=FM13(j+1:j+2)
      fvalu(dy)=float(cvalu(dy),fmiss)
      if (fvalu(dy).ne.fmiss) then
        if (nint(fvalu(dy)).lt.1
     +  .or.nint(fvalu(dy)).gt.31) fvalu(dy)=fmiss
      endif
      if (day.ne.fmiss .and. fvalu(dy).ne.fmiss) then
        if (abs(day-fvalu(dy)).gt.25.) then
          fvalu(mo)=fvalu(mo)+isign(1,day-nint(fvalu(dy)))
          if (nint(fvalu(mo)).lt.1) then
            fvalu(mo)=12.
            fvalu(yr)=fvalu(yr)-1.
          else if (nint(fvalu(mo)).gt.12) then
            fvalu(mo)=1.
            fvalu(yr)=fvalu(yr)+1.
          endif
        endif
      endif
      cvalu(hr)=FM13(j+3:j+4)
      fvalu(hr)=float(cvalu(hr),fmiss)
      if (fvalu(hr).ne.fmiss) then
        if (nint(fvalu(hr)).gt.23) then
          fvalu(hr)=fmiss
        else
          fvalu(ti)=0.
        endif
      endif
      cvalu(wi)=FM13(j+5:j+5)
      if (cvalu(wi).eq.'0'
     +.or.cvalu(wi).eq.'1'
     +.or.cvalu(wi).eq.'3'
     +.or.cvalu(wi).eq.'4')
     +fvalu(wi)=float(cvalu(wi),fmiss)
      j=j+6
C 99L(a)L(a)L(a)
      if (index(FM13(j+1:),' ').ne.6) goto 700
      if (FM13(j+7:j+7).eq.'1'
     +.or.FM13(j+7:j+7).eq.'3'
     +.or.FM13(j+7:j+7).eq.'5'
     +.or.FM13(j+7:j+7).eq.'7') then
        cvalu(lat)=FM13(j+1:j+5)
        if (cvalu(lat)(1:2).eq.'99') then
          fvalu(lat)=float(cvalu(lat),fmiss)
          if (fvalu(lat).ne.fmiss) then
            fvalu(lat)=mod(fvalu(lat)/10.,100.)
            if (nint(fvalu(lat)*10.).gt.900) then
              fvalu(lat)=fmiss
            else
              if (FM13(j+7:j+7).eq.'3'
     +        .or.FM13(j+7:j+7).eq.'5')
     +        fvalu(lat)=-fvalu(lat)
              fvalu(li)=0.
            endif
          endif
        endif
      endif
      j=j+6
C Q(c)L(o)L(o)L(o)L(o)
      if (index(FM13(j+1:),' ').ne.6) goto 700
      cvalu(lon)=FM13(j+1:j+5)
      if (cvalu(lon)(1:1).eq.'1'
     +.or.cvalu(lon)(1:1).eq.'3'
     +.or.cvalu(lon)(1:1).eq.'5'
     +.or.cvalu(lon)(1:1).eq.'7') then
        fvalu(lon)=float(cvalu(lon),fmiss)
        if (fvalu(lon).ne.fmiss) then
          fvalu(lon)=mod(fvalu(lon)/10.,1000.)
          if (nint(fvalu(lon)*10.).gt.1800) then
            fvalu(lon)=fmiss
          else
            if (cvalu(lon)(1:1).eq.'5'
     +      .or.cvalu(lon)(1:1).eq.'7')
     +      fvalu(lon)=mod(360.-fvalu(lon),360.)
            fvalu(li)=0.
          endif
        endif
      endif
      j=j+6
      endif
C
C SECTION 1
C i(R)i(x)hVV
      if (index(FM13(j+1:),' ').ne.6) goto 700
      cvalu(ir)=FM13(j+1:j+1)
      if (cvalu(ir).ge.'0' .and. cvalu(ir).le.'4')
     +fvalu(ir)=float(cvalu(ir),fmiss)
      cvalu(ix)=FM13(j+2:j+2)
      if (cvalu(ix).ge.'1' .and. cvalu(ix).le.'7')
     +fvalu(ix)=float(cvalu(ix),fmiss)
      cvalu(h)=FM13(j+3:j+3)
      fvalu(h)=float(tens(cvalu(h)),fmiss)
      cvalu(vv)=FM13(j+4:j+5)
      if (cvalu(vv)(1:1).eq.'9')
     +fvalu(vv)=float(cvalu(vv),fmiss)
      j=j+6
C Nddff
      if (index(FM13(j+1:),' ').ne.6) goto 700
      cvalu(n)=FM13(j+1:j+1)
      fvalu(n)=float(tens(cvalu(n)),fmiss)
      cvalu(d)=FM13(j+2:j+3)
      if (cvalu(d).eq.'00') then
        fvalu(d)=361.
      else if (cvalu(d).eq.'99') then
        fvalu(d)=362.
      else
        fvalu(d)=float(cvalu(d),fmiss)
        if (fvalu(d).ne.fmiss) then
          fvalu(d)=fvalu(d)*10.
          if (nint(fvalu(d)).gt.360) fvalu(d)=fmiss
        endif
      endif
      if (fvalu(d).ne.fmiss) fvalu(di)=0.
      if (fvalu(wi).ne.fmiss) then
        cvalu(w)=FM13(j+4:j+5)
        fvalu(w)=float(cvalu(w),fmiss)
        if (fvalu(w).ne.fmiss) then
          if (cvalu(wi).eq.'3'
     +    .or.cvalu(wi).eq.'4')
     +    fvalu(w)=fxktms(fvalu(w))
        endif
      endif
      j=j+6
C 00fff
      if (FM13(j:j+1).eq.' 0' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j-2:j+2).eq.'99 00' .and. FM13(j+3:j+5).ne.'///') then
          if (fvalu(wi).ne.fmiss) then
            cvalu(w)=FM13(j+3:j+5)
            fvalu(w)=float(cvalu(w),fmiss)
            if (fvalu(w).ne.fmiss) then
              if (cvalu(wi).eq.'3'
     +        .or.cvalu(wi).eq.'4')
     +        fvalu(w)=fxktms(fvalu(w))
              if (nint(fvalu(w)*10.).gt.999) fvalu(w)=fmiss
            endif
          endif
        endif
        j=j+6
      endif
      if (fvalu(w).ne.fmiss) then
      else
        if (fvalu(wi).ne.fmiss) then
          cvalu(wi)=' '
          fvalu(wi)=fmiss
        endif
      endif
C 1s(n)TTT
      if (FM13(j:j+1).eq.' 1' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(at)=FM13(j+2:j+5)
          if (cvalu(at)(1:1).eq.'0'
     +    .or.cvalu(at)(1:1).eq.'1') then
            fvalu(at)=float(tenths(cvalu(at)),fmiss)
            if (fvalu(at).ne.fmiss) then
              fvalu(at)=mod(fvalu(at)/10.,100.)
              if (cvalu(at)(1:1).eq.'1') fvalu(at)=-fvalu(at)
              call indic(cvalu(at),fvalu(it))
            endif
          endif
        endif
        j=j+6
      endif
C 2s(n)T(d)T(d)T(d)
C 29UUU
      if (FM13(j:j+3).eq.' 222') goto 222
      if (FM13(j:j+1).eq.' 2' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+2:j+2).eq.'9') then
          cvalu(rh)=FM13(j+3:j+5)
          fvalu(rh)=float(cvalu(rh),fmiss)
          if (fvalu(rh).ne.fmiss) then
            if (nint(fvalu(rh)).gt.100) then
              fvalu(rh)=fmiss
            else
              fvalu(rhi)=1.
            endif
          endif
          if (fvalu(at).ne.fmiss) then
            cvalu(dpt)=FM13(j+3:j+5)
            fvalu(dpt)=float(cvalu(dpt),fmiss)
            if (fvalu(dpt).ne.fmiss) then
              if (nint(fvalu(dpt)).gt.100) then
                fvalu(dpt)=fmiss
              else
                fvalu(dpt)=dwpt(fvalu(at),fvalu(dpt))
                fvalu(dpti)=1.
              endif
            endif
          endif
        else
     +  if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(dpt)=FM13(j+2:j+5)
          if (cvalu(dpt)(1:1).eq.'0'
     +    .or.cvalu(dpt)(1:1).eq.'1') then
            fvalu(dpt)=float(tenths(cvalu(dpt)),fmiss)
            if (fvalu(dpt).ne.fmiss) then
              fvalu(dpt)=mod(fvalu(dpt)/10.,100.)
              if (cvalu(dpt)(1:1).eq.'1') fvalu(dpt)=-fvalu(dpt)
              call indic(cvalu(dpt),fvalu(it))
            endif
          endif
        endif
        j=j+6
      endif
C 3P(0)P(0)P(0)P(0)
      if (FM13(j:j+4).eq.' 333 ') goto 333
      if (FM13(j:j+1).eq.' 3' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        j=j+6
      endif
C 4PPPP
C 4a(3)hhh
      if (FM13(j:j+4).eq.' 444 ') goto 444
      if (FM13(j:j+1).eq.' 4' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+2:j+2).eq.'1'
     +  .or.FM13(j+2:j+2).eq.'2'
     +  .or.FM13(j+2:j+2).eq.'5'
     +  .or.FM13(j+2:j+2).eq.'7'
     +  .or.FM13(j+2:j+2).eq.'8') then
        else
     +  if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(slp)=FM13(j+2:j+5)
          fvalu(slp)=float(tenths(cvalu(slp)),fmiss)
          if (fvalu(slp).ne.fmiss) then
            fvalu(slp)=fvalu(slp)/10.
            if (cvalu(slp)(1:1).eq.'0') fvalu(slp)=fvalu(slp)+1000.
          endif
        endif
        j=j+6
      endif
C 5appp
      if (FM13(j:j+4).eq.' 555 ') goto 555
      if (FM13(j:j+1).eq.' 5' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(a)=FM13(j+2:j+2)
          if (cvalu(a).ge.'0' .and. cvalu(a).le.'8')
     +    fvalu(a)=float(cvalu(a),fmiss)
          cvalu(ppp)=FM13(j+3:j+5)
          fvalu(ppp)=float(tenths(cvalu(ppp)),fmiss)
          if (fvalu(ppp).ne.fmiss) fvalu(ppp)=fvalu(ppp)/10.
        endif
        j=j+6
      endif
C 6RRRt(R)
      if (FM13(j:j+1).eq.' 6' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(rrr)=FM13(j+2:j+4)
          fvalu(rrr)=float(cvalu(rrr),fmiss)
          cvalu(tr)=FM13(j+5:j+5)
          fvalu(tr)=float(cvalu(tr),fmiss)
        endif
        j=j+6
      endif
C 7wwW(1)W(2)
C 7w(a)w(a)W(a1)W(a2)
      if (FM13(j:j+1).eq.' 7' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(ww)=FM13(j+2:j+3)
          fvalu(ww)=float(cvalu(ww),fmiss)
          cvalu(w1)=FM13(j+4:j+4)
          fvalu(w1)=float(cvalu(w1),fmiss)
          cvalu(w2)=FM13(j+5:j+5)
          fvalu(w2)=float(cvalu(w2),fmiss)
        endif
        j=j+6
      endif
C 8N(h)C(L)C(M)C(H)
      if (FM13(j:j+1).eq.' 8' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(nh)=FM13(j+2:j+2)
          fvalu(nh)=float(tens(cvalu(nh)),fmiss)
          cvalu(cl)=FM13(j+3:j+3)
          fvalu(cl)=float(tens(cvalu(cl)),fmiss)
          cvalu(cm)=FM13(j+4:j+4)
          fvalu(cm)=float(tens(cvalu(cm)),fmiss)
          cvalu(ch)=FM13(j+5:j+5)
          fvalu(ch)=float(tens(cvalu(ch)),fmiss)
        endif
        j=j+6
      endif
C 9GGgg
      if (FM13(j:j+1).eq.' 9' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+2:j+2).ge.'0' .and. FM13(j+2:j+2).le.'2'
     +.and. FM13(j+3:j+3).ge.'0' .and. FM13(j+3:j+3).le.'9'
     +.and. FM13(j+4:j+4).ge.'0' .and. FM13(j+4:j+4).le.'5'
     +.and. FM13(j+5:j+5).ge.'0' .and. FM13(j+5:j+5).le.'9' .and.
     +.not.(FM13(j+2:j+2).eq.'2' .and. FM13(j+3:j+3).gt.'3')) then
          if (fvalu(yr).ne.fmiss .and. fvalu(mo).ne.fmiss
     +  .and. fvalu(dy).ne.fmiss .and. fvalu(hr).ne.fmiss) then
            fvalu(hr)=float(FM13(j+2:j+3),fmiss)
     +               +float(FM13(j+4:j+5),fmiss)/60.
            if (cvalu(hr).ge.FM13(j+2:j+3)) then
            else
              call rxnddt(ixdtnd(nint(fvalu(dy)),nint(fvalu(mo))
     +        ,nint(fvalu(yr)))-1,day,month,year)
              fvalu(yr)=year
              fvalu(mo)=month
              fvalu(dy)=day
            endif
            fvalu(ti)=2.
          endif
        endif
        j=j+6
      endif
C
C SECTION 2
  222 continue
C 222D(s)v(s)
      if (FM13(j:j+3).eq.' 222') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        cvalu(ds)=FM13(j+4:j+4)
        fvalu(ds)=float(tens(cvalu(ds)),fmiss)
        cvalu(vs)=FM13(j+5:j+5)
        fvalu(vs)=float(tens(cvalu(vs)),fmiss)
        j=j+6
C 0s(s)T(w)T(w)T(w)
        if (FM13(j:j+1).eq.' 0' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(sst)=FM13(j+2:j+5)
            if (cvalu(sst)(1:1).eq.'0'
     +      .or.cvalu(sst)(1:1).eq.'1') then
              fvalu(si)=1.
            else
     +      if (cvalu(sst)(1:1).eq.'2'
     +      .or.cvalu(sst)(1:1).eq.'3') then
              fvalu(si)=0.
            else
     +      if (cvalu(sst)(1:1).eq.'4'
     +      .or.cvalu(sst)(1:1).eq.'5') then
              fvalu(si)=3.
            endif
            if (fvalu(si).ne.fmiss
     +      .or.cvalu(sst)(1:1).eq.'6'
     +      .or.cvalu(sst)(1:1).eq.'7') then
              fvalu(sst)=float(tenths(cvalu(sst)),fmiss)
              if (fvalu(sst).ne.fmiss) then
                fvalu(sst)=mod(fvalu(sst)/10.,100.)
                if (cvalu(sst)(1:1).eq.'1'
     +          .or.cvalu(sst)(1:1).eq.'3'
     +          .or.cvalu(sst)(1:1).eq.'5'
     +          .or.cvalu(sst)(1:1).eq.'7') fvalu(sst)=-fvalu(sst)
                call indic(cvalu(sst),fvalu(it))
              else
                fvalu(si)=fmiss
              endif
            endif
          endif
          j=j+6
        endif
C 1P(wa)P(wa)H(wa)H(wa)
        if (FM13(j:j+1).eq.' 1' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(wp)=FM13(j+2:j+3)
            fvalu(wp)=float(cvalu(wp),fmiss)
            cvalu(wh)=FM13(j+4:j+5)
            fvalu(wh)=float(cvalu(wh),fmiss)
          endif
          j=j+6
        endif
C 2P(w)P(w)H(w)H(w)
        if (FM13(j:j+1).eq.' 2' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(wp)=FM13(j+2:j+3)
            fvalu(wp)=float(cvalu(wp),fmiss)
            cvalu(wh)=FM13(j+4:j+5)
            fvalu(wh)=float(cvalu(wh),fmiss)
          endif
          j=j+6
        endif
C 3d(w1)d(w1)d(w2)d(w2)
        if (FM13(j:j+4).eq.' 333 ') goto 333
        if (FM13(j:j+1).eq.' 3' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(sd)=FM13(j+2:j+3)
            if (cvalu(sd).eq.'99') then
              fvalu(sd)=38.
            else
              fvalu(sd)=float(cvalu(sd),fmiss)
              if (nint(fvalu(sd)).gt.36) fvalu(sd)=fmiss
            endif
            cvalu(sd2)=FM13(j+4:j+5)
            if (cvalu(sd2).eq.'99') then
              fvalu(sd2)=38.
            else
              fvalu(sd2)=float(cvalu(sd2),fmiss)
              if (nint(fvalu(sd2)).gt.36) fvalu(sd2)=fmiss
            endif
          endif
          j=j+6
        endif
C 4P(w1)P(w1)H(w1)H(w1)
        if (FM13(j:j+4).eq.' 444 ') goto 444
        if (FM13(j:j+1).eq.' 4' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(sp)=FM13(j+2:j+3)
            fvalu(sp)=float(cvalu(sp),fmiss)
            cvalu(sh)=FM13(j+4:j+5)
            fvalu(sh)=float(cvalu(sh),fmiss)
          endif
          j=j+6
        endif
C 5P(w2)P(w2)H(w2)H(w2)
        if (FM13(j:j+4).eq.' 555 ') goto 555
        if (FM13(j:j+1).eq.' 5' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(sp2)=FM13(j+2:j+3)
            fvalu(sp2)=float(cvalu(sp2),fmiss)
            cvalu(sh2)=FM13(j+4:j+5)
            fvalu(sh2)=float(cvalu(sh2),fmiss)
          endif
          j=j+6
        endif
C 6I(s)E(s)E(s)R(s)
C ICING
        if (FM13(j:j+1).eq.' 6' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(is)=FM13(j+2:j+2)
            if (cvalu(is).ge.'1' .and. cvalu(is).le.'5')
     +      fvalu(is)=float(cvalu(is),fmiss)
            cvalu(es)=FM13(j+3:j+4)
            fvalu(es)=float(cvalu(es),fmiss)
            cvalu(rs)=FM13(j+5:j+5)
            if (cvalu(rs).ge.'0' .and. cvalu(rs).le.'4')
     +      fvalu(rs)=float(cvalu(rs),fmiss)
          endif
          j=j+6
        else if (FM13(j:j+5).eq.' ICING') then
          goto 700
        endif
C 70H(wa)H(wa)H(wa)
        if (FM13(j:j+2).eq.' 70' .or. FM13(j:j+2).eq.' //'
     +  .or.FM13(j:j+2).eq.' 7/') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+2:j+2).ne.'/' .and. FM13(j+3:j+5).ne.'///') then
cc            cvalu(wh)=FM13(j+3:j+5)
cc            fvalu(wh)=float(cvalu(wh),fmiss)
cc            if (fvalu(wh).ne.fmiss) fvalu(wh)=fvalu(wh)/10.*2.
          endif
          j=j+6
        endif
C 8s(w)T(b)T(b)T(b)
        if (FM13(j:j+1).eq.' 8' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(wbt)=FM13(j+2:j+5)
            if (cvalu(wbt)(1:1).eq.'0'
     +      .or.cvalu(wbt)(1:1).eq.'1') then
              fvalu(wbti)=0.
            else
     +      if (cvalu(wbt)(1:1).eq.'5'
     +      .or.cvalu(wbt)(1:1).eq.'6') then
              fvalu(wbti)=1.
            else
     +      if (cvalu(wbt)(1:1).eq.'2') then
              fvalu(wbti)=2.
            else
     +      if (cvalu(wbt)(1:1).eq.'7') then
              fvalu(wbti)=3.
            endif
            if (fvalu(wbti).ne.fmiss) then
              fvalu(wbt)=float(tenths(cvalu(wbt)),fmiss)
              if (fvalu(wbt).ne.fmiss) then
                fvalu(wbt)=mod(fvalu(wbt)/10.,100.)
c???ice bulb neg
                if (cvalu(wbt)(1:1).eq.'1'
     +          .or.cvalu(wbt)(1:1).eq.'2'
     +          .or.cvalu(wbt)(1:1).eq.'6'
     +          .or.cvalu(wbt)(1:1).eq.'7') fvalu(wbt)=-fvalu(wbt)
                call indic(cvalu(wbt),fvalu(it))
              else
                fvalu(wbti)=fmiss
              endif
            endif
          endif
          j=j+6
        endif
C ICE c(i)S(i)b(i)D(i)z(i)
        if (FM13(j:j+3).eq.' ICE') then
          if (index(FM13(j+1:),' ').eq.9) j=j-1
          j=j+4
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+5).ne.'/////') then
            cvalu(ic1)=FM13(j+1:j+1)
            fvalu(ic1)=float(tens(cvalu(ic1)),fmiss)
            cvalu(ic2)=FM13(j+2:j+2)
            fvalu(ic2)=float(tens(cvalu(ic2)),fmiss)
            cvalu(ic3)=FM13(j+3:j+3)
            fvalu(ic3)=float(tens(cvalu(ic3)),fmiss)
            cvalu(ic4)=FM13(j+4:j+4)
            fvalu(ic4)=float(tens(cvalu(ic4)),fmiss)
            cvalu(ic5)=FM13(j+5:j+5)
            fvalu(ic5)=float(tens(cvalu(ic5)),fmiss)
          endif
          j=j+6
        endif
      endif
C
C SECTION 3
  333 continue
C 333
C 0....
C 1s(n)T(x)T(x)T(x)
C 2s(n)T(n)T(n)T(n)
C 3Ejjj
C 4E'sss
C 5j(1)j(2)j(3)j(4)
C j(5)j(6)j(7)j(8)j(9)
C 6RRRt(R)
C 7R(24)R(24)R(24)R(24)
C 8N(s)Ch(s)h(s)
C 9S(p)S(p)s(p)s(p)
C 8000
C 0....
C 1....
C
C SECTION 4
  444 continue
C 444
C N'C'H'H'C(t)
C
C SECTION 5
  555 continue
C 555
C
  700 imma=core(fvalu,fmiss)//attm1//attm5(fvalu(os-2),fmiss)
cc      if (imma(119:121).eq.'992') imma(35:43)='MASKSTID'
      write(*,'(9a)')imma
     +,'99 0 '//head//FM13(:lentrm(FM13))
      do 790 i=1,imono
        if (cvalu(i).ne.' ' .and. fvalu(i).eq.fmiss)
     +  call savstr(abbr(i)//cvalu(i))
  790 continue
      end
C-----------------------------------------------------------------------3456789
      subroutine indic(cvalu,fvalu)
      character*(*) cvalu(*)
      real fvalu(*)
      integer wbt,it
      parameter(wbt=1,it=1)
c
      if (cvalu(wbt)(4:4).eq.'/') then
        if (nint(fvalu(it)).eq.0
     + .or. nint(fvalu(it)).eq.3) then
          fvalu(it)=3.
        else
          fvalu(it)=2.
        endif
      else
        if (nint(fvalu(it)).eq.2
     + .or. nint(fvalu(it)).eq.3) then
          fvalu(it)=3.
        else
          fvalu(it)=0.
        endif
      endif
      end
C-----------------------------------------------------------------------3456789
      character*8 function tenths(char)
      character*(*) char
      len=lentrm(char)
      if (char(len:len).eq.'/') then
        tenths=char(:len-1)//'0'
      else
        tenths=char
      endif
      end
C-----------------------------------------------------------------------3456789
      character*8 function tens(char)
      character*(*) char
      if (char.eq.'/') then
        tens='10'
      else
        tens=char
      endif
      end
C-----------------------------------------------------------------------3456789
      subroutine buoy(head,FM18)
      implicit integer(a-e,g-z)
      character*(*) head,FM18
     +,abbr*4,cvalu*8,imma*267,core*108,attm1*65,attm5*94,attm6*68
     +,redths*8
      parameter(YR=1,MO=2,DY=3,HR=4,LAT=5,LON=6,IM=7,ATTC=8,TI=9,LI=10
     +,DS=11,VS=12,NID=13,II=14,ID=15,C1=24,DI=26,D=27,WI=28,W=29,VI=30
     +,VV=31,WW=32,W1=33,SLP=34,A=35,PPP=36,IT=37,AT=38,WBTI=39,WBT=40
     +,DPTI=41,DPT=42,SI=43,SST=44,N=45,NH=46,CL=47,HI=48,H=49,CM=50
     +,CH=51,WD=52,WP=53,WH=54,SD=55,SP=56,SH=57,OS=58,OP=59,FM=60
     +,IMMV=61,IX=62,W2=63,WMI=64,SD2=65,SP2=66,SH2=67,IS=68,ES=69,RS=70
     +,IC1=71,IC2=72,IC3=73,IC4=74,IC5=75,IR=76,RRR=77,TR=78,NU=79
     +,QCI=80,QI1=81,QI2=82,QI3=83,QI4=84,QI5=85,QI6=86,QI7=87,QI8=88
     +,QI9=89,QI10=90,QI11=91,QI12=92,QI13=93,QI14=94,QI15=95,QI16=96
     +,QI17=97,QI18=98,QI19=99,QI20=100,QI21=101,HDG=102,COG=103,SOG=104
     +,SLL=105,SLHH=106,RWD=107,RWS=108,QI22=109,QI23=110,QI24=111
     +,QI25=112,QI26=113,QI27=114,QI28=115,QI29=116,RH=117,RHI=118
     +,AWSI=119,IMONO=120)
      dimension abbr(imono),cvalu(imono),fvalu(imono)
      data abbr/'yr','mo','dy','hr','lat','lon','im','attc','ti','li'
     +,'ds','vs','nid','ii','id','id+1','id+2','id+3','id+4','id+5'
     +,'id+6','id+7','id+8','c1','c1+1','di','d','wi','w','vi','vv','ww'
     +,'w1','slp','a','ppp','it','at','wbti','wbt','dpti','dpt','si'
     +,'sst','n','nh','cl','hi','h','cm','ch','wd','wp','wh','sd','sp'
     +,'sh','os','op','fm','immv','ix','w2','wmi','sd2','sp2','sh2','is'
     +,'es','rs','ic1','ic2','ic3','ic4','ic5','ir','rrr','tr','nu'
     +,'qci','qi1','qi2','qi3','qi4','qi5','qi6','qi7','qi8','qi9'
     +,'qi10','qi11','qi12','qi13','qi14','qi15','qi16','qi17','qi18'
     +,'qi19','qi20','qi21','hdg','cog','sog','sll','slhh','rwd','rws'
     +,'qi22','qi23','qi24','qi25','qi26','qi27','qi28','qi29','rh'
     +,'rhi','awsi','imono'/
     +,attm1/' 165'/,attm6/' 668'/,fmiss/-99999./
      external float
      real dwpt
C
      do 9 i=1,imono
        cvalu(i)=' '
        fvalu(i)=fmiss
    9 continue
      fvalu(im)=1.
      fvalu(attc)=3.
C
C SECTION 0
C M(i)M(i)M(j)M(j)
C A(1)b(w)n(b)n(b)n(b)
      j=index(FM18,' ')
      if (j.ne.6) goto 700
      do 90 i=1,min(j-1,9)
        cvalu(id+i-1)=FM18(i:i)
        fvalu(id+i-1)=ichar(cvalu(id+i-1)(1:1))
   90 continue
      attm1(11:16)='994114'
C YYMMJ
      if (index(FM18(j+1:),' ').ne.6) goto 700
      cvalu(dy)=FM18(j+1:j+2)
      fvalu(dy)=float(cvalu(dy),fmiss)
      if (fvalu(dy).ne.fmiss) then
        if (nint(fvalu(dy)).lt.1
     +  .or.nint(fvalu(dy)).gt.31) fvalu(dy)=fmiss
      endif
      cvalu(mo)=FM18(j+3:j+4)
      fvalu(mo)=float(cvalu(mo),fmiss)
      if (fvalu(mo).ne.fmiss) then
        if (nint(fvalu(mo)).lt.1
     +  .or.nint(fvalu(mo)).gt.12) fvalu(mo)=fmiss
      endif
      cvalu(yr)=FM18(j+5:j+5)
      fvalu(yr)=float(cvalu(yr),fmiss)
      if (fvalu(yr).ne.fmiss) then
c        if (nint(fvalu(yr)).ge.7) then
c          fvalu(yr)=1990.+fvalu(yr)
c        else
          fvalu(yr)=2010.+fvalu(yr)
c        endif
      endif
      j=j+6
C GGggi(w)
      if (index(FM18(j+1:),' ').ne.6) goto 700
      cvalu(hr)=FM18(j+1:j+4)
      fvalu(hr)=float(cvalu(hr),fmiss)
      if (fvalu(hr).ne.fmiss) then
        if (nint(fvalu(hr))/100.gt.23 .or.
     +  mod(nint(fvalu(hr)),100).gt.59) then
          fvalu(hr)=fmiss
        else
          fvalu(hr)=nint(fvalu(hr))/100
     +         +mod(nint(fvalu(hr)),100)/60.
          fvalu(ti)=2.
        endif
      endif
      cvalu(wi)=FM18(j+5:j+5)
      if (cvalu(wi).eq.'0'
     +.or.cvalu(wi).eq.'1'
     +.or.cvalu(wi).eq.'3'
     +.or.cvalu(wi).eq.'4')
     +fvalu(wi)=float(cvalu(wi),fmiss)
      j=j+6
C Q(c)L(a)L(a)L(a)L(a)L(a)
      if (index(FM18(j+1:),' ').ne.7) goto 700
      cvalu(lat)=FM18(j+1:j+6)
      if (cvalu(lat)(1:1).eq.'1'
     +.or.cvalu(lat)(1:1).eq.'3'
     +.or.cvalu(lat)(1:1).eq.'5'
     +.or.cvalu(lat)(1:1).eq.'7') then
        fvalu(lat)=float(redths(cvalu(lat)),fmiss)
        if (fvalu(lat).ne.fmiss) then
          fvalu(lat)=mod(fvalu(lat),100000.)/1000.
          if (nint(fvalu(lat)*1000.).gt.90000) then
            fvalu(lat)=fmiss
          else
            if (cvalu(lat)(1:1).eq.'3'
     +      .or.cvalu(lat)(1:1).eq.'5')
     +      fvalu(lat)=-fvalu(lat)
            if (cvalu(lat)(5:6).eq.'//') then
              fvalu(li)=0.
            else
              fvalu(li)=5.
            endif
          endif
        endif
      endif
      j=j+7
C L(o)L(o)L(o)L(o)L(o)L(o)
      if (index(FM18(j+1:),' ').ne.7) goto 700
      if (cvalu(lat)(1:1).eq.'1'
     +.or.cvalu(lat)(1:1).eq.'3'
     +.or.cvalu(lat)(1:1).eq.'5'
     +.or.cvalu(lat)(1:1).eq.'7') then
        cvalu(lon)=FM18(j+1:j+6)
        fvalu(lon)=float(redths(cvalu(lon)),fmiss)
        if (fvalu(lon).ne.fmiss) then
          fvalu(lon)=fvalu(lon)/1000.
          if (nint(fvalu(lon)*1000.).gt.180000) then
            fvalu(lon)=fmiss
          else
            fvalu(lon)=nint(nint(fvalu(lon)*1000.)/10.)*10./1000.
            if (cvalu(lat)(1:1).eq.'5'
     +      .or.cvalu(lat)(1:1).eq.'7')
     +      fvalu(lon)=mod(360.-fvalu(lon),360.)
            if (cvalu(lon)(5:6).eq.'//') then
              if (nint(fvalu(li)).eq.5) then
                fvalu(li)=2.
              else
                fvalu(li)=0.
              endif
            else
              if (nint(fvalu(li)).eq.0) then
                fvalu(li)=2.
              else
                fvalu(li)=5.
              endif
            endif
          endif
        endif
      endif
      j=j+7
C 6Q(l)Q(t)Q(A)/
      if (FM18(j:j+1).eq.' 6' .or. FM18(j:j+1).eq.' /') then
        if (index(FM18(j+1:),' ').ne.6) goto 700
        j=j+6
      endif
C
C SECTION 1
  111 continue
C 111Q(d)Q(x)
      if (FM18(j:j+3).eq.' 111') then
        if (index(FM18(j+1:),' ').ne.6) goto 700
        j=j+6
C 0ddff
        if (FM18(j:j+1).eq.' 0' .or. FM18(j:j+1).eq.' /') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          if (FM18(j+1:j+1).ne.'/' .and. FM18(j+2:j+5).ne.'////') then
            cvalu(d)=FM18(j+2:j+3)
            if (cvalu(d).eq.'00') then
              fvalu(d)=361.
            else if (cvalu(d).eq.'99') then
              fvalu(d)=362.
            else
              fvalu(d)=float(cvalu(d),fmiss)
              if (fvalu(d).ne.fmiss) then
                fvalu(d)=fvalu(d)*10.
                if (nint(fvalu(d)).gt.360) fvalu(d)=fmiss
              endif
            endif
            if (fvalu(d).ne.fmiss) fvalu(di)=0.
            if (fvalu(wi).ne.fmiss) then
              cvalu(w)=FM18(j+4:j+5)
              fvalu(w)=float(cvalu(w),fmiss)
              if (fvalu(w).ne.fmiss) then
                if (cvalu(wi).eq.'3'
     +          .or.cvalu(wi).eq.'4')
     +          fvalu(w)=fxktms(fvalu(w))
              endif
            endif
          endif
          j=j+6
        endif
        if (fvalu(w).ne.fmiss) then
        else
          if (fvalu(wi).ne.fmiss) then
            cvalu(wi)=' '
            fvalu(wi)=fmiss
          endif
        endif
C 1s(n)TTT
        if (FM18(j:j+1).eq.' 1' .or. FM18(j:j+1).eq.' /') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          if (FM18(j+1:j+1).ne.'/' .and. FM18(j+2:j+5).ne.'////') then
            cvalu(at)=FM18(j+2:j+5)
            if (cvalu(at)(1:1).eq.'0'
     +      .or.cvalu(at)(1:1).eq.'1') then
              fvalu(at)=float(cvalu(at),fmiss)
              if (fvalu(at).ne.fmiss) then
                fvalu(at)=mod(fvalu(at)/10.,100.)
                if (cvalu(at)(1:1).eq.'1') fvalu(at)=-fvalu(at)
                fvalu(it)=0.
              endif
            endif
          endif
          j=j+6
        endif
C 2s(n)T(d)T(d)T(d)
C 29UUU
        if (FM18(j:j+3).eq.' 222') goto 222
        if (FM18(j:j+1).eq.' 2' .or. FM18(j:j+1).eq.' /') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          if (FM18(j+2:j+2).eq.'9') then
            cvalu(rh)=FM18(j+3:j+5)
            fvalu(rh)=float(cvalu(rh),fmiss)
            if (fvalu(rh).ne.fmiss) then
              if (nint(fvalu(rh)).gt.100) then
                fvalu(rh)=fmiss
              else
                fvalu(rhi)=1.
              endif
            endif
            if (fvalu(at).ne.fmiss) then
              cvalu(dpt)=FM18(j+3:j+5)
              fvalu(dpt)=float(cvalu(dpt),fmiss)
              if (fvalu(dpt).ne.fmiss) then
                if (nint(fvalu(dpt)).gt.100) then
                  fvalu(dpt)=fmiss
                else
                  fvalu(dpt)=dwpt(fvalu(at),fvalu(dpt))
                  fvalu(dpti)=1.
                endif
              endif
            endif
          else
     +    if (FM18(j+1:j+1).ne.'/' .and. FM18(j+2:j+5).ne.'////') then
            cvalu(dpt)=FM18(j+2:j+5)
            if (cvalu(dpt)(1:1).eq.'0'
     +      .or.cvalu(dpt)(1:1).eq.'1') then
              fvalu(dpt)=float(cvalu(dpt),fmiss)
              if (fvalu(dpt).ne.fmiss) then
                fvalu(dpt)=mod(fvalu(dpt)/10.,100.)
                if (cvalu(dpt)(1:1).eq.'1') fvalu(dpt)=-fvalu(dpt)
                fvalu(it)=0.
              endif
            endif
          endif
          j=j+6
        endif
C 3P(0)P(0)P(0)P(0)
        if (FM18(j:j+3).eq.' 333') goto 333
        if (FM18(j:j+1).eq.' 3' .or. FM18(j:j+1).eq.' /') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          j=j+6
        endif
C 4PPPP
        if (FM18(j:j+4).eq.' 444 ') goto 444
        if (FM18(j:j+1).eq.' 4' .or. FM18(j:j+1).eq.' /') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          if (FM18(j+1:j+1).ne.'/' .and. FM18(j+2:j+5).ne.'////') then
            cvalu(slp)=FM18(j+2:j+5)
            fvalu(slp)=float(cvalu(slp),fmiss)
            if (fvalu(slp).ne.fmiss) then
              fvalu(slp)=fvalu(slp)/10.
              if (cvalu(slp)(1:1).eq.'0') fvalu(slp)=fvalu(slp)+1000.
            endif
          endif
          j=j+6
        endif
C 5appp
        if (FM18(j:j+4).eq.' 555 ') goto 555
        if (FM18(j:j+1).eq.' 5' .or. FM18(j:j+1).eq.' /') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          if (FM18(j+1:j+1).ne.'/' .and. FM18(j+2:j+5).ne.'////') then
            cvalu(a)=FM18(j+2:j+2)
            if (cvalu(a).ge.'0' .and. cvalu(a).le.'8')
     +      fvalu(a)=float(cvalu(a),fmiss)
            cvalu(ppp)=FM18(j+3:j+5)
            fvalu(ppp)=float(cvalu(ppp),fmiss)
            if (fvalu(ppp).ne.fmiss) fvalu(ppp)=fvalu(ppp)/10.
          endif
          j=j+6
        endif
      endif
C
C SECTION 2
  222 continue
C 222Q(d)Q(x)
      if (FM18(j:j+3).eq.' 222') then
        if (index(FM18(j+1:),' ').ne.6) goto 700
        j=j+6
C 0s(n)T(w)T(w)T(w)
        if (FM18(j:j+1).eq.' 0' .or. FM18(j:j+1).eq.' /') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          if (FM18(j+1:j+1).ne.'/' .and. FM18(j+2:j+5).ne.'////') then
            cvalu(sst)=FM18(j+2:j+5)
            if (cvalu(sst)(1:1).eq.'0'
     +      .or.cvalu(sst)(1:1).eq.'1') then
              fvalu(sst)=float(cvalu(sst),fmiss)
              if (fvalu(sst).ne.fmiss) then
                fvalu(sst)=mod(fvalu(sst)/10.,100.)
                if (cvalu(sst)(1:1).eq.'1') fvalu(sst)=-fvalu(sst)
                fvalu(it)=0.
              endif
            endif
          endif
          j=j+6
        endif
C 1P(wa)P(wa)H(wa)H(wa)
        if (FM18(j:j+1).eq.' 1' .or. FM18(j:j+1).eq.' /') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          if (FM18(j+1:j+1).ne.'/' .and. FM18(j+2:j+5).ne.'////') then
            cvalu(wp)=FM18(j+2:j+3)
            fvalu(wp)=float(cvalu(wp),fmiss)
            cvalu(wh)=FM18(j+4:j+5)
            fvalu(wh)=float(cvalu(wh),fmiss)
          endif
          j=j+6
        endif
C 20P(wa)P(wa)P(wa)
        if (FM18(j:j+2).eq.' 20' .or. FM18(j:j+2).eq.' //'
     +  .or.FM18(j:j+2).eq.' 2/') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          if (FM18(j+2:j+2).ne.'/' .and. FM18(j+3:j+5).ne.'///') then
            cvalu(wp)=FM18(j+3:j+5)
            fvalu(wp)=float(cvalu(wp),fmiss)
            if (fvalu(wp).ne.fmiss) fvalu(wp)=fvalu(wp)/10.
          endif
          j=j+6
        endif
C 21H(wa)H(wa)H(wa)
        if (FM18(j:j+2).eq.' 21' .or. FM18(j:j+2).eq.' //'
     +  .or.FM18(j:j+2).eq.' 2/') then
          if (index(FM18(j+1:),' ').ne.6) goto 700
          if (FM18(j+2:j+2).ne.'/' .and. FM18(j+3:j+5).ne.'///') then
            cvalu(wh)=FM18(j+3:j+5)
            fvalu(wh)=float(cvalu(wh),fmiss)
            if (fvalu(wh).ne.fmiss) fvalu(wh)=fvalu(wh)/10.*2.
          endif
          j=j+6
        endif
      endif
C
C SECTION 3
  333 continue
C 333Q(d1)Q(d2)
C 8887k(2)
C 2z(0)z(0)z(0)z(0).....2z(n)z(n)z(n)z(n)
C 3T(0)T(0)T(0)T(0).....3T(n)T(n)T(n)T(n)
C 4S(0)S(0)S(0)S(0).....4S(n)S(n)S(n)S(n)
C 66k(6)9k(3)
C 2z(0)z(0)z(0)z(0).....2z(n)z(n)z(n)z(n)
C d(0)d(0)c(0)c(0)c(0).....d(n)d(n)c(n)c(n)c(n)
C
C SECTION 4
  444 continue
C 444
C 1Q(P)Q(2)Q(TW)Q(4)
C 2Q(N)Q(L)Q(A)/
C Q(c)L(a)L(a)L(a)L(a)L(a)
C L(o)L(o)L(o)L(o)L(o)L(o)
C YYMMJ
C GGgg/
C 7V(B)V(B)d(B)d(B)
C 8V(i)V(i)V(i)V(i)
C 9i(d)Z(d)Z(d)Z(d)
C
C SECTION 5
  555 continue
C 555
C
  700 imma=core(fvalu,fmiss)//attm1//attm5(fvalu(os-2),fmiss)
      write(*,'(9a)')imma
     +,'99 0 '//head//FM18(:lentrm(FM18))
      do 790 i=1,imono
        if (cvalu(i).ne.' ' .and. fvalu(i).eq.fmiss)
     +  call savstr(abbr(i)//cvalu(i))
  790 continue
      end
C-----------------------------------------------------------------------3456789
      real function float(char,fmiss)
      character*(*) char
      float=fmiss
      int=0
      do 190 i=1,lentrm(char)
        if (char(i:i).lt.'0' .or. char(i:i).gt.'9') return
        int=int*10 + ichar(char(i:i))-ichar('0')
  190 continue
      if (char.ne.' ') float=int
      end
C-----------------------------------------------------------------------3456789
      character*8 function redths(char)
      character*(*) char
      len=lentrm(char)
      if (char(len-1:len).eq.'//') then
        redths=char(:len-2)//'00'
      else if (char(len:len).eq.'/') then
        redths=char(:len-1)//'0'
      else
        redths=char
      endif
      end
C-----------------------------------------------------------------------3456789
      subroutine cman(head,FM13)
      implicit integer(a-e,g-z)
      character*(*) head,FM13
     +,abbr*4,cvalu*8,imma*267,core*108,attm1*65,attm5*94,attm6*68
     +,tens*8,tenths*8
      parameter(YR=1,MO=2,DY=3,HR=4,LAT=5,LON=6,IM=7,ATTC=8,TI=9,LI=10
     +,DS=11,VS=12,NID=13,II=14,ID=15,C1=24,DI=26,D=27,WI=28,W=29,VI=30
     +,VV=31,WW=32,W1=33,SLP=34,A=35,PPP=36,IT=37,AT=38,WBTI=39,WBT=40
     +,DPTI=41,DPT=42,SI=43,SST=44,N=45,NH=46,CL=47,HI=48,H=49,CM=50
     +,CH=51,WD=52,WP=53,WH=54,SD=55,SP=56,SH=57,OS=58,OP=59,FM=60
     +,IMMV=61,IX=62,W2=63,WMI=64,SD2=65,SP2=66,SH2=67,IS=68,ES=69,RS=70
     +,IC1=71,IC2=72,IC3=73,IC4=74,IC5=75,IR=76,RRR=77,TR=78,NU=79
     +,QCI=80,QI1=81,QI2=82,QI3=83,QI4=84,QI5=85,QI6=86,QI7=87,QI8=88
     +,QI9=89,QI10=90,QI11=91,QI12=92,QI13=93,QI14=94,QI15=95,QI16=96
     +,QI17=97,QI18=98,QI19=99,QI20=100,QI21=101,HDG=102,COG=103,SOG=104
     +,SLL=105,SLHH=106,RWD=107,RWS=108,QI22=109,QI23=110,QI24=111
     +,QI25=112,QI26=113,QI27=114,QI28=115,QI29=116,RH=117,RHI=118
     +,AWSI=119,IMONO=120)
      dimension abbr(imono),cvalu(imono),fvalu(imono)
      data abbr/'yr','mo','dy','hr','lat','lon','im','attc','ti','li'
     +,'ds','vs','nid','ii','id','id+1','id+2','id+3','id+4','id+5'
     +,'id+6','id+7','id+8','c1','c1+1','di','d','wi','w','vi','vv','ww'
     +,'w1','slp','a','ppp','it','at','wbti','wbt','dpti','dpt','si'
     +,'sst','n','nh','cl','hi','h','cm','ch','wd','wp','wh','sd','sp'
     +,'sh','os','op','fm','immv','ix','w2','wmi','sd2','sp2','sh2','is'
     +,'es','rs','ic1','ic2','ic3','ic4','ic5','ir','rrr','tr','nu'
     +,'qci','qi1','qi2','qi3','qi4','qi5','qi6','qi7','qi8','qi9'
     +,'qi10','qi11','qi12','qi13','qi14','qi15','qi16','qi17','qi18'
     +,'qi19','qi20','qi21','hdg','cog','sog','sll','slhh','rwd','rws'
     +,'qi22','qi23','qi24','qi25','qi26','qi27','qi28','qi29','rh'
     +,'rhi','awsi','imono'/
     +,attm1/' 165'/,attm6/' 668'/,fmiss/-99999./
      external float
      real dwpt
C
      do 9 i=1,imono
        cvalu(i)=' '
        fvalu(i)=fmiss
    9 continue
      cvalu(yr)=head(11:14)
      fvalu(yr)=float(cvalu(yr),fmiss)
      if (fvalu(yr).ne.fmiss) then
        if (nint(fvalu(yr)).lt.1770
     +  .or.nint(fvalu(yr)).gt.2024) fvalu(yr)=fmiss
      endif
      cvalu(mo)=head(7:8)
      fvalu(mo)=float(cvalu(mo),fmiss)
      if (fvalu(mo).ne.fmiss) then
        if (nint(fvalu(mo)).lt.1
     +  .or.nint(fvalu(mo)).gt.12) fvalu(mo)=fmiss
      endif
      cvalu(dy)=head(9:10)
      fvalu(dy)=float(cvalu(dy),fmiss)
      if (fvalu(dy).ne.fmiss) then
        if (nint(fvalu(dy)).lt.1
     +  .or.nint(fvalu(dy)).gt.31) fvalu(dy)=fmiss
      endif
      fvalu(im)=1.
      fvalu(attc)=3.
      j=-6
C YYGGi(w)
      if (index(FM13(j+1:),' ').ne.6) goto 700
      day=fvalu(dy)
      cvalu(dy)=FM13(j+1:j+2)
      fvalu(dy)=float(cvalu(dy),fmiss)
      if (fvalu(dy).ne.fmiss) then
        if (nint(fvalu(dy)).lt.1
     +  .or.nint(fvalu(dy)).gt.31) fvalu(dy)=fmiss
      endif
      if (day.ne.fmiss .and. fvalu(dy).ne.fmiss) then
        if (abs(day-fvalu(dy)).gt.25.) then
          fvalu(mo)=fvalu(mo)+isign(1,day-nint(fvalu(dy)))
          if (nint(fvalu(mo)).lt.1) then
            fvalu(mo)=12.
            fvalu(yr)=fvalu(yr)-1.
          else if (nint(fvalu(mo)).gt.12) then
            fvalu(mo)=1.
            fvalu(yr)=fvalu(yr)+1.
          endif
        endif
      endif
      cvalu(hr)=FM13(j+3:j+4)
      fvalu(hr)=float(cvalu(hr),fmiss)
      if (fvalu(hr).ne.fmiss) then
        if (nint(fvalu(hr)).gt.23) then
          fvalu(hr)=fmiss
        else
          fvalu(ti)=0.
        endif
      endif
      cvalu(wi)=FM13(j+5:j+5)
      if (cvalu(wi).eq.'0'
     +.or.cvalu(wi).eq.'1'
     +.or.cvalu(wi).eq.'3'
     +.or.cvalu(wi).eq.'4')
     +fvalu(wi)=float(cvalu(wi),fmiss)
      j=j+6
C
C SECTION 0
C M(i)M(i)M(j)M(j)
C D....D
C A(1)b(w)n(b)n(b)n(b)
      j=index(FM13,' ')
      if (j.ne.6) goto 700
      do 90 i=1,min(j-1,9)
        cvalu(id+i-1)=FM13(i:i)
        fvalu(id+i-1)=ichar(cvalu(id+i-1)(1:1))
   90 continue
      attm1(11:16)='995114'
      call locate(FM13(:j),fvalu(lat),fvalu(lon))
      if (fvalu(lat).ne.fmiss) fvalu(li)=5.
C
C SECTION 1
C i(R)i(x)hVV
      if (index(FM13(j+1:),' ').ne.6) goto 700
      cvalu(ir)=FM13(j+1:j+1)
      if (cvalu(ir).ge.'0' .and. cvalu(ir).le.'4')
     +fvalu(ir)=float(cvalu(ir),fmiss)
      cvalu(ix)=FM13(j+2:j+2)
      if (cvalu(ix).ge.'1' .and. cvalu(ix).le.'7')
     +fvalu(ix)=float(cvalu(ix),fmiss)
      cvalu(h)=FM13(j+3:j+3)
      fvalu(h)=float(tens(cvalu(h)),fmiss)
      cvalu(vv)=FM13(j+4:j+5)
      if (cvalu(vv)(1:1).eq.'9')
     +fvalu(vv)=float(cvalu(vv),fmiss)
      j=j+6
C Nddff
      if (index(FM13(j+1:),' ').ne.6) goto 700
      cvalu(n)=FM13(j+1:j+1)
      fvalu(n)=float(tens(cvalu(n)),fmiss)
      cvalu(d)=FM13(j+2:j+3)
      if (cvalu(d).eq.'00') then
        fvalu(d)=361.
      else if (cvalu(d).eq.'99') then
        fvalu(d)=362.
      else
        fvalu(d)=float(cvalu(d),fmiss)
        if (fvalu(d).ne.fmiss) then
          fvalu(d)=fvalu(d)*10.
          if (nint(fvalu(d)).gt.360) fvalu(d)=fmiss
        endif
      endif
      if (fvalu(d).ne.fmiss) fvalu(di)=0.
      if (fvalu(wi).ne.fmiss) then
        cvalu(w)=FM13(j+4:j+5)
        fvalu(w)=float(cvalu(w),fmiss)
        if (fvalu(w).ne.fmiss) then
          if (cvalu(wi).eq.'3'
     +    .or.cvalu(wi).eq.'4')
     +    fvalu(w)=fxktms(fvalu(w))
        endif
      endif
      j=j+6
C 00fff
      if (FM13(j:j+1).eq.' 0' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j-2:j+2).eq.'99 00' .and. FM13(j+3:j+5).ne.'///') then
          if (fvalu(wi).ne.fmiss) then
            cvalu(w)=FM13(j+3:j+5)
            fvalu(w)=float(cvalu(w),fmiss)
            if (fvalu(w).ne.fmiss) then
              if (cvalu(wi).eq.'3'
     +        .or.cvalu(wi).eq.'4')
     +        fvalu(w)=fxktms(fvalu(w))
              if (nint(fvalu(w)*10.).gt.999) fvalu(w)=fmiss
            endif
          endif
        endif
        j=j+6
      endif
      if (fvalu(w).ne.fmiss) then
      else
        if (fvalu(wi).ne.fmiss) then
          cvalu(wi)=' '
          fvalu(wi)=fmiss
        endif
      endif
C 1s(n)TTT
      if (FM13(j:j+1).eq.' 1' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(at)=FM13(j+2:j+5)
          if (cvalu(at)(1:1).eq.'0'
     +    .or.cvalu(at)(1:1).eq.'1') then
            fvalu(at)=float(tenths(cvalu(at)),fmiss)
            if (fvalu(at).ne.fmiss) then
              fvalu(at)=mod(fvalu(at)/10.,100.)
              if (cvalu(at)(1:1).eq.'1') fvalu(at)=-fvalu(at)
              call indic(cvalu(at),fvalu(it))
            endif
          endif
        endif
        j=j+6
      endif
C 2s(n)T(d)T(d)T(d)
C 29UUU
      if (FM13(j:j+3).eq.' 222') goto 222
      if (FM13(j:j+1).eq.' 2' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+2:j+2).eq.'9') then
          cvalu(rh)=FM13(j+3:j+5)
          fvalu(rh)=float(cvalu(rh),fmiss)
          if (fvalu(rh).ne.fmiss) then
            if (nint(fvalu(rh)).gt.100) then
              fvalu(rh)=fmiss
            else
              fvalu(rhi)=1.
            endif
          endif
          if (fvalu(at).ne.fmiss) then
            cvalu(dpt)=FM13(j+3:j+5)
            fvalu(dpt)=float(cvalu(dpt),fmiss)
            if (fvalu(dpt).ne.fmiss) then
              if (nint(fvalu(dpt)).gt.100) then
                fvalu(dpt)=fmiss
              else
                fvalu(dpt)=dwpt(fvalu(at),fvalu(dpt))
                fvalu(dpti)=1.
              endif
            endif
          endif
        else
     +  if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(dpt)=FM13(j+2:j+5)
          if (cvalu(dpt)(1:1).eq.'0'
     +    .or.cvalu(dpt)(1:1).eq.'1') then
            fvalu(dpt)=float(tenths(cvalu(dpt)),fmiss)
            if (fvalu(dpt).ne.fmiss) then
              fvalu(dpt)=mod(fvalu(dpt)/10.,100.)
              if (cvalu(dpt)(1:1).eq.'1') fvalu(dpt)=-fvalu(dpt)
              call indic(cvalu(dpt),fvalu(it))
            endif
          endif
        endif
        j=j+6
      endif
C 3P(0)P(0)P(0)P(0)
      if (FM13(j:j+4).eq.' 333 ') goto 333
      if (FM13(j:j+1).eq.' 3' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        j=j+6
      endif
C 4PPPP
C 4a(3)hhh
      if (FM13(j:j+4).eq.' 444 ') goto 444
      if (FM13(j:j+1).eq.' 4' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+2:j+2).eq.'1'
     +  .or.FM13(j+2:j+2).eq.'2'
     +  .or.FM13(j+2:j+2).eq.'5'
     +  .or.FM13(j+2:j+2).eq.'7'
     +  .or.FM13(j+2:j+2).eq.'8') then
        else
     +  if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(slp)=FM13(j+2:j+5)
          fvalu(slp)=float(tenths(cvalu(slp)),fmiss)
          if (fvalu(slp).ne.fmiss) then
            fvalu(slp)=fvalu(slp)/10.
            if (cvalu(slp)(1:1).eq.'0') fvalu(slp)=fvalu(slp)+1000.
          endif
        endif
        j=j+6
      endif
C 5appp
      if (FM13(j:j+4).eq.' 555 ') goto 555
      if (FM13(j:j+1).eq.' 5' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(a)=FM13(j+2:j+2)
          if (cvalu(a).ge.'0' .and. cvalu(a).le.'8')
     +    fvalu(a)=float(cvalu(a),fmiss)
          cvalu(ppp)=FM13(j+3:j+5)
          fvalu(ppp)=float(tenths(cvalu(ppp)),fmiss)
          if (fvalu(ppp).ne.fmiss) fvalu(ppp)=fvalu(ppp)/10.
        endif
        j=j+6
      endif
C 6RRRt(R)
      if (FM13(j:j+1).eq.' 6' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(rrr)=FM13(j+2:j+4)
          fvalu(rrr)=float(cvalu(rrr),fmiss)
          cvalu(tr)=FM13(j+5:j+5)
          fvalu(tr)=float(cvalu(tr),fmiss)
        endif
        j=j+6
      endif
C 7wwW(1)W(2)
C 7w(a)w(a)W(a1)W(a2)
      if (FM13(j:j+1).eq.' 7' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(ww)=FM13(j+2:j+3)
          fvalu(ww)=float(cvalu(ww),fmiss)
          cvalu(w1)=FM13(j+4:j+4)
          fvalu(w1)=float(cvalu(w1),fmiss)
          cvalu(w2)=FM13(j+5:j+5)
          fvalu(w2)=float(cvalu(w2),fmiss)
        endif
        j=j+6
      endif
C 8N(h)C(L)C(M)C(H)
      if (FM13(j:j+1).eq.' 8' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
          cvalu(nh)=FM13(j+2:j+2)
          fvalu(nh)=float(tens(cvalu(nh)),fmiss)
          cvalu(cl)=FM13(j+3:j+3)
          fvalu(cl)=float(tens(cvalu(cl)),fmiss)
          cvalu(cm)=FM13(j+4:j+4)
          fvalu(cm)=float(tens(cvalu(cm)),fmiss)
          cvalu(ch)=FM13(j+5:j+5)
          fvalu(ch)=float(tens(cvalu(ch)),fmiss)
        endif
        j=j+6
      endif
C 9GGgg
      if (FM13(j:j+1).eq.' 9' .or. FM13(j:j+1).eq.' /') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        if (FM13(j+2:j+2).ge.'0' .and. FM13(j+2:j+2).le.'2'
     +.and. FM13(j+3:j+3).ge.'0' .and. FM13(j+3:j+3).le.'9'
     +.and. FM13(j+4:j+4).ge.'0' .and. FM13(j+4:j+4).le.'5'
     +.and. FM13(j+5:j+5).ge.'0' .and. FM13(j+5:j+5).le.'9' .and.
     +.not.(FM13(j+2:j+2).eq.'2' .and. FM13(j+3:j+3).gt.'3')) then
          if (fvalu(yr).ne.fmiss .and. fvalu(mo).ne.fmiss
     +  .and. fvalu(dy).ne.fmiss .and. fvalu(hr).ne.fmiss) then
            fvalu(hr)=float(FM13(j+2:j+3),fmiss)
     +               +float(FM13(j+4:j+5),fmiss)/60.
            if (cvalu(hr).ge.FM13(j+2:j+3)) then
            else
              call rxnddt(ixdtnd(nint(fvalu(dy)),nint(fvalu(mo))
     +        ,nint(fvalu(yr)))-1,day,month,year)
              fvalu(yr)=year
              fvalu(mo)=month
              fvalu(dy)=day
            endif
            fvalu(ti)=2.
          endif
        endif
        j=j+6
      endif
C
C SECTION 2
  222 continue
C 222D(s)v(s)
      if (FM13(j:j+3).eq.' 222') then
        if (index(FM13(j+1:),' ').ne.6) goto 700
        cvalu(ds)=FM13(j+4:j+4)
        fvalu(ds)=float(tens(cvalu(ds)),fmiss)
        cvalu(vs)=FM13(j+5:j+5)
        fvalu(vs)=float(tens(cvalu(vs)),fmiss)
        j=j+6
C 0s(s)T(w)T(w)T(w)
        if (FM13(j:j+1).eq.' 0' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(sst)=FM13(j+2:j+5)
            if (cvalu(sst)(1:1).eq.'0'
     +      .or.cvalu(sst)(1:1).eq.'1') then
              fvalu(si)=1.
            else
     +      if (cvalu(sst)(1:1).eq.'2'
     +      .or.cvalu(sst)(1:1).eq.'3') then
              fvalu(si)=0.
            else
     +      if (cvalu(sst)(1:1).eq.'4'
     +      .or.cvalu(sst)(1:1).eq.'5') then
              fvalu(si)=3.
            endif
            if (fvalu(si).ne.fmiss
     +      .or.cvalu(sst)(1:1).eq.'6'
     +      .or.cvalu(sst)(1:1).eq.'7') then
              fvalu(sst)=float(tenths(cvalu(sst)),fmiss)
              if (fvalu(sst).ne.fmiss) then
                fvalu(sst)=mod(fvalu(sst)/10.,100.)
                if (cvalu(sst)(1:1).eq.'1'
     +          .or.cvalu(sst)(1:1).eq.'3'
     +          .or.cvalu(sst)(1:1).eq.'5'
     +          .or.cvalu(sst)(1:1).eq.'7') fvalu(sst)=-fvalu(sst)
                call indic(cvalu(sst),fvalu(it))
              else
                fvalu(si)=fmiss
              endif
            endif
          endif
          j=j+6
        endif
C 1P(wa)P(wa)H(wa)H(wa)
        if (FM13(j:j+1).eq.' 1' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(wp)=FM13(j+2:j+3)
            fvalu(wp)=float(cvalu(wp),fmiss)
            cvalu(wh)=FM13(j+4:j+5)
            fvalu(wh)=float(cvalu(wh),fmiss)
          endif
          j=j+6
        endif
C 2P(w)P(w)H(w)H(w)
        if (FM13(j:j+1).eq.' 2' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(wp)=FM13(j+2:j+3)
            fvalu(wp)=float(cvalu(wp),fmiss)
            cvalu(wh)=FM13(j+4:j+5)
            fvalu(wh)=float(cvalu(wh),fmiss)
          endif
          j=j+6
        endif
C 3d(w1)d(w1)d(w2)d(w2)
        if (FM13(j:j+4).eq.' 333 ') goto 333
        if (FM13(j:j+1).eq.' 3' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(sd)=FM13(j+2:j+3)
            if (cvalu(sd).eq.'99') then
              fvalu(sd)=38.
            else
              fvalu(sd)=float(cvalu(sd),fmiss)
              if (nint(fvalu(sd)).gt.36) fvalu(sd)=fmiss
            endif
            cvalu(sd2)=FM13(j+4:j+5)
            if (cvalu(sd2).eq.'99') then
              fvalu(sd2)=38.
            else
              fvalu(sd2)=float(cvalu(sd2),fmiss)
              if (nint(fvalu(sd2)).gt.36) fvalu(sd2)=fmiss
            endif
          endif
          j=j+6
        endif
C 4P(w1)P(w1)H(w1)H(w1)
        if (FM13(j:j+4).eq.' 444 ') goto 444
        if (FM13(j:j+1).eq.' 4' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(sp)=FM13(j+2:j+3)
            fvalu(sp)=float(cvalu(sp),fmiss)
            cvalu(sh)=FM13(j+4:j+5)
            fvalu(sh)=float(cvalu(sh),fmiss)
          endif
          j=j+6
        endif
C 5P(w2)P(w2)H(w2)H(w2)
        if (FM13(j:j+4).eq.' 555 ') goto 555
        if (FM13(j:j+1).eq.' 5' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(sp2)=FM13(j+2:j+3)
            fvalu(sp2)=float(cvalu(sp2),fmiss)
            cvalu(sh2)=FM13(j+4:j+5)
            fvalu(sh2)=float(cvalu(sh2),fmiss)
          endif
          j=j+6
        endif
C 6I(s)E(s)E(s)R(s)
C ICING
        if (FM13(j:j+1).eq.' 6' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(is)=FM13(j+2:j+2)
            if (cvalu(is).ge.'1' .and. cvalu(is).le.'5')
     +      fvalu(is)=float(cvalu(is),fmiss)
            cvalu(es)=FM13(j+3:j+4)
            fvalu(es)=float(cvalu(es),fmiss)
            cvalu(rs)=FM13(j+5:j+5)
            if (cvalu(rs).ge.'0' .and. cvalu(rs).le.'4')
     +      fvalu(rs)=float(cvalu(rs),fmiss)
          endif
          j=j+6
        else if (FM13(j:j+5).eq.' ICING') then
          goto 700
        endif
C 70H(wa)H(wa)H(wa)
        if (FM13(j:j+2).eq.' 70' .or. FM13(j:j+2).eq.' //'
     +  .or.FM13(j:j+2).eq.' 7/') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+2:j+2).ne.'/' .and. FM13(j+3:j+5).ne.'///') then
cc            cvalu(wh)=FM13(j+3:j+5)
cc            fvalu(wh)=float(cvalu(wh),fmiss)
cc            if (fvalu(wh).ne.fmiss) fvalu(wh)=fvalu(wh)/10.*2.
          endif
          j=j+6
        endif
C 8s(w)T(b)T(b)T(b)
        if (FM13(j:j+1).eq.' 8' .or. FM13(j:j+1).eq.' /') then
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+1).ne.'/' .and. FM13(j+2:j+5).ne.'////') then
            cvalu(wbt)=FM13(j+2:j+5)
            if (cvalu(wbt)(1:1).eq.'0'
     +      .or.cvalu(wbt)(1:1).eq.'1') then
              fvalu(wbti)=0.
            else
     +      if (cvalu(wbt)(1:1).eq.'5'
     +      .or.cvalu(wbt)(1:1).eq.'6') then
              fvalu(wbti)=1.
            else
     +      if (cvalu(wbt)(1:1).eq.'2') then
              fvalu(wbti)=2.
            else
     +      if (cvalu(wbt)(1:1).eq.'7') then
              fvalu(wbti)=3.
            endif
            if (fvalu(wbti).ne.fmiss) then
              fvalu(wbt)=float(tenths(cvalu(wbt)),fmiss)
              if (fvalu(wbt).ne.fmiss) then
                fvalu(wbt)=mod(fvalu(wbt)/10.,100.)
c???ice bulb neg
                if (cvalu(wbt)(1:1).eq.'1'
     +          .or.cvalu(wbt)(1:1).eq.'2'
     +          .or.cvalu(wbt)(1:1).eq.'6'
     +          .or.cvalu(wbt)(1:1).eq.'7') fvalu(wbt)=-fvalu(wbt)
                call indic(cvalu(wbt),fvalu(it))
              else
                fvalu(wbti)=fmiss
              endif
            endif
          endif
          j=j+6
        endif
C ICE c(i)S(i)b(i)D(i)z(i)
        if (FM13(j:j+3).eq.' ICE') then
          if (index(FM13(j+1:),' ').eq.9) j=j-1
          j=j+4
          if (index(FM13(j+1:),' ').ne.6) goto 700
          if (FM13(j+1:j+5).ne.'/////') then
            cvalu(ic1)=FM13(j+1:j+1)
            fvalu(ic1)=float(tens(cvalu(ic1)),fmiss)
            cvalu(ic2)=FM13(j+2:j+2)
            fvalu(ic2)=float(tens(cvalu(ic2)),fmiss)
            cvalu(ic3)=FM13(j+3:j+3)
            fvalu(ic3)=float(tens(cvalu(ic3)),fmiss)
            cvalu(ic4)=FM13(j+4:j+4)
            fvalu(ic4)=float(tens(cvalu(ic4)),fmiss)
            cvalu(ic5)=FM13(j+5:j+5)
            fvalu(ic5)=float(tens(cvalu(ic5)),fmiss)
          endif
          j=j+6
        endif
      endif
C
C SECTION 3
  333 continue
C 333
C 0....
C 1s(n)T(x)T(x)T(x)
C 2s(n)T(n)T(n)T(n)
C 3Ejjj
C 4E'sss
C 5j(1)j(2)j(3)j(4)
C j(5)j(6)j(7)j(8)j(9)
C 6RRRt(R)
C 7R(24)R(24)R(24)R(24)
C 8N(s)Ch(s)h(s)
C 9S(p)S(p)s(p)s(p)
C 8000
C 0....
C 1....
C
C SECTION 4
  444 continue
C 444
C N'C'H'H'C(t)
C
C SECTION 5
  555 continue
C 555
C
  700 imma=core(fvalu,fmiss)//attm1//attm5(fvalu(os-2),fmiss)
      write(*,'(9a)')imma
     +,'99 0 '//head//FM13(:lentrm(FM13))
      do 790 i=1,imono
        if (cvalu(i).ne.' ' .and. fvalu(i).eq.fmiss)
     +  call savstr(abbr(i)//cvalu(i))
  790 continue
      end
C-----------------------------------------------------------------------3456789
      function lentrm(str)
c length of a string minus trailing blanks
      character str*(*)
      do 190 lentrm=len(str),1,-1
        if (str(lentrm:lentrm).ne.' ') return
  190 continue
      end
c-----------------------------------------------------------------------3456789
      character*108 function core(value,miss)
      implicit integer(a-z)
      real value(*),miss
      parameter(YR=1,MO=2,DY=3,HR=4,LAT=5,LON=6,IM=7,ATTC=8,TI=9,LI=10
     +,DS=11,VS=12,NID=13,II=14,ID=15,C1=24,DI=26,D=27,WI=28,W=29,VI=30
     +,VV=31,WW=32,W1=33,SLP=34,A=35,PPP=36,IT=37,AT=38,WBTI=39,WBT=40
     +,DPTI=41,DPT=42,SI=43,SST=44,N=45,NH=46,CL=47,HI=48,H=49,CM=50
     +,CH=51,WD=52,WP=53,WH=54,SD=55,SP=56,SH=57)
c
      core=' '
c
      if (nint(value(YR)).ge.1770 .and.
     +    nint(value(YR)).le.2024) then
        write(core(1:4),'(i4)')nint(value(YR))
      else
        if (value(YR).ne.miss) call saverr('YR',value(YR))
      endif
c
      if (nint(value(MO)).ge.1 .and.
     +    nint(value(MO)).le.12) then
        write(core(5:6),'(i2)')nint(value(MO))
      else
        if (value(MO).ne.miss) call saverr('MO',value(MO))
      endif
c
      if (nint(value(DY)).ge.1 .and.
     +    nint(value(DY)).le.31) then
        write(core(7:8),'(i2)')nint(value(DY))
      else
        if (value(DY).ne.miss) call saverr('DY',value(DY))
      endif
c
      if (nint(value(HR)/0.01).ge.nint(0.00/0.01) .and.
     +    nint(value(HR)/0.01).le.nint(23.99/0.01)) then
        write(core(9:12),'(i4)')nint(value(HR)/0.01)
      else
        if (value(HR).ne.miss) call saverr('HR',value(HR))
      endif
c
      if (nint(value(LAT)/0.01).ge.nint(-90.00/0.01) .and.
     +    nint(value(LAT)/0.01).le.nint(90.00/0.01)) then
        write(core(13:17),'(i5)')nint(value(LAT)/0.01)
      else
        if (value(LAT).ne.miss) call saverr('LAT',value(LAT))
      endif
c
      if (nint(value(LON)/0.01).ge.nint(-179.99/0.01) .and.
     +    nint(value(LON)/0.01).le.nint(359.99/0.01)) then
        write(core(18:23),'(i6)')nint(value(LON)/0.01)
      else
        if (value(LON).ne.miss) call saverr('LON',value(LON))
      endif
c
      if (nint(value(IM)).ge.0 .and.
     +    nint(value(IM)).le.99) then
        write(core(24:25),'(i2)')nint(value(IM))
      else
        if (value(IM).ne.miss) call saverr('IM',value(IM))
      endif
c
      if (nint(value(ATTC)).ge.0 .and.
     +    nint(value(ATTC)).le.9) then
        write(core(26:26),'(i1)')nint(value(ATTC))
      else
        if (value(ATTC).ne.miss) call saverr('ATTC',value(ATTC))
      endif
c
      if (nint(value(TI)).ge.0 .and.
     +    nint(value(TI)).le.3) then
        write(core(27:27),'(i1)')nint(value(TI))
      else
        if (value(TI).ne.miss) call saverr('TI',value(TI))
      endif
c
      if (nint(value(LI)).ge.0 .and.
     +    nint(value(LI)).le.6) then
        write(core(28:28),'(i1)')nint(value(LI))
      else
        if (value(LI).ne.miss) call saverr('LI',value(LI))
      endif
c
      if (nint(value(DS)).ge.0 .and.
     +    nint(value(DS)).le.9) then
        write(core(29:29),'(i1)')nint(value(DS))
      else
        if (value(DS).ne.miss) call saverr('DS',value(DS))
      endif
c
      if (nint(value(VS)).ge.0 .and.
     +    nint(value(VS)).le.9) then
        write(core(30:30),'(i1)')nint(value(VS))
      else
        if (value(VS).ne.miss) call saverr('VS',value(VS))
      endif
c
      if (nint(value(NID)).ge.0 .and.
     +    nint(value(NID)).le.99) then
        write(core(31:32),'(i2)')nint(value(NID))
      else
        if (value(NID).ne.miss) call saverr('NID',value(NID))
      endif
c
      if (nint(value(II)).ge.0 .and.
     +    nint(value(II)).le.10) then
        write(core(33:34),'(i2)')nint(value(II))
      else
        if (value(II).ne.miss) call saverr('II',value(II))
      endif
c
      if (nint(value(ID)).ge.ichar(' ') .and.
     +    nint(value(ID)).le.ichar('~')) then
        write(core(35:35),'(a1)')char(nint(value(ID)))
      else
        if (value(ID).ne.miss) call saverr('ID',value(ID))
      endif
c
      if (nint(value(ID+1)).ge.ichar(' ') .and.
     +    nint(value(ID+1)).le.ichar('~')) then
        write(core(36:36),'(a1)')char(nint(value(ID+1)))
      else
        if (value(ID+1).ne.miss) call saverr('ID+1',value(ID+1))
      endif
c
      if (nint(value(ID+2)).ge.ichar(' ') .and.
     +    nint(value(ID+2)).le.ichar('~')) then
        write(core(37:37),'(a1)')char(nint(value(ID+2)))
      else
        if (value(ID+2).ne.miss) call saverr('ID+2',value(ID+2))
      endif
c
      if (nint(value(ID+3)).ge.ichar(' ') .and.
     +    nint(value(ID+3)).le.ichar('~')) then
        write(core(38:38),'(a1)')char(nint(value(ID+3)))
      else
        if (value(ID+3).ne.miss) call saverr('ID+3',value(ID+3))
      endif
c
      if (nint(value(ID+4)).ge.ichar(' ') .and.
     +    nint(value(ID+4)).le.ichar('~')) then
        write(core(39:39),'(a1)')char(nint(value(ID+4)))
      else
        if (value(ID+4).ne.miss) call saverr('ID+4',value(ID+4))
      endif
c
      if (nint(value(ID+5)).ge.ichar(' ') .and.
     +    nint(value(ID+5)).le.ichar('~')) then
        write(core(40:40),'(a1)')char(nint(value(ID+5)))
      else
        if (value(ID+5).ne.miss) call saverr('ID+5',value(ID+5))
      endif
c
      if (nint(value(ID+6)).ge.ichar(' ') .and.
     +    nint(value(ID+6)).le.ichar('~')) then
        write(core(41:41),'(a1)')char(nint(value(ID+6)))
      else
        if (value(ID+6).ne.miss) call saverr('ID+6',value(ID+6))
      endif
c
      if (nint(value(ID+7)).ge.ichar(' ') .and.
     +    nint(value(ID+7)).le.ichar('~')) then
        write(core(42:42),'(a1)')char(nint(value(ID+7)))
      else
        if (value(ID+7).ne.miss) call saverr('ID+7',value(ID+7))
      endif
c
      if (nint(value(ID+8)).ge.ichar(' ') .and.
     +    nint(value(ID+8)).le.ichar('~')) then
        write(core(43:43),'(a1)')char(nint(value(ID+8)))
      else
        if (value(ID+8).ne.miss) call saverr('ID+8',value(ID+8))
      endif
c
      if (nint(value(C1)).ge.ichar('A') .and.
     +    nint(value(C1)).le.ichar('Z') .or.
     +    nint(value(C1)).ge.ichar('0') .and.
     +    nint(value(C1)).le.ichar('9')) then
        write(core(44:44),'(a1)')char(nint(value(C1)))
      else
        if (value(C1).ne.miss) call saverr('C1',value(C1))
      endif
c
      if (nint(value(C1+1)).ge.ichar('A') .and.
     +    nint(value(C1+1)).le.ichar('Z') .or.
     +    nint(value(C1+1)).ge.ichar('0') .and.
     +    nint(value(C1+1)).le.ichar('9')) then
        write(core(45:45),'(a1)')char(nint(value(C1+1)))
      else
        if (value(C1+1).ne.miss) call saverr('C1+1',value(C1+1))
      endif
c
      if (nint(value(DI)).ge.0 .and.
     +    nint(value(DI)).le.6) then
        write(core(46:46),'(i1)')nint(value(DI))
      else
        if (value(DI).ne.miss) call saverr('DI',value(DI))
      endif
c
      if (nint(value(D)).ge.1 .and.
     +    nint(value(D)).le.362) then
        write(core(47:49),'(i3)')nint(value(D))
      else
        if (value(D).ne.miss) call saverr('D',value(D))
      endif
c
      if (nint(value(WI)).ge.0 .and.
     +    nint(value(WI)).le.8) then
        write(core(50:50),'(i1)')nint(value(WI))
      else
        if (value(WI).ne.miss) call saverr('WI',value(WI))
      endif
c
      if (nint(value(W)/0.1).ge.nint(0.0/0.1) .and.
     +    nint(value(W)/0.1).le.nint(99.9/0.1)) then
        write(core(51:53),'(i3)')nint(value(W)/0.1)
      else
        if (value(W).ne.miss) call saverr('W',value(W))
      endif
c
      if (nint(value(VI)).ge.0 .and.
     +    nint(value(VI)).le.2) then
        write(core(54:54),'(i1)')nint(value(VI))
      else
        if (value(VI).ne.miss) call saverr('VI',value(VI))
      endif
c
      if (nint(value(VV)).ge.90 .and.
     +    nint(value(VV)).le.99) then
        write(core(55:56),'(i2)')nint(value(VV))
      else
        if (value(VV).ne.miss) call saverr('VV',value(VV))
      endif
c
      if (nint(value(WW)).ge.0 .and.
     +    nint(value(WW)).le.99) then
        write(core(57:58),'(i2)')nint(value(WW))
      else
        if (value(WW).ne.miss) call saverr('WW',value(WW))
      endif
c
      if (nint(value(W1)).ge.0 .and.
     +    nint(value(W1)).le.9) then
        write(core(59:59),'(i1)')nint(value(W1))
      else
        if (value(W1).ne.miss) call saverr('W1',value(W1))
      endif
c
      if (nint(value(SLP)/0.1).ge.nint(870.0/0.1) .and.
     +    nint(value(SLP)/0.1).le.nint(1074.6/0.1)) then
        write(core(60:64),'(i5)')nint(value(SLP)/0.1)
      else
        if (value(SLP).ne.miss) call saverr('SLP',value(SLP))
      endif
c
      if (nint(value(A)).ge.0 .and.
     +    nint(value(A)).le.8) then
        write(core(65:65),'(i1)')nint(value(A))
      else
        if (value(A).ne.miss) call saverr('A',value(A))
      endif
c
      if (nint(value(PPP)/0.1).ge.nint(0.0/0.1) .and.
     +    nint(value(PPP)/0.1).le.nint(51.0/0.1)) then
        write(core(66:68),'(i3)')nint(value(PPP)/0.1)
      else
        if (value(PPP).ne.miss) call saverr('PPP',value(PPP))
      endif
c
      if (nint(value(IT)).ge.0 .and.
     +    nint(value(IT)).le.9) then
        write(core(69:69),'(i1)')nint(value(IT))
      else
        if (value(IT).ne.miss) call saverr('IT',value(IT))
      endif
c
      if (nint(value(AT)/0.1).ge.nint(-99.9/0.1) .and.
     +    nint(value(AT)/0.1).le.nint(99.9/0.1)) then
        write(core(70:73),'(i4)')nint(value(AT)/0.1)
      else
        if (value(AT).ne.miss) call saverr('AT',value(AT))
      endif
c
      if (nint(value(WBTI)).ge.0 .and.
     +    nint(value(WBTI)).le.3) then
        write(core(74:74),'(i1)')nint(value(WBTI))
      else
        if (value(WBTI).ne.miss) call saverr('WBTI',value(WBTI))
      endif
c
      if (nint(value(WBT)/0.1).ge.nint(-99.9/0.1) .and.
     +    nint(value(WBT)/0.1).le.nint(99.9/0.1)) then
        write(core(75:78),'(i4)')nint(value(WBT)/0.1)
      else
        if (value(WBT).ne.miss) call saverr('WBT',value(WBT))
      endif
c
      if (nint(value(DPTI)).ge.0 .and.
     +    nint(value(DPTI)).le.3) then
        write(core(79:79),'(i1)')nint(value(DPTI))
      else
        if (value(DPTI).ne.miss) call saverr('DPTI',value(DPTI))
      endif
c
      if (nint(value(DPT)/0.1).ge.nint(-99.9/0.1) .and.
     +    nint(value(DPT)/0.1).le.nint(99.9/0.1)) then
        write(core(80:83),'(i4)')nint(value(DPT)/0.1)
      else
        if (value(DPT).ne.miss) call saverr('DPT',value(DPT))
      endif
c
      if (nint(value(SI)).ge.0 .and.
     +    nint(value(SI)).le.12) then
        write(core(84:85),'(i2)')nint(value(SI))
      else
        if (value(SI).ne.miss) call saverr('SI',value(SI))
      endif
c
      if (nint(value(SST)/0.1).ge.nint(-99.9/0.1) .and.
     +    nint(value(SST)/0.1).le.nint(99.9/0.1)) then
        write(core(86:89),'(i4)')nint(value(SST)/0.1)
      else
        if (value(SST).ne.miss) call saverr('SST',value(SST))
      endif
c
      if (nint(value(N)).ge.0 .and.
     +    nint(value(N)).le.9) then
        write(core(90:90),'(i1)')nint(value(N))
      else
        if (value(N).ne.miss) call saverr('N',value(N))
      endif
c
      if (nint(value(NH)).ge.0 .and.
     +    nint(value(NH)).le.9) then
        write(core(91:91),'(i1)')nint(value(NH))
      else
        if (value(NH).ne.miss) call saverr('NH',value(NH))
      endif
c
      if (nint(value(CL)).ge.0 .and.
     +    nint(value(CL)).le.10) then
        write(core(92:92),'(z1)')nint(value(CL))
      else
        if (value(CL).ne.miss) call saverr('CL',value(CL))
      endif
c
      if (nint(value(HI)).ge.0 .and.
     +    nint(value(HI)).le.1) then
        write(core(93:93),'(i1)')nint(value(HI))
      else
        if (value(HI).ne.miss) call saverr('HI',value(HI))
      endif
c
      if (nint(value(H)).ge.0 .and.
     +    nint(value(H)).le.10) then
        write(core(94:94),'(z1)')nint(value(H))
      else
        if (value(H).ne.miss) call saverr('H',value(H))
      endif
c
      if (nint(value(CM)).ge.0 .and.
     +    nint(value(CM)).le.10) then
        write(core(95:95),'(z1)')nint(value(CM))
      else
        if (value(CM).ne.miss) call saverr('CM',value(CM))
      endif
c
      if (nint(value(CH)).ge.0 .and.
     +    nint(value(CH)).le.10) then
        write(core(96:96),'(z1)')nint(value(CH))
      else
        if (value(CH).ne.miss) call saverr('CH',value(CH))
      endif
c
      if (nint(value(WD)).ge.0 .and.
     +    nint(value(WD)).le.38) then
        write(core(97:98),'(i2)')nint(value(WD))
      else
        if (value(WD).ne.miss) call saverr('WD',value(WD))
      endif
c
      if (nint(value(WP)).ge.0 .and.
     +    nint(value(WP)).le.30 .or.
     +    nint(value(WP)).eq.99) then
        write(core(99:100),'(i2)')nint(value(WP))
      else
        if (value(WP).ne.miss) call saverr('WP',value(WP))
      endif
c
      if (nint(value(WH)).ge.0 .and.
     +    nint(value(WH)).le.99) then
        write(core(101:102),'(i2)')nint(value(WH))
      else
        if (value(WH).ne.miss) call saverr('WH',value(WH))
      endif
c
      if (nint(value(SD)).ge.0 .and.
     +    nint(value(SD)).le.38) then
        write(core(103:104),'(i2)')nint(value(SD))
      else
        if (value(SD).ne.miss) call saverr('SD',value(SD))
      endif
c
      if (nint(value(SP)).ge.0 .and.
     +    nint(value(SP)).le.30 .or.
     +    nint(value(SP)).eq.99) then
        write(core(105:106),'(i2)')nint(value(SP))
      else
        if (value(SP).ne.miss) call saverr('SP',value(SP))
      endif
c
      if (nint(value(SH)).ge.0 .and.
     +    nint(value(SH)).le.99) then
        write(core(107:108),'(i2)')nint(value(SH))
      else
        if (value(SH).ne.miss) call saverr('SH',value(SH))
      endif
      end
c-----------------------------------------------------------------------3456789
      character*94 function attm5(value,miss)
      implicit integer(a-z)
      real value(*),miss
      parameter(OS=3,OP=4,FM=5,IMMV=6,IX=7,W2=8,WMI=9,SD2=10,SP2=11
     +,SH2=12,IS=13,ES=14,RS=15,IC1=16,IC2=17,IC3=18,IC4=19,IC5=20,IR=21
     +,RRR=22,TR=23,NU=24,QCI=25,QI1=26,QI2=27,QI3=28,QI4=29,QI5=30
     +,QI6=31,QI7=32,QI8=33,QI9=34,QI10=35,QI11=36,QI12=37,QI13=38
     +,QI14=39,QI15=40,QI16=41,QI17=42,QI18=43,QI19=44,QI20=45,QI21=46
     +,HDG=47,COG=48,SOG=49,SLL=50,SLHH=51,RWD=52,RWS=53,QI22=54,QI23=55
     +,QI24=56,QI25=57,QI26=58,QI27=59,QI28=60,QI29=61,RH=62,RHI=63
     +,AWSI=64,IMONO=65)
      character*36 :: base36='0123456789ABCDEFGHIJKLMNOPQRSTUVWXYZ'
c
      attm5=' 594'
c
      if (value(OS).ne.miss) then
        if (nint(value(OS)).ge.0 .and.
     +      nint(value(OS)).le.6) then
          write(attm5(5:5),'(i1)')nint(value(OS))
        else
          call saverr('OS',value(OS))
        endif
      endif
c
      if (value(OP).ne.miss) then
        if (nint(value(OP)).ge.0 .and.
     +      nint(value(OP)).le.9) then
          write(attm5(6:6),'(i1)')nint(value(OP))
        else
          call saverr('OP',value(OP))
        endif
      endif
c
      if (value(FM).ne.miss) then
        if (nint(value(FM)).ge.0 .and.
     +      nint(value(FM)).le.35) then
          write(attm5(7:7),'(a1)')base36(nint(value(FM))+1:)
        else
          call saverr('FM',value(FM))
        endif
      endif
c
      if (value(IMMV).ne.miss) then
        if (nint(value(IMMV)).ge.0 .and.
     +      nint(value(IMMV)).le.35) then
          write(attm5(8:8),'(a1)')base36(nint(value(IMMV))+1:)
        else
          call saverr('IMMV',value(IMMV))
        endif
      endif
c
      if (value(IX).ne.miss) then
        if (nint(value(IX)).ge.1 .and.
     +      nint(value(IX)).le.7) then
          write(attm5(9:9),'(i1)')nint(value(IX))
        else
          call saverr('IX',value(IX))
        endif
      endif
c
      if (value(W2).ne.miss) then
        if (nint(value(W2)).ge.0 .and.
     +      nint(value(W2)).le.9) then
          write(attm5(10:10),'(i1)')nint(value(W2))
        else
          call saverr('W2',value(W2))
        endif
      endif
c
      if (value(WMI).ne.miss) then
        if (nint(value(WMI)).ge.0 .and.
     +      nint(value(WMI)).le.9) then
          write(attm5(11:11),'(i1)')nint(value(WMI))
        else
          call saverr('WMI',value(WMI))
        endif
      endif
c
      if (value(SD2).ne.miss) then
        if (nint(value(SD2)).ge.0 .and.
     +      nint(value(SD2)).le.38) then
          write(attm5(12:13),'(i2)')nint(value(SD2))
        else
          call saverr('SD2',value(SD2))
        endif
      endif
c
      if (value(SP2).ne.miss) then
        if (nint(value(SP2)).ge.0 .and.
     +      nint(value(SP2)).le.30 .or.
     +      nint(value(SP2)).eq.99) then
          write(attm5(14:15),'(i2)')nint(value(SP2))
        else
          call saverr('SP2',value(SP2))
        endif
      endif
c
      if (value(SH2).ne.miss) then
        if (nint(value(SH2)).ge.0 .and.
     +      nint(value(SH2)).le.99) then
          write(attm5(16:17),'(i2)')nint(value(SH2))
        else
          call saverr('SH2',value(SH2))
        endif
      endif
c
      if (value(IS).ne.miss) then
        if (nint(value(IS)).ge.1 .and.
     +      nint(value(IS)).le.5) then
          write(attm5(18:18),'(i1)')nint(value(IS))
        else
          call saverr('IS',value(IS))
        endif
      endif
c
      if (value(ES).ne.miss) then
        if (nint(value(ES)).ge.0 .and.
     +      nint(value(ES)).le.99) then
          write(attm5(19:20),'(i2)')nint(value(ES))
        else
          call saverr('ES',value(ES))
        endif
      endif
c
      if (value(RS).ne.miss) then
        if (nint(value(RS)).ge.0 .and.
     +      nint(value(RS)).le.4) then
          write(attm5(21:21),'(i1)')nint(value(RS))
        else
          call saverr('RS',value(RS))
        endif
      endif
c
      if (value(IC1).ne.miss) then
        if (nint(value(IC1)).ge.0 .and.
     +      nint(value(IC1)).le.10) then
          write(attm5(22:22),'(z1)')nint(value(IC1))
        else
          call saverr('IC1',value(IC1))
        endif
      endif
c
      if (value(IC2).ne.miss) then
        if (nint(value(IC2)).ge.0 .and.
     +      nint(value(IC2)).le.10) then
          write(attm5(23:23),'(z1)')nint(value(IC2))
        else
          call saverr('IC2',value(IC2))
        endif
      endif
c
      if (value(IC3).ne.miss) then
        if (nint(value(IC3)).ge.0 .and.
     +      nint(value(IC3)).le.10) then
          write(attm5(24:24),'(z1)')nint(value(IC3))
        else
          call saverr('IC3',value(IC3))
        endif
      endif
c
      if (value(IC4).ne.miss) then
        if (nint(value(IC4)).ge.0 .and.
     +      nint(value(IC4)).le.10) then
          write(attm5(25:25),'(z1)')nint(value(IC4))
        else
          call saverr('IC4',value(IC4))
        endif
      endif
c
      if (value(IC5).ne.miss) then
        if (nint(value(IC5)).ge.0 .and.
     +      nint(value(IC5)).le.10) then
          write(attm5(26:26),'(z1)')nint(value(IC5))
        else
          call saverr('IC5',value(IC5))
        endif
      endif
c
      if (value(IR).ne.miss) then
        if (nint(value(IR)).ge.0 .and.
     +      nint(value(IR)).le.4) then
          write(attm5(27:27),'(i1)')nint(value(IR))
        else
          call saverr('IR',value(IR))
        endif
      endif
c
      if (value(RRR).ne.miss) then
        if (nint(value(RRR)).ge.0 .and.
     +      nint(value(RRR)).le.999) then
          write(attm5(28:30),'(i3)')nint(value(RRR))
        else
          call saverr('RRR',value(RRR))
        endif
      endif
c
      if (value(TR).ne.miss) then
        if (nint(value(TR)).ge.1 .and.
     +      nint(value(TR)).le.9) then
          write(attm5(31:31),'(i1)')nint(value(TR))
        else
          call saverr('TR',value(TR))
        endif
      endif
c
      if (value(NU).ne.miss) then
        if (nint(value(NU)).ge.ichar(' ') .and.
     +      nint(value(NU)).le.ichar('~')) then
          write(attm5(32:32),'(a1)')char(nint(value(NU)))
        else
          call saverr('NU',value(NU))
        endif
      endif
c
      if (value(QCI).ne.miss) then
        if (nint(value(QCI)).ge.0 .and.
     +      nint(value(QCI)).le.9) then
          write(attm5(33:33),'(i1)')nint(value(QCI))
        else
          call saverr('QCI',value(QCI))
        endif
      endif
c
      if (value(QI1).ne.miss) then
        if (nint(value(QI1)).ge.0 .and.
     +      nint(value(QI1)).le.9) then
          write(attm5(34:34),'(i1)')nint(value(QI1))
        else
          call saverr('QI1',value(QI1))
        endif
      endif
c
      if (value(QI2).ne.miss) then
        if (nint(value(QI2)).ge.0 .and.
     +      nint(value(QI2)).le.9) then
          write(attm5(35:35),'(i1)')nint(value(QI2))
        else
          call saverr('QI2',value(QI2))
        endif
      endif
c
      if (value(QI3).ne.miss) then
        if (nint(value(QI3)).ge.0 .and.
     +      nint(value(QI3)).le.9) then
          write(attm5(36:36),'(i1)')nint(value(QI3))
        else
          call saverr('QI3',value(QI3))
        endif
      endif
c
      if (value(QI4).ne.miss) then
        if (nint(value(QI4)).ge.0 .and.
     +      nint(value(QI4)).le.9) then
          write(attm5(37:37),'(i1)')nint(value(QI4))
        else
          call saverr('QI4',value(QI4))
        endif
      endif
c
      if (value(QI5).ne.miss) then
        if (nint(value(QI5)).ge.0 .and.
     +      nint(value(QI5)).le.9) then
          write(attm5(38:38),'(i1)')nint(value(QI5))
        else
          call saverr('QI5',value(QI5))
        endif
      endif
c
      if (value(QI6).ne.miss) then
        if (nint(value(QI6)).ge.0 .and.
     +      nint(value(QI6)).le.9) then
          write(attm5(39:39),'(i1)')nint(value(QI6))
        else
          call saverr('QI6',value(QI6))
        endif
      endif
c
      if (value(QI7).ne.miss) then
        if (nint(value(QI7)).ge.0 .and.
     +      nint(value(QI7)).le.9) then
          write(attm5(40:40),'(i1)')nint(value(QI7))
        else
          call saverr('QI7',value(QI7))
        endif
      endif
c
      if (value(QI8).ne.miss) then
        if (nint(value(QI8)).ge.0 .and.
     +      nint(value(QI8)).le.9) then
          write(attm5(41:41),'(i1)')nint(value(QI8))
        else
          call saverr('QI8',value(QI8))
        endif
      endif
c
      if (value(QI9).ne.miss) then
        if (nint(value(QI9)).ge.0 .and.
     +      nint(value(QI9)).le.9) then
          write(attm5(42:42),'(i1)')nint(value(QI9))
        else
          call saverr('QI9',value(QI9))
        endif
      endif
c
      if (value(QI10).ne.miss) then
        if (nint(value(QI10)).ge.0 .and.
     +      nint(value(QI10)).le.9) then
          write(attm5(43:43),'(i1)')nint(value(QI10))
        else
          call saverr('QI10',value(QI10))
        endif
      endif
c
      if (value(QI11).ne.miss) then
        if (nint(value(QI11)).ge.0 .and.
     +      nint(value(QI11)).le.9) then
          write(attm5(44:44),'(i1)')nint(value(QI11))
        else
          call saverr('QI11',value(QI11))
        endif
      endif
c
      if (value(QI12).ne.miss) then
        if (nint(value(QI12)).ge.0 .and.
     +      nint(value(QI12)).le.9) then
          write(attm5(45:45),'(i1)')nint(value(QI12))
        else
          call saverr('QI12',value(QI12))
        endif
      endif
c
      if (value(QI13).ne.miss) then
        if (nint(value(QI13)).ge.0 .and.
     +      nint(value(QI13)).le.9) then
          write(attm5(46:46),'(i1)')nint(value(QI13))
        else
          call saverr('QI13',value(QI13))
        endif
      endif
c
      if (value(QI14).ne.miss) then
        if (nint(value(QI14)).ge.0 .and.
     +      nint(value(QI14)).le.9) then
          write(attm5(47:47),'(i1)')nint(value(QI14))
        else
          call saverr('QI14',value(QI14))
        endif
      endif
c
      if (value(QI15).ne.miss) then
        if (nint(value(QI15)).ge.0 .and.
     +      nint(value(QI15)).le.9) then
          write(attm5(48:48),'(i1)')nint(value(QI15))
        else
          call saverr('QI15',value(QI15))
        endif
      endif
c
      if (value(QI16).ne.miss) then
        if (nint(value(QI16)).ge.0 .and.
     +      nint(value(QI16)).le.9) then
          write(attm5(49:49),'(i1)')nint(value(QI16))
        else
          call saverr('QI16',value(QI16))
        endif
      endif
c
      if (value(QI17).ne.miss) then
        if (nint(value(QI17)).ge.0 .and.
     +      nint(value(QI17)).le.9) then
          write(attm5(50:50),'(i1)')nint(value(QI17))
        else
          call saverr('QI17',value(QI17))
        endif
      endif
c
      if (value(QI18).ne.miss) then
        if (nint(value(QI18)).ge.0 .and.
     +      nint(value(QI18)).le.9) then
          write(attm5(51:51),'(i1)')nint(value(QI18))
        else
          call saverr('QI18',value(QI18))
        endif
      endif
c
      if (value(QI19).ne.miss) then
        if (nint(value(QI19)).ge.0 .and.
     +      nint(value(QI19)).le.9) then
          write(attm5(52:52),'(i1)')nint(value(QI19))
        else
          call saverr('QI19',value(QI19))
        endif
      endif
c
      if (value(QI20).ne.miss) then
        if (nint(value(QI20)).ge.0 .and.
     +      nint(value(QI20)).le.9) then
          write(attm5(53:53),'(i1)')nint(value(QI20))
        else
          call saverr('QI20',value(QI20))
        endif
      endif
c
      if (value(QI21).ne.miss) then
        if (nint(value(QI21)).ge.0 .and.
     +      nint(value(QI21)).le.9) then
          write(attm5(54:54),'(i1)')nint(value(QI21))
        else
          call saverr('QI21',value(QI21))
        endif
      endif
c
      if (value(HDG).ne.miss) then
        if (nint(value(HDG)).ge.0 .and.
     +      nint(value(HDG)).le.360) then
          write(attm5(55:57),'(i3)')nint(value(HDG))
        else
          call saverr('HDG',value(HDG))
        endif
      endif
c
      if (value(COG).ne.miss) then
        if (nint(value(COG)).ge.0 .and.
     +      nint(value(COG)).le.360) then
          write(attm5(58:60),'(i3)')nint(value(COG))
        else
          call saverr('COG',value(COG))
        endif
      endif
c
      if (value(SOG).ne.miss) then
        if (nint(value(SOG)).ge.0 .and.
     +      nint(value(SOG)).le.99) then
          write(attm5(61:62),'(i2)')nint(value(SOG))
        else
          call saverr('SOG',value(SOG))
        endif
      endif
c
      if (value(SLL).ne.miss) then
        if (nint(value(SLL)).ge.0 .and.
     +      nint(value(SLL)).le.99) then
          write(attm5(63:64),'(i2)')nint(value(SLL))
        else
          call saverr('SLL',value(SLL))
        endif
      endif
c
      if (value(SLHH).ne.miss) then
        if (nint(value(SLHH)).ge.-99 .and.
     +      nint(value(SLHH)).le.99) then
          write(attm5(65:67),'(i3)')nint(value(SLHH))
        else
          call saverr('SLHH',value(SLHH))
        endif
      endif
c
      if (value(RWD).ne.miss) then
        if (nint(value(RWD)).ge.1 .and.
     +      nint(value(RWD)).le.362) then
          write(attm5(68:70),'(i3)')nint(value(RWD))
        else
          call saverr('RWD',value(RWD))
        endif
      endif
c
      if (value(RWS).ne.miss) then
        if (nint(value(RWS)/0.1).ge.nint(0.0/0.1) .and.
     +      nint(value(RWS)/0.1).le.nint(99.9/0.1)) then
          write(attm5(71:73),'(i3)')nint(value(RWS)/0.1)
        else
          call saverr('RWS',value(RWS))
        endif
      endif
c
      if (value(QI22).ne.miss) then
        if (nint(value(QI22)).ge.0 .and.
     +      nint(value(QI22)).le.9) then
          write(attm5(74:74),'(i1)')nint(value(QI22))
        else
          call saverr('QI22',value(QI22))
        endif
      endif
c
      if (value(QI23).ne.miss) then
        if (nint(value(QI23)).ge.0 .and.
     +      nint(value(QI23)).le.9) then
          write(attm5(75:75),'(i1)')nint(value(QI23))
        else
          call saverr('QI23',value(QI23))
        endif
      endif
c
      if (value(QI24).ne.miss) then
        if (nint(value(QI24)).ge.0 .and.
     +      nint(value(QI24)).le.9) then
          write(attm5(76:76),'(i1)')nint(value(QI24))
        else
          call saverr('QI24',value(QI24))
        endif
      endif
c
      if (value(QI25).ne.miss) then
        if (nint(value(QI25)).ge.0 .and.
     +      nint(value(QI25)).le.9) then
          write(attm5(77:77),'(i1)')nint(value(QI25))
        else
          call saverr('QI25',value(QI25))
        endif
      endif
c
      if (value(QI26).ne.miss) then
        if (nint(value(QI26)).ge.0 .and.
     +      nint(value(QI26)).le.9) then
          write(attm5(78:78),'(i1)')nint(value(QI26))
        else
          call saverr('QI26',value(QI26))
        endif
      endif
c
      if (value(QI27).ne.miss) then
        if (nint(value(QI27)).ge.0 .and.
     +      nint(value(QI27)).le.9) then
          write(attm5(79:79),'(i1)')nint(value(QI27))
        else
          call saverr('QI27',value(QI27))
        endif
      endif
c
      if (value(QI28).ne.miss) then
        if (nint(value(QI28)).ge.0 .and.
     +      nint(value(QI28)).le.9) then
          write(attm5(80:80),'(i1)')nint(value(QI28))
        else
          call saverr('QI28',value(QI28))
        endif
      endif
c
      if (value(QI29).ne.miss) then
        if (nint(value(QI29)).ge.0 .and.
     +      nint(value(QI29)).le.9) then
          write(attm5(81:81),'(i1)')nint(value(QI29))
        else
          call saverr('QI29',value(QI29))
        endif
      endif
c
      if (value(RH).ne.miss) then
        if (nint(value(RH)/0.1).ge.nint(0.0/0.1) .and.
     +      nint(value(RH)/0.1).le.nint(100.0/0.1)) then
          write(attm5(82:85),'(i4)')nint(value(RH)/0.1)
        else
          call saverr('RH',value(RH))
        endif
      endif
c
      if (value(RHI).ne.miss) then
        if (nint(value(RHI)).ge.0 .and.
     +      nint(value(RHI)).le.4) then
          write(attm5(86:86),'(i1)')nint(value(RHI))
        else
          call saverr('RHI',value(RHI))
        endif
      endif
c
      if (value(AWSI).ne.miss) then
        if (nint(value(AWSI)).ge.0 .and.
     +      nint(value(AWSI)).le.2) then
          write(attm5(87:87),'(i1)')nint(value(AWSI))
        else
          call saverr('AWSI',value(AWSI))
        endif
      endif
c
      if (value(IMONO).ne.miss) then
        if (nint(value(IMONO)).ge.0 .and.
     +      nint(value(IMONO)).le.9999999) then
          write(attm5(88:94),'(i7)')nint(value(IMONO))
        else
          call saverr('IMONO',value(IMONO))
        endif
      endif
      end
c-----------------------------------------------------------------------3456789
      subroutine saverr(abbr,value)
      character*(*) abbr
      real value
      character*40 str
c
      write(str,'(a8,f16.3)')abbr//'        ',value
      call savstr(str)
      return
c
      entry prnerr()
      call prnstr
      end
c-----------------------------------------------------------------------3456789
      subroutine savstr(str)
      implicit integer(a-e,g-z)
      parameter(nmax=9999)
      character arr1(nmax)*40,str*(*)
      dimension arr2(nmax)
      data n/0/
      save arr1,arr2,n
c
      do 190 j=1,n
        if(str.ne.arr1(j))goto 190
        arr2(j)=arr2(j)+1
        return
  190 continue
      if(n+1.gt.nmax)stop 'nmax'
      n=n+1
      arr1(n)=str
      arr2(n)=1
      return
c
      entry prnstr
      do 290 j=1,n
        jmin=j
        do 280 j2=j+1,n
          if(llt(arr1(j2),arr1(jmin)))jmin=j2
  280   continue
        write(8,'(1x,a,t31,i9)')arr1(jmin),arr2(jmin)
        if(jmin.ne.j)then
          arr1(jmin)=arr1(j)
          arr2(jmin)=arr2(j)
        endif
  290 continue
      n=0
      end
c-----------------------------------------------------------------------
        FUNCTION DWPT(T,RH)

C       INCLUDE 'LIB_DEV:[GUDOC]EDFVAXBOX.FOR/LIST'
C       Baker, Schlatter  17-MAY-1982     Original version.

C   THIS FUNCTION RETURNS THE DEW POINT (CELSIUS) GIVEN THE TEMPERATURE
C   (CELSIUS) AND RELATIVE HUMIDITY (%). THE FORMULA IS USED IN THE
C   PROCESSING OF U.S. RAWINSONDE DATA AND IS REFERENCED IN PARRY, H.
C   DEAN, 1969: "THE SEMIAUTOMATIC COMPUTATION OF RAWINSONDES,"
C   TECHNICAL MEMORANDUM WBTM EDL 10, U.S. DEPARTMENT OF COMMERCE,
C   ENVIRONMENTAL SCIENCE SERVICES ADMINISTRATION, WEATHER BUREAU,
C   OFFICE OF SYSTEMS DEVELOPMENT, EQUIPMENT DEVELOPMENT LABORATORY,
C   SILVER SPRING, MD (OCTOBER), PAGE 9 AND PAGE II-4, LINE 460.

        X = 1.-0.01*RH

C   COMPUTE DEW POINT DEPRESSION.

        DPD =(14.55+0.114*T)*X+((2.5+0.007*T)*X)**3+(15.9+0.117*T)*X**14
        DWPT = T-DPD
        RETURN
        END
EOR
rm a.out
gfortran p.f lmrlib.o
./a.out