| 1 | module conc |
|---|
| 2 | use dimphy |
|---|
| 3 | IMPLICIT NONE |
|---|
| 4 | |
|---|
| 5 | |
|---|
| 6 | REAL, ALLOCATABLE, SAVE :: mmean(:,:) |
|---|
| 7 | !$OMP THREADPRIVATE(mmean) |
|---|
| 8 | REAL, ALLOCATABLE, SAVE :: rho(:,:) |
|---|
| 9 | !$OMP THREADPRIVATE(rho) |
|---|
| 10 | REAL, ALLOCATABLE, SAVE :: Akknew(:,:) |
|---|
| 11 | !$OMP THREADPRIVATE(Akknew) |
|---|
| 12 | REAL, ALLOCATABLE, SAVE :: rnew(:,:) |
|---|
| 13 | !$OMP THREADPRIVATE(rnew) |
|---|
| 14 | REAL, ALLOCATABLE, SAVE :: cpnew(:,:) |
|---|
| 15 | !$OMP THREADPRIVATE(cpnew) |
|---|
| 16 | REAL, ALLOCATABLE, SAVE :: jfotsout(:,:,:) |
|---|
| 17 | !$OMP THREADPRIVATE(jfotsout) |
|---|
| 18 | REAL, ALLOCATABLE, SAVE :: coefit4(:,:) |
|---|
| 19 | !$OMP THREADPRIVATE(coefit4) |
|---|
| 20 | REAL, ALLOCATABLE, SAVE :: coefit3(:,:) |
|---|
| 21 | !$OMP THREADPRIVATE(coefit3) |
|---|
| 22 | REAL, ALLOCATABLE, SAVE :: coefit2(:,:) |
|---|
| 23 | !$OMP THREADPRIVATE(coefit2) |
|---|
| 24 | REAL, ALLOCATABLE, SAVE :: coefit1(:,:) |
|---|
| 25 | !$OMP THREADPRIVATE(coefit1) |
|---|
| 26 | REAL, ALLOCATABLE, SAVE :: coefit0(:,:) |
|---|
| 27 | !$OMP THREADPRIVATE(coefit0) |
|---|
| 28 | REAL, ALLOCATABLE, SAVE :: fluxtop(:) |
|---|
| 29 | !$OMP THREADPRIVATE(fluxtop) |
|---|
| 30 | ! REAL, ALLOCATABLE, SAVE :: freccen(:) |
|---|
| 31 | !$OMP THREADPRIVATE(freccen) |
|---|
| 32 | REAL, ALLOCATABLE, SAVE :: t0(:) |
|---|
| 33 | !$OMP THREADPRIVATE(t0) |
|---|
| 34 | ! REAL, ALLOCATABLE, SAVE :: crscabsi2(:,:) |
|---|
| 35 | !$OMP THREADPRIVATE(crscabsi2) |
|---|
| 36 | REAL, ALLOCATABLE, SAVE :: jabsifotsintpar(:,:,:) |
|---|
| 37 | !$OMP THREADPRIVATE(crscabsi2) |
|---|
| 38 | REAL, ALLOCATABLE, SAVE :: c1_16(:,:) |
|---|
| 39 | !$OMP THREADPRIVATE(c1_16) |
|---|
| 40 | REAL, ALLOCATABLE, SAVE :: c17_24(:) |
|---|
| 41 | !$OMP THREADPRIVATE(c17_24) |
|---|
| 42 | REAL, ALLOCATABLE, SAVE :: c25_29(:) |
|---|
| 43 | !$OMP THREADPRIVATE(c25_29) |
|---|
| 44 | REAL, ALLOCATABLE, SAVE :: c30_31(:) |
|---|
| 45 | !$OMP THREADPRIVATE(c30_31) |
|---|
| 46 | REAL, ALLOCATABLE, SAVE :: c32(:) |
|---|
| 47 | !$OMP THREADPRIVATE(c32) |
|---|
| 48 | REAL, ALLOCATABLE, SAVE :: c33(:) |
|---|
| 49 | !$OMP THREADPRIVATE(c33) |
|---|
| 50 | REAL, ALLOCATABLE, SAVE :: c34(:) |
|---|
| 51 | !$OMP THREADPRIVATE(c34) |
|---|
| 52 | REAL, ALLOCATABLE, SAVE :: c35(:) |
|---|
| 53 | !$OMP THREADPRIVATE(c35) |
|---|
| 54 | REAL, ALLOCATABLE, SAVE :: c36(:) |
|---|
| 55 | !$OMP THREADPRIVATE(c36) |
|---|
| 56 | REAL, ALLOCATABLE, SAVE :: ct1(:) |
|---|
| 57 | !$OMP THREADPRIVATE(ct1) |
|---|
| 58 | REAL, ALLOCATABLE, SAVE :: ct2(:) |
|---|
| 59 | !$OMP THREADPRIVATE(ct2) |
|---|
| 60 | REAL, ALLOCATABLE, SAVE :: p1(:) |
|---|
| 61 | !$OMP THREADPRIVATE(p1) |
|---|
| 62 | REAL, ALLOCATABLE, SAVE :: p2(:) |
|---|
| 63 | !$OMP THREADPRIVATE(p2) |
|---|
| 64 | |
|---|
| 65 | |
|---|
| 66 | CONTAINS |
|---|
| 67 | !================================================================= |
|---|
| 68 | subroutine conc_init |
|---|
| 69 | USE dimphy |
|---|
| 70 | implicit none |
|---|
| 71 | #include"param.h" |
|---|
| 72 | |
|---|
| 73 | ALLOCATE(mmean(klon,klev)) |
|---|
| 74 | ALLOCATE(rho(klon,klev)) |
|---|
| 75 | ALLOCATE(Akknew(klon,klev)) |
|---|
| 76 | ALLOCATE(rnew(klon,klev)) |
|---|
| 77 | ALLOCATE(cpnew(klon,klev)) |
|---|
| 78 | ALLOCATE(jfotsout(ninter,nabs,klev)) |
|---|
| 79 | ALLOCATE(coefit0(ninter,nabs)) |
|---|
| 80 | ALLOCATE(coefit1(ninter,nabs)) |
|---|
| 81 | ALLOCATE(coefit2(ninter,nabs)) |
|---|
| 82 | ALLOCATE(coefit3(ninter,nabs)) |
|---|
| 83 | ALLOCATE(coefit4(ninter,nabs)) |
|---|
| 84 | ALLOCATE(fluxtop(ninter)) |
|---|
| 85 | ! ALLOCATE(freccen(ninter)) |
|---|
| 86 | ALLOCATE(t0(nz2)) |
|---|
| 87 | ! ALLOCATE(crscabsi2(nabs,16)) |
|---|
| 88 | ALLOCATE(c1_16(nz2,16)) |
|---|
| 89 | ALLOCATE(c17_24(nz2) ) |
|---|
| 90 | ALLOCATE(c25_29(nz2) ) |
|---|
| 91 | ALLOCATE(c30_31(nz2) ) |
|---|
| 92 | ALLOCATE(c32(nz2) ) |
|---|
| 93 | ALLOCATE(c33(nz2) ) |
|---|
| 94 | ALLOCATE(c34(nz2) ) |
|---|
| 95 | ALLOCATE(c35(nz2) ) |
|---|
| 96 | ALLOCATE(c36(nz2) ) |
|---|
| 97 | ALLOCATE(ct2(ninter) ) |
|---|
| 98 | ALLOCATE(ct1(ninter) ) |
|---|
| 99 | ALLOCATE(p1(ninter) ) |
|---|
| 100 | ALLOCATE(p2(ninter) ) |
|---|
| 101 | ALLOCATE(jabsifotsintpar(nz2,nabs,ninter)) |
|---|
| 102 | end subroutine conc_init |
|---|
| 103 | |
|---|
| 104 | |
|---|
| 105 | end module conc |
|---|
| 106 | |
|---|
