/[lmdze]/trunk/libf/phylmd/physiq.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/physiq.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 54 by guez, Tue Dec 6 15:07:04 2011 UTC revision 61 by guez, Fri Apr 20 14:58:43 2012 UTC
# Line 12  contains Line 12  contains
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
# Line 34  contains Line 35  contains
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
# Line 904  contains Line 905  contains
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
# Line 1658  contains Line 1655  contains
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'
Legend:
Removed from v.54  
changed lines
  Added in v.61
  ViewVC Help
Powered by ViewVC 1.1.21