/[lmdze]/trunk/Sources/phylmd/phytrac.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/phytrac.f

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

revision 224 by guez, Tue Apr 25 15:31:48 2017 UTC revision 225 by guez, Mon Oct 16 12:35:41 2017 UTC
# Line 78  contains Line 78  contains
78      REAL pde_u(klon, llm) ! flux detraine dans le panache montant      REAL pde_u(klon, llm) ! flux detraine dans le panache montant
79      REAL pen_d(klon, llm) ! flux entraine dans le panache descendant      REAL pen_d(klon, llm) ! flux entraine dans le panache descendant
80      REAL coefh(klon, llm) ! coeff melange couche limite      REAL coefh(klon, llm) ! coeff melange couche limite
81        real fm_therm(klon, llm+1), entr_therm(klon, llm) ! thermiques
82      ! thermiques:      REAL, intent(in):: yu1(klon), yv1(klon) ! vent au premier niveau
     real fm_therm(klon, llm+1), entr_therm(klon, llm)  
   
     ! Couche limite:  
     REAL yu1(klon) ! vents au premier niveau  
     REAL yv1(klon) ! vents au premier niveau  
83    
84      ! Arguments n\'ecessaires pour les sources et puits de traceur :      ! Arguments n\'ecessaires pour les sources et puits de traceur :
85      real, intent(in):: ftsol(:, :) ! (klon, nbsrf) surface temperature (K)      real, intent(in):: ftsol(:, :) ! (klon, nbsrf) surface temperature (K)
# Line 140  contains Line 135  contains
135    
136      CHARACTER itn      CHARACTER itn
137    
138      ! nature du traceur      logical, save:: aerosol(nqmx - 2) ! Nature du traceur
   
     logical aerosol(nqmx - 2) ! Nature du traceur  
139      ! ! aerosol(it) = true => aerosol      ! ! aerosol(it) = true => aerosol
140      ! ! aerosol(it) = false => gaz      ! ! aerosol(it) = false => gaz
141      logical clsol(nqmx - 2) ! couche limite sol calcul\'ee  
142      logical radio(nqmx - 2) ! d\'ecroisssance radioactive      logical, save:: clsol(nqmx - 2) ! couche limite sol flux
143      save aerosol, clsol, radio      ! calcul\'ee, sinon prescrit
144        logical, save:: radio(nqmx - 2) ! d\'ecroisssance radioactive
145    
146      ! convection tiedtke      ! convection tiedtke
147      INTEGER i, k, it      INTEGER i, k, it
# Line 179  contains Line 173  contains
173      integer isplit, varid      integer isplit, varid
174    
175      ! Controls:      ! Controls:
     logical:: couchelimite = .true.  
176      logical:: convection = .true.      logical:: convection = .true.
     logical, save:: inirnpb  
177    
178      !--------------------------------------      !--------------------------------------
179    
# Line 189  contains Line 181  contains
181      call assert(shape(tr_seri) == (/klon, llm, nqmx - 2/), "phytrac tr_seri")      call assert(shape(tr_seri) == (/klon, llm, nqmx - 2/), "phytrac tr_seri")
182    
183      if (firstcal) then      if (firstcal) then
        inirnpb = .true.  
   
184         ! Initialisation de certaines variables pour le radon et le plomb         ! Initialisation de certaines variables pour le radon et le plomb
185         ! Initialisation du traceur dans le sol (couche limite radonique)         ! Initialisation du traceur dans le sol (couche limite radonique)
186         trs(:, 2:) = 0.         trs(:, 2:) = 0.
# Line 207  contains Line 197  contains
197    
198         ! Initialisation de la nature des traceurs         ! Initialisation de la nature des traceurs
199    
200         DO it = 1, nqmx - 2         aerosol = .FALSE. ! Tous les traceurs sont des gaz par defaut
201            aerosol(it) = .FALSE. ! Tous les traceurs sont des gaz par defaut         radio = .FALSE. ! par d\'efaut pas de passage par "radiornpb"
           radio(it) = .FALSE. ! par d\'efaut pas de passage par "radiornpb"  
           clsol(it) = .FALSE. ! Par defaut couche limite avec flux prescrit  
        ENDDO  
202    
203         if (nqmx >= 5) then         if (nqmx >= 5) then
204            call press_coefoz ! read input pressure levels for ozone coefficients            call press_coefoz ! read input pressure levels for ozone coefficients
205         end if         end if
     ENDIF  
206    
     if (inirnpb) THEN  
207         ! Initialisation du traceur dans le sol (couche limite radonique)         ! Initialisation du traceur dans le sol (couche limite radonique)
208         radio(1)= .true.         radio(1)= .true.
209         radio(2)= .true.         radio(2)= .true.
210         clsol(1)= .true.         clsol(:2)= .true.
211         clsol(2)= .true.         clsol(3:)= .false.
212         aerosol(2) = .TRUE. ! le Pb est un aerosol         aerosol(2) = .TRUE. ! le Pb est un aerosol
213         call initrrnpb(pctsrf, masktr, fshtr, hsoltr, tautr, vdeptr, scavtr)         call initrrnpb(pctsrf, masktr, fshtr, hsoltr, tautr, vdeptr, scavtr)
        inirnpb=.false.  
214      endif      endif
215    
216      if (convection) then      if (convection) then
# Line 285  contains Line 269  contains
269    
270      ! Calcul de l'effet de la couche limite      ! Calcul de l'effet de la couche limite
271    
272      if (couchelimite) then      DO k = 1, llm
273         DO k = 1, llm         DO i = 1, klon
274            DO i = 1, klon            delp(i, k) = paprs(i, k)-paprs(i, k+1)
              delp(i, k) = paprs(i, k)-paprs(i, k+1)  
           ENDDO  
275         ENDDO         ENDDO
276        ENDDO
277    
278         ! MAF modif pour tenir compte du cas traceur      ! MAF modif pour tenir compte du cas traceur
279         DO it=1, nqmx - 2      DO it=1, nqmx - 2
280            if (clsol(it)) then         if (clsol(it)) then
281               ! couche limite avec quantite dans le sol calculee            ! couche limite avec quantite dans le sol calculee
282               CALL cltracrn(it, pdtphys, yu1, yv1, coefh, t_seri, ftsol, &            CALL cltracrn(it, pdtphys, yu1, yv1, coefh, t_seri, ftsol, &
283                    pctsrf, tr_seri(:, :, it), trs(:, it), paprs, pplay, delp, &                 pctsrf, tr_seri(:, :, it), trs(:, it), paprs, pplay, delp, &
284                    masktr(1, it), fshtr(1, it), hsoltr(it), tautr(it), &                 masktr(1, it), fshtr(1, it), hsoltr(it), tautr(it), &
285                    vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs)                 vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs)
286               DO k = 1, llm            DO k = 1, llm
287                  DO i = 1, klon               DO i = 1, klon
288                     tr_seri(i, k, it) = tr_seri(i, k, it) + d_tr_cl(i, k, it)                  tr_seri(i, k, it) = tr_seri(i, k, it) + d_tr_cl(i, k, it)
                 ENDDO  
289               ENDDO               ENDDO
290              ENDDO
291    
292               trs(:, it) = trs(:, it) + d_trs            trs(:, it) = trs(:, it) + d_trs
293            else         else
294               ! couche limite avec flux prescrit            ! couche limite avec flux prescrit
295               !MAF provisoire source / traceur a creer            !MAF provisoire source / traceur a creer
296               DO i=1, klon            DO i=1, klon
297                  source(i) = 0. ! pas de source, pour l'instant               source(i) = 0. ! pas de source, pour l'instant
298               ENDDO            ENDDO
299    
300               CALL cltrac(pdtphys, coefh, t_seri, tr_seri(:, :, it), source, &            CALL cltrac(pdtphys, coefh, t_seri, tr_seri(:, :, it), source, &
301                    paprs, pplay, delp, d_tr_cl(1, 1, it))                 paprs, pplay, delp, d_tr_cl(1, 1, it))
302               DO k = 1, llm            DO k = 1, llm
303                  DO i = 1, klon               DO i = 1, klon
304                     tr_seri(i, k, it) = tr_seri(i, k, it) + d_tr_cl(i, k, it)                  tr_seri(i, k, it) = tr_seri(i, k, it) + d_tr_cl(i, k, it)
                 ENDDO  
305               ENDDO               ENDDO
306            endif            ENDDO
307         ENDDO         endif
308      endif      ENDDO
309    
310      ! Calcul de l'effet du puits radioactif      ! Calcul de l'effet du puits radioactif
311    
Legend:
Removed from v.224  
changed lines
  Added in v.225
  ViewVC Help
Powered by ViewVC 1.1.21