SUBROUTINE SETSPV(IPRINT) INTEGER NSPECV,NSPC1V PARAMETER (NSPECV=24,NSPC1V=25) REAL BWNV(NSPC1V),WNOV(NSPECV),DWNV(NSPECV) REAL WLNV(NSPECV) INTEGER NTERM(NSPECV) REAL SOLARF(NSPECV),PEXPON(NSPECV) REAL ATERM(4,NSPECV),BTERM(4,NSPECV) COMMON /SPECTV/ BWNV,WNOV,DWNV,WLNV COMMON /VISGAS/SOLARF,NTERM,PEXPON,ATERM,BTERM DATA WLNV/ & 0.325, 0.375, 0.425, 0.475, 0.525, 0.575, & 0.640, 0.715, 0.789, 0.850, 0.891, 0.935, & 0.998, 1.073, 1.144, 1.213, 1.292, 1.381, & 1.484, 1.603, 1.742, 1.909, 2.111, 2.361/ data solarf/ & 525.430, 680.440, 1051.900, 1173.300, 1082.000, 1056.500, & 1495.100, 1123.100, 1053.200, 541.820, 397.630, 484.210, & 624.540, 528.390, 383.250, 372.260, 360.090, 339.230, & 315.600, 292.580, 267.600, 239.700, 208.040, 180.550/ data NTERM/1,1,4,4,4,3,4,4,3,4,3,3,3,4,4,3,3,4,3,4,4,3,4,4/ data pexpon/ & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, & 0.000, 0.000, 0.000, 0.149, 0.156, 0.186, & 0.302, 0.097, 1.150, 1.040, 1.030, 1.040, & 1.080, 1.070, 1.090, 1.050, 1.050, 0.959/ data ATERM/ & 1.000000, 0.000000, 0.000000, 0.000000, 1.000000, 0.000000, & 0.000000, 0.000000, 0.700000, 0.093900, 0.015000, 0.191100, & 0.300030, 0.135014, 0.460546, 0.104410, 0.164416, 0.375038, & 0.295630, 0.164917, 0.444755, 0.355864, 0.199380, 0.000000, & 0.198120, 0.331633, 0.335834, 0.134413, 0.327300, 0.335500, & 0.146800, 0.190400, 0.277055, 0.333367, 0.389578, 0.000000, & 0.126666, 0.416016, 0.364590, 0.092728, 0.286216, 0.492476, & 0.221308, 0.000000, 0.445402, 0.453717, 0.100882, 0.000000, & 0.351017, 0.371854, 0.277129, 0.000000, 0.434400, 0.477200, & 0.077600, 0.010800, 0.292343, 0.374023, 0.242032, 0.091602, & 0.501604, 0.385625, 0.112771, 0.000000, 0.597075, 0.308155, & 0.094771, 0.000000, 0.116916, 0.447207, 0.338013, 0.097864, & 0.541475, 0.367862, 0.090663, 0.000000, 0.468164, 0.212875, & 0.213276, 0.105685, 0.245440, 0.416617, 0.242734, 0.095209, & 0.330423, 0.503512, 0.166065, 0.000000, 0.307538, 0.361097, & 0.232456, 0.098909, 0.115609, 0.353355, 0.413319, 0.117718/ data BTERM/ &0.0000E+00,0.0000E+00,0.0000E+00,0.0000E+00,0.0000E+00,0.0000E+00, &0.0000E+00,0.0000E+00,0.0000E+00,7.8260E-04,2.3210E-03,5.5120E-03, &0.0000E+00,9.4390E-04,4.5750E-03,4.2690E-02,1.9500E-04,2.5340E-03, &1.5020E-02,8.9550E-02,2.5430E-03,1.5420E-02,3.1900E-02,0.0000E+00, &7.1740E-03,2.6010E-02,9.7160E-02,4.4030E-01,4.4550E-02,1.8750E-01, &7.9000E-01,2.8660E+00,5.6920E-02,1.9440E-01,7.9630E-01,0.0000E+00, &3.1190E-02,4.9600E-01,2.2970E+00,2.8660E+01,1.2740E+00,1.4250E+01, &1.4070E+02,0.0000E+00,9.4170E-02,6.1850E-01,6.9190E+00,0.0000E+00, &1.8690E+00,1.2560E+01,1.4730E+02,0.0000E+00,1.3600E-01,9.3830E-01, &4.2070E+00,1.4080E+02,6.1300E-02,1.3220E+00,2.2710E+01,8.5640E+02, &0.0000E+00,6.0630E-01,3.9880E+01,0.0000E+00,0.0000E+00,6.0850E-01, &4.5140E+01,0.0000E+00,2.7190E-01,1.9420E+00,3.1990E+01,1.5370E+03, &0.0000E+00,5.8440E-01,5.0650E+01,0.0000E+00,6.2920E-03,7.8350E-01, &3.5630E+01,1.2870E+03,1.7680E+00,3.0850E+01,4.6500E+02,1.3940E+04, &0.0000E+00,6.6970E-01,4.6030E+01,0.0000E+00,0.0000E+00,7.3460E-01, &3.9680E+01,3.1920E+03,9.5930E+00,3.7750E+02,5.1130E+03,2.7940E+05/ C ****** C SET UP SPECRAL INTERVALS THESE ARE BASED ON KATHY RAGES PROGRAM C CONVERT PER KM-AMAGATS TO PER GM CM-3 (SEE NOTES FOR CONSTANT) C&& FUV=2.5 c fuv=1. C&& SOLARF(1) = FUV*SOLARF(1) C&& SOLARF(2) = FUV*SOLARF(2) C DO 101 K=1,NSPECV DO 102 NT=1,4 BTERM(NT,K)=BTERM(NT,K)/(1.E5 * 16./22.4E3) 102 CONTINUE 101 CONTINUE C C SET UP MEAN WAVENUMBERS AND DELTAS C UNITS ON WAVENUMBER ARE CM-1 C UNITS ON WAVELN ARE MICRONS BWNV(1)=1.E4/.3 DO 100 K=2,NSPC1V BWLN=1.E4/BWNV(K-1) EWLN=2.*WLNV(K-1)-BWLN BWNV(K)=1.E4/EWLN 100 CONTINUE DO 160 K=1,NSPECV WNOV(K)=1.E4/WLNV(K) DWNV(K)=BWNV(K)-BWNV(K+1) 160 CONTINUE C IF THERE IS ONLY ONE SPECTRAL INTERVAL THEN TOTAL E IS USED AND IF (NSPECV .EQ. 1) DWNV(1) =1.0 C PRINT OUT SPECTRAL INTERVALS IF (IPRINT .GT. 1) THEN WRITE (6,190) DO 200 K=1,NSPECV WRITE (6,210)K,WLNV(K),WNOV(K),BWNV(K) & ,BWNV(K)+DWNV(K),DWNV(K) 200 CONTINUE WRITE(6,320) 320 FORMAT (///' J WAVELN SOLAR FLUX N NP [', & 16X,'A TERMS',14X,'] [',16X,'B TERMS',14X,']'/) DO 300 J=1,24 K=1 SUM=ATERM(K,J)+ATERM(K+1,J)+ATERM(K+2,J)+ATERM(K+3,J) WRITE(6,20)J,WLNV(J),SOLARF(J),NTERM(J),PEXPON(J) &,(ATERM(K,J),K=1,4),(BTERM(K,J),K=1,4) 300 CONTINUE c20 FORMAT(1X,I2,F7.4,1PE14.5,0PI2,F6.3,3X,1P8E10.3) 20 FORMAT(1X,I2,F7.4,1PE14.5,I2,F6.3,3X,1P8E10.3) END IF 210 FORMAT(1X,I3,F10.3,F10.2,F10.2,'-',F8.2,F10.3) 190 FORMAT(///' SPECTRAL INTERVALS'// & ' SNUM MICRONS WAVENU INTERVAL',11X,'DELTA-WN') C ****** END SPECTRAL INTERVAL SET UP ************* RETURN END