Commit bc41363f authored by Oliver Henke-Seemann's avatar Oliver Henke-Seemann
Browse files

Fix eines Fehlers mit ini_melt_curves

parent 66abd816
......@@ -598,14 +598,14 @@ program CHIC
! if (pI%C_U.ne.0.0_dp) call initialize_heat(pC,pI,pF,pG,pX,pE%H_from_surf,pM%Lambda,cell) ! needs to be before crust_init
if (pH%use_magn_H.ne.0.0_dp) call initialize_magn_heat(pF,pH,mesh%rc,mesh%nr,pI%prof_nm,0.0_dp,cell)
! if (pM%use_melt.eq.1) then ! always load melt curves since cut_T_solidus otherwise doesn't work
if (pM%use_melt.eq.1) then ! always load melt curves since cut_T_solidus otherwise doesn't work
call init_melt_curves(mesh,mc,pF,pM,pG%dim,cell,pI%X_FeM)
if (pG%dim.eq.2) then
mc%DeltaF(0:mesh%nl+1,1,0:mesh%nr+1) = 0.0_dp
else if (pG%dim.eq.3) then
mc%DeltaF(0:mesh%nl+1,0:mesh%ny+1,0:mesh%nr+1) = 0.0_dp
endif
! endif
endif
if ((pT%linear.ge.7).or.(pT%CutSolidus)) call cut_T_solidus(field,mc,pI,pM)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment