| 1 | ! |
|---|
| 2 | ! $Header$ |
|---|
| 3 | ! |
|---|
| 4 | SUBROUTINE cltrac(dtime,coef,t,tr,flux,paprs,pplay,delp, |
|---|
| 5 | s d_tr) |
|---|
| 6 | IMPLICIT none |
|---|
| 7 | c====================================================================== |
|---|
| 8 | c Auteur(s): O. Boucher (LOA/LMD) date: 19961127 |
|---|
| 9 | c inspire de clvent |
|---|
| 10 | c Objet: diffusion verticale de traceurs avec flux fixe a la surface |
|---|
| 11 | c ou/et flux du type c-drag |
|---|
| 12 | c====================================================================== |
|---|
| 13 | c Arguments: |
|---|
| 14 | c dtime----input-R- intervalle du temps (en second) |
|---|
| 15 | c coef-----input-R- le coefficient d'echange (m**2/s) l>1 |
|---|
| 16 | c t--------input-R- temperature (K) |
|---|
| 17 | c tr-------input-R- la q. de traceurs |
|---|
| 18 | c flux-----input-R- le flux de traceurs a la surface |
|---|
| 19 | c paprs----input-R- pression a inter-couche (Pa) |
|---|
| 20 | c pplay----input-R- pression au milieu de couche (Pa) |
|---|
| 21 | c delp-----input-R- epaisseur de couche (Pa) |
|---|
| 22 | c cdrag----input-R- cdrag pour le flux de surface (non active) |
|---|
| 23 | c tr0------input-R- traceurs a la surface ou dans l'ocean (non active) |
|---|
| 24 | c d_tr-----output-R- le changement de tr |
|---|
| 25 | c flux_tr--output-R- flux de tr |
|---|
| 26 | c====================================================================== |
|---|
| 27 | #include "dimensions.h" |
|---|
| 28 | #include "dimphy.h" |
|---|
| 29 | REAL dtime |
|---|
| 30 | REAL coef(klon,klev) |
|---|
| 31 | REAL t(klon,klev), tr(klon,klev) |
|---|
| 32 | REAL paprs(klon,klev+1), pplay(klon,klev), delp(klon,klev) |
|---|
| 33 | REAL d_tr(klon,klev) |
|---|
| 34 | REAL flux(klon), cdrag(klon), tr0(klon) |
|---|
| 35 | c REAL flux_tr(klon,klev) |
|---|
| 36 | c====================================================================== |
|---|
| 37 | #include "YOMCST.h" |
|---|
| 38 | c====================================================================== |
|---|
| 39 | INTEGER i, k |
|---|
| 40 | REAL zx_ctr(klon,2:klev) |
|---|
| 41 | REAL zx_dtr(klon,2:klev) |
|---|
| 42 | REAL zx_buf(klon) |
|---|
| 43 | REAL zx_coef(klon,klev) |
|---|
| 44 | REAL local_tr(klon,klev) |
|---|
| 45 | REAL zx_alf1(klon), zx_alf2(klon), zx_flux(klon) |
|---|
| 46 | c====================================================================== |
|---|
| 47 | DO k = 1, klev |
|---|
| 48 | DO i = 1, klon |
|---|
| 49 | local_tr(i,k) = tr(i,k) |
|---|
| 50 | ENDDO |
|---|
| 51 | ENDDO |
|---|
| 52 | c |
|---|
| 53 | |
|---|
| 54 | c====================================================================== |
|---|
| 55 | DO i = 1, klon |
|---|
| 56 | zx_alf1(i) = (paprs(i,1)-pplay(i,2))/(pplay(i,1)-pplay(i,2)) |
|---|
| 57 | zx_alf2(i) = 1.0 - zx_alf1(i) |
|---|
| 58 | zx_flux(i) = -flux(i)*dtime*RG |
|---|
| 59 | c--pour le moment le flux est prescrit |
|---|
| 60 | c--cdrag et zx_coef(1) vaut 0 |
|---|
| 61 | cdrag(i) = 0.0 |
|---|
| 62 | tr0(i) = 0.0 |
|---|
| 63 | zx_coef(i,1) = cdrag(i)*dtime*RG |
|---|
| 64 | ENDDO |
|---|
| 65 | c====================================================================== |
|---|
| 66 | DO k = 2, klev |
|---|
| 67 | DO i = 1, klon |
|---|
| 68 | zx_coef(i,k) = coef(i,k)*RG/(pplay(i,k-1)-pplay(i,k)) |
|---|
| 69 | . *(paprs(i,k)*2/(t(i,k)+t(i,k-1))/RD)**2 |
|---|
| 70 | zx_coef(i,k) = zx_coef(i,k)*dtime*RG |
|---|
| 71 | ENDDO |
|---|
| 72 | ENDDO |
|---|
| 73 | c====================================================================== |
|---|
| 74 | DO i = 1, klon |
|---|
| 75 | zx_buf(i) = delp(i,1) + zx_coef(i,1)*zx_alf1(i) + zx_coef(i,2) |
|---|
| 76 | zx_ctr(i,2) = (local_tr(i,1)*delp(i,1)+ |
|---|
| 77 | . zx_coef(i,1)*tr0(i)-zx_flux(i))/zx_buf(i) |
|---|
| 78 | zx_dtr(i,2) = (zx_coef(i,2)-zx_alf2(i)*zx_coef(i,1)) / |
|---|
| 79 | . zx_buf(i) |
|---|
| 80 | ENDDO |
|---|
| 81 | c |
|---|
| 82 | DO k = 3, klev |
|---|
| 83 | DO i = 1, klon |
|---|
| 84 | zx_buf(i) = delp(i,k-1) + zx_coef(i,k) |
|---|
| 85 | . + zx_coef(i,k-1)*(1.-zx_dtr(i,k-1)) |
|---|
| 86 | zx_ctr(i,k) = (local_tr(i,k-1)*delp(i,k-1) |
|---|
| 87 | . +zx_coef(i,k-1)*zx_ctr(i,k-1) )/zx_buf(i) |
|---|
| 88 | zx_dtr(i,k) = zx_coef(i,k)/zx_buf(i) |
|---|
| 89 | ENDDO |
|---|
| 90 | ENDDO |
|---|
| 91 | DO i = 1, klon |
|---|
| 92 | local_tr(i,klev) = ( local_tr(i,klev)*delp(i,klev) |
|---|
| 93 | . +zx_coef(i,klev)*zx_ctr(i,klev) ) |
|---|
| 94 | . / ( delp(i,klev) + zx_coef(i,klev) |
|---|
| 95 | . -zx_coef(i,klev)*zx_dtr(i,klev) ) |
|---|
| 96 | ENDDO |
|---|
| 97 | DO k = klev-1, 1, -1 |
|---|
| 98 | DO i = 1, klon |
|---|
| 99 | local_tr(i,k) = zx_ctr(i,k+1) + zx_dtr(i,k+1)*local_tr(i,k+1) |
|---|
| 100 | ENDDO |
|---|
| 101 | ENDDO |
|---|
| 102 | c====================================================================== |
|---|
| 103 | c== flux_tr est le flux de traceur (positif vers bas) |
|---|
| 104 | c DO i = 1, klon |
|---|
| 105 | c flux_tr(i,1) = zx_coef(i,1)/(RG*dtime) |
|---|
| 106 | c ENDDO |
|---|
| 107 | c DO k = 2, klev |
|---|
| 108 | c DO i = 1, klon |
|---|
| 109 | c flux_tr(i,k) = zx_coef(i,k)/(RG*dtime) |
|---|
| 110 | c . * (local_tr(i,k)-local_tr(i,k-1)) |
|---|
| 111 | c ENDDO |
|---|
| 112 | c ENDDO |
|---|
| 113 | c====================================================================== |
|---|
| 114 | DO k = 1, klev |
|---|
| 115 | DO i = 1, klon |
|---|
| 116 | d_tr(i,k) = local_tr(i,k) - tr(i,k) |
|---|
| 117 | ENDDO |
|---|
| 118 | ENDDO |
|---|
| 119 | c |
|---|
| 120 | RETURN |
|---|
| 121 | END |
|---|
