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 |
35 |
use drag_noro_m, only: drag_noro |
use drag_noro_m, only: drag_noro |
36 |
USE fcttre, ONLY: foeew, qsatl, qsats, thermcep |
USE fcttre, ONLY: foeew, qsatl, qsats, thermcep |
37 |
USE hgardfou_m, ONLY: hgardfou |
USE hgardfou_m, ONLY: hgardfou |
38 |
USE histcom, ONLY: histsync |
USE histsync_m, ONLY: histsync |
39 |
USE histwrite_m, ONLY: histwrite |
USE histwrite_m, ONLY: histwrite |
40 |
USE indicesol, ONLY: clnsurf, epsfra, is_lic, is_oce, is_sic, is_ter, & |
USE indicesol, ONLY: clnsurf, epsfra, is_lic, is_oce, is_sic, is_ter, & |
41 |
nbsrf |
nbsrf |
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' |