| 5 | |
| 6 | New scaling parameters converting Phantom Units to astrobear CU |
| 7 | {{{ |
| 8 | ! a(:,:,:,5)=a(:,:,:,5)*a(:,:,:,1)*1.907e15*efact / pscale !< rescale internal energy for simulation gamma |
| 9 | a(:,:,:,1)=a(:,:,:,1)*5.901/rscale |
| 10 | a(:,:,:,5)=a(:,:,:,5)*1.907e15*a(:,:,:,1)*efact / pscale !< rescale internal energy for simulation gamma !! moved from top by B.Liu, 1.12.2022 |
| 11 | a(:,:,:,2:4)=a(:,:,:,2:4)*4.367e7/velscale * spread(a(:,:,:,1),4,3) |
| 12 | a(:,:,:,5)=a(:,:,:,5)+.5*sum(a(:,:,:,2:4)**2,4)/a(:,:,:,1) |
| 13 | |
| 14 | ! Calculate overlap |
| 15 | mb(:,2)=min(ibox(:,2), info%mglobal(:,2)) |
| 16 | |
| 17 | }}} |