[MITgcm-devel] LAYERS_THERMODYNAMICS

[Hajoon Song]

Dear all,

I am trying to use the LAYERS_THERMODYNAMICS option in the layers package.
I think some of the units in "layers_diagnostics_init.F" may be incorrect.

For example, in case "layers_num=3", the line 211says that "LTha3RHO" has a unit of 'm deg./s'.
This diagnostic variable contains "layers_TtendAdvh" in "layers_calc.F", and 
"TtendAdvh" in "layers_fluxcalc.F" seems to have a unit of 'm/s'. 

I could be wrong, so please correct me if so.

Best,
Hajoon