Changeset 911 for codes/icosagcm/devel/src
- Timestamp:
- 06/17/19 00:31:27 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/initial/etat0_venus.f90
r741 r911 168 168 DO ij=1,ngrid 169 169 X_ij = 1./B(ij,1) 170 C(ij, 2)= -A(ij,2) * X_ij170 C(ij,1) = -A(ij,2) * X_ij 171 171 xu(ij,1) = Ru(ij,1) * X_ij 172 172 xv(ij,1) = Rv(ij,1) * X_ij … … 175 175 !DIR$ SIMD 176 176 DO ij=1,ngrid 177 X_ij = 1./( B(ij,l) + A(ij,l)*C(ij,l ) )178 C(ij,l +1)= -A(ij,l+1) * X_ij ! zero for l=llm177 X_ij = 1./( B(ij,l) + A(ij,l)*C(ij,l-1) ) 178 C(ij,l) = -A(ij,l+1) * X_ij ! zero for l=llm 179 179 xu(ij,l) = (Ru(ij,l)+A(ij,l)*xu(ij,l-1)) * X_ij 180 180 xv(ij,l) = (Rv(ij,l)+A(ij,l)*xv(ij,l-1)) * X_ij
Note: See TracChangeset
for help on using the changeset viewer.