Hello, Are there any routines that could properly set geostrophically balanced initial conditions ( from a given density field calculate U,V,Eta ) ? I made it myself but initial conditions are not in exact balance, and I'm getting ageostrophic flow. Thanks in advance for help, Georgy.