Ocean Modeling Discussion

In the ana_dqdsst file, I can modify the timescale of SST relaxation as seen below. #ifdef MY_APPLICATION fac=day2sec/30*1.0_r8 ! 30 day relaxation scale 1/s/decC DO j=JstrT,JendT DO i=IstrT,IendT dqdt(i,j)=fac*Hz(i,j,N(ng)) END DO END DO #else If I want a one-day relaxation ( strong ), then I just ...

If I understand correctly, then to relax SST, we need to do flux correction by switching on QCORRECTION. This is done by using dQdsst which is kinematic surface net heat flux sensitivity to SST. How good it is to use monthly climatological dqdsst computed from other climatological variables like ...

Thanks for the reply. I have recompiled the model using the following cppdefs and forced the model with wind stress ( not wind speed ). #undef TS_MPDATA #undef INLINE_2DIO #define SPONGE #define UV_QDRAG #define VISC_GRID #define UV_ADV #define TS_A6VADVECTION #define TS_A4HADVECTION #define UV ...

Normally I force the model with wind speed at 10m. What changes do I need to do if I want to force the model with wind stress ? Actually I am developing LETKF data assimilation for ROMS. It is using ensemble atmospheric fluxes and I am provided with wind stress rather than wind speed. I am not ...

Thanks for the mail Kate. I am however a little confused. Do you mean to say that I need to get into the code and modify ? If I switch on QCORRECTION, where does the model get SST from ( roms_frc.nc ?? ). Please correct me if I am wrong. In order to relax the model SST to some SST, I need to do the ...

Suppose I want to relax the model SST to a pre-arranged SST netcdf file over a relaxation time scale of T days. How do I do it in ROMS ?

We are trying to completely eliminate wave reflections from the boundaries which are closed. We have used SPONGE in the *.h file and made changes in ana_hmixcoef.h. We have only modified cff2 to 100 and Iwrk to be 8. But still we are unable to cancel the reflections from the boundaries. Are we doing ...

Thanks. Actually we are trying to run a series of controlled experiments. For example, I was wondering if I put a constant f ( Coriolis parameter ) in the grid file, will this imply that there will be no beta-effect ? Or are there more simpler way to introduce f = constant in governing equations ?

Thanks for the reply Kate. Actually, I wanted to know if ROMS can be run only in the baroclinic mode thereby avoiding all barotropic contributions ?

If we put NDTFAST > DT, does that mean that the ocean state will be purely baroclinic in nature and all barotropic contributions will be filtered out ?

We are trying to run roms in the barotropic mode for the Indian Ocean wherein there is open boundary in South, East & West. I was wondering how to generate daily realistic boundary condition for ubar, vbar and zeta for any particular year ( say 2011 ).

Thank you for your response. I got the equations. However if I want to run the model in the barotropic mode dropping all the non-linear terms from the equations for some controlled experiments, how do I do it ?
Also, how is it possible to generate barotropic boundary conditions ?

Thanks Kate for the reply.
But what's the forcing term look like ? Is it curl of the wind stress/H where H is the bottom topography ? Does it have any bottom friction term ? Is there no dependence on beta ?

When ROMS is run in the barotropic mode, what is the explicit form of the barotropic equation for SSH that is being solved ?

I am planning to assimilate ( using Ensemble Kalman Filter technique ) sea level anomaly using satellite observations. I am new to ROMS and have gone through the documentation ( not in a detailed manner though ). I would like to ask how zeta is computed in ROMS. Is there a dynamical equation for ...