| 1 | subroutine init_chimie_B |
|---|
| 2 | c |
|---|
| 3 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
|---|
| 4 | c Scheme B |
|---|
| 5 | c |
|---|
| 6 | c subject: |
|---|
| 7 | c -------- |
|---|
| 8 | c |
|---|
| 9 | c Initialization of the photochemical scheme |
|---|
| 10 | c |
|---|
| 11 | c VERSION: 30/04/2003 |
|---|
| 12 | c |
|---|
| 13 | c This is dependent of the chemistry version: |
|---|
| 14 | c for scheme B its only purpose is reading the |
|---|
| 15 | c photodissociation rates lookup table |
|---|
| 16 | c |
|---|
| 17 | c Author: Franck Lefevre for scheme B |
|---|
| 18 | c |
|---|
| 19 | c Arguments: |
|---|
| 20 | c ---------- |
|---|
| 21 | c |
|---|
| 22 | c The output variable is jphot and is put in common chimiedata. |
|---|
| 23 | c |
|---|
| 24 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
|---|
| 25 | c |
|---|
| 26 | implicit none |
|---|
| 27 | c |
|---|
| 28 | #include "dimensions.h" |
|---|
| 29 | #include "chimiedata.h" |
|---|
| 30 | #include "datafile.h" |
|---|
| 31 | c |
|---|
| 32 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
|---|
| 33 | c local: |
|---|
| 34 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
|---|
| 35 | c |
|---|
| 36 | integer fic, ij, iozo, isza, itemp, iz, ierr |
|---|
| 37 | c |
|---|
| 38 | real xsza |
|---|
| 39 | integer lnblnk |
|---|
| 40 | external lnblnk |
|---|
| 41 | c |
|---|
| 42 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
|---|
| 43 | c |
|---|
| 44 | fic = 81 |
|---|
| 45 | |
|---|
| 46 | open(fic, form = 'formatted', status = 'old', |
|---|
| 47 | $file =datafile(1:lnblnk(datafile))//'/jmars.20030707',iostat=ierr) |
|---|
| 48 | |
|---|
| 49 | IF (ierr.NE.0) THEN |
|---|
| 50 | write(*,*)'Error : cannot open file jmars.20030707 ' |
|---|
| 51 | write(*,*)'(in aeronomars/init_chimie_B.F)' |
|---|
| 52 | write(*,*)'It should be in :',datafile(1:lnblnk(datafile)),'/' |
|---|
| 53 | write(*,*)'1) You can change this directory address in ' |
|---|
| 54 | write(*,*)' file phymars/datafile.h' |
|---|
| 55 | write(*,*)'2) If necessary jmars.20030707 (and other datafiles)' |
|---|
| 56 | write(*,*)' can be obtained online on:' |
|---|
| 57 | write(*,*)' http://www.lmd.jussieu.fr/~forget/datagcm/datafile' |
|---|
| 58 | STOP |
|---|
| 59 | ENDIF |
|---|
| 60 | |
|---|
| 61 | c |
|---|
| 62 | print*, 'lecture table des J...' |
|---|
| 63 | c |
|---|
| 64 | do itemp = 1,ntemp |
|---|
| 65 | do iozo = 1,nozo |
|---|
| 66 | do isza = 1,nsza |
|---|
| 67 | do iz = 200,0,-1 |
|---|
| 68 | read(fic,*) colairtab(iz), xsza, table_ozo(iozo) |
|---|
| 69 | read(fic,*) (jphot(itemp,isza,iz,iozo,ij), ij=1,7) |
|---|
| 70 | read(fic,*) (jphot(itemp,isza,iz,iozo,ij), ij=8,11) |
|---|
| 71 | do ij = 1,nd |
|---|
| 72 | if (jphot(itemp,isza,iz,iozo,ij) .eq. 1.e-30) then |
|---|
| 73 | jphot(itemp,isza,iz,iozo,ij) = 0. |
|---|
| 74 | end if |
|---|
| 75 | end do |
|---|
| 76 | end do |
|---|
| 77 | end do |
|---|
| 78 | end do |
|---|
| 79 | end do |
|---|
| 80 | c |
|---|
| 81 | print*, 'lecture table des J...ok' |
|---|
| 82 | close(fic) |
|---|
| 83 | c |
|---|
| 84 | return |
|---|
| 85 | end |
|---|
