106 |
rmuu = darwin_rmuu |
rmuu = darwin_rmuu |
107 |
|
|
108 |
c find deepest wet layer |
c find deepest wet layer |
109 |
kbot = kmax |
kbot = MIN(kmax, Nr) |
110 |
DO WHILE (H(kbot).EQ.0 .AND. kbot.GT.0) |
DO WHILE (H(kbot).EQ.0 .AND. kbot.GT.1) |
111 |
kbot = kbot - 1 |
kbot = kbot - 1 |
112 |
ENDDO |
ENDDO |
113 |
|
IF (H(kbot).EQ.0) kbot = kbot - 1 |
114 |
|
|
115 |
DO nl = 1,tlam |
DO nl = 1,tlam |
116 |
DO k=1,Nr |
DO k=1,Nr |
201 |
amp1(nl,k) = c1 |
amp1(nl,k) = c1 |
202 |
amp2(nl,k) = c2 |
amp2(nl,k) = c2 |
203 |
ENDDO |
ENDDO |
204 |
|
IF (kbot .LT. Nr) THEN |
205 |
|
Estop(nl,kbot+1) = Esbot(nl,kbot) |
206 |
|
Eutop(nl,kbot+1) = Eubot(nl,kbot) |
207 |
|
ENDIF |
208 |
|
|
209 |
C endif thresh |
C endif thresh |
210 |
ENDIF |
ENDIF |