Hi,
I am a new user in ROMS, and I try to model an overflow,
I need a constant salt forcing at the northern side of the domain,
Is there a way to relax some grids of the domain into specified number,
just like regular relaxation, but not only at boundaries also in the domain,
or anything works for this, such as source or something like that?
Regional relaxation
1) define
#define TCLIMATOLOGY
#define TCLM_NUDGING
2) create a climatology file with your salinity fields
3) go into set_nudgcoef.F and define the region in which the relaxation has to take place, setting wrk (see below) different from zero ONLY where you want to relax. Units of wrk are [1/s].
! Initialize.
!
DO j=JstrR,JendR
DO i=IstrR,IendR
wrk(i,j)=0.0_r8
END DO
END DO
# ifdef ZCLM_NUDGING
DO j=JstrR,JendR
DO i=IstrR,IendR
CLIMA(ng)%Znudgcof(i,j)=wrk(i,j)
END DO
END DO
# endif
# ifdef TCLM_NUDGING
DO j=JstrR,JendR
DO i=IstrR,IendR
CLIMA(ng)%Tnudgcof(i,j,itemp)=wrk(i,j)
CLIMA(ng)%Tnudgcof(i,j,isalt)=wrk(i,j)
END DO
END DO
# endif
Hope this help,
Jacopo
#define TCLIMATOLOGY
#define TCLM_NUDGING
2) create a climatology file with your salinity fields
3) go into set_nudgcoef.F and define the region in which the relaxation has to take place, setting wrk (see below) different from zero ONLY where you want to relax. Units of wrk are [1/s].
! Initialize.
!
DO j=JstrR,JendR
DO i=IstrR,IendR
wrk(i,j)=0.0_r8
END DO
END DO
# ifdef ZCLM_NUDGING
DO j=JstrR,JendR
DO i=IstrR,IendR
CLIMA(ng)%Znudgcof(i,j)=wrk(i,j)
END DO
END DO
# endif
# ifdef TCLM_NUDGING
DO j=JstrR,JendR
DO i=IstrR,IendR
CLIMA(ng)%Tnudgcof(i,j,itemp)=wrk(i,j)
CLIMA(ng)%Tnudgcof(i,j,isalt)=wrk(i,j)
END DO
END DO
# endif
Hope this help,
Jacopo