Bodyforce

[thomas.roc]

Hi everyone,

It's probably a dummy question but can anyone explain me the physical meaning of "wrk" in rhs3d.F in the part reserved for the bodyforce expression.
!
! Apply bottom stress as a bodyforce: determine the thickness (m)
! of the bottom layer; then add in bottom stress as a bodyfoce.
!

Code: Select all

      DO j=JstrV-1,Jend
        DO i=IstrU-1,Iend
         [color=#008040] wrk(i,j)=0.0_r8[/color]
        END DO
      END DO
      DO k=1,levbfrc(ng)
        DO j=JstrV-1,Jend
          DO i=IstrU-1,Iend
            wrk(i,j)=wrk(i,j)+Hz(i,j,k)
          END DO
        END DO
      END DO

I thought it should represent a work per unit density but, as far as I understood by reading the code, it does fit with the unit since "wrk"=[m].

Thanks


-T-