| | 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 | }}} |