source: LMDZ.3.3/branches/rel-LF/libf/phylmd/cltrac.F @ 404

Last change on this file since 404 was 2, checked in by lmdz, 25 years ago

Initial revision

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