Humidity Boundary condition Atmospheric module

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
SimBu
Posts: 30
Joined: Tue Sep 13, 2022 4:13 pm

Humidity Boundary condition Atmospheric module

Post by SimBu »

Dear all,

While using the humid atmospheric module, I noticed that if the meteo file wasn't provided, it is not possible to define humidity for inlet or inletOutlet :

I tried several things, including neumann, diriclet conditions, but not a single one worked.
Is it mandatory to provide meteo file if we want to define humidity ?

Code: Select all

!< [try_1]
call getfbr('north',nlelt,lstelt)
do ilelt = 1, nlelt
      ifac = lstelt(ilelt)
 ! Found that : " nscaus user-defined scalars and nscapp scalars related to a "specific physics". nscal=nscaus+nscapp "
 
  if ( nscasp.gt.0 ) then    
    do ii = (nscal-nscapp + 1), nscal
      h_spec = 1.d0
      rcodcl(ifac,isca(ii),1) = h_spec
    enddo
  end if 
enddo
 ! [try_1]> 
 
  !< [try_2]
  ! More brutal, I know that it is not a proper way to do it but it doesn't work anyway
call getfbr('north',nlelt,lstelt)
do ilelt = 1, nlelt
      ifac = lstelt(ilelt)
      call field_get_val_s(ihumid, cpro_humid)!FIXME
      h_spec = 1.d0
      cpro_humid(ifac) = h_spec
enddo
 ! [try_2]>
 
I came to the conclusion that it was impossible, am I right ?
So it is not possible to use humid module with idealized meteo profile or standard inlet ?
Post Reply