95 |
vst = 0. |
vst = 0. |
96 |
ssq = 0. |
ssq = 0. |
97 |
|
|
98 |
CALL BULKF_FORMULA_LANL( |
IF ( blk_nIter.EQ.0 ) THEN |
99 |
|
CALL BULKF_FORMULA_LANL( |
100 |
I uwind(i,j,bi,bj), vwind(i,j,bi,bj), wspeed(i,j,bi,bj), |
I uwind(i,j,bi,bj), vwind(i,j,bi,bj), wspeed(i,j,bi,bj), |
101 |
I Tair(i,j,bi,bj), Qair(i,j,bi,bj), cloud(i,j,bi,bj), Tsf, |
I Tair(i,j,bi,bj), Qair(i,j,bi,bj), cloud(i,j,bi,bj), Tsf, |
102 |
O flwup, flh, fsh, df0dT, ust, vst, evap, ssq, dEvdT, |
O flwup, flh, fsh, df0dT, ust, vst, evap, ssq, dEvdT, |
103 |
I iceornot, myThid ) |
I iceornot, myThid ) |
104 |
|
ELSE |
105 |
|
CALL BULKF_FORMULA_LAY( |
106 |
|
I uwind(i,j,bi,bj), vwind(i,j,bi,bj), wspeed(i,j,bi,bj), |
107 |
|
I Tair(i,j,bi,bj), Qair(i,j,bi,bj), Tsf, |
108 |
|
O flwup, flh, fsh, df0dT, ust, vst, evap, ssq, dEvdT, |
109 |
|
I iceornot, i,j,bi,bj,myThid ) |
110 |
|
ENDIF |
111 |
|
|
112 |
flwNet_dwn = flwdwn(i,j,bi,bj) - flwup |
flwNet_dwn = flwdwn(i,j,bi,bj) - flwup |
113 |
flxExceptSw = flwNet_dwn + fsh + flh |
flxExceptSw = flwNet_dwn + fsh + flh |