Simulation_initBlok.F90

  if (NDIM == 1 .AND. meshGeom == SPHERICAL) then

     call Grid_getCellCoords(IAXIS,blockID, CENTER, .true.,xCenter,iSizeGC)
     call Grid_getCellCoords(IAXIS,blockID,LEFT_EDGE, .true.,xLeft, iSizeGC)
     call Grid_getCellCoords(IAXIS,blockID,RIGHT_EDGE,.true.,xRight,iSizeGC)

     do k = blkLimits(LOW,KAXIS),blkLimits(HIGH,KAXIS)
        do j = blkLimits(LOW,JAXIS),blkLimits(HIGH,JAXIS)
           do i = blkLimits(LOW,IAXIS),blkLimits(HIGH,IAXIS)

              ! the interpolation will be done using the parabolic interpolation
              ! routine
              if (xCenter(i) <= xzn(n1d_total)) then
                 do ivar = 1, NUNK_VARS
                    call parabolic_interp(xzn, model_1d(:,ivar), &
                         n1d_total, xLeft(i), xCenter(i), xRight(i), var_interp)
                    solnData(ivar,i,j,k) = var_interp
                 end do
              else
                 call Driver_abortFlash("Beyond the 1D model data")
              endif
              ! multiply initial radial velocity, if desired
              solnData(VELX_VAR,i,j,k) = sim_velMult*solnData(VELX_VAR,i,j,k)

              sumY = 1.0

#ifdef FLASH_MULTISPECIES
              sum = 0.e0
              do n = SPECIES_BEGIN,SPECIES_END
                 solnData(n,i,j,k) = &
                      max(sim_smallx, &
                      min(1.e0,solnData(n,i,j,k)))
                 sum = sum + solnData(n,i,j,k)
              enddo
              suminv = 1.e0 / sum
              do n = SPECIES_BEGIN, SPECIES_END
                 solnData(n,i,j,k) =  &
                      max(sim_smallx, min(1.e0,suminv*&
                      solnData(n,i,j,k)))
              enddo
              call Multispecies_getSumInv(A,sumY,solnData(SPECIES_BEGIN:SPECIES_END,i,j,k))
#endif
#ifdef SUMY_MSCALAR
              solnData(SUMY_MSCALAR,i,j,k) = sumY
#endif
           axis(IAXIS) = i
           axis(JAXIS) = j
           axis(KAXIS) = k

           call RadTrans_mgdEFromT(blockId, axis, 11604.55, tradActual)
           call Grid_putPointData(blockId, CENTER, TRAD_VAR, EXTERIOR, axis,tradActual)

           enddo
        enddo
     enddo

  endif