source: trunk/LMDZ.TITAN.old/libf/phytitan/setpia.F @ 3094

Last change on this file since 3094 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 SETPIA(IPRINT,IFLAG)
2      PARAMETER (NSPECI=46,NSPC1I=47, NTEMPS=14)
3      REAL PIANN,PIACC,PIACN,PIAHN,T,W
4      COMMON /PIAC/ PIANN(NSPECI,NTEMPS),PIACC(NSPECI,NTEMPS)
5     & ,PIACN(NSPECI,NTEMPS),PIAHN(NSPECI,NTEMPS),TMIN,TMAX
6      COMMON /SPECTI/ BWNI(NSPC1I),WNOI(NSPECI),DWNI(NSPECI)
7     &           ,WLNI(NSPECI)
8      DIMENSION T(NTEMPS),W(NSPECI)
9      DATA TMAX,TMIN/190.,60./
10      DO 11 K=1,NSPECI
11      DO 11 NT=1,NTEMPS
12      PIANN(K,NT)=0.0
13      PIACC(K,NT)=0.0
14      PIACN(K,NT)=0.0
15      PIAHN(K,NT)=0.0
16 11   CONTINUE
17      NWAVES=NSPECI
18      DO 99 INU=1,NSPECI
19      W(INU)=WNOI(INU)
20      IF (WNOI(INU) .LT. 1000.) GOTO 99
21      NWAVES=INU-1
22      GO TO 991
23 99   CONTINUE
24 991  CONTINUE
25      DT=(TMAX-TMIN)/(NTEMPS-1)
26      DO 12 I=1,NTEMPS
27      T(I)=TMIN + (I-1)*DT
28 12   CONTINUE
29      CALL REGIS (W,NSPECI,NWAVES,T,NTEMPS,NTEMPS,
30     &  PIANN,PIACC,PIACN,PIAHN)
31      RETURN
32CC TO BE DELETED...
33 1234 NT=NTEMPS
34      NNU=NSPECI
35      WRITE(6,101)
36101   FORMAT(/55X,'N2-N2 ABSORPTION'//)
37      WRITE(6,100) (T(IT),IT=1,NT)
38100   FORMAT(2X,'NU\T',13F9.0/)
39      DO 3 INU=1,NNU
403     WRITE(6,110) WNOI(INU),(PIANN(INU,IT),IT=1,NT)
41      WRITE(6,103)
42103   FORMAT(/55X,'CH4-CH4 ABSORPTION'//)
43      WRITE(6,100) (T(IT),IT=1,NT)
44      DO 5 INU=1,NNU
455     WRITE(6,110) WNOI(INU),(PIACC(INU,IT),IT=1,NT)
46      WRITE(6,105)
47105   FORMAT(/50X,'N2-CH4 + CH4-N2 ABSORPTION'//)
48      DO 7 INU=1,NNU
497     WRITE(6,110) WNOI(INU),(PIACN(INU,IT),IT=1,NT)
50      WRITE(6,102)
51102   FORMAT(/55X,'N2-H2 + H2-N2 ABSORPTION'//)
52      WRITE(6,100) (T(IT),IT=1,NT)
53      DO 4 INU=1,NNU
544     WRITE(6,110) WNOI(INU),(PIAHN(INU,IT),IT=1,NT)
55110   FORMAT(1X,F5.0,2X,13(1PD9.2))
56      RETURN
57      END
Note: See TracBrowser for help on using the repository browser.