- Timestamp:
- 01/29/15 23:18:18 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/theta_rhodz.f90
r322 r323 163 163 REAL(rstd),INTENT(OUT) :: theta_rhodz(iim*jjm,llm) 164 164 INTEGER,INTENT(IN) :: offset 165 INTEGER :: i,j,ij,l 166 167 !$OMP BARRIER 168 DO l = ll_begin, ll_end 169 DO j=jj_begin-offset,jj_end+offset 170 DO i=ii_begin-offset,ii_end+offset 171 ij=(j-1)*iim+i 172 theta_rhodz(ij,l) = theta(ij,l) * ( (ap(l)-ap(l+1)) + ( bp(l)- bp(l+1))* ps(ij) )/g 165 REAL(rstd) :: rhodz 166 INTEGER :: i,j,ij,l 167 168 !$OMP BARRIER 169 DO l = ll_begin, ll_end 170 DO j=jj_begin-offset,jj_end+offset 171 DO i=ii_begin-offset,ii_end+offset 172 ij=(j-1)*iim+i 173 rhodz = ( ap(l)-ap(l+1) + (bp(l)-bp(l+1))*ps(ij) )/g 174 theta_rhodz(ij,l) = theta(ij,l) * rhodz 173 175 ENDDO 174 176 ENDDO … … 188 190 REAL(rstd),INTENT(OUT) :: theta(iim*jjm,llm) 189 191 INTEGER,INTENT(IN) :: offset 190 INTEGER :: i,j,ij,l 191 192 !$OMP BARRIER 193 DO l = ll_begin, ll_end 194 DO j=jj_begin-offset,jj_end+offset 195 DO i=ii_begin-offset,ii_end+offset 196 ij=(j-1)*iim+i 197 theta(ij,l) = theta_rhodz(ij,l) / ( (ap(l)-ap(l+1)) + ( bp(l)- bp(l+1))* ps(ij) )/g 192 REAL(rstd) :: rhodz 193 INTEGER :: i,j,ij,l 194 195 !$OMP BARRIER 196 DO l = ll_begin, ll_end 197 DO j=jj_begin-offset,jj_end+offset 198 DO i=ii_begin-offset,ii_end+offset 199 ij=(j-1)*iim+i 200 rhodz = ( ap(l)-ap(l+1) + (bp(l)-bp(l+1))*ps(ij) )/g 201 theta(ij,l) = theta_rhodz(ij,l) / rhodz 198 202 ENDDO 199 203 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.