calculate salinity in PSU according to
subroutine update_salinity() !! calculate salinity in PSU according to $$S = 35 + !! 1000*S_{UVic}$$ integer i,j,k ! local loop counters if (model%variable_needs_values(id_practical_salinity)) then do j=2,jmt-1 do k=1,km do i=2,imt-1 if (kmt(i,j) > 0) then salt(i,k,j) = 35._rk+1000._rke*t(i,k,j,isalt,0) end if end do end do end do end if end subroutine update_salinity