[MITgcm-support] Uniform diffkrfile

[Dimitris Menemenlis]

Geof, this is beyond my area of expertise.  Your results do not make sense.
The value of diffKzS does not appear to impact S in any way, only T.
Occasionally there have been packages, e.g., KPP that assume that 
diffKzT=diffKzS.  What is content of your data.pkg file?

I don't know what else quick to suggest.


If you know how to use pkg/diagnostics, I would try dumping
>   93 |ADVr_TH | 23 |WM      LR      |degC.m^3/s      |Vertical   Advective Flux of Pot.Temperature
>   96 |DFrE_TH | 23 |WM      LR      |degC.m^3/s      |Vertical Diffusive Flux of Pot.Temperature (Explicit part)
>   96 |DFrE_TH | 23 |WM      LR      |degC.m^3/s      |Vertical Diffusive Flux of Pot.Temperature (Explicit part)
>  100 |ADVr_SLT| 23 |WM      LR      |psu.m^3/s       |Vertical   Advective Flux of Salinity
>  103 |DFrE_SLT| 23 |WM      LR      |psu.m^3/s       |Vertical Diffusive Flux of Salinity    (Explicit part)
>  106 |DFrI_SLT| 23 |WM      LR      |psu.m^3/s       |Vertical Diffusive Flux of Salinity    (Implicit part)

Another thing one could do is to turn off some of the model terms, e.g.,

> C     momAdvection  :: Flag which turns advection of momentum on and off.                                                                     
> C     momForcing    :: Flag which turns external forcing of momentum on                                                                       
> C                     and off.                                                                                                                
> C     momPressureForcing :: Flag which turns pressure term in momentum equation                                                               
> C                          on and off.                                                                                                        
> C     metricTerms   :: Flag which turns metric terms on or off.                                                                               
> C     useNHMTerms   :: If TRUE use non-hydrostatic metric terms.                                                                              
> C     useCoriolis   :: Flag which turns the coriolis terms on and off.                                                                        
> C     tempAdvection :: Flag which turns advection of temperature on                                                                           
> C                     and off.                                                                                                                
> C     tempForcing   :: Flag which turns external forcing of temperature on                                                                    
> C                     and off.                                                                                                                
> C     saltAdvection :: Flag which turns advection of salinity on                                                                              
> C                     and off.                                                                                                                
> C     saltForcing   :: Flag which turns external forcing of salinity on                                                                       
> C                     and off.