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