Changeset 1985

Timestamp:

Feb 19, 2014, 6:53:24 AM (10 years ago)

Author:

fhourdin

Message:

Modification de du schéma implicite des thermique pour que la concentration

soit inchangée en cas de flux de masse uniformément nul.

Implicit scheme for thermal plume modified to keep unchanged the concentration

when the mass flux is zero.

File:

• LMDZ5/trunk/libf/phylmd/thermcell_dq.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/thermcell_dq.F90

@@ -40 +40 @@
 
 ! Old explicite scheme
-      if (impl==-1) then
-         call thermcell_dq_o(ngrid,nlay,ptimestep,fm,entr,  &
+      if (impl<=-1) then
+         call thermcell_dq_o(ngrid,nlay,impl,ptimestep,fm,entr,  &
      &           masse,q,dq,qa,lev_out)
          return
@@ -124 +124 @@
    else
       do k=nlay-1,1,-1
-         q(:,k)=(masse(:,k)*q(:,k)/ptimestep+fqa(:,k)-fqa(:,k+1)+fm(:,k+1)*q(:,k+1)) &
-     &               /(fm(:,k)+masse(:,k)/ptimestep)
+! FH debut de modif : le calcul ci dessous modifiait numériquement
+! la concentration quand le flux de masse etait nul car on divisait
+! puis multipliait par masse/ptimestep.
+!        q(:,k)=(masse(:,k)*q(:,k)/ptimestep+fqa(:,k)-fqa(:,k+1)+fm(:,k+1)*q(:,k+1)) &
+!    &               /(fm(:,k)+masse(:,k)/ptimestep)
+         q(:,k)=(q(:,k)+ptimestep/masse(:,k)*(fqa(:,k)-fqa(:,k+1)+fm(:,k+1)*q(:,k+1))) &
+      &               /(1.+fm(:,k)*ptimestep/masse(:,k))
+! FH fin de modif.
       enddo
    endif
@@ -144 +150 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-      subroutine thermcell_dq_o(ngrid,nlay,ptimestep,fm,entr,  &
+      subroutine thermcell_dq_o(ngrid,nlay,impl,ptimestep,fm,entr,  &
      &           masse,q,dq,qa,lev_out)
       implicit none
@@ -157 +163 @@
 !=======================================================================
 
-      integer ngrid,nlay
+      integer ngrid,nlay,impl
 
       real ptimestep