Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv.f90
- Timestamp:
- Aug 2, 2024, 9:58:25 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv.f90
r5134 r5159 2 2 3 3 SUBROUTINE gradiv(klevel, xcov, ycov, ld, gdx, gdy) 4 ! 4 5 5 ! Auteur : P. Le Van 6 ! 6 7 7 ! *************************************************************** 8 ! 8 9 9 ! ld 10 10 ! calcul de (grad (div) ) du vect. v .... 11 ! 11 12 12 ! xcov et ycov etant les composant.covariantes de v 13 13 ! **************************************************************** 14 14 ! xcov , ycov et ld sont des arguments d'entree pour le s-prog 15 15 ! gdx et gdy sont des arguments de sortie pour le s-prog 16 ! 17 ! 16 17 18 18 USE lmdz_filtreg, ONLY: filtreg 19 19 USE lmdz_ssum_scopy, ONLY: scopy 20 20 USE lmdz_comdissipn, ONLY: tetaudiv, tetaurot, tetah, cdivu, crot, cdivh 21 21 22 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 23 USE lmdz_paramet 22 24 IMPLICIT NONE 23 25 ! 24 INCLUDE "dimensions.h" 25 INCLUDE "paramet.h" 26 27 26 28 27 29 INTEGER :: klevel 28 ! 30 29 31 REAL :: xcov(ip1jmp1, klevel), ycov(ip1jm, klevel) 30 32 REAL :: gdx(ip1jmp1, klevel), gdy(ip1jm, klevel) … … 33 35 34 36 INTEGER :: l, ij, iter, ld 35 ! 36 ! 37 ! 37 38 39 38 40 CALL SCOPY(ip1jmp1 * klevel, xcov, 1, gdx, 1) 39 41 CALL SCOPY(ip1jm * klevel, ycov, 1, gdy, 1) 40 ! 42 41 43 DO iter = 1, ld 42 ! 44 43 45 CALL diverg(klevel, gdx, gdy, div) 44 46 CALL filtreg(div, jjp1, klevel, 2, 1, .TRUE., 2) 45 47 CALL grad(klevel, div, gdx, gdy) 46 ! 48 47 49 DO l = 1, klevel 48 50 DO ij = 1, ip1jmp1 … … 53 55 END DO 54 56 END DO 55 ! 57 56 58 END DO 57 59 RETURN
Note: See TracChangeset
for help on using the changeset viewer.