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]>
So it is not possible to use humid module with idealized meteo profile or standard inlet ?