Directory of image this file is from
This file as a plain text file

C		MKARRY is used to scan through a given data file until the
C		the given time is found for the POT, then it fills the FLUX
C		arrray with the A/D values for the 37 calorimeters and the
C		the values from the POT and the PYRHELIOMETER channels. Then
C		these values are converted to engineering units with the con-
C		version  ENG=ADVAL*COEF(N) where COEF(N)  is a conversion 
C		factor for channel N.



	SUBROUTINE MKARRY(SCTIME,INVALD,INPUT,FMAX,BMAX)

	INTEGER SCREEN,OUTPUT,INPUT

	REAL FLUX(39,73),DATE(13),POSIT(13),STIME(13),HOUR
	REAL MIN,SEC,DAY,SCTIME,BMAX,FMAX,COEF(37),LINE(80)
	REAL PHI

	LOGICAL INVALD(128)

	COMMON /FLUCKS/FLUX

	DATA SCREEN/4/

C		The values of the coeficients for the calorimeters are 
C		entered here...

	DATA COEF/.0863647,.0801392,.0899048,.1039430,.0831299,.0852661
     $            .0831299,.0791626,.0791626,.0820930,.0863640,.0801392
     $            .0875244,.0791626,.0831299,.0923462,.0886841,.0831299
     $            .0886841,.0852661,.0875244,.0899048,.0831299,.0811157
     $            .0820923,.0831299,.0875244,.0791626,.0831299,.0899048
     $            .0841675,.0852661,.0831299,.0886841,.0806841,.0820923/

C		Now go get the flux array filled from the data file.

	CALL GETFLX(FLUX,INVALD,SCTIME,INPUT)

C		Begin processing the array
C		First convert to watts/centimeter-squred and find the maximum flux
C		over the entire plane FMAX, and the maximum flux around the border
C		BMAX.
	FMAX=0.0
	BMAX=0.0
	DO 20 M=1,73
	   FLUX(39,M)=FLUX(39,M)*0.6531
	   IF (FLUX(39,M).EQ.0.0) FLUX(39,M)=FLUX(39,M-1)
	   DO 30 N=1,37
	      IF (FLUX(N,M).GT.2048.0) FLUX(N,M)=0.0
	      PHI=FLUX(N,M)*COEF(N)
	      IF (PHI.GE.FMAX) FMAX=PHI
	      IF (((M.EQ.1).OR.(M.EQ.73).OR.(N.EQ.1).OR.(N.EQ.37))
     1  	   .AND.(PHI.GE.BMAX))  BMAX=PHI
	      FLUX(N,M)=PHI
30	   CONTINUE
	   FLUX(4,M)=(FLUX(3,M)+FLUX(5,M))/2
	   CALL CURPOS(14,M)
	   CALL CHROUT(32)
20	CONTINUE
	RETURN
	END

	   DO 30 N=1,37
	      IF (FLUX(N,M).GT.2048.0) FLUX(N,M)=0.0
	      PHI=FLUX(N,M)*COEF(N)
	      IF (PHI.GE.FMAX) FMAX=PHI
	      IF (((M.EQ.1).OR.(M.EQ.73).OR.(N.EQ.1).OR.(N.EQ.37))
     1  	   .AND.(PHI.GE.BMAX))  BMAX=PHI
	      FLUX(N,M)=PHI
30	   CONTINUE
	   FLUX(4,M)=(FLUX(3,M)+FLUX(5,M))/2340   255   212   162   118    87    70    54    38    26    16    16    12    11     4     1
    32   628   594   523   372   280   231   175   127    94    74    56    38    26    16    16    12     9     4     1
    33   628   618   550   397   300   249   187   135   100    78    60    38    26    18    18    14    15     4     3
    34   627   642   576   422   318   265   199   144   107    82    60    40    26    18    18    18     7     4     0
    35   628   662   597   442   335   277   208   150   114    86    62    40    26    16    18    15    17     4     1
    36   627   675   613   456   347   287   211   157   119    91    62    40    26    16    18    17    15     6     1
    37   627   679   619   462   353   289   216   161   127    96    66    40    26    16    16    17    15     6     3
    38   628   683   624   467   358   291   216   165   132   100    68    40    26    14    16    17    15     6     1
    39   629   675   614   464   358   289   216   165   132   102    70    42    23    14    16    17    15     6     1
    40   629   655   595   449   348   276   209   159   127   102    70    40    23    12    14    14    13     4     1
    41   629   628   569   431   334   265   199   152   121    97    68    40    21    12    12    14    10     4     1
    42   631   596   542   411   318   251   188   142   112    90    62    38    21    11    12    12     7     4     1
    43   632   566   514   389   299   235   174   130   102    81    56    34    19