1 | SUBROUTINE cltrac_spl(dtime,coef,yu1,yv1,t,tr, & |
---|
2 | flux,paprs,pplay,d_tr) |
---|
3 | |
---|
4 | USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO & |
---|
5 | , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA & |
---|
6 | , R_ecc, R_peri, R_incl & |
---|
7 | , RA, RG, R1SA & |
---|
8 | , RSIGMA & |
---|
9 | , R, RMD, RMV, RD, RV, RCPD & |
---|
10 | , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12 & |
---|
11 | , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w & |
---|
12 | , RCW, RCS & |
---|
13 | , RLVTT, RLSTT, RLMLT, RTT, RATM & |
---|
14 | , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS & |
---|
15 | , RALPD, RBETD, RGAMD |
---|
16 | USE dimphy |
---|
17 | USE dimensions_mod, ONLY: iim, jjm, llm, ndm |
---|
18 | IMPLICIT none |
---|
19 | !====================================================================== |
---|
20 | ! Auteur(s): O. Boucher (LOA/LMD) date: 19961127 |
---|
21 | ! inspire de clvent |
---|
22 | ! Objet: diffusion verticale de traceurs avec flux fixe a la surface |
---|
23 | ! ou/et flux du type c-drag |
---|
24 | !====================================================================== |
---|
25 | ! Arguments: |
---|
26 | ! dtime----input-R- intervalle du temps (en second) |
---|
27 | ! coef-----input-R- le coefficient d'echange (m**2/s) l>1 |
---|
28 | ! yu1------input-R- le vent dans le 1iere couche |
---|
29 | ! yv1------input-R- le vent dans le 1iere couche |
---|
30 | ! t--------input-R- temperature (K) |
---|
31 | ! tr-------input-R- la q. de traceurs |
---|
32 | ! flux-----input-R- le flux de traceurs a la surface |
---|
33 | ! paprs----input-R- pression a inter-couche (Pa) |
---|
34 | ! pplay----input-R- pression au milieu de couche (Pa) |
---|
35 | ! delp-----input-R- epaisseur de couche (Pa) |
---|
36 | ! cdrag----input-R- cdrag pour le flux de surface (non active) |
---|
37 | ! tr0------input-R- traceurs a la surface ou dans l'ocean (non active) |
---|
38 | ! d_tr-----output-R- le changement de tr |
---|
39 | ! flux_tr--output-R- flux de tr |
---|
40 | !====================================================================== |
---|
41 | |
---|
42 | REAL :: dtime |
---|
43 | REAL :: coef(klon,klev) |
---|
44 | REAL :: yu1(klon), yv1(klon) |
---|
45 | REAL :: t(klon,klev), tr(klon,klev) |
---|
46 | REAL :: paprs(klon,klev+1), pplay(klon,klev), delp(klon,klev) |
---|
47 | REAL :: d_tr(klon,klev) |
---|
48 | REAL :: flux(klon), cdrag(klon), tr0(klon) |
---|
49 | ! REAL flux_tr(klon,klev) |
---|
50 | !====================================================================== |
---|
51 | !====================================================================== |
---|
52 | INTEGER :: i, k |
---|
53 | REAL :: zx_ctr(klon,2:klev) |
---|
54 | REAL :: zx_dtr(klon,2:klev) |
---|
55 | REAL :: zx_buf(klon) |
---|
56 | REAL :: zx_coef(klon,klev) |
---|
57 | REAL :: local_tr(klon,klev) |
---|
58 | REAL :: zx_alf1(klon), zx_alf2(klon), zx_flux(klon) |
---|
59 | !====================================================================== |
---|
60 | ! CHECKING VALUES |
---|
61 | ! print *,'CHECKING VALUES IN CLTRAC (INI)' |
---|
62 | ! print *,'d_tr = ',sum(d_tr),MINVAL(d_tr),MAXVAL(d_tr) |
---|
63 | ! print *,'flux = ',sum(flux),MINVAL(flux),MAXVAL(flux) |
---|
64 | ! print *,'tr = ',sum(tr),MINVAL(tr),MAXVAL(tr) |
---|
65 | !====================================================================== |
---|
66 | DO k = 1, klev |
---|
67 | DO i = 1, klon |
---|
68 | local_tr(i,k) = tr(i,k) |
---|
69 | delp(i,k) = paprs(i,k)-paprs(i,k+1) |
---|
70 | ENDDO |
---|
71 | ENDDO |
---|
72 | !====================================================================== |
---|
73 | DO i = 1, klon |
---|
74 | zx_alf1(i) = (paprs(i,1)-pplay(i,2))/(pplay(i,1)-pplay(i,2)) |
---|
75 | zx_alf2(i) = 1.0 - zx_alf1(i) |
---|
76 | zx_flux(i) = -flux(i)*dtime*RG |
---|
77 | !--pour le moment le flux est prescrit |
---|
78 | cdrag(i) = 0.0 |
---|
79 | ! cdrag(i) = coef(i,1) * (1.0+SQRT(yu1(i)**2+yv1(i)**2)) |
---|
80 | ! . * pplay(i,1)/(RD*t(i,1)) |
---|
81 | tr0(i) = 0.0 |
---|
82 | zx_coef(i,1) = cdrag(i)*dtime*RG |
---|
83 | ENDDO |
---|
84 | !====================================================================== |
---|
85 | DO k = 2, klev |
---|
86 | DO i = 1, klon |
---|
87 | zx_coef(i,k) = coef(i,k)*RG/(pplay(i,k-1)-pplay(i,k)) & |
---|
88 | *(paprs(i,k)*2/(t(i,k)+t(i,k-1))/RD)**2 |
---|
89 | zx_coef(i,k) = zx_coef(i,k)*dtime*RG |
---|
90 | ENDDO |
---|
91 | ENDDO |
---|
92 | !====================================================================== |
---|
93 | DO i = 1, klon |
---|
94 | zx_buf(i) = delp(i,1) + zx_coef(i,1)*zx_alf1(i) + zx_coef(i,2) |
---|
95 | zx_ctr(i,2) = (local_tr(i,1)*delp(i,1)+ & |
---|
96 | zx_coef(i,1)*tr0(i)-zx_flux(i))/zx_buf(i) |
---|
97 | zx_dtr(i,2) = (zx_coef(i,2)-zx_alf2(i)*zx_coef(i,1)) / & |
---|
98 | zx_buf(i) |
---|
99 | ENDDO |
---|
100 | ! |
---|
101 | DO k = 3, klev |
---|
102 | DO i = 1, klon |
---|
103 | zx_buf(i) = delp(i,k-1) + zx_coef(i,k) & |
---|
104 | + zx_coef(i,k-1)*(1.-zx_dtr(i,k-1)) |
---|
105 | zx_ctr(i,k) = (local_tr(i,k-1)*delp(i,k-1) & |
---|
106 | +zx_coef(i,k-1)*zx_ctr(i,k-1) )/zx_buf(i) |
---|
107 | zx_dtr(i,k) = zx_coef(i,k)/zx_buf(i) |
---|
108 | ENDDO |
---|
109 | ENDDO |
---|
110 | DO i = 1, klon |
---|
111 | local_tr(i,klev) = ( local_tr(i,klev)*delp(i,klev) & |
---|
112 | +zx_coef(i,klev)*zx_ctr(i,klev) ) & |
---|
113 | / ( delp(i,klev) + zx_coef(i,klev) & |
---|
114 | -zx_coef(i,klev)*zx_dtr(i,klev) ) |
---|
115 | ENDDO |
---|
116 | DO k = klev-1, 1, -1 |
---|
117 | DO i = 1, klon |
---|
118 | local_tr(i,k) = zx_ctr(i,k+1) + zx_dtr(i,k+1)*local_tr(i,k+1) |
---|
119 | ENDDO |
---|
120 | ENDDO |
---|
121 | !====================================================================== |
---|
122 | ! print *,'CHECKING VALUES IN CLTRAC (FIN)' |
---|
123 | ! print *,'local_tr = ',sum(local_tr),MINVAL(local_tr), |
---|
124 | ! . MAXVAL(local_tr) |
---|
125 | ! print *,'zx_ctr = ',sum(zx_ctr),MINVAL(zx_ctr),MAXVAL(zx_ctr) |
---|
126 | ! print *,'zx_dtr = ',sum(zx_dtr),MINVAL(zx_dtr),MAXVAL(zx_dtr) |
---|
127 | ! print *,'tr = ',sum(tr),MINVAL(tr),MAXVAL(tr) |
---|
128 | !====================================================================== |
---|
129 | !== flux_tr est le flux de traceur (positif vers bas) |
---|
130 | ! DO i = 1, klon |
---|
131 | ! flux_tr(i,1) = zx_coef(i,1)/(RG*dtime) |
---|
132 | ! ENDDO |
---|
133 | ! DO k = 2, klev |
---|
134 | ! DO i = 1, klon |
---|
135 | ! flux_tr(i,k) = zx_coef(i,k)/(RG*dtime) |
---|
136 | ! . * (local_tr(i,k)-local_tr(i,k-1)) |
---|
137 | ! ENDDO |
---|
138 | ! ENDDO |
---|
139 | !====================================================================== |
---|
140 | DO k = 1, klev |
---|
141 | DO i = 1, klon |
---|
142 | d_tr(i,k) = local_tr(i,k) - tr(i,k) |
---|
143 | ENDDO |
---|
144 | ENDDO |
---|
145 | ! print *,'CHECKING VALUES IN CLTRAC (END)' |
---|
146 | ! print *,'d_tr = ',sum(d_tr),MINVAL(d_tr),MAXVAL(d_tr) |
---|
147 | ! |
---|
148 | RETURN |
---|
149 | END SUBROUTINE cltrac_spl |
---|