from_sol_to_temp_mesh2D Subroutine

public subroutine from_sol_to_temp_mesh2D(me, sim)

Arguments

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

Called by

proc~~from_sol_to_temp_mesh2d~~CalledByGraph proc~from_sol_to_temp_mesh2d from_sol_to_temp_mesh2D proc~main_loop main_loop proc~main_loop->proc~from_sol_to_temp_mesh2d program~reims_p reims_p program~reims_p->proc~main_loop

Source Code

subroutine from_sol_to_temp_mesh2D(me,sim)    
    type(mesh2D_t), intent(inout) :: me
    type(simulation_t), intent(in) :: sim

    integer :: nb_e

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

    do nb_e=1,me%M2D_Prop%nb_elements
        me%StVar%temp(nb_e)=me%StVar%temp(nb_e)+me%big%bvx(nb_e) ! update of mesh2D temperature by the solution
        
        ! For step management
        if(.not. sim%explicit) then
            me%err = me%err + ( (1.0_dp/6.0_dp) * (sim%dtPrev1+sim%dt) * ( &
                (me%StVar    %temp(nb_e)-me%StVarOld %temp(nb_e))               / sim%dt - &
                (me%StVarOld %temp(nb_e)-me%StVarOld2%temp(nb_e)) * (sim%dt+sim%dtPrev1) / sim%dtPrev1**2 + &
                (me%StVarOld2%temp(nb_e)-me%StVarOld3%temp(nb_e)) *  sim%dt/sim%dtPrev1  / sim%dtPrev2 ))**2
            me%err_den = me%err_den + (me%StVarOld%temp(nb_e))**2
        endif
    enddo
      
end subroutine from_sol_to_temp_mesh2D