[MITgcm-support] Negative Salinities
Sergio Jaramillo
sju at eos.ubc.ca
Tue May 25 12:31:21 EDT 2004
Hi Martin,
Thanks a lot for your comments. I noticed that the CPP-flag does not
improve anything for me. The problem does occur after the first time
step, but I wouldn't think that it has anything to do with my forcing
fields as I haven't changed this from my old runs, and I don't use any
salinity filed as an input. I force the model with a body force that I
wrote into external_forcing.F. I will put here part of my data file in
case anyone sees something wrong and could give me any advice:
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
Tref = 14.1283, 11.9377, 9.89651, 8.154, 6.71021,
5.6647, 4.91792, 4.32048, 3.72305, 3.02605,
2.32905, 1.7814, 1.18397, 0.486969, -0.110463,
-0.558536, -1.0564, -1.70361, -2.35083, -2.84869,
-3.29676, -3.69505, -4.19291, -4.69077, -5.18863,
-5.58691, -5.9852, -6.28392, -6.6822, -7.18006,
-7.92685, -8.92257,
sRef= 32*35.,
viscAz=1.E-6,
viscAh=1.E-6,
no_slip_sides=.TRUE.,
no_slip_bottom=.TRUE.,
bottomDragLinear=0.,
viscA4=0.E12,
diffKhT=1.E-6,
diffKzT=1.5E-9,
diffKhS=1.E-6,
diffKzS=1.5E-9,
tempAdvScheme= 33,
saltAdvScheme= 33,
saltForcing= .FALSE.,
tempForcing=.FALSE.,
staggerTimeStep=.TRUE.,
exactConserv=.TRUE.,
f0=0.52,
beta=0.E-11,
tAlpha=2.E-3,
rhonil=1.0043E3,
sBeta =0.E-4,
gravity=9.81,
gBaro=9.81,
rigidLid=.FALSE.,
implicitFreeSurface=.TRUE.,
eosType='LINEAR',
hFacMin=0.2,
nonHydrostatic=.FALSE.,
readBinaryPrec=64,
globalFiles=.FALSE.,
vectorInvariantMomentum=.FALSE.,
&
Sergio Jaramillo
UBC
>
>Sergio,
>
>no use in defining this CPP-flag, it will not solve your problem. This
>routine raises a flag/warning, if it finds one or more negative
>salinity values in the current horizontal slab and resets the negative
>values to zero.
>But if you have salinities of +/-1.5e+139, your problem is likely to be
>something completely different. If this happens after one or two time
>steps, I am guessing, that there are some memory issues, such as a
>double precision field is read into a single precision variable or vice
>versa (for example, your forcing fields).
>
>Martin
>
>On Saturday, May 22, 2004, at 02:47 AM, Sergio Jaramillo wrote:
>
>
>
>>Dear MITgcm support,
>>
>>I just noticed that the salinity field that I get from my runs is
>>going negative like -1.5e+139!. Why could this be? I only input
>>sRef=32*35, and I do all the stratification using the temperature
>>profile. I also have set sBeta=0. Also I didn't have this problem
>>before , when I used an older version of the model (like a year ago,
>>v1.15 according to find_rho.F ).
>>
>>Finally I find this in find_rho.F:
>>
>>#ifdef CHECK_SALINITY_FOR_NEGATIVE_VALUES
>> CALL LOOK_FOR_NEG_SALINITY( bi, bj, iMin, iMax, jMin, jMax, k,
>> & sFld, myThid )
>>#endif
>>
>>Should I define this in CPP_OPTIONS.h?
>>
>>
>>Thanks a lot,
>>
>>Sergio Jaramillo
>>UBC
>>
>>_______________________________________________
>>MITgcm-support mailing list
>>MITgcm-support at mitgcm.org
>>http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
>>
>>
>>
>Martin Losch // mailto:mlosch at awi-bremerhaven.de
>Alfred-Wegener-Institut für Polar- und Meeresforschung
>Postfach 120161, 27515 Bremerhaven, Germany
>Tel./Fax: ++49(471)4831-1872/1797
>http://www.awi-bremerhaven.de/People/show?mlosch
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20040525/440552e1/attachment.htm>
More information about the MITgcm-support
mailing list