Changeset 5246 for LMDZ6/trunk/libf/dyn3d_common/gr_v_scal.f90
- Timestamp:
- Oct 21, 2024, 2:58:45 PM (23 hours ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d_common/gr_v_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_v_scal(nx,x_v,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_v(ip1jm,nx),x_scal(ip1jmp1,nx)39 INTEGER :: nx 40 REAL :: x_v(ip1jm,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 (airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l))53 s/(airev(ij-iip1)+airev(ij))54 55 56 57 58 59 60 61 49 DO l=1,nx 50 DO ij=iip2,ip1jm 51 x_scal(ij,l)= & 52 (airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l)) & 53 /(airev(ij-iip1)+airev(ij)) 54 ENDDO 55 DO ij=1,iip1 56 x_scal(ij,l)=0. 57 ENDDO 58 DO ij=ip1jm+1,ip1jmp1 59 x_scal(ij,l)=0. 60 ENDDO 61 ENDDO 62 62 63 64 END 63 RETURN 64 END SUBROUTINE gr_v_scal
Note: See TracChangeset
for help on using the changeset viewer.