/[lmdze]/trunk/phylmd/clvent.f

Diff of /trunk/phylmd/clvent.f

Parent Directory | Revision Log | View Patch Patch

-trunk/libf/phylmd/clvent.f
revision 3 by guez,
Wed Feb 27 13:16:39 2008 UTC
+trunk/phylmd/clvent.f
revision 254 by guez,
Mon Feb  5 10:39:38 2018 UTC
 Line 1
-       SUBROUTINE clvent(knon,dtime, u1lay,v1lay,coef,t,ven,
+ module clvent_m
-      e                  paprs,pplay,delp,
-      s                  d_ven,flux_v)
+   IMPLICIT none
-       use dimens_m
-       use dimphy
+ contains
-       use iniprint
-       use YOMCST
+   SUBROUTINE clvent(dtime, u1lay, v1lay, coef, cdrag, t, ven, paprs, pplay, &
-       IMPLICIT none
+        delp, d_ven, flux_v)
- c======================================================================
- c Auteur(s): Z.X. Li (LMD/CNRS) date: 19930818
+     ! Author: Z. X. Li (LMD/CNRS)
- c Objet: diffusion vertical de la vitesse "ven"
+     ! Date: 1993/08/18
- c======================================================================
+     ! Objet : diffusion verticale de la vitesse
- c Arguments:
- c dtime----input-R- intervalle du temps (en second)
+     USE dimphy, ONLY: klev
- c u1lay----input-R- vent u de la premiere couche (m/s)
+     use nr_util, only: assert
- c v1lay----input-R- vent v de la premiere couche (m/s)
+     USE suphec_m, ONLY: rd, rg
- c coef-----input-R- le coefficient d'echange (m**2/s) multiplie par
- c                   le cisaillement du vent (dV/dz); la premiere
+     REAL, intent(in):: dtime ! intervalle de temps (en s)
- c                   valeur indique la valeur de Cdrag (sans unite)
- c t--------input-R- temperature (K)
+     REAL, intent(in):: u1lay(:), v1lay(:) ! (knon)
- c ven------input-R- vitesse horizontale (m/s)
+     ! vent de la premiere couche (m / s)
- c paprs----input-R- pression a inter-couche (Pa)
- c pplay----input-R- pression au milieu de couche (Pa)
+     REAL, intent(in):: coef(:, 2:) ! (knon, 2:klev)
- c delp-----input-R- epaisseur de couche (Pa)
+     ! Coefficient d'echange (m**2 / s) multiplié par le cisaillement du
- c
+     ! vent (dV / dz)
- c
- c d_ven----output-R- le changement de "ven"
+     REAL, intent(in):: cdrag(:) ! (knon) sans unité
- c flux_v---output-R- (diagnostic) flux du vent: (kg m/s)/(m**2 s)
+     REAL, intent(in):: t(:, :) ! (knon, klev) ! temperature (K)
- c======================================================================
+     REAL, intent(in):: ven(:, :) ! (knon, klev) vitesse horizontale (m / s)
-       INTEGER knon
+     REAL, intent(in):: paprs(:, :) ! (knon, klev + 1) pression a
-       REAL dtime
+     ! inter-couche (Pa)
-       REAL u1lay(klon), v1lay(klon)
+     real, intent(in):: pplay(:, :) ! (knon, klev) pression au milieu
-       REAL coef(klon,klev)
+     ! de couche (Pa)
-       REAL t(klon,klev), ven(klon,klev)
+     real, intent(in):: delp(:, :) ! (knon, klev) epaisseur de couche (Pa)
-       REAL paprs(klon,klev+1), pplay(klon,klev), delp(klon,klev)
+     REAL, intent(out):: d_ven(:, :) ! (knon, klev) ! le changement de "ven"
-       REAL d_ven(klon,klev)
-       REAL flux_v(klon,klev)
+     REAL, intent(out):: flux_v(:) ! (knon)
- c======================================================================
+     ! (diagnostic) flux du vent à la surface, en (kg m / s) / (m**2 s)
- c======================================================================
+     ! flux_v est le flux de moment angulaire (positif vers bas)
-       INTEGER i, k
-       REAL zx_cv(klon,2:klev)
+     ! Local:
-       REAL zx_dv(klon,2:klev)
+     INTEGER k
-       REAL zx_buf(klon)
+     REAL zx_cv(size(u1lay), 2:klev) ! (knon, 2:klev)
-       REAL zx_coef(klon,klev)
+     REAL zx_dv(size(u1lay), 2:klev) ! (knon, 2:klev)
-       REAL local_ven(klon,klev)
+     REAL zx_buf(size(u1lay)) ! (knon)
-       REAL zx_alf1(klon), zx_alf2(klon)
+     REAL zx_coef(size(u1lay), klev) ! (knon, klev)
- c======================================================================
+     REAL local_ven(size(u1lay), klev) ! (knon, klev)
-       DO k = 1, klev
-       DO i = 1, knon
+     !------------------------------------------------------------------
-          local_ven(i,k) = ven(i,k)
-       ENDDO
+     call assert(size(u1lay) == [size(v1lay), size(coef, 1), size(t, 1), &
-       ENDDO
+          size(ven, 1), size(paprs, 1), size(pplay, 1), size(delp, 1), &
- c======================================================================
+          size(d_ven, 1), size(flux_v)], "clvent knon")
-       DO i = 1, knon
- ccc         zx_alf1(i) = (paprs(i,1)-pplay(i,2))/(pplay(i,1)-pplay(i,2))
+     zx_coef(:, 1) = cdrag * (1. + SQRT(u1lay**2 + v1lay**2)) * pplay(:, 1) &
-          zx_alf1(i) = 1.0
+          / (RD * t(:, 1)) * dtime * RG
-          zx_alf2(i) = 1.0 - zx_alf1(i)
-          zx_coef(i,1) = coef(i,1)
+     DO k = 2, klev
-      .                 * (1.0+SQRT(u1lay(i)**2+v1lay(i)**2))
+        zx_coef(:, k) = coef(:, k) * RG / (pplay(:, k - 1) - pplay(:, k)) &
-      .                 * pplay(i,1)/(RD*t(i,1))
+             * (paprs(:, k) * 2 / (t(:, k) + t(:, k - 1)) / RD)**2 * dtime * RG
-          zx_coef(i,1) = zx_coef(i,1) * dtime*RG
+     ENDDO
-       ENDDO
- c======================================================================
+     zx_buf = delp(:, 1) + zx_coef(:, 1) + zx_coef(:, 2)
-       DO k = 2, klev
+     zx_cv(:, 2) = ven(:, 1) * delp(:, 1) / zx_buf
-       DO i = 1, knon
+     zx_dv(:, 2) = zx_coef(:, 2) / zx_buf
-          zx_coef(i,k) = coef(i,k)*RG/(pplay(i,k-1)-pplay(i,k))
-      .                  *(paprs(i,k)*2/(t(i,k)+t(i,k-1))/RD)**2
+     DO k = 3, klev
-          zx_coef(i,k) = zx_coef(i,k) * dtime*RG
+        zx_buf = delp(:, k - 1) + zx_coef(:, k) &
-       ENDDO
+             + zx_coef(:, k - 1) * (1. - zx_dv(:, k - 1))
-       ENDDO
+        zx_cv(:, k) = (ven(:, k - 1) * delp(:, k - 1) &
- c======================================================================
+             + zx_coef(:, k - 1) * zx_cv(:, k - 1)) / zx_buf
-       DO i = 1, knon
+        zx_dv(:, k) = zx_coef(:, k) / zx_buf
-          zx_buf(i) = delp(i,1) + zx_coef(i,1)*zx_alf1(i)+zx_coef(i,2)
+     ENDDO
-          zx_cv(i,2) = local_ven(i,1)*delp(i,1) / zx_buf(i)
-          zx_dv(i,2) = (zx_coef(i,2)-zx_alf2(i)*zx_coef(i,1))
+     local_ven(:, klev) = (ven(:, klev) * delp(:, klev) &
-      .                /zx_buf(i)
+          + zx_coef(:, klev) * zx_cv(:, klev)) &
-       ENDDO
+          / (delp(:, klev) + zx_coef(:, klev) &
-       DO k = 3, klev
+          - zx_coef(:, klev) * zx_dv(:, klev))
-       DO i = 1, knon
-          zx_buf(i) = delp(i,k-1) + zx_coef(i,k)
+     DO k = klev - 1, 1, - 1
-      .                         + zx_coef(i,k-1)*(1.-zx_dv(i,k-1))
+        local_ven(:, k) = zx_cv(:, k + 1) + zx_dv(:, k + 1) * local_ven(:, k + 1)
-          zx_cv(i,k) = (local_ven(i,k-1)*delp(i,k-1)
+     ENDDO
-      .                  +zx_coef(i,k-1)*zx_cv(i,k-1) )/zx_buf(i)
-          zx_dv(i,k) = zx_coef(i,k)/zx_buf(i)
+     flux_v = zx_coef(:, 1) / (RG * dtime) * local_ven(:, 1)
-       ENDDO
+     d_ven = local_ven - ven
-       ENDDO
-       DO i = 1, knon
+   END SUBROUTINE clvent
-          local_ven(i,klev) = ( local_ven(i,klev)*delp(i,klev)
-      .                        +zx_coef(i,klev)*zx_cv(i,klev) )
+ end module clvent_m
-      .                   / ( delp(i,klev) + zx_coef(i,klev)
-      .                       -zx_coef(i,klev)*zx_dv(i,klev) )
-       ENDDO
-       DO k = klev-1, 1, -1
-       DO i = 1, knon
-          local_ven(i,k) = zx_cv(i,k+1) + zx_dv(i,k+1)*local_ven(i,k+1)
-       ENDDO
-       ENDDO
- c======================================================================
- c== flux_v est le flux de moment angulaire (positif vers bas)
- c== dont l'unite est: (kg m/s)/(m**2 s)
-       DO i = 1, knon
-          flux_v(i,1) = zx_coef(i,1)/(RG*dtime)
-      .                 *(local_ven(i,1)*zx_alf1(i)
-      .                  +local_ven(i,2)*zx_alf2(i))
-       ENDDO
-       DO k = 2, klev
-       DO i = 1, knon
-          flux_v(i,k) = zx_coef(i,k)/(RG*dtime)
-      .               * (local_ven(i,k)-local_ven(i,k-1))
-       ENDDO
-       ENDDO
- c
-       DO k = 1, klev
-       DO i = 1, knon
-          d_ven(i,k) = local_ven(i,k) - ven(i,k)
-       ENDDO
-       ENDDO
- c
-       RETURN
-       END

 Legend:



Removed from v.3
 


changed lines


 
Added in v.254
 Legend:



Removed from v.3
 


changed lines


 
Added in v.254
-Removed from v.3
+Added in v.254

	ViewVC Help
Powered by ViewVC 1.1.21