98 |
USE yomcst, ONLY : rd, rg, rkappa |
USE yomcst, ONLY : rd, rg, rkappa |
99 |
USE conf_phys_m, ONLY : iflag_pbl |
USE conf_phys_m, ONLY : iflag_pbl |
100 |
USE gath_cpl, ONLY : gath2cpl |
USE gath_cpl, ONLY : gath2cpl |
101 |
|
use hbtm_m, only: hbtm |
102 |
|
|
103 |
IMPLICIT NONE |
IMPLICIT NONE |
104 |
|
|
251 |
! -- LOOP |
! -- LOOP |
252 |
|
|
253 |
REAL yt10m(klon), yq10m(klon) |
REAL yt10m(klon), yq10m(klon) |
254 |
!IM cf. AM : pbl, hbtm2 (Comme les autres diagnostics on cumule ds |
!IM cf. AM : pbl, hbtm (Comme les autres diagnostics on cumule ds |
255 |
! physiq ce qui permet de sortir les grdeurs par sous surface) |
! physiq ce qui permet de sortir les grdeurs par sous surface) |
256 |
REAL pblh(klon, nbsrf) |
REAL pblh(klon, nbsrf) |
257 |
REAL plcl(klon, nbsrf) |
REAL plcl(klon, nbsrf) |
290 |
PARAMETER (t_coup=273.15) |
PARAMETER (t_coup=273.15) |
291 |
|
|
292 |
CHARACTER (len=20) :: modname = 'clmain' |
CHARACTER (len=20) :: modname = 'clmain' |
|
LOGICAL check |
|
|
PARAMETER (check=.FALSE.) |
|
293 |
|
|
294 |
!------------------------------------------------------------ |
!------------------------------------------------------------ |
295 |
|
|
296 |
! initialisation Anne |
! initialisation Anne |
297 |
ytherm = 0. |
ytherm = 0. |
298 |
|
|
|
IF (check) THEN |
|
|
PRINT *, modname, ' klon=', klon |
|
|
END IF |
|
|
|
|
299 |
IF (debugindex .AND. first_appel) THEN |
IF (debugindex .AND. first_appel) THEN |
300 |
first_appel = .FALSE. |
first_appel = .FALSE. |
301 |
|
|
444 |
END IF |
END IF |
445 |
END DO |
END DO |
446 |
|
|
|
IF (check) THEN |
|
|
PRINT *, 'CLMAIN, nsrf, knon =', nsrf, knon |
|
|
END IF |
|
|
|
|
447 |
! variables pour avoir une sortie IOIPSL des INDEX |
! variables pour avoir une sortie IOIPSL des INDEX |
448 |
IF (debugindex) THEN |
IF (debugindex) THEN |
449 |
tabindx = 0. |
tabindx = 0. |
805 |
|
|
806 |
END DO |
END DO |
807 |
|
|
|
!IM cf AM : pbl, HBTM |
|
808 |
DO i = 1, knon |
DO i = 1, knon |
809 |
y_cd_h(i) = ycoefh(i, 1) |
y_cd_h(i) = ycoefh(i, 1) |
810 |
y_cd_m(i) = ycoefm(i, 1) |
y_cd_m(i) = ycoefm(i, 1) |
811 |
END DO |
END DO |
812 |
! print*, 'appel hbtm2' |
CALL hbtm(knon, ypaprs, ypplay, yt2m, yt10m, yq2m, yq10m, yustar, & |
813 |
CALL hbtm(knon, ypaprs, ypplay, yt2m, yt10m, yq2m, yq10m, yustar, y_flux_t, & |
y_flux_t, y_flux_q, yu, yv, yt, yq, ypblh, ycapcl, yoliqcl, & |
814 |
y_flux_q, yu, yv, yt, yq, ypblh, ycapcl, yoliqcl, ycteicl, ypblt, ytherm, & |
ycteicl, ypblt, ytherm, ytrmb1, ytrmb2, ytrmb3, ylcl) |
|
ytrmb1, ytrmb2, ytrmb3, ylcl) |
|
|
! print*, 'fin hbtm2' |
|
815 |
|
|
816 |
DO j = 1, knon |
DO j = 1, knon |
817 |
i = ni(j) |
i = ni(j) |