- Timestamp:
- Jul 23, 2024, 7:14:34 PM (8 weeks ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/laplacien_gam.f90
r5104 r5105 2 2 ! $Header$ 3 3 4 SUBROUTINE laplacien_gam ( klevel, cuvsga, cvusga, unsaigam , 5 *unsapolnga, unsapolsga, teta, divgra )4 SUBROUTINE laplacien_gam ( klevel, cuvsga, cvusga, unsaigam , & 5 unsapolnga, unsapolsga, teta, divgra ) 6 6 7 cP. Le Van8 c 9 c************************************************************10 c 11 c.... calcul de (div( grad )) de teta .....12 c************************************************************13 cklevel et teta sont des arguments d'entree pour le s-prog14 cdivgra est un argument de sortie pour le s-prog15 c 16 17 c 18 19 20 7 ! P. Le Van 8 ! 9 ! ************************************************************ 10 ! 11 ! .... calcul de (div( grad )) de teta ..... 12 ! ************************************************************ 13 ! klevel et teta sont des arguments d'entree pour le s-prog 14 ! divgra est un argument de sortie pour le s-prog 15 ! 16 IMPLICIT NONE 17 ! 18 INCLUDE "dimensions.h" 19 INCLUDE "paramet.h" 20 INCLUDE "comgeom.h" 21 21 22 c 23 c............ variables en arguments ..........24 c 25 INTEGERklevel26 REALteta( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )27 REAL cuvsga(ip1jm) , cvusga( ip1jmp1 ),unsaigam(ip1jmp1),28 *unsapolnga, unsapolsga29 c 30 c........... variables locales .................31 c 32 REALghy(ip1jm,llm), ghx(ip1jmp1,llm)33 c......................................................22 ! 23 ! ............ variables en arguments .......... 24 ! 25 INTEGER :: klevel 26 REAL :: teta( ip1jmp1,klevel ), divgra( ip1jmp1,klevel ) 27 REAL :: cuvsga(ip1jm) , cvusga( ip1jmp1 ),unsaigam(ip1jmp1), & 28 unsapolnga, unsapolsga 29 ! 30 ! ........... variables locales ................. 31 ! 32 REAL :: ghy(ip1jm,llm), ghx(ip1jmp1,llm) 33 ! ...................................................... 34 34 35 c 36 c 37 c... cvuscugam = ( cvu/ cu ) ** (- gamdissip )38 c... cuvscvgam = ( cuv/ cv ) ** (- gamdissip ) calcules dans inigeom ..39 c... unsairegam = 1. / aire ** (- gamdissip )40 c 35 ! 36 ! 37 ! ... cvuscugam = ( cvu/ cu ) ** (- gamdissip ) 38 ! ... cuvscvgam = ( cuv/ cv ) ** (- gamdissip ) calcules dans inigeom .. 39 ! ... unsairegam = 1. / aire ** (- gamdissip ) 40 ! 41 41 42 43 c 44 45 c 46 CALL diverg_gam ( klevel, cuvsga, cvusga, unsaigam ,47 *unsapolnga, unsapolsga, ghx , ghy , divgra )42 CALL SCOPY ( ip1jmp1 * klevel, teta, 1, divgra, 1 ) 43 ! 44 CALL grad ( klevel, divgra, ghx, ghy ) 45 ! 46 CALL diverg_gam ( klevel, cuvsga, cvusga, unsaigam , & 47 unsapolnga, unsapolsga, ghx , ghy , divgra ) 48 48 49 c 49 ! 50 50 51 51 52 53 END 52 RETURN 53 END SUBROUTINE laplacien_gam
Note: See TracChangeset
for help on using the changeset viewer.