| |
| | [No title] (Site not responding. Last check: 2007-11-02) |
 | | !=============================================================================== if (firstCall) then MyPE = dBasePropertyInteger('MyProcessor') MasterPE = dBasePropertyInteger('MasterProcessor') idens = dBaseKeyNumber('dens') iener = dBaseKeyNumber('ener') ivelx = dBaseKeyNumber('velx') ively = dBaseKeyNumber('vely') ivelz = dBaseKeyNumber('velz') ipres = dBaseKeyNumber('pres') igamc = dBaseKeyNumber('gamc') igame = dBaseKeyNumber('game') inuc_begin = dBaseSpecies(1) iPt = dBaseKeyNumber('Point') call get_parm_from_context(GLOBAL_PARM_CONTEXT, 'gamma', gamma) call get_parm_from_context(GLOBAL_PARM_CONTEXT, 'smallp', smallp) call get_parm_from_context(GLOBAL_PARM_CONTEXT, 'smallx', smallx) ! |
 | | Compute the gas energy and set the gamma-values needed for ! |
 | | ekin_zone = 0.5 * (velx_zone**2 + vely_zone**2 + velz_zone**2) ener_zone = pres_zone / (gamma-1.) ener_zone = ener_zone / rho_zone ener_zone = ener_zone + ekin_zone ener_zone = max(ener_zone, smallp) call dBasePutData(inuc_begin,iPt, i,j,k, block_no, & & 1.0-(ionmax-1)*smallx, 1) do n = 2, ionmax call dBasePutData(inuc_begin-1+n,iPt, i,j,k, block_no, smallx, 1) enddo ! |
| flash.uchicago.edu /maillists_public/flash-users/att-0682/02-init_block.F90 (283 words) |
|