from_sol_to_temp_solid Subroutine

public subroutine from_sol_to_temp_solid(me, sim)

Arguments

Type IntentOptional Attributes Name
type(solid_t), intent(inout) :: me
type(simulation_t), intent(in) :: sim

Called by

proc~~from_sol_to_temp_solid~~CalledByGraph proc~from_sol_to_temp_solid from_sol_to_temp_solid proc~main_loop main_loop proc~main_loop->proc~from_sol_to_temp_solid program~reims_p reims_p program~reims_p->proc~main_loop

Source Code

subroutine from_sol_to_temp_solid(me,sim)    
    type(solid_t), intent(inout) :: me
    type(simulation_t), intent(in) :: sim

    integer :: jj

    me%err = 0_dp
    me%err_den = 0_dp

    do jj=1,me%MC_Prop%NbCells
        me%StVar%MCtemp(jj)=me%StVar%MCtemp(jj)+me%big%bvS(jj) ! update of solid temperature by the solution
        
        ! Write to HDF5
        me%hdf%data(jj,1) = me%StVar%MCtemp(jj)

        ! For step management
        if(.not. sim%explicit) then
            me%err = me%err + ( (1.0_dp/6.0_dp) * (sim%dtPrev1+sim%dt) * ( &
                (me%StVar%MCtemp(jj)-me%StVarOld%MCtemp(jj))               / sim%dt - &
                (me%StVarOld%MCtemp(jj)-me%StVarOld2%MCtemp(jj))  * (sim%dt+sim%dtPrev1) / sim%dtPrev1**2 + &
                (me%StVarOld2%MCtemp(jj)-me%StVarOld3%MCtemp(jj)) *  sim%dt/sim%dtPrev1  / sim%dtPrev2 ))**2
            me%err_den = me%err_den + (me%StVarOld%MCtemp(jj))**2
        endif
    enddo
      
end subroutine from_sol_to_temp_solid