7 |
C !ROUTINE: OPPS_CALC |
C !ROUTINE: OPPS_CALC |
8 |
|
|
9 |
C !INTERFACE: ====================================================== |
C !INTERFACE: ====================================================== |
10 |
subroutine OPPS_CALC( |
SUBROUTINE OPPS_CALC( |
11 |
U tracerEnv, |
U tracerEnv, |
12 |
I wVel, |
I wVel, |
13 |
I kMax, nTracer, nTracerInuse, |
I kMax, nTracer, nTracerInuse, |
920 |
c and eq 5 from the paper |
c and eq 5 from the paper |
921 |
c |
c |
922 |
if(debug) then |
if(debug) then |
923 |
write(0,*)"Dd,De,entr,k ",Dd(k2),De(k2),entrainrate,k2 |
write(0,*)"Dd,De,entr,k ",Dd(k2),De(k2),entrainrate,k2 |
924 |
endif |
endif |
925 |
c |
c |
926 |
wsqr = wsqr - wsqr*abs(entrainrate)+ 9.81* |
wsqr = wsqr - wsqr*abs(entrainrate)+ 9.81* |