SUBROUTINE CROSEC(IJC,IKC) #include "com2d.h" C O2 DISSOCIATION CROSS SECTIONS DO 400 I=5,18 XSECT(1,I)=XSCHRUN(I,IJC,IKC) C NEED TO ASSIGN CROSS SECTIONS ONLY IN SCHUMANN-RUNGE BANDS. C THE CROSS SECTIONS OUTSIDE THE SCHUMANN-RUNGE BANDS HAVE ALREADY C BEEN ASSIGNED. 400 CONTINUE C NO DISSOCIATION CROSS SECTIONS DO 520 I=8,9 XSECT(16,I)=NOSIG(1,IJC,IKC) 520 CONTINUE DO 500 I=13,14 XSECT(16,I)=NOSIG(2,IJC,IKC) 500 CONTINUE C C PUT IN TEMPERATURE DEPENDENCE FOR CROSS SECTIONS FOR NO2 GE 2667A - 9/3/93 - JPL-94 DO 521 IJ=26,39 521 XSECT(7,IJ) = XNO2(IJ) + ( 1.E-22*NO2T(IJ)*(TS - 273.15) ) C C C J COEFFICIENTS FOR F113, F114, F115 FROM SIMON ET AL. (1988) C DONE 12/30/88 C C FIRST FOR F113 IS=33 DO 901 ILF113=10,22 LMN=1.E-1*WVL(ILF113) SIGLOG10=A0(1) + A1(1)*LMN + A2(1)*LMN*LMN + A3(1)*LMN*LMN*LMN * + A4(1)*LMN*LMN*LMN*LMN + * (TS-273.)*(B0(1) + B1(1)*LMN + B2(1)*LMN*LMN + * B3(1)*LMN*LMN*LMN + B4(1)*LMN*LMN*LMN*LMN) SIGF113=10.**SIGLOG10 XSECT(IS,ILF113)=SIGF113 901 CONTINUE C SECOND FOR F114 IS=34 DO 903 ILF114=3,21 LMN=1.E-1*WVL(ILF114) SIGLOG10=A0(2) + A1(2)*LMN + A2(2)*LMN*LMN + A3(2)*LMN*LMN*LMN * + A4(2)*LMN*LMN*LMN*LMN + * (TS-273.)*(B0(2) + B1(2)*LMN + B2(2)*LMN*LMN + * B3(2)*LMN*LMN*LMN + B4(2)*LMN*LMN*LMN*LMN) SIGF114=10.**SIGLOG10 XSECT(IS,ILF114)=SIGF114 903 CONTINUE C THIRD FOR F115 IS=35 DO 905 ILF115=3,20 LMN=1.E-1*WVL(ILF115) SIGLOG10=A0(3) + A1(3)*LMN + A2(3)*LMN*LMN + A3(3)*LMN*LMN*LMN * + A4(3)*LMN*LMN*LMN*LMN SIGF115=10.**SIGLOG10 XSECT(IS,ILF115)=SIGF115 905 CONTINUE c HOBr cross section installed 950216, David B. Considine is=62 xsecthobr=0. do 906 ilhobr=23,37 LMN=1.E-1*WVL(ilhobr) do 907 iilhobr=1,3 xsecthobr=xsecthobr+hobra(iilhobr)*exp(-1.*hobrb(iilhobr) c *(hobrlam(iilhobr)-lmn)**2) 907 continue xsect(is,ilhobr)=xsecthobr*1.e-18 906 continue C C FIND TEMP INDEX FOR TEMP DEPENDENT XSECTIONS C IJTS = IFIX(TS + 0.5) - 119 IF (TS .LT. 120.) IJTS = 1 IF (TS .GT. 320.) IJTS = 201 C C Load in all temp dependent X-sections C DO 701 I = 1,IL$ XSECT(2,I) = XSECTTD(1,I,IJTS) XSECT(3,I) = XSECTTD(2,I,IJTS) XO3(I) = XSECT(2,I) + XSECT(3,I) XSECT(6,I) = XSECTTD(17,I,IJTS) XSECT(8,I) = XSECTTD(3,I,IJTS) XSECT(9,I) = XSECTTD(4,I,IJTS) XSECT(10,I) = XSECTTD(5,I,IJTS) XSECT(11,I) = XSECTTD(6,I,IJTS) XSECT(14,I) = XSECTTD(7,I,IJTS) XSECT(15,I) = XSECTTD(8,I,IJTS) XSECT(20,I) = XSECTTD(9,I,IJTS) XSECT(21,I) = XSECTTD(10,I,IJTS) XSECT(22,I) = XSECTTD(11,I,IJTS) XSECT(26,I) = XSECTTD(12,I,IJTS) XSECT(32,I) = XSECTTD(13,I,IJTS) XSECT(37,I) = XSECTTD(14,I,IJTS) XSECT(38,I) = XSECTTD(18,I,IJTS) XSECT(45,I) = XSECTTD(15,I,IJTS) XSECT(56,I) = XSECTTD(19,I,IJTS) XSECT(58,I) = XSECTTD(16,I,IJTS) 701 CONTINUE C C PRESSURE-DEPENDENT QUANTUM YIELD FOR CH2O=H2+CO FOR C WAVELENGTHS GREATER THAN 329 NM. GROUND PRESSURE=1013MB. C REMOVE CONTRIBUTION INCLUDED IN INPUT (QUANTUM YIELD=1). C PRATIO=PRP/1013.0 DO 680 I=34,35 XCH2O = XSECTTD(6,I,IJTS) LMN=1.e-1*WVL(I) R=(1.0-EXP(112.80-0.3470*LMN))/(1.0+PRATIO . *(LMN-329.0)/(364.0-LMN)) 680 XSECT(11,I)=R*XCH2O C C C NOW WRITE OUT CROSS SECTIONS TO FORT.81 C C IF (IYR .GE. 19 .AND. DAY360 .GE. 350. C * .AND. IJC .EQ. 8 .AND. IKC .EQ. 17) THEN C C DO 140 I=1,PH$ C WRITE (81,210) C WRITE (81,20) TS, IYR, DAY360, IPHOT(I), PHOTCHAR(I) C C WRITE (81,13) (XSECT(I,JL),JL=1,IL$) C140 CONTINUE C C END IF 13 FORMAT(1P8E10.3) 20 FORMAT(F7.1, 3X, I3, 2X, F5.1, 5X, I2,3X,A37) 210 FORMAT(A10) SAVE RETURN END