Calculating density

The following functions will be useful:

gsw.p_from_z(z, lat) (i.e. if needed)
gsw.SA_from_SP(SP, p, lon, lat)
gsw.CT_from_pt(SA, pt)
gsw.rho(SA, CT, p)

SP = practical salinity, SA = absolute salinity, p = pressure (dbar), CT = conservative temperature, pt = potential temperature.



In [1]:

    
import numpy
import gsw

Broadcasting and gsw



In [6]:

    
z = [-10., -50., -125., -250., -600., -1000.]
lat = 4.
gsw.p_from_z(z, lat)









    Out[6]:





array([   10.05572704,    50.28354425,   125.73185732,   251.54028663,
         604.2099135 ,  1007.9900587 ])



In [7]:

    
z = [-10., -50., -125., -250., -600., -1000.]
lat = [4., 4., 4., 4., 4., 4.]
gsw.p_from_z(z, lat)









    Out[7]:





array([   10.05572704,    50.28354425,   125.73185732,   251.54028663,
         604.2099135 ,  1007.9900587 ])



In [11]:

    
z = numpy.ones((3,4))
lat = numpy.ones((3,4))
gsw.p_from_z(z, lat)









    Out[11]:





array([[-1.00552148, -1.00552148, -1.00552148, -1.00552148],
       [-1.00552148, -1.00552148, -1.00552148, -1.00552148],
       [-1.00552148, -1.00552148, -1.00552148, -1.00552148]])



In [ ]: