Changeset 1750 for LMDZ5/branches/testing/libf/phylmd/cltrac.F90

Timestamp:

Apr 25, 2013, 5:27:27 PM (11 years ago)

Author:

Laurent Fairhead

Message:

Version testing basée sur r1745

---

Testing release based on r1745

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/cltrac.F90 (modified) (5 diffs)

LMDZ5/branches/testing/libf/phylmd/cltrac.F90

@@ -2 +2 @@
 ! $Id $
 !
-SUBROUTINE cltrac(dtime,coef,t,tr,flux,paprs,pplay,delp,d_tr)
+SUBROUTINE cltrac(dtime,coef,t,tr,flux,paprs,pplay,delp, &
+                  d_tr,d_tr_dry,flux_tr_dry)                    !jyg
+
   USE dimphy
   IMPLICIT NONE
@@ -13 +15 @@
 ! Arguments:
 !-----------
-! dtime....input-R- intervalle du temps (en secondes)
-! coef.....input-R- le coefficient d'echange (m**2/s) l>1
-! t........input-R- temperature (K)
-! tr.......input-R- la q. de traceurs
-! flux.....input-R- le flux de traceurs a la surface
-! paprs....input-R- pression a inter-couche (Pa)
-! pplay....input-R- pression au milieu de couche (Pa)
-! delp.....input-R- epaisseur de couche (Pa)
-! cdrag....input-R- cdrag pour le flux de surface (non active)
-! tr0......input-R- traceurs a la surface ou dans l'ocean (non active)
-! d_tr.....output-R- le changement de tr
-! flux_tr..output-R- flux de tr
+! dtime.......input-R- intervalle du temps (en secondes)
+! coef........input-R- le coefficient d'echange (m**2/s) l>1
+! t...........input-R- temperature (K)
+! tr..........input-R- la q. de traceurs
+! flux........input-R- le flux de traceurs a la surface
+! paprs.......input-R- pression a inter-couche (Pa)
+! pplay.......input-R- pression au milieu de couche (Pa)
+! delp........input-R- epaisseur de couche (Pa)
+! cdrag.......input-R- cdrag pour le flux de surface (non active)
+! tr0.........input-R- traceurs a la surface ou dans l'ocean (non active)
+! d_tr........output-R- le changement de tr
+! d_tr_dry....output-R- le changement de tr du au depot sec (1st layer)
+! flux_tr_dry.output-R- depot sec
+!!! flux_tr..output-R- flux de tr
 !======================================================================
   include "YOMCST.h"
@@ -40 +44 @@
 !
   REAL ,DIMENSION(klon,klev),INTENT(OUT) :: d_tr
+  REAL ,DIMENSION(klon),INTENT(OUT)       :: d_tr_dry          !jyg
+  REAL ,DIMENSION(klon),INTENT(OUT)       :: flux_tr_dry       !jyg
 !  REAL ,DIMENSION(klon,klev),INTENT(OUT) :: flux_tr
 !
@@ -66 +72 @@
      zx_alf1(i) = (paprs(i,1)-pplay(i,2))/(pplay(i,1)-pplay(i,2))
      zx_alf2(i) = 1.0 - zx_alf1(i)
-     zx_flux(i) =  -flux(i)*dtime*RG
+     flux_tr_dry(i) = -flux(i)*dtime                              !jyg
+     zx_flux(i) =  flux_tr_dry(i)*RG                              !jyg
+!!     zx_flux(i) =  -flux(i)*dtime*RG                            !jyg
 ! Pour le moment le flux est prescrit cdrag et zx_coef(1) vaut 0
      cdrag(i) = 0.0 
@@ -95 +103 @@
      zx_dtr(i,2) = (zx_coef(i,2)-zx_alf2(i)*zx_coef(i,1)) /   & 
           zx_buf(i)
+     d_tr_dry(i) = -zx_flux(i)/zx_buf(i)                          !jyg
   ENDDO