12 |
|
|
13 |
! This is the main procedure for the "physics" part of the program. |
! This is the main procedure for the "physics" part of the program. |
14 |
|
|
15 |
|
use aaam_bud_m, only: aaam_bud |
16 |
USE abort_gcm_m, ONLY: abort_gcm |
USE abort_gcm_m, ONLY: abort_gcm |
17 |
use ajsec_m, only: ajsec |
use ajsec_m, only: ajsec |
18 |
USE calendar, ONLY: ymds2ju |
USE calendar, ONLY: ymds2ju |
905 |
! Mettre en action les conditions aux limites (albedo, sst, etc.). |
! Mettre en action les conditions aux limites (albedo, sst, etc.). |
906 |
|
|
907 |
! Prescrire l'ozone et calculer l'albedo sur l'ocean. |
! Prescrire l'ozone et calculer l'albedo sur l'ocean. |
908 |
if (nqmx >= 5) then |
wo = ozonecm(REAL(julien), paprs) |
|
wo = qx(:, :, 5) * zmasse / dobson_u / 1e3 |
|
|
else IF (MOD(itap - 1, lmt_pas) == 0) THEN |
|
|
wo = ozonecm(REAL(julien), paprs) |
|
|
ENDIF |
|
909 |
|
|
910 |
! Évaporation de l'eau liquide nuageuse : |
! Évaporation de l'eau liquide nuageuse : |
911 |
DO k = 1, llm |
DO k = 1, llm |
1655 |
ENDDO |
ENDDO |
1656 |
ENDDO |
ENDDO |
1657 |
|
|
1658 |
!IM calcul composantes axiales du moment angulaire et couple des montagnes |
CALL aaam_bud(ra, rg, romega, rlat, rlon, pphis, zustrdr, zustrli, & |
1659 |
|
zustrph, zvstrdr, zvstrli, zvstrph, paprs, u, v, aam, torsfc) |
|
CALL aaam_bud(27, klon, llm, time, ra, rg, romega, rlat, rlon, pphis, & |
|
|
zustrdr, zustrli, zustrph, zvstrdr, zvstrli, zvstrph, paprs, u, v, & |
|
|
aam, torsfc) |
|
1660 |
|
|
1661 |
IF (if_ebil >= 2) THEN |
IF (if_ebil >= 2) THEN |
1662 |
ztit = 'after orography' |
ztit = 'after orography' |