--- trunk/libf/dyn3d/laplacien_gam.f 2008/02/27 13:16:39 3 +++ trunk/dyn3d/laplacien_gam.f 2013/11/15 18:45:49 76 @@ -13,6 +13,7 @@ c klevel et teta sont des arguments d'entree pour le s-prog c divgra est un argument de sortie pour le s-prog c + use grad_m, only: grad use dimens_m use paramet_m use comgeom @@ -22,7 +23,7 @@ c c ............ variables en arguments .......... c - INTEGER klevel + INTEGER, intent(in):: klevel REAL teta( ip1jmp1,klevel ), divgra( ip1jmp1,klevel ) REAL cuvsga(ip1jm) , cvusga( ip1jmp1 ),unsaigam(ip1jmp1), * unsapolnga, unsapolsga @@ -32,13 +33,6 @@ REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm) c ...................................................... -c -c -c ... cvuscugam = ( cvu/ cu ) ** (- gamdissip ) -c ... cuvscvgam = ( cuv/ cv ) ** (- gamdissip ) calcules dans inigeom .. -c ... unsairegam = 1. / aire ** (- gamdissip ) -c - CALL SCOPY ( ip1jmp1 * klevel, teta, 1, divgra, 1 ) c CALL grad ( klevel, divgra, ghx, ghy )