[MITgcm-support] Relaxation constant for salinity, tauSaltClimRelax in Global Ocean Simulation at 4deg resolution.
Krishnakumar Rajagopalan
krishna_raj_2010 at yahoo.com
Thu Sep 23 02:57:04 EDT 2010
Hi David,
Thanks a lot for your comments.
As you pointed out the units of ds/dt and (1/tauSaltClimRelax)*(theta - theta*)
are inconsistent. I only wanted to suggest a relationship between
(theta-theta*) and salinity.
Regards
Krishnakumar
________________________________
From: David Ferreira <dfer at mit.edu>
To: mitgcm-support at mitgcm.org
Sent: Tue, September 21, 2010 5:38:20 PM
Subject: Re: [MITgcm-support] Relaxation constant for salinity, tauSaltClimRelax
in Global Ocean Simulation at 4deg resolution.
Hi Krishnakumar,
My question is regarding the relaxation constant for salinity in the
tutorial (tutorial_global_oce_latlon). In the case of potential
temperature, we calculated the tauThetaClimRelax according to the
publication by Barnier et. al(Barnier et. al : Thermal forcing for a
global ocean circulation model using a three-year climatology of
ECMWF analyses, Journal of Marine Systems, 1995.) and are getting
values in the range of 60days, which is also the value chosen in
tutorial.
>For salinity, however, I am unable to get hold of any
>reference/publication which discusses in detail the calculation of
>tauSaltClimRelax. Please let me know how this quantity
>(tauSaltClimRelax) could be calculated.
>
There is no physical justification to tauSaltClimRelax. The non-physical
justification is that it "works", i.e it keeps the solution close to S*
(usually some observed climatology).
>Finally, I feel that the net fresh water surface flux, EmPmR, would
>depend primarily on temperature than on salinity. If true, the
>relaxation term in the salinity equation could look like,
>
>dS/dt= (1/tauSaltClimRelax)*(theta - theta*) + flux (T*,S*)
>
This looks very strange. It doesn't have the right dimension. Also, even if
evaporation could be approximated as a function of SST, it is less clear for
precipitation.
rather than
>
>dS/dt=(1/tauSaltClimRelax)*(S-S*) + flux(S*)
>
>But then, I am not sure whether the salinity solution would converge
>to S*. And of course the determination of tauSaltClimRelax would still
>require some reference.
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20100922/f41b33ae/attachment.htm>
More information about the MITgcm-support
mailing list