132 |
C ------ Augment the bin values |
C ------ Augment the bin values |
133 |
layers_UFlux(i,j,kgu(i,j),bi,bj) = |
layers_UFlux(i,j,kgu(i,j),bi,bj) = |
134 |
& layers_UFlux(i,j,kgu(i,j),bi,bj) + |
& layers_UFlux(i,j,kgu(i,j),bi,bj) + |
135 |
& dZZ * uVel(i,j,kci,bi,bj) * hFacW(i,j,kci,bi,bj) |
& dZZf(kk) * uVel(i,j,kci,bi,bj) * hFacW(i,j,kci,bi,bj) |
136 |
|
|
137 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
138 |
layers_HU(i,j,kgu(i,j),bi,bj) = layers_HU(i,j,kgu(i,j),bi,bj) |
layers_HU(i,j,kgu(i,j),bi,bj) = layers_HU(i,j,kgu(i,j),bi,bj) |
139 |
& + dZZ * hFacW(i,j,kci,bi,bj) |
& + dZZf(kk) * hFacW(i,j,kci,bi,bj) |
140 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
141 |
|
|
142 |
#endif /* LAYERS_UFLUX */ |
#endif /* LAYERS_UFLUX */ |
182 |
C ------ Augment the bin values |
C ------ Augment the bin values |
183 |
layers_VFlux(i,j,kgv(i,j),bi,bj) = |
layers_VFlux(i,j,kgv(i,j),bi,bj) = |
184 |
& layers_VFlux(i,j,kgv(i,j),bi,bj) |
& layers_VFlux(i,j,kgv(i,j),bi,bj) |
185 |
& + dZZ * vVel(i,j,kci,bi,bj) * hFacS(i,j,kci,bi,bj) |
& + dZZf(kk) * vVel(i,j,kci,bi,bj) * hFacS(i,j,kci,bi,bj) |
186 |
|
|
187 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
188 |
layers_HV(i,j,kgv(i,j),bi,bj) = layers_HV(i,j,kgv(i,j),bi,bj) |
layers_HV(i,j,kgv(i,j),bi,bj) = layers_HV(i,j,kgv(i,j),bi,bj) |
189 |
& + dZZ * hFacS(i,j,kci,bi,bj) |
& + dZZf(kk) * hFacS(i,j,kci,bi,bj) |
190 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
191 |
|
|
192 |
#endif /* LAYERS_VFLUX */ |
#endif /* LAYERS_VFLUX */ |
193 |
|
|
|
C k loop |
|
194 |
ENDDO |
ENDDO |
|
|
|
195 |
ENDDO |
ENDDO |
196 |
ENDDO |
ENDDO |
197 |
|
|