Writing an equantion in Compaq Visual fortran 6.1

[collins-esi]

hi guys , i am in trouble , i am in Phd degree and i need to write this equation in Compaq Visual fortran 6.1 but nothing work
i need your help, please

thank you very much

 

[salgerman]

I suggest you do a few things:
- break down the equation into several lines to make it easier to read and make it look as much as possible to the textbook
- leave the "sign(0.8,-0.2)" operation out of the double loop, its value does not change with the indices...it will always be -0.8
- BECAUSE the exponent "(am-1)" turns out to be negative, I would suggest you move everything to the denominator so that you can raise it to a positive power

do i = 1, L1
  do j = 1, M1
    A = ( (u(i,j)-u(i-1,j))/xdif(i) )**2 + ( (v(i,j)-v(i,j-1))/ydif(j) )**2
    B = ( (v(i,j)-v(i-1,j))/xdif(i)      +   (u(i,j)-u(i,j-1))/ydif(j) )**2
    
    Viscapp(i,j) = 1.0 / ( ( sqrt( 2*A + B ) )**(1.8) )
  end do
end do