I am still a novice at ocean modeling, and I am having an issue with cold water (< 5º in shallow water) forming at my boundaries during August in the Gulf of Mexico. I have checked and double checked my climatology, boundary, and forcing files, but can't see anything that should cause this. I have my boundary conditions sets as:
LBC(isFsur) == Che Che Che Che ! free-surface
LBC(isUbar) == Shc Shc Shc Shc ! 2D U-momentum
LBC(isVbar) == Shc Shc Shc Shc ! 2D V-momentum
LBC(isUvel) == Rad Rad Rad Rad ! 3D U-momentum
LBC(isVvel) == Rad Rad Rad Rad ! 3D V-momentum
LBC(isMtke) == Per Per Per Per ! mixing TKE
LBC(isTvar) == RadNud RadNud RadNud RadNud \ ! temperature
RadNud RadNud RadNud RadNud ! salinity
The problem begins right at the beginning of the model run. Does anyone know what might cause this? Any help would be greatly appreciated. I would attach some short mp4 animations, but it doesn't seem that I am allowed.
The cold water is forming at the ocean boundaries, not at the land/sea boundaries. I have TNUDG set to 2*0.0d0, which now that I look at it, seems to imply zero. Could that be my problem? As far as my grid goes, it is the same resolution as my climatology and boundary files, although finer than my forcing file.
So I figured out you can upload animated gifs, so I am attaching a couple here. I am only doing short model runs at this point. It is the bottom layer that the cold is creeping in, and the surface has the opposite problem. I did try setting my nudging timescales, but the problem still remains:
! Nudging/relaxation time scales, inverse scales will be computed
! internally, [1:Ngrids].
TNUDG == 2*10.0d0 ! days
ZNUDG == 0.0d0 ! days
M2NUDG == 0.0d0 ! days
M3NUDG == 10.0d0 ! days
! Factor between passive (outflow) and active (inflow) open boundary
! conditions, [1:Ngrids]. If OBCFAC > 1, nudging on inflow is stronger
! than on outflow (recommended).
OBCFAC == 1.0d0 ! nondimensional
I'm hoping this equates to an inflow timescale of 1 day and outflow timescale of 10 days, as in Fennel et al. (2009), but I might be wrong.
Attachments
ms_sound_3d_model_surface.gif (410.7 KiB) Viewed 5909 times
ms_sound_3d_model_bottom.gif (415.25 KiB) Viewed 5909 times
My boundaries were derived from my climatology file using SeaPy. I created an animation of the cross-section of the southern boundary and everything looked reasonable. So how would I set my inflow/outflow to 1 and 10 days? Thanks.
! Nudging/relaxation time scales, inverse scales will be computed
! internally, [1:Ngrids].
TNUDG == 2*10.0d0 ! days
ZNUDG == 0.0d0 ! days
M2NUDG == 0.0d0 ! days
M3NUDG == 0.0d0 ! days
! Factor between passive (outflow) and active (inflow) open boundary
! conditions, [1:Ngrids]. If OBCFAC > 1, nudging on inflow is stronger
! than on outflow (recommended).
OBCFAC == 10.0d0 ! nondimensional
Actually, looking over my ocean.in file it seems that the TCLM_NUDGING macro (among others) is obsolete. So I will try to turn on the logical switches for LTracerCLM and LnudgeTCLM and defining ANA_NUDGCOEF in the header file.
Well the error above was cause by me not having set temp:time in my climatology file. So I managed to get the model running again, but I still have the same problem. Any insights would be greatly appreciated.
I think I've figured out my problem. I think the data arrays in my boundary files are backward. At first I thought the lines at the edges of my animations were artifacts of my laziness. But that is the actually output in my _his file. It seems the script I used to create my boundary file was expecting the climatology data in the inverse order. I guess there's only one way to find out at this point. At least I learned a lot about nudging