! ! Define UMaine CoSiNE ecosystem parameters. ! Vinfo( 1)='BioIter' Vinfo( 2)='number of iterations to achieve convergence' status=def_var(ng, model, ncid, varid, nf90_int, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='reg1' Vinfo( 2)='microzooplankton excretion rate to ammonium' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='reg2' Vinfo( 2)='mesozooplankton excretion rate to ammonium' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='gmaxs1' Vinfo( 2)='maximum specific growth rate of small phytoplankton' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='gmaxs2' Vinfo( 2)='maximum specific growth rate of diatom' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='gmaxs3' Vinfo( 2)='maximum specific growth rate of coccolithophores' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='beta1' Vinfo( 2)='microzooplankton maximum grazing rate' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='beta2' Vinfo( 2)='mesozooplankton maximum grazing rate' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akz1' Vinfo( 2)='half saturation for microzooplankton grazing' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akz2' Vinfo( 2)='half saturation for mesozooplankton grazing' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='PARfrac' Vinfo( 2)='photosynthetically available radiation fraction' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='alphachl_s1' Vinfo( 2)='slope chl-specific P-I curve of small phytoplankton' Vinfo( 3)='meter2 watts-1 day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='alphachl_s2' Vinfo( 2)='slope chl-specific P-I curve of diatom' Vinfo( 3)='meter2 watts-1 day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='alphachl_s3' Vinfo( 2)='slope chl-specific P-I curve of coccolithophores' Vinfo( 3)='meter2 watts-1 day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='pis1' Vinfo( 2)='ammonium inhibition parameter for small phytoplankton' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='pis2' Vinfo( 2)='ammonium inhibition parameter for diatom' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='pis3' Vinfo( 2)='ammonium inhibition parameter for coccolithophores' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akno3s1' Vinfo( 2)='half saturation of nitrate uptake by iSphy' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akno3s2' Vinfo( 2)='half saturation of nitrate uptake by diatom' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akno3s3' Vinfo( 2)='half saturation of nitrate uptake by coccolithophores' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='aknh4s1' Vinfo( 2)='half saturation of ammonium uptake by iSphy' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='aknh4s2' Vinfo( 2)='half saturation of ammonium uptake by iLphy' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='aknh4s3' Vinfo( 2)='half saturation of ammonium uptake by coccolithophores' Vinfo( 3)='millimole_N meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akpo4s1' Vinfo( 2)='half saturation of phosphate uptake by iSphy' Vinfo( 3)='millimole_P meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akpo4s2' Vinfo( 2)='half saturation of phosphate uptake by iLphy' Vinfo( 3)='millimole_P meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akpo4s3' Vinfo( 2)='half saturation of phosphate uptake by coccolithophores' Vinfo( 3)='millimole_P meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akco2s1' Vinfo( 2)='half saturation of co2 uptake by iSphy' Vinfo( 3)='millimole_C meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akco2s2' Vinfo( 2)='half saturation of co2 uptake by iLphy' Vinfo( 3)='millimole_C meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='akco2s3' Vinfo( 2)='half saturation of co2 uptake by coccolithophores' Vinfo( 3)='millimole_C meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='aksio4s2' Vinfo( 2)='half saturation of silicate uptake by iLphy' Vinfo( 3)='millimole_Si meter-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ES1' Vinfo( 2)='Phytoplankton exudation parameter for small phytoplankton' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ES2' Vinfo( 2)='Phytoplankton exudation parameter for diatom' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ES3' Vinfo( 2)='Phytoplankton exudation parameter for coccolithophores' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ak1' Vinfo( 2)='light attenuation coefficient of water' Vinfo( 3)='meter-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ak2' Vinfo( 2)='phytoplankton specific light attenuation rate' Vinfo( 3)='meter2 millimole_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='Qmax' Vinfo( 2)='Maximum phytoplankton N:C ratio' Vinfo( 3)='mol_N mol_C-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='Qmin' Vinfo( 2)='Minimum phytoplankton N:C ratio' Vinfo( 3)='mol_N mol_C-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='lambdano3_s1' Vinfo( 2)='Cost of biosynthesis for small phytoplankton' Vinfo( 3)='mol_C mol_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='lambdano3_s2' Vinfo( 2)='Cost of biosynthesis for diatom' Vinfo( 3)='mol_C mol_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='lambdano3_s3' Vinfo( 2)='Cost of biosynthesis for coccolithophores' Vinfo( 3)='mol_C mol_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='thetaNmax_s1' Vinfo( 2)='Maximum Chl:N for small phytoplankton' Vinfo( 3)='g_Chl mol_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='thetaNmax_s2' Vinfo( 2)='Maximum Chl:N for small phytoplankton' Vinfo( 3)='g_Chl mol_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='thetaNmax_s3' Vinfo( 2)='Maximum Chl:N for small phytoplankton' Vinfo( 3)='g_Chl mol_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma' Vinfo( 2)='mesozooplankton specific mortality rate' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma1' Vinfo( 2)='grazing efficiency of microzooplankton' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma2' Vinfo( 2)='grazing efficiency of mesozooplankton for N' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma22' Vinfo( 2)='grazing efficiency of mesozooplankton for C' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma3' Vinfo( 2)='death rate of small phytoplankton' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma4' Vinfo( 2)='death rate of large phytoplankton' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma10' Vinfo( 2)='death rate of coccolithophores' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma12' Vinfo( 2)='death rate of bacteria' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma5' Vinfo( 2)='decay rate of detritus' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma7' Vinfo( 2)='nitrafication rate' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma11' Vinfo( 2)='Maximum ammonium uptake rate by bacteria' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='bgamma13' Vinfo( 2)='Maximum semi-labile hydrolysis' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='mtos1' Vinfo( 2)='Mortality to dissolved pool of small phytoplankton' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='mtos2' Vinfo( 2)='Mortality to dissolved pool of diatom' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='mtos3' Vinfo( 2)='Mortality to dissolved pool of coccolithophores' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='flz1' Vinfo( 2)='Feeding loss by small zooplankton' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='flz2' Vinfo( 2)='Feeding loss by large zooplankton' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='lk1' Vinfo( 2)='Phytoplankton leakage fraction of small phytoplankton' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='lk2' Vinfo( 2)='Phytoplankton leakage fraction of diatom' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='lk3' Vinfo( 2)='Phytoplankton leakage fraction of coccolithophores' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ratiol1' Vinfo( 2)='Labile fraction' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ratiol2' Vinfo( 2)='Labile fraction for phytoplankton' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='wsdn' Vinfo( 2)='sinking velocity of detritus N' Vinfo( 3)='meter day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='wsdc' Vinfo( 2)='sinking velocity of detritus C' Vinfo( 3)='meter day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='wsdsi' Vinfo( 2)='sinking velocity of detritus silicate' Vinfo( 3)='meter day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='wsdca' Vinfo( 2)='sinking velocity of particulate inorganic carbon' Vinfo( 3)='meter day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='wsp1' Vinfo( 2)='sinking velocity of small phytoplankton' Vinfo( 3)='meter day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='wsp2' Vinfo( 2)='sinking velocity of diatom' Vinfo( 3)='meter day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='wsp3' Vinfo( 2)='sinking velocity of coccolithophores' Vinfo( 3)='meter day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='pco2a' Vinfo( 2)='air CO2 partial pressure' Vinfo( 3)='parts per million by volume' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='p2n' Vinfo( 2)='Phosphorus to nitrogen ratio' Vinfo( 3)='mole_P mole_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='o2no' Vinfo( 2)='oxygen to nitrate ratio' Vinfo( 3)='mole_O2 mole_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='o2nh' Vinfo( 2)='oxygen to ammonium ratio' Vinfo( 3)='mole_O2 mole_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='cnb' Vinfo( 2)='C:N in bacteria' Vinfo( 3)='mole_C mole_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='apsilon' Vinfo( 2)='Ratio of PIC to organic carbon in coccolithophores' Vinfo( 3)='mole_C mole_N-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ro5' Vinfo( 2)='grazing preference for diatom' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ro6' Vinfo( 2)='grazing preference for mesozooplankton' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ro7' Vinfo( 2)='grazing preference for detritus' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ro10' Vinfo( 2)='grazing preference for coccolithophores' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='rop' Vinfo( 2)='grazing preference for small phytoplankton' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='rob' Vinfo( 2)='grazing preference for bacteria' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='kabac' Vinfo( 2)='Half saturation for ammonium uptake by bacteria' Vinfo( 3)='mmol_N m-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='klbac' Vinfo( 2)='Half saturation for labile DOC uptake' Vinfo( 3)='mmol_C m-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ksdoc' Vinfo( 2)='Half saturation for semi-labile DOC uptake' Vinfo( 3)='mmol_C m-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ksdon' Vinfo( 2)='Half saturation for semi-labile DON uptake' Vinfo( 3)='mmol_N m-3' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ratiob' Vinfo( 2)='Bacteria growth loss fraction' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='ratiobc' Vinfo( 2)='Color fraction of Bacteria loss' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='RtUVLDOC' Vinfo( 2)='Rate of conversion of colored labile DOC to labile DOC' Vinfo( 3)='mmol_C m-2 d-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='RtUVSDOC' Vinfo( 2)='Conversion of colored semi-labile DOC to labile DOC' Vinfo( 3)='mmol_C m-2 d-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='RtUVLDIC' Vinfo( 2)='Rate of conversion of colored labile DOC to DIC' Vinfo( 3)='mmol_C m-2 d-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='RtUVSDIC' Vinfo( 2)='Rate of conversion of colored semi-labile DOC to DIC' Vinfo( 3)='mmol_C m-2 d-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='colorFR1' Vinfo( 2)='Color fraction for labile DOC' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='colorFR2' Vinfo( 2)='Color fraction for semi-labile DOC' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN #ifdef IRON_LIMIT Vinfo( 1)='T_Fe' Vinfo( 2)='Iron updake timescale' Vinfo( 3)='day' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='A_Fe' Vinfo( 2)='Empirical FE:C power' Vinfo( 3)='nondimensional' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='B_Fe' Vinfo( 2)='Empirical FE:C coefficient' Vinfo( 3)='meter-1 C' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='S1_FeC' Vinfo( 2)='Small phytoplankton Fe:C at F=0.5' Vinfo( 3)='muM-Fe/M-C' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='S2_FeC' Vinfo( 2)='Large phytoplankton Fe:C at F=0.5' Vinfo( 3)='muM-Fe/M-C' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='S3_FeC' Vinfo( 2)='Coccolithophores Fe:C at F=0.5' Vinfo( 3)='muM-Fe/M-C' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN Vinfo( 1)='FeRR' Vinfo( 2)='Fe remineralization rate' Vinfo( 3)='day-1' status=def_var(ng, model, ncid, varid, NF_TYPE, & & 1, (/0/), Aval, Vinfo, ncname, & & SetParAccess = .FALSE.) IF (exit_flag.ne.NoError) RETURN #endif �