6 |
d_tr,d_trs ) |
d_tr,d_trs ) |
7 |
|
|
8 |
! This procedure is clean: no C preprocessor directive, no include line. |
! This procedure is clean: no C preprocessor directive, no include line. |
|
|
|
9 |
! From phylmd/cltracrn.F,v 1.2 2005/05/25 13:10:09 |
! From phylmd/cltracrn.F,v 1.2 2005/05/25 13:10:09 |
10 |
|
|
11 |
use indicesol, only: nbsrf |
use indicesol, only: nbsrf |
12 |
use dimphy, only: klon, klev |
use dimphy, only: klon, klev |
13 |
use YOMCST, only: RD, rg |
use SUPHEC_M, only: RD, rg |
14 |
|
|
15 |
IMPLICIT none |
IMPLICIT none |
16 |
!====================================================================== |
!====================================================================== |
45 |
! d_trs----output-R- le changement de "trs" |
! d_trs----output-R- le changement de "trs" |
46 |
!====================================================================== |
!====================================================================== |
47 |
!====================================================================== |
!====================================================================== |
48 |
REAL dtime |
REAL, intent(in):: dtime |
49 |
REAL u1lay(klon), v1lay(klon) |
REAL u1lay(klon), v1lay(klon) |
50 |
REAL coef(klon,klev) |
REAL coef(klon,klev) |
51 |
REAL, intent(in):: t(klon,klev) ! temperature (K) |
REAL, intent(in):: t(klon,klev) ! temperature (K) |
52 |
real ftsol(klon,nbsrf), pctsrf(klon,nbsrf) |
real ftsol(klon,nbsrf), pctsrf(klon,nbsrf) |
53 |
REAL tr(klon,klev), trs(klon) |
REAL tr(klon,klev), trs(klon) |
54 |
REAL, intent(in):: paprs(klon,klev+1) |
REAL, intent(in):: paprs(klon,klev+1) |
55 |
real pplay(klon,klev), delp(klon,klev) |
real, intent(in):: pplay(klon,klev) |
56 |
|
real delp(klon,klev) |
57 |
REAL masktr(klon) |
REAL masktr(klon) |
58 |
REAL fshtr(klon) |
REAL fshtr(klon) |
59 |
REAL hsoltr |
REAL hsoltr |
112 |
!AA Il doit y avoir coherence (dc la meme chose ici) |
!AA Il doit y avoir coherence (dc la meme chose ici) |
113 |
|
|
114 |
DO i = 1, klon |
DO i = 1, klon |
|
!AA zx_alpha1(i) = (paprs(i,1)-pplay(i,2))/(pplay(i,1)-pplay(i,2)) |
|
115 |
zx_alpha1(i) = 1.0 |
zx_alpha1(i) = 1.0 |
116 |
zx_alpha2(i) = 1.0 - zx_alpha1(i) |
zx_alpha2(i) = 1.0 - zx_alpha1(i) |
117 |
ENDDO |
ENDDO |