--- trunk/libf/phylmd/physiq.f90 2011/12/06 15:07:04 54 +++ trunk/libf/phylmd/physiq.f90 2012/01/30 12:54:02 57 @@ -12,6 +12,7 @@ ! This is the main procedure for the "physics" part of the program. + use aaam_bud_m, only: aaam_bud USE abort_gcm_m, ONLY: abort_gcm use ajsec_m, only: ajsec USE calendar, ONLY: ymds2ju @@ -904,11 +905,7 @@ ! Mettre en action les conditions aux limites (albedo, sst, etc.). ! Prescrire l'ozone et calculer l'albedo sur l'ocean. - if (nqmx >= 5) then - wo = qx(:, :, 5) * zmasse / dobson_u / 1e3 - else IF (MOD(itap - 1, lmt_pas) == 0) THEN - wo = ozonecm(REAL(julien), paprs) - ENDIF + wo = ozonecm(REAL(julien), paprs) ! Évaporation de l'eau liquide nuageuse : DO k = 1, llm @@ -1658,11 +1655,8 @@ ENDDO ENDDO - !IM calcul composantes axiales du moment angulaire et couple des montagnes - - CALL aaam_bud(27, klon, llm, time, ra, rg, romega, rlat, rlon, pphis, & - zustrdr, zustrli, zustrph, zvstrdr, zvstrli, zvstrph, paprs, u, v, & - aam, torsfc) + CALL aaam_bud(ra, rg, romega, rlat, rlon, pphis, zustrdr, zustrli, & + zustrph, zvstrdr, zvstrli, zvstrph, paprs, u, v, aam, torsfc) IF (if_ebil >= 2) THEN ztit = 'after orography'