SUBROUTINE CONTROL

#include "com2d.h"

c ITRANS:  no of transported species
c INDAYS:  no of days for this calculation
C IPF:  photolysis calculation frequency
C IRK:  parameter to read in tropospheric diffusion
C IOUTP:  interval in days of printout to unit 33
	READ(7,1111)ITRANS,INDAYS,IPF,IRK,IOUTP

	type *,itrans,indays,ipf,irk
C READ IN LSTSTATE (LOGICAL DEFINING IF RUN IS STEADY-STATE OR TIME-
C DEPENDENT.  READ IN IYEARSS 
C (YEAR DESIRED FOR STEADY-STATE SCENARIO).
	READ(7,1900)LSTSTATE,IYEARSS
1900	FORMAT(1X,L1,1X,I4)
	type *,lsstate,iyearss

	READ(7,1111)(INTTRA(I1),I1=1,ITRANS)
C   INTTRA IS THE ARRAY OF SPECIES WHICH WILL BE 
c	TRANSPORTED IN THIS RUN 
	type *,inttra
1111       FORMAT(16I5)

	READ(7,1113)DP,DTH,EKZZI,EKYYI
1113       FORMAT(8e10.2)
C   DP=0.2844 IS THE PRESSURE INTERVAL,  LOG(P(I+1))-LOG(P(I))
C   DPHI=10. IS THE LATITUDE INTERVAL,  IT GOES FROM -85. TO +85.
C  EKZZI:  lowest vertical eddy diffusion
C  EKYYI:  lowest horizontal eddy diffusion
	type *,dp,dth,ekzzi,ekyyi

	READ(7,1113)(LAT(I),I=1,L$)
C   LAT IS THE ARRAY OF LATITUDE VALUES IN THE MODEL.
	type *,'LAT',lat
C  READ CONTROL LOGICALS
	READ(7,1114)LDAV,YSMSCAT,YSCORSR,YSLFPRNT,
     c		YSO3FIX,YSHSCT
1114	FORMAT(30L1)
C     LDAV=TRUE   THEN DIURNAL AVERAGE APPROXIMATION IS ON.
C     LDAV= FALSE   THEN DAYTIME (NOON)  CALCULATION IS DONE.
C    YSMSCAT= TRUE   THEN MULTIPLE SCATTERING IS USED TO CALCULATE
C     THE RADIATION FIELD.  BOTH RAYLEIGH SCATTERING AND A GROUND
C     ALBEDO OF 0.3 ARE USED.
C    YSCORSR=TRUE   THEN CORRECT SCHUMANN RUNGE BAND CROSS SECTIONS
C     IN THE MANNER SUGGESTED BY J. FREDERICK ON 5/26/83
C    YSLFPRNT Inital value should be false; set true when IYR changes
C    YSO3FIX= TRUE  THEN O3 IS FIXED
C    YSHSCT= TRUE   THEN stratospheric planes fly!!
	type *,ldav,ysmscat,yscorsr,yslfprnt,yso3fix,yshsct
c  set up tropospheric diffusion
	READ(7,1111)(IZZ(I1),I1=1,IRK)
	READ(7,1113)(FZZ(I1),I1=1,IRK)
C    IZZ IS THE ARRAY OF ALTITUDE INDICES THAT THE CORRECTION FACTORS
C       CONTAINED IN ARRAY FZZ ARE APPLIED.
	READ(7,1111)(IYY(I1),I1=1,IRK)
 	READ(7,1113)(FYY(I1),I1=1,IRK)
C    IYY IS THE ARRAY OF ALTITUDE INDICES THAT THE CORRECTION FACTORS
C       CONTAINED IN ARRAY FYY ARE APPLIED.
	type *,izz,fzz,kyy,fyy
	save

	return

	end