Changeset 1153 for palm/trunk/SOURCE/buoyancy.f90
- Timestamp:
- May 10, 2013 2:33:08 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/buoyancy.f90
r1132 r1153 20 20 ! Currrent revisions: 21 21 ! ----------------- 22 ! 22 ! code adjustment of accelerator version for PGI 12.3 / CUDA 5.0 23 23 ! 24 24 ! Former revisions: … … 238 238 DO i = i_left, i_right 239 239 DO j = j_south, j_north 240 !$acc loop vector(32) 241 DO k = 1, nzt-1 242 IF ( k > nzb_s_inner(j,i) ) THEN 243 tend(k,j,i) = tend(k,j,i) + atmos_ocean_sign * g * 0.5 * & 244 ( & 245 ( var(k,j,i) - hom(k,1,pr,0) ) / hom(k,1,pr,0) + & 246 ( var(k+1,j,i) - hom(k+1,1,pr,0) ) / hom(k+1,1,pr,0) & 247 ) 248 ENDIF 240 !$acc loop independent vector(32) 241 DO k = nzb_s_inner(j,i)+1, nzt-1 242 tend(k,j,i) = tend(k,j,i) + atmos_ocean_sign * g * 0.5 * & 243 ( & 244 ( var(k,j,i) - hom(k,1,pr,0) ) / hom(k,1,pr,0) + & 245 ( var(k+1,j,i) - hom(k+1,1,pr,0) ) / hom(k+1,1,pr,0) & 246 ) 249 247 ENDDO 250 248 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.