mocsy  2.0
Fortran 95 routines to model ocean carbonate system thermodynamics
 All Classes Namespaces Files Functions Variables Pages
Variables
mocsy_GLODAP Namespace Reference

Variables

string usage
 
tuple nb_arguments = len(sys.argv)
 
tuple preindustrial_option = int(sys.argv[1])
 
string OCMIP5_DAT = 'http://dods.ipsl.jussieu.fr/cgi-bin/nph-dods/ocmip/phase5/DATA'
 
string GLODAPsubdir = 'gridded/glodap'
 
string WOA2009subdir = 'gridded/WOA2009'
 
 DIC_dsubdir = GLODAPsubdir
 
string DIC_dfile = 'TCO2.nc'
 
string DIC_dvar = 'TCO2'
 
 Alk_dsubdir = GLODAPsubdir
 
string Alk_dfile = 'TALK.nc'
 
string Alk_dvar = 'TALK'
 
 Temp_dsubdir = WOA2009subdir
 
string Temp_dfile = "temperature_annual_1deg.nc"
 
string Temp_dvar = 't_an'
 
 Salt_dsubdir = WOA2009subdir
 
string Salt_dfile = "salinity_annual_1deg.nc"
 
string Salt_dvar = 's_an'
 
 PO4_dsubdir = WOA2009subdir
 
string PO4_dfile = "phosphate_annual_1deg.nc"
 
string PO4_dvar = 'p_an'
 
 Si_dsubdir = WOA2009subdir
 
string Si_dfile = "silicate_annual_1deg.nc"
 
string Si_dvar = 'i_an'
 
string DIC_datfile = OCMIP5_DAT+"/"
 
string Alk_datfile = OCMIP5_DAT+"/"
 
string Temp_datfile = OCMIP5_DAT+"/"
 
string Salt_datfile = OCMIP5_DAT+"/"
 
string PO4_datfile = OCMIP5_DAT+"/"
 
string Si_datfile = OCMIP5_DAT+"/"
 
tuple j = cdms.open(DIC_datfile)
 
tuple dicdat = j(DIC_dvar)
 
tuple depth1D = dicdat.getLevel()
 
tuple lat2D = dicdat.getLatitude()
 
tuple f = cdms.open(DIC_datfile)
 
tuple dicd = f(DIC_dvar)
 
tuple dic_axes = dicd.getAxisList()
 
tuple dic_grid = dicd.getGrid()
 
string aDIC_dsubdir = 'gridded/glodap'
 
string aDIC_dfile = 'ANTCO2.nc'
 
string aDIC_dvar = 'ANTCO2'
 
string aDIC_datfile = OCMIP5_DAT+"/"
 
tuple h = cdms.open(aDIC_datfile)
 
tuple anth_dic = h(aDIC_dvar)
 
tuple alkd = f(Alk_dvar)
 
tuple td = f(Temp_dvar)
 
tuple sd = f(Salt_dvar)
 
tuple po4d = f(PO4_dvar)
 
tuple sio2d = f(Si_dvar)
 
 depth = depth1D
 
 lat = lat2D
 
tuple common_mask = np.zeros(dicd.shape)
 
tuple v = cdms.MV2.reshape(v, dicd.shape)
 
 m = v.mask
 
tuple pfd = mocsy.mdepth2press(depth,lat)
 
tuple rhois = mocsy.mrhoinsitu(sd, td, pfd)
 
list patmd = td[:]
 
string optcon = 'mol/kg'
 
tuple pH = np.reshape(pH, dicd.shape)
 
tuple pco2 = np.reshape(pco2, dicd.shape)
 
tuple fco2 = np.reshape(fco2, dicd.shape)
 
tuple co2 = np.reshape(co2, dicd.shape)
 
tuple hco3 = np.reshape(hco3, dicd.shape)
 
tuple co3 = np.reshape(co3, dicd.shape)
 
tuple OmegaA = np.reshape(OmegaA, dicd.shape)
 
tuple OmegaC = np.reshape(OmegaC, dicd.shape)
 
tuple BetaD = np.reshape(BetaD, dicd.shape)
 
tuple Alk = np.reshape(alkd, dicd.shape)
 
tuple DIC = np.reshape(dicd.ravel(), dicd.shape)
 
tuple DENis = np.reshape(DENis, dicd.shape)
 
tuple p = np.reshape(p, dicd.shape)
 
tuple Tis = np.reshape(Tis, dicd.shape)
 
string dirnameo = '/prodigfs/OCMIP5/DATA/gridded/glodap/derived2'
 
tuple version = os.popen("date +'v20%y%m%d'")
 
string filenameo = "GLODAP_pH_Carb_1765_3Da_"
 
string dirfilenameo = dirnameo+"/"
 
tuple dirfileo = cdms.open(dirfilenameo, "w")
 

Variable Documentation

string mocsy_GLODAP.aDIC_datfile = OCMIP5_DAT+"/"

Definition at line 164 of file mocsy_GLODAP.py.

string mocsy_GLODAP.aDIC_dfile = 'ANTCO2.nc'

Definition at line 162 of file mocsy_GLODAP.py.

string mocsy_GLODAP.aDIC_dsubdir = 'gridded/glodap'

Definition at line 161 of file mocsy_GLODAP.py.

string mocsy_GLODAP.aDIC_dvar = 'ANTCO2'

Definition at line 163 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.Alk = np.reshape(alkd, dicd.shape)

Definition at line 298 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Alk_datfile = OCMIP5_DAT+"/"

Definition at line 141 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Alk_dfile = 'TALK.nc'

Definition at line 120 of file mocsy_GLODAP.py.

mocsy_GLODAP.Alk_dsubdir = GLODAPsubdir

Definition at line 119 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Alk_dvar = 'TALK'

Definition at line 121 of file mocsy_GLODAP.py.

float mocsy_GLODAP.alkd = f(Alk_dvar)

Definition at line 176 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.anth_dic = h(aDIC_dvar)

Definition at line 166 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.BetaD = np.reshape(BetaD, dicd.shape)

Definition at line 297 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.co2 = np.reshape(co2, dicd.shape)

Definition at line 292 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.co3 = np.reshape(co3, dicd.shape)

Definition at line 294 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.common_mask = np.zeros(dicd.shape)

Definition at line 236 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.DENis = np.reshape(DENis, dicd.shape)

Definition at line 301 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.depth = depth1D

Definition at line 208 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.depth1D = dicdat.getLevel()

Definition at line 151 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.DIC = np.reshape(dicd.ravel(), dicd.shape)

Definition at line 299 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.dic_axes = dicd.getAxisList()

Definition at line 156 of file mocsy_GLODAP.py.

string mocsy_GLODAP.DIC_datfile = OCMIP5_DAT+"/"

Definition at line 140 of file mocsy_GLODAP.py.

string mocsy_GLODAP.DIC_dfile = 'TCO2.nc'

Definition at line 116 of file mocsy_GLODAP.py.

mocsy_GLODAP.DIC_dsubdir = GLODAPsubdir

Definition at line 115 of file mocsy_GLODAP.py.

string mocsy_GLODAP.DIC_dvar = 'TCO2'

Definition at line 117 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.dic_grid = dicd.getGrid()

Definition at line 157 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.dicd = f(DIC_dvar)

Definition at line 155 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.dicdat = j(DIC_dvar)

Definition at line 150 of file mocsy_GLODAP.py.

string mocsy_GLODAP.dirfilenameo = dirnameo+"/"

Definition at line 416 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.dirfileo = cdms.open(dirfilenameo, "w")

Definition at line 417 of file mocsy_GLODAP.py.

string mocsy_GLODAP.dirnameo = '/prodigfs/OCMIP5/DATA/gridded/glodap/derived2'

Definition at line 407 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.f = cdms.open(DIC_datfile)

Definition at line 154 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.fco2 = np.reshape(fco2, dicd.shape)

Definition at line 291 of file mocsy_GLODAP.py.

string mocsy_GLODAP.filenameo = "GLODAP_pH_Carb_1765_3Da_"

Definition at line 412 of file mocsy_GLODAP.py.

string mocsy_GLODAP.GLODAPsubdir = 'gridded/glodap'

Definition at line 112 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.h = cdms.open(aDIC_datfile)

Definition at line 165 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.hco3 = np.reshape(hco3, dicd.shape)

Definition at line 293 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.j = cdms.open(DIC_datfile)

Definition at line 149 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.lat = lat2D

Definition at line 209 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.lat2D = dicdat.getLatitude()

Definition at line 152 of file mocsy_GLODAP.py.

mocsy_GLODAP.m = v.mask

Definition at line 239 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.nb_arguments = len(sys.argv)

Definition at line 91 of file mocsy_GLODAP.py.

string mocsy_GLODAP.OCMIP5_DAT = 'http://dods.ipsl.jussieu.fr/cgi-bin/nph-dods/ocmip/phase5/DATA'

Definition at line 110 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.OmegaA = np.reshape(OmegaA, dicd.shape)

Definition at line 295 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.OmegaC = np.reshape(OmegaC, dicd.shape)

Definition at line 296 of file mocsy_GLODAP.py.

string mocsy_GLODAP.optcon = 'mol/kg'

Definition at line 275 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.p = np.reshape(p, dicd.shape)

Definition at line 302 of file mocsy_GLODAP.py.

list mocsy_GLODAP.patmd = td[:]

Definition at line 264 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.pco2 = np.reshape(pco2, dicd.shape)

Definition at line 290 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.pfd = mocsy.mdepth2press(depth,lat)

Definition at line 258 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.pH = np.reshape(pH, dicd.shape)

Definition at line 289 of file mocsy_GLODAP.py.

string mocsy_GLODAP.PO4_datfile = OCMIP5_DAT+"/"

Definition at line 144 of file mocsy_GLODAP.py.

string mocsy_GLODAP.PO4_dfile = "phosphate_annual_1deg.nc"

Definition at line 133 of file mocsy_GLODAP.py.

mocsy_GLODAP.PO4_dsubdir = WOA2009subdir

Definition at line 132 of file mocsy_GLODAP.py.

string mocsy_GLODAP.PO4_dvar = 'p_an'

Definition at line 134 of file mocsy_GLODAP.py.

list mocsy_GLODAP.po4d = f(PO4_dvar)

Definition at line 194 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.preindustrial_option = int(sys.argv[1])

Definition at line 95 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.rhois = mocsy.mrhoinsitu(sd, td, pfd)

Definition at line 259 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Salt_datfile = OCMIP5_DAT+"/"

Definition at line 143 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Salt_dfile = "salinity_annual_1deg.nc"

Definition at line 129 of file mocsy_GLODAP.py.

mocsy_GLODAP.Salt_dsubdir = WOA2009subdir

Definition at line 128 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Salt_dvar = 's_an'

Definition at line 130 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.sd = f(Salt_dvar)

Definition at line 188 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Si_datfile = OCMIP5_DAT+"/"

Definition at line 145 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Si_dfile = "silicate_annual_1deg.nc"

Definition at line 137 of file mocsy_GLODAP.py.

mocsy_GLODAP.Si_dsubdir = WOA2009subdir

Definition at line 136 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Si_dvar = 'i_an'

Definition at line 138 of file mocsy_GLODAP.py.

list mocsy_GLODAP.sio2d = f(Si_dvar)

Definition at line 200 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.td = f(Temp_dvar)

Definition at line 182 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Temp_datfile = OCMIP5_DAT+"/"

Definition at line 142 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Temp_dfile = "temperature_annual_1deg.nc"

Definition at line 125 of file mocsy_GLODAP.py.

mocsy_GLODAP.Temp_dsubdir = WOA2009subdir

Definition at line 124 of file mocsy_GLODAP.py.

string mocsy_GLODAP.Temp_dvar = 't_an'

Definition at line 126 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.Tis = np.reshape(Tis, dicd.shape)

Definition at line 303 of file mocsy_GLODAP.py.

string mocsy_GLODAP.usage
Initial value:
1 = """
2  Usage: python mocsy_GLODAP.py [preindustrial_option]
3 
4  preindus-option: Option for computing from preindustrial DIC or total DIC
5  0 = total (1994), i.e., modern era, this is the default
6  1 = preindustrial (1765), after removing anthropogenic DIC from total DIC
7 
8  From observed DIC & Alk (GLODAP), T, S, PO4, and SiO2 (WOA2009)
9  compute following 3D variables:
10  - pH,
11  - CO2*, HCO3- and CO3-- concentrations,
12  - OmegaA & OmegaC (saturation states for aragonite and calcite),
13  - Revelle factor (BetaD or homogeneous buffer factor)
14  Save them in appropriate netCDF files
15  preindustrial_option = 0
16  """

Definition at line 74 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.v = cdms.MV2.reshape(v, dicd.shape)

Definition at line 238 of file mocsy_GLODAP.py.

tuple mocsy_GLODAP.version = os.popen("date +'v20%y%m%d'")

Definition at line 410 of file mocsy_GLODAP.py.

string mocsy_GLODAP.WOA2009subdir = 'gridded/WOA2009'

Definition at line 113 of file mocsy_GLODAP.py.