SUBROUTINE INPUT C C Base 8A version has P and PRESS arrays (which are in COMMON) defined here C #include "com2d.h" DATA AREA/0.0152,0.0451,0.0737,0.1000,0.1233, * 0.1428,0.1580,0.1683,0.1736,0.1736,0.1683, * 0.1580,0.1428,0.1233,0.1000,0.0737,0.0451, * 0.0152/ DATA IDAYMN/5,35,65,95,125,155,185,215,245,275, * 305,335/ C INITIALIZE PHYSICAL CONSTANTS. DP=0.2844 DZSTAR=7.E5*DP RD=2.87e6 WTA=28.97e0 RE=6.371e3 RCGS=RE*1.e5 RINV=1.0/RCGS PI=4.*ATAN(1.) RPI=1./PI DTR=PI/180. RTD=1./DTR BK=1.380622e-16 AMU=1.660531e-24 CP=3.5 H=6.6262e-27 CL=3.e10 OMEGAD=15. HBAR=6.79 ZBAR=30. R0=6371. YRL=360. DT0=86400. A=R0*1.e5 DY=.1745329252 DY2=2.*DY DPHI=DTH*DTR DY1=A*DPHI DP2=DP*DP DPHI2=DPHI*DPHI DO 104 IK=1,Z$ P(IK) = -DP*(IK-0.5) 104 PRESS(IK) = 1013.0E0*EXP(P(IK)) c T is the number of the run in days T=0. C FIX TIME OF DAY AT NOON. USED IN ZENITH ANGLE CALCULATION ONLY. CLOCK=12. C POSITIONS OF O2 AND O3 IN THE COLUMN DENSITY ARRAY. BOTH ARE C ASSUMED TO BE DEFINED. ISCO2=1 ISCO3=2 JLRT1=1 JURT1=IL$ DO 335 IJ=1,L$ PHI(IJ)=LAT(IJ)*DTR COSIN=COS(PHI(IJ)) TANGENT(IJ)=SIN(PHI(IJ))/COSIN RCSINV(IJ)=1./(RCGS*COSIN) c zero column and terms for lifetime calculations DO 335 IU=1,30 COLUMN(IU,IJ)=0.0 CLOSS(IU,IJ)=0.0 CLOSSST(IU,IJ)=0.0 CLOSSTR(IU,IJ)=0.0 335 CONTINUE do 230 iu=1,30 do 230 ij=1,l$ do 230 ik=1,z$ rloss(iu,ij,ik)=0.0 230 continue do 100 IJ=1,l$+1 phih(IJ)=(90.-10.*(IJ-1))*pi/180. 100 CONTINUE C set surface area of NAT aerosols equal to zero do 300 ij=1,l$ do 300 ik=1,z$ nataer(ij,ik)=0.0e0 300 continue c SET YSLFPRNT FOR LIFETIME CALCULATION YSLFPRNT=.FALSE. LIFECHAR(1)='N2O ' LIFECHAR(2)='CFCL3 ' LIFECHAR(3)='CF2CL2 ' LIFECHAR(4)='CCL4 ' LIFECHAR(5)='CH3CL ' LIFECHAR(6)='CH4 ' LIFECHAR(7)='CH3CCL3 ' LIFECHAR(8)='CO ' LIFECHAR(9)='H2 ' LIFECHAR(10)='CH3BR ' LIFECHAR(11)='CHCLF2 ' LIFECHAR(12)='C2CL3F3 ' LIFECHAR(13)='C2CL2F4 ' LIFECHAR(14)='C2CLF5 ' LIFECHAR(15)='CBRCLF2 ' LIFECHAR(16)='CBRF3 ' DO 9200 II=17,30 LIFECHAR(II)='XXX ' 9200 CONTINUE SAVE RETURN END