In set_vbc.F, you'll find the following code helps a lot:
# ifdef SALINITY ! ! ! Multiply fresh water flux with surface salinity. If appropriate, ! apply correction. ! ! DO j=JstrR,JendR DO i=IstrR,IendR # if defined SCORRECTION stflx(i,j,isalt)=stflx(i,j,isalt)*t(i,j,N(ng),nrhs,isalt) & & Tnudg(isalt,ng)*Hz(i,j,N(ng))* & & (t(i,j,N(ng),nrhs,isalt)sss(i,j)) # elif defined SRELAXATION stflx(i,j,isalt)=Tnudg(isalt,ng)*Hz(i,j,N(ng))* & & (t(i,j,N(ng),nrhs,isalt)sss(i,j)) # else stflx(i,j,isalt)=stflx(i,j,isalt)*t(i,j,N(ng),nrhs,isalt) # endif btflx(i,j,isalt)=btflx(i,j,isalt)*t(i,j,1,nrhs,isalt) END DO END DO # endif
