Changeset 1742 for LMDZ5/trunk/libf/phylmd/traclmdz_mod.F90

Timestamp:

Apr 5, 2013, 1:49:35 PM (11 years ago)

Author:

idelkadi

Message:

1- Inclusion des developpements de la these de Romain Pilon sur le
lessivage des aerosols :

a/ par les pluies convectives (modifs cv30_routines et cv3_routines pour

sortir les champs nécessaires au calcul off-line ; modif cvltr)

b/ par les pluies stratiformes (modifs phytrac et introduction

lsc_scav).

2- Choix entre plusieurs schemas pour les pluies stratiformes, commande
par iflag_lscav.

3- Quelques corrections dans la convection "Nouvelle Physique" pour
assurer la conservation des traceurs (cv3p1_mixing et cva_driver) (travail
de Robin Locatelli).

File:

• LMDZ5/trunk/libf/phylmd/traclmdz_mod.F90 (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/traclmdz_mod.F90

@@ -2 +2 @@
 !
 MODULE traclmdz_mod
+
 ! 
 ! In this module all tracers specific to LMDZ are treated. This module is used 
@@ -117 +118 @@
     REAL, DIMENSION(klev)          :: mintmp, maxtmp
     LOGICAL                        :: zero
-
+! RomP >>> profil initial Be7
+      integer ilesfil
+      parameter (ilesfil=1)
+      integer  irr,kradio
+      real     beryllium(klon,klev)
+! profil initial Pb210
+      integer ilesfil2
+      parameter (ilesfil2=1)
+      integer  irr2,kradio2
+      real     plomb(klon,klev)
+!! RomP <<<
 ! --------------------------------------------
 ! Allocation
@@ -148 +159 @@
 
     lessivage  = .TRUE.
+!!jyg(20130206) : le choix d activation du lessivage est fait dans phytrac avec iflag_lscav
+!!    call getin('lessivage',lessivage)
+!!    if(lessivage) then
+!!     print*,'lessivage lsc ON'
+!!    else
+!!     print*,'lessivage lsc OFF'
+!!    endif
     aerosol(:) = .FALSE.  ! Tous les traceurs sont des gaz par defaut
     
@@ -161 +179 @@
        ELSE IF ( tname(iiq) == "PB") THEN
           id_pb=it ! plomb
+! RomP >>> profil initial de PB210 
+     open (ilesfil2,file='prof.pb210',status='old',iostat=irr2)
+     IF (irr2 == 0) THEN 
+      read(ilesfil2,*) kradio2
+      print*,'number of levels for pb210 profile ',kradio2
+      do k=kradio2,1,-1
+       read (ilesfil2,*) plomb(:,k)
+      enddo
+      close(ilesfil2)
+      do k=1,klev
+       do i=1,klon
+         tr_seri(i,k,id_pb)=plomb(i,k)
+!!        print*, 'tr_seri',i,k,tr_seri(i,k,id_pb)
+        enddo
+      enddo
+     ELSE
+       print *, 'Prof.pb210 does not exist: use restart values'
+     ENDIF
+! RomP <<<
        ELSE IF ( tname(iiq) == "Aga" .OR. tname(iiq) == "AGA" ) THEN
           ! Age of stratospheric air
@@ -183 +220 @@
           radio(id_be) = .TRUE.
           aerosol(id_be) = .TRUE. ! le Be est un aerosol
-          CALL init_be(pctsrf,masktr(:,id_be),tautr(id_be),vdeptr(id_be),scavtr(id_be),srcbe)
+!jyg le 13/03/2013 ; ajout de pplay en argument de init_be
+!!!          CALL init_be(pctsrf,masktr(:,id_be),tautr(id_be),vdeptr(id_be),scavtr(id_be),srcbe)
+          CALL init_be(pctsrf,pplay,masktr(:,id_be),tautr(id_be),vdeptr(id_be),scavtr(id_be),srcbe)
           WRITE(lunout,*) 'Initialisation srcBe: OK'
+! RomP >>> profil initial de Be7
+      open (ilesfil,file='prof.be7',status='old',iostat=irr)
+      IF (irr == 0) THEN 
+       read(ilesfil,*) kradio
+       print*,'number of levels for Be7 profile ',kradio
+       do k=kradio,1,-1
+        read (ilesfil,*) beryllium(:,k)
+       enddo
+       close(ilesfil)
+       do k=1,klev
+         do i=1,klon
+          tr_seri(i,k,id_be)=beryllium(i,k)
+!!        print*, 'tr_seri',i,k,tr_seri(i,k,id_be)
+         enddo
+       enddo
+     ELSE
+       print *, 'Prof.Be7 does not exist: use restart values'
+     ENDIF
+! RomP <<<
        ELSE IF (tname(iiq)=="O3" .OR. tname(iiq)=="o3") THEN
           ! Recherche de l'ozone : parametrization de la chimie par Cariolle
@@ -280 +338 @@
        cdragh,  coefh, yu1, yv1, ftsol, pctsrf, xlat, xlon, couchelimite, sh, &
        rh, pphi, ustar, zu10m, zv10m, &
-       tr_seri, source, solsym, d_tr_cl, zmasse)
+!!          tr_seri, source, solsym, d_tr_cl, zmasse)                      !RomP
+          tr_seri, source, solsym, d_tr_cl,d_tr_dec, zmasse)               !RomP
     
     USE dimphy
@@ -316 +375 @@
 !--------------
 !
-    REAL,DIMENSION(klon),INTENT(IN)      :: cdragh  ! coeff drag pour T et Q
-    REAL,DIMENSION(klon,klev),INTENT(IN) :: coefh   ! diffusivite turb (m**2/s)
-    REAL,DIMENSION(klon),INTENT(IN)      :: yu1     ! vents au premier niveau
-    REAL,DIMENSION(klon),INTENT(IN)      :: yv1     ! vents au premier niveau
+    REAL,DIMENSION(klon),INTENT(IN)      :: cdragh     ! coeff drag pour T et Q
+    REAL,DIMENSION(klon,klev),INTENT(IN) :: coefh      ! coeff melange CL (m**2/s)
+    REAL,DIMENSION(klon),INTENT(IN)      :: yu1        ! vents au premier niveau
+    REAL,DIMENSION(klon),INTENT(IN)      :: yv1        ! vents au premier niveau
     LOGICAL,INTENT(IN)                   :: couchelimite
-    REAL,DIMENSION(klon,klev),INTENT(IN) :: sh      ! humidite specifique
+    REAL,DIMENSION(klon,klev),INTENT(IN) :: sh         ! humidite specifique
     REAL,DIMENSION(klon,klev),INTENT(IN) :: rh      ! Humidite relative
     REAL,DIMENSION(klon,klev),INTENT(IN) :: pphi    ! geopotentie