Hello All,
I am running an Indian Ocean model on mom4 with two inland seas, the
Persian Gulf and the Read Sea and no xland mixing. The run came down after
2.5 years and after a considerable amount of cussing and fretting, I
isolated the grid
point where the instability set in (the temp went negative and salinity
went to 80 g/kg). It happens in the Persian Gulf at
a corner point which is only two levels (20m) deep. Adjacent to it are
other points also two levels deep. When I put the following code at the
bottom of ocean_tracer.F90,
!Swathi..A dumb check to see if the instability at x=52.2, y=24.84 will go
away with this
!Remove this after you fix it
if (iscomp <= 36 .and. iecomp >= 36) then
if (jscomp <=87 .and. jecomp >= 87) then
tprog(36,87,1:2,1,taup1) = 26.6
tprog(36,87,1:2,2,taup1) = 39.1
endif
endif
!
I went past the instability and I was able to integrate for two more
years and the model came down again. I havent located the point where the
instablity set in (it is a pain to do) but my guess it probably
happened again at a corner point.
I am reasonalbly sure that the problem is with kpp as I have a parallel
run with pp mixing and full convection and it is continuing fine.
I wonder if anyone has noticed something similar and has a fix.
Thank you.
Swathi
-------------------------------------------------------
This SF.Net email sponsored by: Free pre-built ASP.NET sites including
Data Reports, E-commerce, Portals, and Forums are available now.
Download today and enter to win an XBOX or Visual Studio .NET.
http://aspnet.click-url.com/go/psa00100003ave/direct;at.aspnet_072303_01/01
_______________________________________________
FMS-mom4 mailing list
[EMAIL PROTECTED]
https://lists.sourceforge.net/lists/listinfo/fms-mom4
