program readm c c read cac summary of month data (format AB) c character nb*120 lout=6 nr=0 print 1000 1000 format( * ' Variables numbered as in format ',//, * ' 1 2 3 4 5 6 7 8 9 10 ' * ' 11 12 13 14 15 16 17 18 19 20 21 22 23' * ,/, * '---- ----- ---- ------ ------- ------ ------ -- ------ -- --' * '---- -- ------ -- ------ -- -------- -- -- -- --- --- --------') 5 read(5,'(a)',err=50,end=90) nb nce=0 nr=nr+1 c c*** set bad elevation to missing if(nb(28:33).eq.'******') nb(28:33)=' 9999.' c read(nb,10)iyrmo,idwmo,call,xlat,xlon,elev,tmen,nump,tmax, 2 maxf,tmin,minf,tmxmx,mxdy,tmnmn,mndy,pcpt,npcpdy,npcptr,ndly, 3 nsyn,nwpcp,pcpsum 10 format(i4,1x,i5,a4,f6.2,f7.2,f6.0,f6.1,i2,f6.1,i2,f6.1, 2 i2,f6.1,1x,i2,f6.1,1x,i2f8.1,1x,3i2,2i3,f8.1) c select data records of interest c c if(iyr.lt.80) go to 5 c if(idwmo.lt.67774 .or. idwmo.gt.67775) go to 5 c if(xlat.lt.6. or.xlat.gt.10.) go to 5 c if(xlon.lt. 75. .or. xlon.gt. 85.) go to 5 c c print variables c c iyr YEAR Characters 1-2 I2 c imo MONTH Characters 3-4 I2 c FILLER Character 5 1X c idwmo WMO-STATION-ID Characters 6-10 I5 c call AIRWAYS-CALLSIGN Characters 11-14 A4 c xlat LATITUDE-CODE Characters 15-20 F6.2 c xlon LONGITUDE-CODE Characters 21-27 F7.2 c elev ELEVATION Characters 28-33 F6.0 c tmen MEAN-MONTHLY-TEMP Characters 34-39 F6.1 c nump NMBR-PAIRED-MAX/MIN TEMPS Characters 40-41 I2 c tmax MEAN-MAX-TEMP Characters 42-47 F6.1 c maxf NMBR-MAX-TEMPS-REPORTED Characters 48-49 I2 c tmin MEAN-MIN-TEMP Characters 50-55 F6.1 c minf NMBR-MIN-TEMPS-REPORTED Characters 56-57 I2 c tmxmx HIGHEST-MAX-TEMP Characters 58-63 F6.1 c FILLER Character 64 1X c mxdy DAY-HIGH-MAX-TEMP-RPTD Characters 65-66 I2 c tmnmn LOWEST-MIN-TEMP Characters 67-72 F6.1 c FILLER Character 73 1X c mndy DAY-LOW-MIN-TEMP-RPTD Characters 74-75 I2 c pcpt SUM-OF-DAILY-PRECIP Characters 76-83 F8.1 c FILLER Character 84 1X c npcpdy NMBR-DAYS-PRECIP-REPORTED Characters 85-86 I2 c npcptr NMBR-OF-TRACE-FLAGS-SET-TO-T Characters 87-88 I2 c ndly NMBR-DAILY-SUMMARY-REPORTS Characters 89-90 I2 c nsyn NMBR-SYNOPTIC-REPORTS-RECEIVED Characters 91-93 I3 c nwpcp NMBR-REPORTS-PRECIP-RPTD-PAST/PRSNT-WEA Characters 94-96 I3 c pcpsum SUM-DAILY-PRECIP-AMOUNTS Characters 97-104 F8.1 c c complete explanations of variables is available in format file c print 1003,iyrmo,idwmo,call,xlat,xlon,elev,tmen,nump,tmax, 2 maxf,tmin,minf,tmxmx,mxdy,tmnmn,mndy,pcpt,npcpdy,npcptr,ndly, 3 nsyn,nwpcp,pcpsum 1003 format(i4,i6,1x,a4,f7.2,f8.2,f7.0,f7.1,i3,f7.1,i3,f7.1, 2 i3,f7.1,i3,f7.1,1x,i2,f9.1,3i3,2i4,f9.1) go to 5 50 continue print 51,nr 51 format(' error, record',i10) go to 5 90 continue nce=nce+1 if(nce.lt.2) go to 5 print 1002,nr 1002 format(' eof - recs read =',i8) end