- Timestamp:
- 09/05/08 11:49:42 (16 years ago)
- Branches:
- master
- Children:
- fafe4d6
- Parents:
- 2ebb8b6
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
bgs.f
r2ebb8b6 rbd2278d 156 156 bv(nph,3)=xiv(nph,3)-zat(iN(icurraa)) 157 157 ab(nph)=bv(nph,1)*bv(nph,1)+bv(nph,2)*bv(nph,2) 158 #+bv(nph,3)*bv(nph,3)158 & +bv(nph,3)*bv(nph,3) 159 159 iph(nph)=iphi(icurraa) 160 160 endif … … 168 168 bv(nph,3)=xiv(nph,3)-zat(iCa(icurraa)) 169 169 ab(nph)=bv(nph,1)*bv(nph,1)+bv(nph,2)*bv(nph,2) 170 #+bv(nph,3)*bv(nph,3)170 & +bv(nph,3)*bv(nph,3) 171 171 iph(nph)=ipsi(icurraa) 172 172 endif … … 185 185 do j=1,nph 186 186 dv(i,j,1)=(1.0/ab(j))*(bv(j,2)*(rv(i,3)-xiv(j,3))- 187 cbv(j,3)*(rv(i,2)-xiv(j,2)))187 & bv(j,3)*(rv(i,2)-xiv(j,2))) 188 188 dv(i,j,2)=(-1.0/ab(j))*(bv(j,1)*(rv(i,3)-xiv(j,3))- 189 cbv(j,3)*(rv(i,1)-xiv(j,1)))189 & bv(j,3)*(rv(i,1)-xiv(j,1))) 190 190 dv(i,j,3)=(1.0/ab(j))*(bv(j,1)*(rv(i,2)-xiv(j,2))- 191 cbv(j,2)*(rv(i,1)-xiv(j,1)))191 & bv(j,2)*(rv(i,1)-xiv(j,1))) 192 192 enddo 193 193 enddo … … 273 273 bv(nph,3)=xiv(nph,3)-zat(iN(icurraa)) 274 274 ab(nph)=bv(nph,1)*bv(nph,1)+bv(nph,2)*bv(nph,2) 275 #+bv(nph,3)*bv(nph,3)275 & +bv(nph,3)*bv(nph,3) 276 276 iph(nph)=iphi(icurraa) 277 277 endif … … 285 285 bv(nph,3)=xiv(nph,3)-zat(iCa(icurraa)) 286 286 ab(nph)=bv(nph,1)*bv(nph,1)+bv(nph,2)*bv(nph,2) 287 #+bv(nph,3)*bv(nph,3)287 & +bv(nph,3)*bv(nph,3) 288 288 iph(nph)=ipsi(icurraa) 289 289 endif … … 302 302 do j=1,nph 303 303 dv(i,j,1)=(1.0/ab(j))*(bv(j,2)*(rv(i,3)-xiv(j,3))- 304 cbv(j,3)*(rv(i,2)-xiv(j,2)))304 & bv(j,3)*(rv(i,2)-xiv(j,2))) 305 305 dv(i,j,2)=(-1.0/ab(j))*(bv(j,1)*(rv(i,3)-xiv(j,3))- 306 cbv(j,3)*(rv(i,1)-xiv(j,1)))306 & bv(j,3)*(rv(i,1)-xiv(j,1))) 307 307 dv(i,j,3)=(1.0/ab(j))*(bv(j,1)*(rv(i,2)-xiv(j,2))- 308 cbv(j,2)*(rv(i,1)-xiv(j,1)))308 & bv(j,2)*(rv(i,1)-xiv(j,1))) 309 309 enddo 310 310 enddo
Note:
See TracChangeset
for help on using the changeset viewer.