search

eguil / Density_bining

Density bining code
2 stars 5 forks source link

Compute overturning stream function #59

Open eguil opened 7 years ago

eguil commented 7 years ago

bin velocity to compute MSF in density space per basin

mignotj commented 6 years ago

Ongoing work at LOCEAN: add variable transport (h*v) to the general binning procedure. Vertical integral will have to be added separately: online or offline?

eguil commented 6 years ago

To check code after changes, verify that value computed from: print ' Total volume in rho coordinates source grid (ref = 1.33 e+18) : ', voltot print ' Mean Temp./Salinity in rho coordinates source grid : ', temtot, saltot remain the same

Using ciclad:/home/nillod/Dev/Density_bin/Density_bining/bak/tmp/quickTest.py Nicolas gets:

Total volume in z coordinates source grid (ref = 1.33 e+18) :  1.3761059837e+18
  Mean Temp./Salinity in z coordinates source grid            :  2.952495425 34.7336599091
eguil commented 6 years ago

First version coded May 2018. Still issues with z_zt vs. z_zw. May require to use only z_zw in interp mapping and integrate T and S as well as V. We should not have to correct for a half level.

eguil commented 5 years ago

Progress made but still a computing issue as result is wrong

eguil commented 5 years ago

Should look like this: ORLNH_deepMSFs_y2100.pdf

eguil commented 5 years ago

Trying new zonal mean calculation - see issue #71 - but no change