source: trunk/LMDZ.TITAN.old/libf/phytitan/setspi.F @ 3303

Last change on this file since 3303 was 3, checked in by slebonnois, 14 years ago

Creation de repertoires:

  • chantiers : pour communiquer sur nos projets de modifs
  • documentation : pour stocker les docs

Ajout de:

  • libf/phytitan : physique de Titan
  • libf/chimtitan: chimie de Titan
  • libf/phyvenus : physique de Venus
File size: 1.7 KB
Line 
1      SUBROUTINE SETSPI(IPRINT)
2      PARAMETER (NSPECI=46,NSPC1I=47)
3      COMMON /SPECTI/ BWNI(NSPC1I),WNOI(NSPECI),DWNI(NSPECI)
4     &           ,WLNI(NSPECI)
5C SET UP SPECRAL INTERVALS
6      BWNI(1)=8.
7      BWNI(2)=15.
8      BWNI(3)=25.
9      BWNI(4)=37.5
10      DO 100 K=5,27
11      BWNI(K)=BWNI(K-1)+25.
12 100  CONTINUE
13      BWNI(28)=645.
14      BWNI(29)=645.+35.6760
15      BWNI(30)=BWNI(29)+4.324/2.
16      BWNI(31)=BWNI(30)+4.324/2.
17      BWNI(32)=BWNI(31)+12.327/4.
18      BWNI(33)=BWNI(32)+12.327/4.
19      BWNI(34)=BWNI(33)+12.327/4.
20      BWNI(35)=BWNI(34)+12.327/4.
21      BWNI(36)=BWNI(35)+35.6290
22      BWNI(37)=BWNI(36)+7.044/4.
23      BWNI(38)=BWNI(37)+7.044/4.
24      BWNI(39)=BWNI(38)+7.044/4.
25      BWNI(40)=BWNI(39)+7.044/4.
26      BWNI(41)=BWNI(40)+16.32/3.
27      BWNI(42)=BWNI(41)+16.32/3.
28      BWNI(43)=BWNI(42)+16.32/3.
29      BWNI(44)=BWNI(43)+156.48
30      BWNI(45)=BWNI(44)+27.2/3.
31      BWNI(46)=BWNI(45)+27.2/3.
32      BWNI(47)=BWNI(46)+27.2/3.
33C
34C SET UP MEAN WAVENUMBERS AND DELTAS
35C UNITS ON WAVENUMBER ARE CM-1
36C UNITS ON WAVELN ARE MICRONS
37      DO 160 K=1,NSPECI
38      WNOI(K)=( BWNI(K+1)+BWNI(K) )*0.5
39      DWNI(K)=BWNI(K+1)-BWNI(K)
40      WLNI(K)=1.E+4/WNOI(K)
41  160 CONTINUE
42C IF THERE IS ONLY ONE INTERVAL THEN TOTAL ENERGY IS USED AND
43      IF (NSPECI .EQ. 1) DWNI(1) = 1.0
44C PRINT OUT SPECTRAL INTERVALS
45      IF (IPRINT .GT. 1) THEN
46          WRITE (6,190)
47          DO 200 K=1,NSPECI
48          WRITE (6,210)K,WLNI(K),WNOI(K),BWNI(K)
49     &    ,BWNI(K+1),DWNI(K)
50  200     CONTINUE
51      END IF
52  210 FORMAT(1X,I3,F10.3,F10.2,F10.2,'-',F8.2,F10.3)
53  190 FORMAT(///'  SPECTRAL INTERVALS'//
54     &       ' SNUM  MICRONS   WAVENU      INTERVAL',11X,'DELTA-WN')
55C ***** END SPECTRAL INTERVAL SET UP *************
56      RETURN
57      END
Note: See TracBrowser for help on using the repository browser.