Changeset 1067 for LMDZ4/trunk/libf/phylmd/cltracrn.F

Timestamp:

Dec 17, 2008, 2:30:13 PM (16 years ago)

Author:

Laurent Fairhead

Message:

• Modifications lie au premier niveau du modele pour la diffusion turbulent

du vent.

• Preparation pour un couplage des courrant oceaniques.

JG

File:

• LMDZ4/trunk/libf/phylmd/cltracrn.F (modified) (4 diffs)

LMDZ4/trunk/libf/phylmd/cltracrn.F

@@ -1 @@
-!
-! $Header$
 !
       SUBROUTINE cltracrn( itr, dtime,u1lay, v1lay,
-     e              coef,t,ftsol,pctsrf,
+     e              cdrag,coef,t,ftsol,pctsrf,
      e              tr,trs,paprs,pplay,delp,
      e              masktr,fshtr,hsoltr,tautr,vdeptr,
@@ -26 +24 @@
 c u1lay----input-R- vent u de la premiere couche (m/s)
 c v1lay----input-R- vent v de la premiere couche (m/s)
-c coef-----input-R- le coefficient d'echange (m**2/s) l>1
+c cdrag----input-R- cdrag
+c coef-----input-R- le coefficient d'echange (m**2/s) l>1, valable uniquement pour k entre 2 et klev
 c t--------input-R- temperature (K)
 c paprs----input-R- pression a inter-couche (Pa)
@@ -50 +49 @@
       REAL dtime
       REAL u1lay(klon), v1lay(klon)
+      REAL cdrag(klon)
       REAL coef(klon,klev)
       REAL t(klon,klev), ftsol(klon,nbsrf), pctsrf(klon,nbsrf) 
@@ -119 +119 @@
 c======================================================================
       DO i = 1, klon
-         zx_coef(i,1) = coef(i,1)
+         zx_coef(i,1) = cdrag(i)
      .                 * (1.0+SQRT(u1lay(i)**2+v1lay(i)**2))
      .                 * pplay(i,1)/(RD*t(i,1))