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

Diff of /trunk/phylmd/Interface_surf/clvent.f

Parent Directory | Revision Log | View Patch Patch

-trunk/libf/phylmd/clvent.f90
revision 70 by guez,
Mon Jun 24 15:39:52 2013 UTC
+trunk/phylmd/Interface_surf/clvent.f
revision 311 by guez,
Mon Dec  3 17:52:21 2018 UTC
 Line 4 
 module clvent_m
  contains
-   SUBROUTINE clvent(knon, dtime, u1lay, v1lay, coef, t, ven, paprs, pplay, &
+   SUBROUTINE clvent(u1lay, v1lay, coef, cdrag, t, ven, paprs, pplay, delp, &
-        delp, d_ven, flux_v)
+        d_ven, flux_v)
      ! Author: Z. X. Li (LMD/CNRS)
      ! Date: 1993/08/18
      ! Objet : diffusion verticale de la vitesse
-     USE dimphy, ONLY: klev, klon
+     ! Library:
+     use nr_util, only: assert
+     use comconst, only: dtphys
+     USE dimphy, ONLY: klev
      USE suphec_m, ONLY: rd, rg
-     INTEGER knon
+     REAL, intent(in):: u1lay(:), v1lay(:) ! (knon)
-     REAL, intent(in):: dtime
+     ! vent de la premiere couche (m / s)
-     ! dtime----input-R- intervalle du temps (en second)
+     REAL, intent(in):: coef(:, 2:) ! (knon, 2:klev)
-     REAL u1lay(klon), v1lay(klon)
+     ! Coefficient d'echange (m**2 / s) multiplié par le cisaillement du
-     ! u1lay----input-R- vent u de la premiere couche (m/s)
+     ! vent (dV / dz)
-     ! v1lay----input-R- vent v de la premiere couche (m/s)
+     REAL, intent(in):: cdrag(:) ! (knon) sans dimension
-     REAL, intent(in):: coef(:, :) ! (knon, klev)
+     REAL, intent(in):: t(:, :) ! (knon, klev) ! temperature (K)
-     ! Coefficient d'echange (m**2/s) multiplié par le cisaillement du
+     REAL, intent(in):: ven(:, :) ! (knon, klev) vitesse horizontale (m / s)
-     ! vent (dV/dz). La première valeur indique la valeur de Cdrag (sans
+     REAL, intent(in):: paprs(:, :) ! (knon, klev + 1) pression a
-     ! unité).
+     ! inter-couche (Pa)
+     real, intent(in):: pplay(:, :) ! (knon, klev) pression au milieu
-     REAL t(klon, klev), ven(klon, klev)
+     ! de couche (Pa)
-     ! t--------input-R- temperature (K)
+     real, intent(in):: delp(:, :) ! (knon, klev) epaisseur de couche (Pa)
-     ! ven------input-R- vitesse horizontale (m/s)
-     REAL paprs(klon, klev+1), pplay(klon, klev), delp(klon, klev)
+     REAL, intent(out):: d_ven(:, :) ! (knon, klev)
-     ! paprs----input-R- pression a inter-couche (Pa)
+     ! le changement de "ven", en m s-1
-     ! pplay----input-R- pression au milieu de couche (Pa)
-     ! delp-----input-R- epaisseur de couche (Pa)
+     REAL, intent(out):: flux_v(:) ! (knon)
-     REAL d_ven(klon, klev)
+     ! (diagnostic) stress du vent à la surface, en Pa
-     ! d_ven----output-R- le changement de "ven"
+     ! flux_v est le flux de quantité de mouvement (positif vers bas)
-     REAL flux_v(klon, klev)
-     ! flux_v---output-R- (diagnostic) flux du vent: (kg m/s)/(m**2 s)
      ! Local:
-     INTEGER i, k
+     INTEGER k
-     REAL zx_cv(klon, 2:klev)
+     REAL zx_cv(size(u1lay), 2:klev) ! (knon, 2:klev)
-     REAL zx_dv(klon, 2:klev)
+     REAL zx_dv(size(u1lay), 2:klev) ! (knon, 2:klev)
-     REAL zx_buf(klon)
+     REAL zx_buf(size(u1lay)) ! (knon)
-     REAL zx_coef(klon, klev)
+     REAL zx_coef(size(u1lay), klev) ! (knon, klev) in Pa
-     REAL local_ven(klon, klev)
+     REAL local_ven(size(u1lay), klev) ! (knon, klev) in m s-1
-     REAL zx_alf1(klon), zx_alf2(klon)
      !------------------------------------------------------------------
-     DO k = 1, klev
+     call assert(size(u1lay) == [size(v1lay), size(coef, 1), size(t, 1), &
-        DO i = 1, knon
+          size(ven, 1), size(paprs, 1), size(pplay, 1), size(delp, 1), &
-           local_ven(i, k) = ven(i, k)
+          size(d_ven, 1), size(flux_v)], "clvent knon")
-        ENDDO
-     ENDDO
+     zx_coef(:, 1) = cdrag * (1. + SQRT(u1lay**2 + v1lay**2)) * pplay(:, 1) &
+          / (RD * t(:, 1)) * dtphys * RG
-     DO i = 1, knon
-        zx_alf1(i) = 1.0
+     forall (k = 2:klev) zx_coef(:, k) = coef(:, k) * RG / (pplay(:, k - 1) &
-        zx_alf2(i) = 1.0 - zx_alf1(i)
+          - pplay(:, k)) * (paprs(:, k) * 2 / (t(:, k) + t(:, k - 1)) / RD)**2 &
-        zx_coef(i, 1) = coef(i, 1) * (1. + SQRT(u1lay(i)**2 + v1lay(i)**2)) &
+          * dtphys * RG
-             * pplay(i, 1) / (RD * t(i, 1))
-        zx_coef(i, 1) = zx_coef(i, 1) * dtime * RG
+     zx_buf = delp(:, 1) + zx_coef(:, 1) + zx_coef(:, 2)
-     ENDDO
+     zx_cv(:, 2) = ven(:, 1) * delp(:, 1) / zx_buf
+     zx_dv(:, 2) = zx_coef(:, 2) / zx_buf
-     DO k = 2, klev
-        DO i = 1, knon
-           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
-           zx_coef(i, k) = zx_coef(i, k) * dtime * RG
-        ENDDO
-     ENDDO
-     DO i = 1, knon
-        zx_buf(i) = delp(i, 1) + zx_coef(i, 1)*zx_alf1(i)+zx_coef(i, 2)
-        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)) &
-             /zx_buf(i)
-     ENDDO
      DO k = 3, klev
-        DO i = 1, knon
+        zx_buf = delp(:, k - 1) + zx_coef(:, k) &
-           zx_buf(i) = delp(i, k-1) + zx_coef(i, k) &
+             + zx_coef(:, k - 1) * (1. - zx_dv(:, k - 1))
-                + zx_coef(i, k-1)*(1.-zx_dv(i, k-1))
+        zx_cv(:, k) = (ven(:, k - 1) * delp(:, k - 1) &
-           zx_cv(i, k) = (local_ven(i, k-1)*delp(i, k-1) &
+             + zx_coef(:, k - 1) * zx_cv(:, k - 1)) / zx_buf
-                +zx_coef(i, k-1)*zx_cv(i, k-1) )/zx_buf(i)
+        zx_dv(:, k) = zx_coef(:, k) / zx_buf
-           zx_dv(i, k) = zx_coef(i, k)/zx_buf(i)
-        ENDDO
-     ENDDO
-     DO i = 1, knon
-        local_ven(i, klev) = ( local_ven(i, klev)*delp(i, klev) &
-             +zx_coef(i, klev)*zx_cv(i, klev) ) &
-             / ( 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
-     ! flux_v est le flux de moment angulaire (positif vers bas) dont
+     local_ven(:, klev) = (ven(:, klev) * delp(:, klev) &
-     ! l'unite est: (kg m/s)/(m**2 s)
+          + zx_coef(:, klev) * zx_cv(:, klev)) &
-     DO i = 1, knon
+          / (delp(:, klev) + zx_coef(:, klev) &
-        flux_v(i, 1) = zx_coef(i, 1)/(RG*dtime) &
+          - zx_coef(:, klev) * zx_dv(:, klev))
-             *(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
-     DO k = 1, klev
+     DO k = klev - 1, 1, - 1
-        DO i = 1, knon
+        local_ven(:, k) = zx_cv(:, k + 1) + zx_dv(:, k + 1) * local_ven(:, k + 1)
-           d_ven(i, k) = local_ven(i, k) - ven(i, k)
-        ENDDO
      ENDDO
+     flux_v = zx_coef(:, 1) / (RG * dtphys) * local_ven(:, 1)
+     d_ven = local_ven - ven
    END SUBROUTINE clvent
  end module clvent_m

 Legend:



Removed from v.70
 


changed lines


 
Added in v.311
 Legend:



Removed from v.70
 


changed lines


 
Added in v.311
-Removed from v.70
+Added in v.311

	ViewVC Help
Powered by ViewVC 1.1.21