Changeset 5246 for LMDZ6/trunk/libf/dyn3d_common/gr_u_scal.f90
- Timestamp:
- Oct 21, 2024, 2:58:45 PM (22 hours ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d_common/gr_u_scal.f90
r5245 r5246 2 2 ! $Header$ 3 3 ! 4 5 c%W% %G%6 c=======================================================================7 c 8 cAuthor: Frederic Hourdin original: 11/11/929 c-------10 c 11 cSubject:12 c------13 c 14 cMethod:15 c--------16 c 17 cInterface:18 c----------19 c 20 cInput:21 c------22 c 23 cOutput:24 c-------25 c 26 c=======================================================================27 28 c-----------------------------------------------------------------------29 cDeclararations:30 c---------------4 SUBROUTINE gr_u_scal(nx,x_u,x_scal) 5 !%W% %G% 6 !======================================================================= 7 ! 8 ! Author: Frederic Hourdin original: 11/11/92 9 ! ------- 10 ! 11 ! Subject: 12 ! ------ 13 ! 14 ! Method: 15 ! -------- 16 ! 17 ! Interface: 18 ! ---------- 19 ! 20 ! Input: 21 ! ------ 22 ! 23 ! Output: 24 ! ------- 25 ! 26 !======================================================================= 27 IMPLICIT NONE 28 !----------------------------------------------------------------------- 29 ! Declararations: 30 ! --------------- 31 31 32 33 34 32 INCLUDE "dimensions.h" 33 INCLUDE "paramet.h" 34 INCLUDE "comgeom.h" 35 35 36 cArguments:37 c----------36 ! Arguments: 37 ! ---------- 38 38 39 INTEGERnx40 REALx_u(ip1jmp1,nx),x_scal(ip1jmp1,nx)39 INTEGER :: nx 40 REAL :: x_u(ip1jmp1,nx),x_scal(ip1jmp1,nx) 41 41 42 cLocal:43 c------42 ! Local: 43 ! ------ 44 44 45 INTEGERl,ij45 INTEGER :: l,ij 46 46 47 c-----------------------------------------------------------------------47 !----------------------------------------------------------------------- 48 48 49 50 51 x_scal(ij,l)=52 s (aireu(ij)*x_u(ij,l)+aireu(ij-1)*x_u(ij-1,l))53 s/(aireu(ij)+aireu(ij-1))54 55 49 DO l=1,nx 50 DO ij=ip1jmp1,2,-1 51 x_scal(ij,l)= & 52 (aireu(ij)*x_u(ij,l)+aireu(ij-1)*x_u(ij-1,l)) & 53 /(aireu(ij)+aireu(ij-1)) 54 ENDDO 55 ENDDO 56 56 57 57 CALL SCOPY(nx*jjp1,x_scal(iip1,1),iip1,x_scal(1,1),iip1) 58 58 59 60 END 59 RETURN 60 END SUBROUTINE gr_u_scal
Note: See TracChangeset
for help on using the changeset viewer.