656 |
INTEGER ip_ebil ! PRINT level for energy conserv. diag. |
INTEGER ip_ebil ! PRINT level for energy conserv. diag. |
657 |
SAVE ip_ebil |
SAVE ip_ebil |
658 |
DATA ip_ebil/0/ |
DATA ip_ebil/0/ |
659 |
INTEGER if_ebil ! level for energy conserv. dignostics |
INTEGER, SAVE:: if_ebil ! level for energy conservation diagnostics |
|
SAVE if_ebil |
|
660 |
!+jld ec_conser |
!+jld ec_conser |
661 |
REAL d_t_ec(klon, llm) ! tendance du a la conersion Ec -> E thermique |
REAL d_t_ec(klon, llm) ! tendance du a la conersion Ec -> E thermique |
662 |
REAL ZRCPD |
REAL ZRCPD |
991 |
! Mettre en action les conditions aux limites (albedo, sst, etc.). |
! Mettre en action les conditions aux limites (albedo, sst, etc.). |
992 |
! Prescrire l'ozone et calculer l'albedo sur l'ocean. |
! Prescrire l'ozone et calculer l'albedo sur l'ocean. |
993 |
|
|
994 |
!!$ if (nq >= 5) then |
if (nq >= 5) then |
995 |
!!$ wo = qx(:, :, 5) * zmasse / dobson_u / 1e3 |
wo = qx(:, :, 5) * zmasse / dobson_u / 1e3 |
996 |
!!$ else IF (MOD(itap - 1, lmt_pas) == 0) THEN |
else IF (MOD(itap - 1, lmt_pas) == 0) THEN |
|
IF (MOD(itap - 1, lmt_pas) == 0) THEN |
|
997 |
CALL ozonecm(REAL(julien), rlat, paprs, wo) |
CALL ozonecm(REAL(julien), rlat, paprs, wo) |
998 |
ENDIF |
ENDIF |
999 |
|
|
1282 |
! (driver commun aux versions 3 et 4) |
! (driver commun aux versions 3 et 4) |
1283 |
|
|
1284 |
IF (ok_cvl) THEN ! new driver for convectL |
IF (ok_cvl) THEN ! new driver for convectL |
1285 |
CALL concvl (iflag_con, pdtphys, paprs, pplay, t_seri, q_seri, & |
CALL concvl(iflag_con, pdtphys, paprs, pplay, t_seri, q_seri, & |
1286 |
u_seri, v_seri, tr_seri, ntra, & |
u_seri, v_seri, tr_seri, ntra, & |
1287 |
ema_work1, ema_work2, & |
ema_work1, ema_work2, & |
1288 |
d_t_con, d_q_con, d_u_con, d_v_con, d_tr, & |
d_t_con, d_q_con, d_u_con, d_v_con, d_tr, & |