/[lmdze]/trunk/phylmd/phytrac.f

Diff of /trunk/phylmd/phytrac.f

Parent Directory | Revision Log | View Patch Patch

-revision 35 by guez,
Tue Jun  8 15:37:21 2010 UTC
+revision 62 by guez,
Thu Jul 26 14:37:37 2012 UTC
 Line 8 
 module phytrac_m
  contains
    SUBROUTINE phytrac(rnpb, itap, lmt_pas, julien, gmtime, firstcal, lafin, &
-        nq_phys, pdtphys, u, t_seri, paprs, pplay, pmfu, pmfd, pen_u, &
+        nq_phys, pdtphys, u, t_seri, paprs, pplay, pmfu, pmfd, pen_u, pde_u, &
-        pde_u, pen_d, pde_d, coefh, fm_therm, entr_therm, yu1, yv1, ftsol, &
+        pen_d, pde_d, coefh, fm_therm, entr_therm, yu1, yv1, ftsol, pctsrf, &
-        pctsrf, frac_impa, frac_nucl, pphis, pphi, albsol, rh, cldfra, rneb, &
+        frac_impa, frac_nucl, pphis, albsol, rh, cldfra, rneb, diafra, cldliq, &
-        diafra, cldliq, pmflxr, pmflxs, prfl, psfl, da, &
+        pmflxr, pmflxs, prfl, psfl, da, phi, mp, upwd, dnwd, tr_seri, zmasse)
-        phi, mp, upwd, dnwd, tr_seri, zmasse)
-     ! From phylmd/phytrac.F, version 1.15 2006/02/21 08:08:30
+     ! From phylmd/phytrac.F, version 1.15 2006/02/21 08:08:30 (SVN revision 679)
      ! Authors: Frédéric Hourdin, Abderrahmane Idelkadi, Marie-Alice
      ! Foujols, Olivia
-Line 24 
 contains
+Line 23 
 contains
      ! les vrais traceurs (en nombre "nbtr", sans la vapeur d'eau ni l'eau
      ! liquide) dans "phytrac".
+     ! Modifications pour les traceurs :
+     ! - uniformisation des parametrisations dans phytrac
+     ! - stockage des moyennes des champs nécessaires en mode traceur off-line
      use dimens_m, only: llm
      use indicesol, only: nbsrf
      use dimphy, only: klon, nbtr
      use clesphys, only: ecrit_tra
      use clesphys2, only: iflag_con
      use abort_gcm_m, only: abort_gcm
-     use YOMCST, only: rg
+     use SUPHEC_M, only: rg
      use ctherm, only: iflag_thermals
      use regr_pr_comb_coefoz_m, only: regr_pr_comb_coefoz
      use phyetat0_m, only: rlat
-Line 38 
 contains
+Line 41 
 contains
      use ini_histrac_m, only: ini_histrac
      use radiornpb_m, only: radiornpb
      use minmaxqfi_m, only: minmaxqfi
-     use numer_rec, only: assert
+     use nr_util, only: assert
      use press_coefoz_m, only: press_coefoz
      logical, intent(in):: rnpb
-Line 56 
 contains
+Line 59 
 contains
      real, intent(inout):: tr_seri(:, :, :) ! (klon, llm, nbtr)
      ! (mass fractions of tracers, excluding water, at mid-layers)
-     real u(klon, llm)
+     real, intent(in):: u(klon, llm)
      real rh(klon, llm)     ! humidite relative
      real cldliq(klon, llm) ! eau liquide nuageuse
      real cldfra(klon, llm) ! fraction nuageuse (tous les nuages)
-Line 73 
 contains
+Line 76 
 contains
      real, intent(in):: pplay(klon, llm)
      ! (pression pour le mileu de chaque couche, en Pa)
-     real pphi(klon, llm) ! geopotentiel
+     real, intent(in):: pphis(klon)
-     real pphis(klon)
      logical, intent(in):: firstcal ! first call to "calfis"
      logical, intent(in):: lafin ! fin de la physique
-Line 82 
 contains
+Line 84 
 contains
      REAL prfl(klon, llm+1),   psfl(klon, llm+1)     !--lessivage large-scale
      !   convection:
-     REAL pmfu(klon, llm)  ! flux de masse dans le panache montant
+     REAL, intent(in):: pmfu(klon, llm) ! flux de masse dans le panache montant
      REAL pmfd(klon, llm)  ! flux de masse dans le panache descendant
      REAL pen_u(klon, llm) ! flux entraine dans le panache montant
-Line 93 
 contains
+Line 95 
 contains
      REAL pde_u(klon, llm) ! flux detraine dans le panache montant
      REAL pen_d(klon, llm) ! flux entraine dans le panache descendant
      REAL pde_d(klon, llm) ! flux detraine dans le panache descendant
-     ! KE
+     ! Kerry Emanuel
      real da(klon, llm), phi(klon, llm, llm), mp(klon, llm)
-     REAL upwd(klon, llm)      ! saturated updraft mass flux
+     REAL, intent(in):: upwd(klon, llm) ! saturated updraft mass flux
-     REAL dnwd(klon, llm)      ! saturated downdraft mass flux
+     REAL, intent(in):: dnwd(klon, llm) ! saturated downdraft mass flux
      !   Couche limite:
-Line 234 
 contains
+Line 236 
 contains
         end if
      ENDIF
-     ! Initialisation du traceur dans le sol (couche limite radonique)
      if (inirnpb) THEN
+        ! Initialisation du traceur dans le sol (couche limite radonique)
         radio(1)= .true.
         radio(2)= .true.
         clsol(1)= .true.
         clsol(2)= .true.
         aerosol(2) = .TRUE. ! le Pb est un aerosol
         call initrrnpb(ftsol, pctsrf, masktr, fshtr, hsoltr, tautr, vdeptr, &
              scavtr)
         inirnpb=.false.
      endif
-     ! Calcul de l'effet de la convection
      if (convection) then
+        ! Calcul de l'effet de la convection
         DO it=1, nq_phys
-           if (iflag_con.eq.2) then
+           if (iflag_con == 2) then
-              ! tiedke
+              ! Tiedke
               CALL nflxtr(pdtphys, pmfu, pmfd, pen_u, pde_u, pen_d, pde_d, &
                    paprs, tr_seri(1, 1, it), d_tr_cv(1, 1, it))
-           else if (iflag_con.eq.3) then
+           else if (iflag_con == 3) then
-              ! KE
+              ! Emanuel
-              call cvltr(pdtphys, da, phi, mp, paprs, &
+              call cvltr(pdtphys, da, phi, mp, paprs, tr_seri(1, 1, it), upwd, &
-                   tr_seri(1, 1, it), upwd, dnwd, d_tr_cv(1, 1, it))
+                   dnwd, d_tr_cv(1, 1, it))
            endif
            DO k = 1, llm
-Line 453 
 contains
+Line 452 
 contains
        ! From phylmd/write_histrac.h, version 1.9 2006/02/21 08:08:30
        use dimens_m, only: iim, jjm, llm
-       use histcom, only: histsync
+       use histsync_m, only: histsync
        use histwrite_m, only: histwrite
        use temps, only: itau_phy
        use iniadvtrac_m, only: tnom

 Legend:



Removed from v.35
 


changed lines


 
Added in v.62
 Legend:



Removed from v.35
 


changed lines


 
Added in v.62
-Removed from v.35
+Added in v.62

	ViewVC Help
Powered by ViewVC 1.1.21