Force temperature in a given region

Force temperature in a given region

This is an example of cs_user_extra_operations which sets temperature to 20 in a given region starting at t = 12s

Local variables to be added

integer          iel
integer          iscal
double precision, dimension(:), pointer :: cvar_scal

Body

Do this with precaution... The user is responsible for the validity of results.

Here is the corresponding code:

iscal = iscalt
 
if (ttcabs .ge. 12.d0) then
 
  if (iscal.gt.0 .and. iscal.le.nscal) then
    call field_get_val_s(ivarfl(isca(iscal)), cvar_scal)
    do iel = 1, ncel
      cvar_scal(iel) = 20.d0
    enddo
  endif
 
  write(nfecra,3000)
 
endif