Boussinesq coefficients for Buoyancy Calculations

[johnsonb]

Hello. I am trying to complete buoyancy calculations on POP output and am trying to figure out how to account for the thermal and haline Boussinesq coefficients. On page 95, the POP reference manual for CESM 1.0 (my integration is using CESM 1.2) mentions the expansion coefficients when discussing the density ratio Rρ = αT∂zΘ/(βS∂zS) and provides a reference to LMD (Large, McWilliams & Doney, 1994).

Page 389 of LMD mentions that α and β are evaluated at local values of T and S.

My question is: If they are evaluated locally, is there a subroutine somewhere in the POP code that I can look at to complete my calculation?

Thank you for your help!

 

[klindsay]

Hi,
alpha and beta are indeed computed locally.They are computed in subroutine state of state_mod.F90, called in subroutine bldepth in vmix_kpp.F90.alpha is state's argument DRHODT and beta is state's argument DRHODS.Keith 

 

[klindsay]

Hi,
alpha and beta are indeed computed locally.They are computed in subroutine state of state_mod.F90, called in subroutine bldepth in vmix_kpp.F90.alpha is state's argument DRHODT and beta is state's argument DRHODS.Keith 

 

[klindsay]

Hi,
alpha and beta are indeed computed locally.They are computed in subroutine state of state_mod.F90, called in subroutine bldepth in vmix_kpp.F90.alpha is state's argument DRHODT and beta is state's argument DRHODS.Keith 

 

[klindsay]

Hi,
alpha and beta are indeed computed locally.They are computed in subroutine state of state_mod.F90, called in subroutine bldepth in vmix_kpp.F90.alpha is state's argument DRHODT and beta is state's argument DRHODS.Keith 

 

[klindsay]

Hi,
alpha and beta are indeed computed locally.They are computed in subroutine state of state_mod.F90, called in subroutine bldepth in vmix_kpp.F90.alpha is state's argument DRHODT and beta is state's argument DRHODS.Keith 

 

[johnsonb]

Hi Keith,

Thank you for your help! I was able to find those files in pop2/source. Thank you also for also highlighting the fact that state is called by vmix_kpp.F90. It is very helpful to see when the subroutine is called.

Hope you're well,

 

[johnsonb]

Hi Keith,

Thank you for your help! I was able to find those files in pop2/source. Thank you also for also highlighting the fact that state is called by vmix_kpp.F90. It is very helpful to see when the subroutine is called.

Hope you're well,

 

[johnsonb]

Hi Keith,

Thank you for your help! I was able to find those files in pop2/source. Thank you also for also highlighting the fact that state is called by vmix_kpp.F90. It is very helpful to see when the subroutine is called.

Hope you're well,

 

[johnsonb]

Hi Keith,

Thank you for your help! I was able to find those files in pop2/source. Thank you also for also highlighting the fact that state is called by vmix_kpp.F90. It is very helpful to see when the subroutine is called.

Hope you're well,

 

[johnsonb]

Hi Keith,

Thank you for your help! I was able to find those files in pop2/source. Thank you also for also highlighting the fact that state is called by vmix_kpp.F90. It is very helpful to see when the subroutine is called.

Hope you're well,