Changeset 5079 for LMDZ6/trunk/libf/filtrez
- Timestamp:
- Jul 19, 2024, 11:28:59 AM (5 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/filtrez/inifgn.F
r2598 r5079 28 28 pi = 2.* ASIN(1.) 29 29 C 30 DO 5i=1,iim30 DO i=1,iim 31 31 dlonu(i)= xprimu( i ) 32 32 dlonv(i)= xprimv( i ) 33 5 CONTINUE33 END DO 34 34 35 DO 12i=1,iim35 DO i=1,iim 36 36 sddv(i) = SQRT(dlonv(i)) 37 37 sddu(i) = SQRT(dlonu(i)) 38 38 unsddu(i) = 1./sddu(i) 39 39 unsddv(i) = 1./sddv(i) 40 12 CONTINUE40 END DO 41 41 C 42 DO 17j=1,iim43 DO 17i=1,iim42 DO j=1,iim 43 DO i=1,iim 44 44 vec(i,j) = 0. 45 45 vec1(i,j) = 0. 46 46 eignfnv(i,j) = 0. 47 47 eignfnu(i,j) = 0. 48 17 CONTINUE 48 END DO 49 END DO 49 50 c 50 51 c 51 52 eignfnv(1,1) = -1. 52 53 eignfnv(iim,1) = 1. 53 DO 20i=1,imm154 DO i=1,imm1 54 55 eignfnv(i+1,i+1)= -1. 55 56 eignfnv(i,i+1) = 1. 56 20 CONTINUE57 DO 25j=1,iim58 DO 25i=1,iim57 END DO 58 DO j=1,iim 59 DO i=1,iim 59 60 eignfnv(i,j) = eignfnv(i,j)/(sddu(i)*sddv(j)) 60 25 CONTINUE 61 DO 30 j=1,iim 62 DO 30 i=1,iim 61 END DO 62 END DO 63 DO j=1,iim 64 DO i=1,iim 63 65 eignfnu(i,j) = -eignfnv(j,i) 64 30 CONTINUE 66 END DO 67 END DO 65 68 c 66 69 #ifdef CRAY
Note: See TracChangeset
for help on using the changeset viewer.