Changeset 4396 for LMDZ6/trunk/libf/phylmd/thermcell_down.F90

Timestamp:

Jan 24, 2023, 2:16:46 PM (17 months ago)

Author:

evignon

Message:

Attention les yeux, commission qui enterine la premiere version
operationnelle des downdrafts dans LMDZ. Speciale dedicace a Fleur

File:

• LMDZ6/trunk/libf/phylmd/thermcell_down.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/thermcell_down.F90

@@ -1 @@
-   SUBROUTINE thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,eup,dup,edn,ddn,masse,trac)
+   SUBROUTINE thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,eup,dup,edn,ddn,masse,trac,dtrac)
 
 !-----------------------------------------------------------------
@@ -25 +25 @@
    real,intent(in), dimension(ngrid,nlay) :: ddn ! detrainment from downdrafts * dz [same unit as flux]
    real,intent(in), dimension(ngrid,nlay) :: masse ! mass of layers = rho dz 
-   real,intent(inout), dimension(ngrid,nlay) :: trac ! tracer 
+   real,intent(in), dimension(ngrid,nlay) :: trac ! tracer 
    integer, intent(in), dimension(ngrid) :: lmax ! max level index at which downdraft are present
+   real,intent(out),dimension(ngrid,nlay) ::dtrac ! tendance du traceur
 
    
@@ -32 +33 @@
 
    real, dimension(ngrid,nlay+1) :: fup,fdn,fthu,fthd,fthe,fthtot
-   real, dimension(ngrid,nlay) :: tracu,tracd,dtrac
+   real, dimension(ngrid,nlay) :: tracu,tracd
    real :: www
    integer ig,ilay
@@ -68 +69 @@
    do ilay=1,nlay,1
       do ig=1,ngrid
-         if (ilay.le.lmax(ig) .and. lmax(ig)>1) then
+         if (ilay.lt.lmax(ig) .and. lmax(ig)>1) then
             fup(ig,ilay+1)=fup(ig,ilay)+eup(ig,ilay)-dup(ig,ilay)
             if (ilay == 1 ) then
@@ -83 +84 @@
       enddo 
    enddo
-   !Boucle pour calculer le flux up
+   !Boucle pour calculer le flux du traceur flux updraft, flux downdraft, flux compensatoire
    do ilay=2,nlay,1
      do ig=1,ngrid
@@ -95 +96 @@
        endif
        fthe(ig,ilay)=-(fup(ig,ilay)-fdn(ig,ilay))*trac(ig,ilay)
+       !! si on voulait le prendre en compte on
+       !fthe(ig,ilay)=-(fup(ig,ilay)-fdn(ig,ilay))*trac(ig,ilay-1)
        fthtot(ig,ilay)=fthu(ig,ilay)+fthd(ig,ilay)+fthe(ig,ilay)
      enddo
@@ -101 +104 @@
    do ilay=1,nlay,1
      do ig=1,ngrid
-       dtrac(ig,ilay)=(fthtot(ig,ilay)-fthtot(ig,ilay+1))*(ptimestep/masse(ig,ilay))
+       dtrac(ig,ilay)=(fthtot(ig,ilay)-fthtot(ig,ilay+1))*(1./masse(ig,ilay))
 !       trac(ig,ilay)=trac(ig,ilay) + (fthtot(ig,ilay)-fthtot(ig,ilay+1))*(ptimestep/masse(ig,ilay))
      enddo
    enddo
-   if (1==0) then
-    do ilay=1,nlay,1
-     do ig=1,ngrid
-       trac(ig,ilay)=trac(ig,ilay) + (fup(ig,ilay)*tracu(ig,ilay-1)-fup(ig,ilay+1)*tracu(ig,ilay) + &
-       & (fup(ig,ilay+1)+fdn(ig,ilay+1))*trac(ig,ilay+1) - (fup(ig,ilay)+fdn(ig,ilay))*trac(ig,ilay) + &
-       & fdn(ig,ilay+1)*tracd(ig,ilay+1)-fdn(ig,ilay)*tracd(ig,ilay))*(ptimestep/masse(ig,ilay))
-     enddo
-    enddo
-   endif
 ! Il reste a coder : 
 ! d(rho trac)/dt = - d/dz(rho w'trac')