Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/psextbar.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/psextbar.f90
r5136 r5159 5 5 USE lmdz_comgeom 6 6 7 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 8 USE lmdz_paramet 7 9 IMPLICIT NONE 8 10 9 11 !======================================================================= 10 ! 12 11 13 ! Auteur: P. Le Van 12 14 ! ------- 13 ! 15 14 16 ! Objet: 15 17 ! ------ 16 ! 18 17 19 ! ********************************************************************** 18 20 ! calcul des moyennes en x et en y de (pression au sol*aire variable) .. 19 21 ! ********************************************************************** 20 ! 22 21 23 ! ps est un argum. d'entree pour le s-pg .. 22 24 ! psexbarxy est un argum. de sortie pour le s-pg .. 23 ! 25 24 26 ! Methode: 25 27 ! -------- 26 ! 28 27 29 ! A chaque point scalaire P (i,j) est affecte 4 coefficients d'aires 28 30 ! alpha1(i,j) calcule au point ( i+1/4,j-1/4 ) … … 30 32 ! alpha3(i,j) calcule au point ( i-1/4,j+1/4 ) 31 33 ! alpha4(i,j) calcule au point ( i-1/4,j-1/4 ) 32 ! 34 33 35 ! Avec alpha1(i,j) = aire(i+1/4,j-1/4)/ aire(i,j) 34 ! 36 35 37 ! N.B . Pour plus de details, voir s-pg ... iniconst ... 36 ! 37 ! 38 ! 38 39 40 39 41 ! alpha4 . . alpha1 . alpha4 40 42 ! (i,j) (i,j) (i+1,j) 41 ! 43 42 44 ! P . U . . P 43 45 ! (i,j) (i,j) (i+1,j) 44 ! 46 45 47 ! alpha3 . . alpha2 .alpha3 46 48 ! (i,j) (i,j) (i+1,j) 47 ! 49 48 50 ! V . Z . . V 49 51 ! (i,j) 50 ! 52 51 53 ! alpha4 . . alpha1 .alpha4 52 54 ! (i,j+1) (i,j+1) (i+1,j+1) 53 ! 55 54 56 ! P . U . . P 55 57 ! (i,j+1) (i+1,j+1) 56 ! 57 ! 58 ! 59 ! 58 59 60 61 60 62 ! On a : 61 ! 63 62 64 ! pbarx(i,j) = Pext(i ,j) * ( alpha1(i ,j) + alpha2(i,j)) + 63 65 ! Pext(i+1,j) * ( alpha3(i+1,j) + alpha4(i+1,j) ) 64 66 ! localise au point ... U (i,j) ... 65 ! 67 66 68 ! pbary(i,j) = Pext(i,j ) * ( alpha2(i,j ) + alpha3(i,j ) + 67 69 ! Pext(i,j+1) * ( alpha1(i,j+1) + alpha4(i,j+1) 68 70 ! localise au point ... V (i,j) ... 69 ! 71 70 72 ! pbarxy(i,j)= Pext(i,j) *alpha2(i,j) + Pext(i+1,j) *alpha3(i+1,j) + 71 73 ! Pext(i,j+1)*alpha1(i,j+1)+ Pext(i+1,j+1)*alpha4(i+1,j+1) 72 74 ! localise au point ... Z (i,j) ... 73 ! 74 ! 75 ! 75 76 77 76 78 !======================================================================= 77 79 78 80 79 INCLUDE "dimensions.h" 80 INCLUDE "paramet.h" 81 82 81 83 82 84 REAL :: ps( ip1jmp1 ), psexbarxy ( ip1jm ), pext( ip1jmp1 )
Note: See TracChangeset
for help on using the changeset viewer.