Input_sounding + Sounding_type=ConstantDensity does not work!

[Adiram97]

Hello,

I am trying to run DNS with Boussinesq equations by providing a profile for velocity and temperature in the form of input_sounding. But when I try sounding_type=constantDensity, it gives me an error saying "Arrays are hard-wired upto 256 high only"

I would just like to give a profile. So if there is issue with non-ideal sounding_type, I can even shift to some other format to provide a profile to start the simulation. But I am not sure if there is any other format available to provide a profile

I can even shift to init_type=ConstantDensity instead of Input_sounding. But is there a way to provide a velocity profile if the init_type is constant density?

Thank you very much

[asalmgren]

init_type is used to define the base state -- so I would suggest keeping that as ConstantDensity but then use the files in Source/Prob -- specifically you can use ERF*UserDefined.H -- to set the velocity and theta initial values. Then set erf.prob_name = UserDefined in your inputs file.

[Adiram97]

hi,
Thank you for your response
I am still trying to understand how to define Velocities in ERF_InitCustomPertVels_UserDefined.H. Also because, I need the functionalities of abl problem type. So I am still trying to understand how to define profiles.

I'll update you once I have something