Changeset 325 for trunk/LMDZ.MARS/libf/phymars/vdif_kc.F

Timestamp:

Oct 24, 2011, 3:56:49 PM (13 years ago)

Author:

acolaitis

Message:

Included variation of mean molar mass in potential temperature definition: this modification realistically accounts for mixing in the polar night when the atm is enriched in Ar. This follows present computations in convadj for the polar night case.

File:

• trunk/LMDZ.MARS/libf/phymars/vdif_kc.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/vdif_kc.F

@@ -1 @@
-      SUBROUTINE vdif_kc(dt,g,zlev,zlay,u,v,teta,cd,q2,km,kn)
+      SUBROUTINE vdif_kc(dt,g,zlev,zlay,u,v,teta,cd,q2,km,kn,zq)
       IMPLICIT NONE
 c.......................................................................
 #include "dimensions.h"
 #include "dimphys.h"
+#include "tracer.h"
+#include "callkeys.h"
 c.......................................................................
 c
@@ -37 +39 @@
       REAL km(ngridmx,nlayermx+1)
       REAL kn(ngridmx,nlayermx+1)
+      REAL zq(ngridmx,nlayermx,nqmx)
 c.......................................................................
 c
@@ -204 +207 @@
      &  cm4=-9.E+0 *a1*a2
      &          )
+
+c AC: variables for theta_m computation
+
+      INTEGER ico2,iq
+      SAVE ico2
+      REAL m_co2, m_noco2, A , B
+      SAVE A, B
+      LOGICAL firstcall
+      save firstcall
+      data firstcall/.true./
+      REAL zhc(ngridmx,nlayermx)
+c.......................................................................
+c  Initialization
+c.......................................................................
+
+      if(firstcall) then
+        ico2=0
+        if (tracer) then
+!     Prepare Special treatment if one of the tracers is CO2 gas
+           do iq=1,nqmx
+             if (noms(iq).eq."co2") then
+                ico2=iq
+                m_co2 = 44.01E-3  ! CO2 molecular mass (kg/mol)
+                m_noco2 = 33.37E-3  ! Non condensible mol mass (kg/mol)
+!               Compute A and B coefficient use to compute
+!               mean molecular mass Mair defined by
+!               1/Mair = q(ico2)/m_co2 + (1-q(ico2))/m_noco2
+!               1/Mair = A*q(ico2) + B
+                A =(1/m_co2 - 1/m_noco2)
+                B=1/m_noco2
+             end if
+           enddo
+        endif
+
+      firstcall=.false.
+      endif !of if firstcall
+
+c.......................................................................
+c  Special treatment for co2
+c.......................................................................
+
+      if (ico2.ne.0) then
+!     Special case if one of the tracers is CO2 gas
+         DO ilay=1,nlay
+           DO igrid=1,ngrid
+            zhc(igrid,ilay) = teta(igrid,ilay)*(A*zq(igrid,ilay,ico2)+B)
+           ENDDO
+         ENDDO
+       else
+          zhc(:,:)=teta(:,:)
+       end if
+
 c.......................................................................
 c  traitment des valeur de q2 en entree
@@ -272 +327 @@
 c
         n2(igrid,ilev)=g*unsdzdec(igrid,ilev)
-     &                   *(teta(igrid,ilev)-teta(igrid,ilev-1))
-     &                   /(teta(igrid,ilev)+teta(igrid,ilev-1)) *2.E+0
+     &                   *(zhc(igrid,ilev)-zhc(igrid,ilev-1))
+     &                   /(zhc(igrid,ilev)+zhc(igrid,ilev-1)) *2.E+0
 c
 c --->