/[lmdze]/trunk/libf/dyn3d/dissip.f90

Diff of /trunk/libf/dyn3d/dissip.f90

Parent Directory | Revision Log | View Patch Patch

-revision 44 by guez,
Wed Apr 13 12:29:18 2011 UTC
+revision 54 by guez,
Tue Dec  6 15:07:04 2011 UTC
 Line 1
- SUBROUTINE dissip(vcov, ucov, teta, p, dv, du, dh)
+ module dissip_m
-   ! From dyn3d/dissip.F, version 1.1.1.1 2004/05/19 12:53:05
-   ! Avec nouveaux operateurs star : gradiv2 , divgrad2, nxgraro2
-   ! Author: P. Le Van
-   ! Objet : dissipation horizontale
-   USE dimens_m, ONLY : llm
-   USE paramet_m, ONLY : iip1, iip2, ip1jm, ip1jmp1, llmp1
-   USE comdissnew, ONLY : lstardis, nitergdiv, nitergrot, niterh
-   USE inidissip_m, ONLY : dtdiss, tetah, tetaudiv, tetaurot
    IMPLICIT NONE
-   ! Arguments:
+ contains
-   REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm), teta(ip1jmp1, llm)
-   REAL, INTENT (IN) :: p(ip1jmp1, llmp1)
+   SUBROUTINE dissip(vcov, ucov, teta, p, dv, du, dh)
-   REAL dv(ip1jm, llm), du(ip1jmp1, llm), dh(ip1jmp1, llm)
+     ! From dyn3d/dissip.F, version 1.1.1.1 2004/05/19 12:53:05
-   ! Local:
+     ! Avec nouveaux opérateurs star : gradiv2, divgrad2, nxgraro2
-   REAL gdx(ip1jmp1, llm), gdy(ip1jm, llm)
+     ! Author: P. Le Van
-   REAL grx(ip1jmp1, llm), gry(ip1jm, llm)
+     ! Objet : dissipation horizontale
-   REAL te1dt(llm), te2dt(llm), te3dt(llm)
-   REAL deltapres(ip1jmp1, llm)
+     USE dimens_m, ONLY : iim, jjm, llm
+     USE paramet_m, ONLY : iip1, iip2, ip1jmp1, llmp1
-   INTEGER l, ij
+     USE comdissnew, ONLY : lstardis, nitergdiv, nitergrot, niterh
+     USE inidissip_m, ONLY : dtdiss, tetah, tetaudiv, tetaurot, cdivu, crot, &
-   !-----------------------------------------------------------------------
+          cdivh
+     use gradiv2_m, only: gradiv2
-   ! initialisations:
+     ! Arguments:
-   DO l = 1, llm
+     REAL vcov((iim + 1) * jjm, llm), ucov(ip1jmp1, llm), teta(ip1jmp1, llm)
-      te1dt(l) = tetaudiv(l)*dtdiss
+     REAL, INTENT (IN) :: p(ip1jmp1, llmp1)
-      te2dt(l) = tetaurot(l)*dtdiss
+     REAL dv((iim + 1) * jjm, llm), du(ip1jmp1, llm), dh(ip1jmp1, llm)
-      te3dt(l) = tetah(l)*dtdiss
-   END DO
+     ! Local:
-   du = 0.
+     REAL gdx(ip1jmp1, llm), gdy((iim + 1) * jjm, llm)
-   dv = 0.
+     REAL grx(ip1jmp1, llm), gry((iim + 1) * jjm, llm)
-   dh = 0.
+     REAL te1dt(llm), te2dt(llm), te3dt(llm)
+     REAL deltapres(ip1jmp1, llm)
-   ! Calcul de la dissipation:
+     INTEGER l, ij
-   ! Calcul de la partie grad (div) :
+     !-----------------------------------------------------------------------
-   IF (lstardis) THEN
-      CALL gradiv2(llm, ucov, vcov, nitergdiv, gdx, gdy)
+     ! Initializations:
-   ELSE
+     te1dt = tetaudiv * dtdiss
-      CALL gradiv(llm, ucov, vcov, nitergdiv, gdx, gdy)
+     te2dt = tetaurot * dtdiss
-   END IF
+     te3dt = tetah * dtdiss
+     du = 0.
-   DO l = 1, llm
+     dv = 0.
-      DO ij = 1, iip1
+     dh = 0.
-         gdx(ij, l) = 0.
-         gdx(ij+ip1jm, l) = 0.
+     ! Calcul de la dissipation:
-      END DO
+     ! Calcul de la partie grad (div) :
-      DO ij = iip2, ip1jm
-         du(ij, l) = du(ij, l) - te1dt(l)*gdx(ij, l)
+     IF (lstardis) THEN
-      END DO
+        CALL gradiv2(llm, ucov, vcov, nitergdiv, gdx, gdy, cdivu)
-      DO ij = 1, ip1jm
+     ELSE
-         dv(ij, l) = dv(ij, l) - te1dt(l)*gdy(ij, l)
+        CALL gradiv(llm, ucov, vcov, nitergdiv, gdx, gdy, cdivu)
-      END DO
+     END IF
-   END DO
+     DO l = 1, llm
-   ! calcul de la partie n X grad (rot) :
+        DO ij = 1, iip1
+           gdx(ij, l) = 0.
-   IF (lstardis) THEN
+           gdx(ij+(iim + 1) * jjm, l) = 0.
-      CALL nxgraro2(llm, ucov, vcov, nitergrot, grx, gry)
+        END DO
-   ELSE
-      CALL nxgrarot(llm, ucov, vcov, nitergrot, grx, gry)
+        DO ij = iip2, (iim + 1) * jjm
-   END IF
+           du(ij, l) = du(ij, l) - te1dt(l) * gdx(ij, l)
+        END DO
+        DO ij = 1, (iim + 1) * jjm
-   DO l = 1, llm
+           dv(ij, l) = dv(ij, l) - te1dt(l) * gdy(ij, l)
-      DO ij = 1, iip1
+        END DO
-         grx(ij, l) = 0.
+     END DO
-      END DO
+     ! calcul de la partie n X grad (rot) :
-      DO ij = iip2, ip1jm
-         du(ij, l) = du(ij, l) - te2dt(l)*grx(ij, l)
+     IF (lstardis) THEN
-      END DO
+        CALL nxgraro2(llm, ucov, vcov, nitergrot, grx, gry, crot)
-      DO ij = 1, ip1jm
+     ELSE
-         dv(ij, l) = dv(ij, l) - te2dt(l)*gry(ij, l)
+        CALL nxgrarot(llm, ucov, vcov, nitergrot, grx, gry, crot)
-      END DO
+     END IF
-   END DO
-   ! calcul de la partie div (grad) :
+     DO l = 1, llm
+        DO ij = 1, iip1
-   IF (lstardis) THEN
+           grx(ij, l) = 0.
-      DO l = 1, llm
+        END DO
-         DO ij = 1, ip1jmp1
-            deltapres(ij, l) = amax1(0., p(ij, l)-p(ij, l+1))
+        DO ij = iip2, (iim + 1) * jjm
-         END DO
+           du(ij, l) = du(ij, l) - te2dt(l) * grx(ij, l)
-      END DO
+        END DO
+        DO ij = 1, (iim + 1) * jjm
-      CALL divgrad2(llm, teta, deltapres, niterh, gdx)
+           dv(ij, l) = dv(ij, l) - te2dt(l) * gry(ij, l)
-   ELSE
+        END DO
-      CALL divgrad(llm, teta, niterh, gdx)
+     END DO
-   END IF
+     ! calcul de la partie div (grad) :
-   DO l = 1, llm
-      DO ij = 1, ip1jmp1
+     IF (lstardis) THEN
-         dh(ij, l) = dh(ij, l) - te3dt(l)*gdx(ij, l)
+        DO l = 1, llm
-      END DO
+           DO ij = 1, ip1jmp1
-   END DO
+              deltapres(ij, l) = max(0., p(ij, l) - p(ij, l + 1))
+           END DO
+        END DO
+        CALL divgrad2(llm, teta, deltapres, niterh, gdx, cdivh)
+     ELSE
+        CALL divgrad(llm, teta, niterh, gdx, cdivh)
+     END IF
+     forall (l = 1: llm) dh(:, l) = dh(:, l) - te3dt(l) * gdx(:, l)
+   END SUBROUTINE dissip
- END SUBROUTINE dissip
+ end module dissip_m

 Legend:



Removed from v.44
 


changed lines


 
Added in v.54
 Legend:



Removed from v.44
 


changed lines


 
Added in v.54
-Removed from v.44
+Added in v.54

	ViewVC Help
Powered by ViewVC 1.1.21