[Mitgcm-support] help setting up tank experiement

[mitgcm-support at dev.mitgcm.org]

jean-michel,

 could you spare a minute to look a problem sai is having.
 story so far is this:

 1. we are setting up a tank experiment 60cmx60cmx15cm
 2. domain is a cylindrical tank with a coffee can of ice in
    the middle.
 3. coordinate system is cartesian.
 4. we setup a depth mask "Bathy.bin"
 5. starting from global_ocean.90x40x15
    we edit SIZE.h etc... for a dx=0.005 and dz=0.01 i.e
    Nx=Ny=124,Nr=17 (couple of layers of land at the bottom!!),
    visc, kapp and deltat.
 6. we initialise theta to 20C +/- 0.001*randn().
 7. we run the model with very small deltat.
 8. it blows up.
 9. cfl for kappa, visc etc... seem ok.
10. there is a weird initial values in Ro_surf of 2.7 x 10^-17
11. maybe there is some threshold in the absolute hfac or something.
12. Ro_surf should be zero everywhere - it is in global_ocean.

 i'm kind of busy today, otherwise i could take a look.

 thanks,

chris