[MITgcm-devel] rho 3-d in common block

[Jean-Michel Campin]

Hi Patrick,

I am ready to check-in the modifications to put the sensity (anomaly)
rho in common block (in DYNVARS.h) that will save 1 computation of rho.

The list of files to modify is:
M inc/DYNVARS.h
M src/calc_phi_hyd.F
M src/diags_rho.F
M src/do_oceanic_phys.F
M src/ini_dynvars.F
M src/ini_pressure.F
+ 2 S/R in pkg/down_slope:
M dwnslp_calc_flow.F
M dwnslp_calc_rho.F

Unfortunatly, it's making a big mess in the adjoint test:
lab_sea & global_ocean.cs32x15 , 
even keeping src/calc_phi_hyd.F unchanged does not help,
but global_ocean.90x40x15 is ok (did not try other ad test).

Let me know when you have time to look at this; then I will probably
check-in those changes, and let you fix the broken dependencies.

Cheers,
Jean-Michel