Changeset 899 for codes/icosagcm/trunk/src/transport
- Timestamp:
- 06/13/19 16:45:41 (5 years ago)
- Location:
- codes/icosagcm/trunk/src/transport
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/transport/advect.F90
r896 r899 59 59 REAL(rstd),INTENT(OUT) :: gradq3d(iim*jjm,llm,3) 60 60 REAL(rstd) :: maxq,minq,minq_c,maxq_c 61 REAL(rstd) :: alphamx,alphami,alpha ,maggrd 62 REAL(rstd) :: leng1,leng2 61 REAL(rstd) :: alphamx,alphami,alpha,maggrd 63 62 REAL(rstd) :: arr(2*iim*jjm) 64 63 REAL(rstd) :: ar(iim*jjm) 65 64 REAL(rstd) :: gradtri(2*iim*jjm,llm,3) 66 INTEGER :: ij,k, ind,l65 INTEGER :: ij,k,l 67 66 REAL(rstd) :: detx,dety,detz,det 68 67 REAL(rstd) :: A(3,3), a11,a12,a13,a21,a22,a23,a31,a32,a33 … … 243 242 !DIR$ SIMD 244 243 DO ij=ij_begin,ij_end 245 ! maggrd = dot_product (gradq3d(ij,l,:),gradq3d(ij,l,:))244 ! maggrd = dot_product_3d(gradq3d(ij,l,:),gradq3d(ij,l,:)) 246 245 maggrd = gradq3d(ij,l,1)*gradq3d(ij,l,1) + gradq3d(ij,l,2)*gradq3d(ij,l,2) + gradq3d(ij,l,3)*gradq3d(ij,l,3) 247 246 maggrd = sqrt(maggrd) … … 275 274 REAL(rstd), INTENT(OUT) :: dq1,dq2,dq3,det 276 275 REAL(rstd) :: dq(3) 277 278 REAL(rstd) ::detx,dety,detz279 INTEGER :: info280 INTEGER :: IPIV(3)281 276 282 277 REAL(rstd) :: A(3,3) 283 REAL(rstd) :: B(3)284 278 285 279 ! TODO : replace A by A1,A2,A3 … … 345 339 REAL(rstd),INTENT(IN) :: tau 346 340 347 REAL(rstd) :: v_e(3), up_e , qe, ed(3)341 REAL(rstd) :: v_e(3), up_e 348 342 INTEGER :: ij,l 349 343 … … 422 416 REAL(rstd), INTENT(INOUT) :: qfluxt(3*iim*jjm,MERGE(llm,1,diagflux_on)) ! time-integrated tracer flux 423 417 424 REAL(rstd) :: dq,dmass,qe, ed(3),newmass418 REAL(rstd) :: dq,dmass,qe, newmass 425 419 REAL(rstd) :: qflux(3*iim*jjm,llm) 426 INTEGER :: ij, k,l420 INTEGER :: ij,l 427 421 428 422 CALL trace_start("compute_advect_horiz") -
codes/icosagcm/trunk/src/transport/advect_tracer.f90
r599 r899 242 242 243 243 REAL(rstd) :: dzqmax, newmass, sigw, qq, w 244 INTEGER :: i ,ij,l,j,ijb,ije244 INTEGER :: ij,l,ijb,ije 245 245 246 246 CALL trace_start("vlz")
Note: See TracChangeset
for help on using the changeset viewer.