C..... F _ P T A B L E C..... C..... Programmer: Daniel C. Wilkinson C..... C..... Date: June 11, 1982 C..... C..... Computer: Data General M600 C..... C..... Purpose: To create a table of prompt flare for C..... publication in SGD. The table is to have C..... column headings, and flare reports for different C..... days separated by a blank line. C..... C..... Abstract: The FLFP.SORT proceedure prepares the flare C..... data by sorting chronologically and zero filling C..... blank columns within the date, start time, max C..... time, end time, latitude, and CMD. If any of C..... these fields starts as all blanks then it is C..... is left alone. This allows these fields to be C..... easily output to the table with character C..... formats. The Central Meridian passage date is C..... calculated via a subroutine. Flare duration is C..... calculated from start and end times in minutes. C..... C..... H-ALPHA S O L A R F L A R E REPORT FORMAT C..... C..... Column Fmt Description C ---------------------------------------------------------------------- C..... 1- 2 I2 Data code; always 31. C..... 3- 5 I3 Station code; numeric; used within center. C..... 6- 7 I2 Year C..... 8- 9 I2 Month C..... 10-11 I2 Day C..... 12-13 A2 Asterisks mark record with unconfirmed change. C..... 14-17 I4 Start time. C..... 18 A1 Qualifier: D=after,E=before,U=uncertain. C..... 19-22 I4 End time. C..... 23 A1 Qualifier: D=after,E=before,U=uncertain. C..... 24-27 I4 Max time; time of maximum brightness. C..... 28 A1 Qualifier: D=after,E=before,U=uncertain. C..... 29 A1 N or S for north or south latitude. C..... 30-31 I2 Latitude C..... 32 A1 E or W for east or west central meridian distance. C..... 33-34 I2 Central meridian distance. C..... 35 A1 Importance based on flare area = S,1,2 or 3. C..... 36 A1 Brightness: F=faint, N=normal, B=bright. C..... 37 A1 Completeness = C,P,V or S; indicates the kind of C..... observation and the completeness of it. C..... 38-41 F4.0 Time flare area was measured. C..... 42-46 F5.0 Area; apparent area in millionths of solar disk. C..... 47-50 F4.1 Area; corrected area in square degrees; implied C..... decimal point between columns 49 and 50. C..... 51 A1 Blank C..... 52-56 F5.2 Line width; width in Angstroms of H-alpha line. C..... 57-59 A3 Intensity; brightness of H-alpha emission line C..... expressed as percentage above continuum. C..... 60 A1 X-ray class: C,M,X code the maximum power of 10 C..... the 1-8 Angstrom flux attains. C..... 61-63 F3.1 X-ray intensity: a number from 1.0 to 9.9 that C..... multiplies the X-ray class. C..... 64-67 I4 Calcium plage region in which flare occurred. C..... 68-71 A4 Station name abbreviation. C..... 72 I1 Seeing: atmospheric stability during observations. C..... 73-80 A8 Remarks: USAF flare reports contain letter-coded C..... remarks in columns 73-75, status codes in 76, C..... and NOAA/USAF region numbers in 77-80. C..... 81-85 I5 NOAA/USAF sunspot region number. C..... 86 A1 Blank; may be used to add a letter to a region. C..... 87-88 I2 Central meridian passage year. C..... 89-90 I2 Central meridian passage month. C..... 91-94 F4.1 Central meridain passage day. C..... 95 A1 Blank C..... 96-100 I5 Grouping number assigned by WDC-A Boulder. C----------------------------------------------------------------------- IMPLICIT INTEGER (A-Z) REAL - SHOUR,SMIN, - EHOUR,EMIN, - DTIME, - CAREA,XRAYMULT CHARACTER - YEAR*2,MONTH*2,DAY*2, - STIME*4,SQUAL*1, - ETIME*4,EQUAL*1, - MTIME*4,MQUAL*1, - NS*1,LAT*2, - EASTORWEST*1,CMDISTANCE*2, - IMPORTANCE*1,BRIGHTNESS*1,OBSTYPE*1, - TIMEOFMEASURE*4,AAREA*5,ACAREA*5, - XRAYCLASS*1,AXRAYMULT*4,PLAGE*4, - STATIONNAME*4,SEEING*1,REMARKS*8, - NOAAREGION*6,CMPYEAR*2,CMPMONTH*2,CMPDAY*4, - - DQUAL*1, - - HOLDDAY*2, - - HEADING1*120,HEADING2*120,HEADING3*120,HEADING4*120, - HEADING5*120,HEADING6*120,HEADING7*120,HEADING8*120, - HEADING9*120, - FOOT*120, - INFILE*30,OUTFILE*30, - SORTKEY*65, SORTDES*65, - MONTHNAMES*10,SGDMONTH*3, - PAGE*3 DIMENSION - MONTHNAMES(12),SGDMONTH(12) PARAMETER - (INUNIT=20, - TABLEUNIT=21, - HEADUNIT=22, - FOOTUNIT=23) DATA - PAGE/'.PA'/ - INCOUNT/0/, - PAGECOUNT/1/, - LINESPERPAGE/68/, - LINECOUNT/0/, - HOLDDAY/'XX'/, - (MONTHNAMES(I),I=1,12)/'JANUARY ', - 'FEBRUARY ', - 'MARCH ', - 'APRIL ', - 'MAY ', - 'JUNE ', - 'JULY ', - 'AUGUST ', - 'SEPTEMBER ', - 'OCTOBER ', - 'NOVEMBER ', - 'DECEMBER '/, - (SGDMONTH(I),I=1,12)/'Jan','Feb','Mar','Apr','May', - 'Jun','Jul','Aug','Sep','Oct', - 'Nov','Dec'/ C..... Write a message to the user and request that he/she enter C..... the name of the input FL file to be published. WRITE(*,1500) 1500 FORMAT(////, - /,' ********* F_PTABLE *********', - /,' ******* *******', - /,' ***** *****', - /,' ***** F.PTABLE sorts the incoming file on *****', - /,' ***** date and time and creates a Prompt *****', - /,' ***** Flare Table. *****', - /,' ***** *****', - /,' ******* *******', - /,' ********* *********') 99 PRINT *,' ' PRINT *,' Input filename : ' READ(*,1520)INFILE 1520 FORMAT(A30) CALL FEXIST(INFILE,EXIST) IF (EXIST.NE.0)THEN PRINT *,' ERROR ! File does not exist: ',INFILE GO TO 99 ENDIF PRINT *,' Output filename : ' READ(*,1520)OUTFILE PRINT *,' Starting SGD page # : ' READ(*,*) SGDPAGE SORTDES = ' DATE START UT MAX UT END UT' SORTKEY = 'SORT(6,6,C,A, 14,4,C,A, 24,4,C,A 19,4,C,A)' PRINT *,' S O R T I N G' PRINT *,' ---------------------------------------------------' PRINT *,' Input file : ',INFILE PRINT *,' Output file : ',INFILE PRINT *,' ' PRINT *,' KEYS: ',SORTDES PRINT *,' : ',SORTKEY CALL SORTPF(INFILE,INFILE,SORTKEY,SORTNUM,SORTERR) PRINT *,' ' PRINT *,' Nunber of Records Sorted : ',SORTNUM PRINT *,' Sort Completion Code : ',SORTERR PRINT *,' ---------------------------------------------------' PRINT *,' ' PRINT *,' Please Standby ... ! ' C..... Open the F file and the Heading file. OPEN( UNIT=INUNIT, FILE= INFILE, STATUS='OLD' ) OPEN( UNIT=TABLEUNIT, FILE= OUTFILE, STATUS='UNKNOWN' ) OPEN( UNIT=HEADUNIT, FILE= 'F_PHEAD', STATUS='OLD' ) OPEN( UNIT=FOOTUNIT, FILE= 'F_GFOOT', STATUS='OLD' ) C..... Read in the heading lines. READ(HEADUNIT,1525)HEADING1,HEADING2,HEADING3, - HEADING4,HEADING5,HEADING6,HEADING7,HEADING8, - HEADING9 1525 FORMAT(10(A110/)) C..... This is the top of the input loop. C..... Read in a record from the F file. C..... L OOOOOOOO OOOOOOOO PPPPPPPP VV C..... L O O O O P P VV C..... L O O O O P P VV C..... L O O O O PPPPPPPP VV C..... L O O O O P V VV V C..... L O O O O P V VV V C..... LLLLLLLL OOOOOOOO OOOOOOOO P VV 100 READ(INUNIT,1530,END=300)YEAR,MONTH,DAY, - STIME,SQUAL, - ETIME,EQUAL, - MTIME,MQUAL, - NS,LAT, - EASTORWEST,CMDISTANCE, - IMPORTANCE,BRIGHTNESS,OBSTYPE, - TIMEOFMEASURE,AAREA,CAREA, - XRAYCLASS,XRAYMULT,PLAGE, - STATIONNAME,SEEING,REMARKS, - NOAAREGION,CMPYEAR,CMPMONTH,CMPDAY 1530 FORMAT(2X,3X,A2,A2,A2,2X, - A4,A1, - A4,A1, - A4,A1, - A1,A2, - A1,A2, - A1,A1,A1, - A4,A5,F4.1,1X, - 5X,3X,A1,F3.1,A4, - A4,A1,A8, - A6,A2,A2,A4) INCOUNT=INCOUNT+1 IF(INCOUNT.NE.1)GO TO 700 C..... Write a heading to the first PAGE. READ(MONTH,1510)IMONTH 1510 FORMAT(I2) PAGETYPE=0 IF(MOD(SGDPAGE,2).NE.0)PAGETYPE=1 IF(PAGETYPE.EQ.0)WRITE(TABLEUNIT,1515)SGDPAGE, - SGDMONTH(IMONTH),YEAR,HEADING1(7:70), - HEADING2,HEADING3, - MONTHNAMES(IMONTH),YEAR, - HEADING5,HEADING6,HEADING7,HEADING8,HEADING9 1515 FORMAT(1X,I3,/, - A3,1X,A2,3X,A63,/, - 2(3X,A70/), - 3X,25X,A10,'19',A2,/, - 4(3X,A110/), - 3X,A110) IF(PAGETYPE.EQ.1)WRITE(TABLEUNIT,1516)SGDPAGE, - HEADING1,SGDMONTH(IMONTH),YEAR, - HEADING2,HEADING3, - MONTHNAMES(IMONTH),YEAR, - HEADING5,HEADING6, - HEADING7,HEADING8,HEADING9 1516 FORMAT(3X,62X,I3,/, - 3X,A61,A3,1X,A2,/, - 2(3X,A70/), - 3X,25X,A10,'19',A2,/, - 4(3X,A110/), - 3X,A110) 700 CONTINUE C..... Calculate the DTIME from the Begin and End times. C..... The calculations are made by first converting these times C..... into minutes. C..... When the event ends on the next day,whenever the C..... end time is less than the begin time,the duration is C..... calculated by adding 1440 minutes (one day)to the end time C..... before substracting the begin time. READ(STIME,1540)SHOUR,SMIN 1540 FORMAT(F2.0,F2.0) READ(ETIME,1540)EHOUR,EMIN SMIN=SHOUR*60.0+SMIN EMIN=EHOUR*60.0+EMIN IF(EMIN.LT.SMIN) - DTIME=(EMIN+1440.0)-SMIN IF(EMIN.GT.SMIN) - DTIME=EMIN-SMIN C..... Round the duration to the units place. The function C..... ANINT rounds to the nearest integral value. DTIME=ANINT(DTIME) C..... Select an appropriate duration qualifier. If the start C..... qualifier is E or the end qualifier is D the duration C..... qualifier should be D. If either the start or the end C..... qualifiers are U then the duration qualifier should be U. DQUAL=' ' IF(SQUAL.EQ.'E'.OR.EQUAL.EQ.'D')DQUAL='D' IF(SQUAL.EQ.'U'.OR.EQUAL.EQ.'U')DQUAL='U' IF(CAREA.EQ.0)THEN ACAREA = ' ' ELSE WRITE(ACAREA,1550)CAREA 1550 FORMAT(F5.1) ENDIF IF(XRAYMULT.EQ.0)THEN AXRAYMULT = ' ' ELSE WRITE(AXRAYMULT,1560)XRAYMULT 1560 FORMAT(F4.1) ENDIF C....................................................................... C....................................................................... C..... ..... C..... P A G E C O N T R O L S E C T I O N ..... C..... ..... C....................................................................... C....................................................................... IF(LINECOUNT.LT.LINESPERPAGE)GO TO 800 C..... Place a dashed line at the bottom of page. WRITE(TABLEUNIT,1576)HEADING9 1576 FORMAT(3X,A110) C..... Place the page number at the bottom of page. Also insert C..... a page eject character '<014>'. WRITE(TABLEUNIT,1577)PAGE 1577 FORMAT(A3) C..... Increment the PAGE counter. PAGECOUNT=PAGECOUNT+1 C..... Place a heading at the top of a new PAGE. SGDPAGE=SGDPAGE+1 PAGETYPE=0 IF(MOD(SGDPAGE,2).NE.0)PAGETYPE=1 IF(PAGETYPE.EQ.0)WRITE(TABLEUNIT,1515)SGDPAGE, - SGDMONTH(IMONTH),YEAR,HEADING1(7:110), - HEADING2,HEADING3, - MONTHNAMES(IMONTH),YEAR, - HEADING5,HEADING6,HEADING7,HEADING8,HEADING9 IF(PAGETYPE.EQ.1)WRITE(TABLEUNIT,1516)SGDPAGE, - HEADING1,SGDMONTH(IMONTH),YEAR, - HEADING2,HEADING3, - MONTHNAMES(IMONTH),YEAR, - HEADING5,HEADING6,HEADING7,HEADING8,HEADING9 C..... Reset the line count to zero. LINECOUNT=0 800 CONTINUE C..... If the DAY has changed and this isn't the first line C..... in a PAGE, skip a line. IF(HOLDDAY.NE.DAY.AND.LINECOUNT.NE.0)THEN WRITE(TABLEUNIT,1590) 1590 FORMAT(1H ) LINECOUNT=LINECOUNT+1 ENDIF IF(HOLDDAY.EQ.DAY.AND.LINECOUNT.NE.0)DAY=' ' IF(DAY.NE.' ')HOLDDAY=DAY C....................................................................... C....................................................................... C....................................................................... C..... Write out a reformatted record to the table. WRITE(TABLEUNIT,1600)STATIONNAME,DAY, - STIME,SQUAL, - MTIME,MQUAL, - ETIME,EQUAL, - NS,LAT,EASTORWEST,CMDISTANCE, - NOAAREGION,CMPMONTH,CMPDAY, - INT(DTIME),DQUAL, - IMPORTANCE,BRIGHTNESS,XRAYCLASS,AXRAYMULT, - SEEING,OBSTYPE, - TIMEOFMEASURE,AAREA,ACAREA, - REMARKS 1600 FORMAT(3X,A4,1X,A2, - 1X,A4,A1, - 1X,A4,A1, - 1X,A4,A1, - 1X,A1,A2,1X,A1,A2, - A6,2X,A2,1X,A4, - 1X,I4.0,A1, - 2X,A1,A1,1X,A1,A4, - 2X,A1,3X,A1,4X, - A4,4X,A5,5X,A5, - 5X,A8) C..... Increment the output line count. LINECOUNT=LINECOUNT+1 GO TO 100 C..... L OOOOOOOO OOOOOOOO PPPPPPPP VV C..... L O O O O P P V VV V C..... L O O O O P P V VV V C..... L O O O O PPPPPPPP VV C..... L O O O O P VV C..... L O O O O P VV C..... LLLLLLLL OOOOOOOO OOOOOOOO P VV C..... Error in the read from FLARE file. 200 WRITE(*,1610) 1610 FORMAT(' Error in the read.') C..... Input file has been depleted. Place a dashed line at the C..... bottom of the last page of output. Close all files. 300 WRITE(TABLEUNIT,1576)HEADING9 PRINT *,' ' PRINT *,' Pages output : ',PAGECOUNT PRINT *,' Ending SGD page # : ',SGDPAGE CALL FOOTNOTE( TABLEUNIT, FOOTUNIT, LINECOUNT ) STOP END