209 |
WVSurf(I2,myThid) = 0. |
WVSurf(I2,myThid) = 0. |
210 |
K = kGrd(I2) |
K = kGrd(I2) |
211 |
IF (K.GT.1) THEN |
IF (K.GT.1) THEN |
212 |
WVSurf(I2,myThid) = (LOG(SIGH(K))-SIGL(K))*WVI(K-1,2) |
C- full cell version of Franco Molteni formula: |
213 |
|
c WVSurf(I2,myThid) = (LOG(SIGH(K))-SIGL(K))*WVI(K-1,2) |
214 |
|
C- partial cell version using true log-P extrapolation: |
215 |
|
WVSurf(I2,myThid) = (LOG(PSA(I2))-SIGL(K))*WVI(K-1,1) |
216 |
C- like in the old code: |
C- like in the old code: |
217 |
c WVSurf(I2,myThid) = WVI(K,2) |
c WVSurf(I2,myThid) = WVI(K,2) |
218 |
ENDIF |
ENDIF |