115 |
tmpx=tmpfld2d(i,j,bi,bj) |
tmpx=tmpfld2d(i,j,bi,bj) |
116 |
if (maskC(i,j,kk,bi,bj) .ne. 0.) then |
if (maskC(i,j,kk,bi,bj) .ne. 0.) then |
117 |
fctiletot = fctiletot |
fctiletot = fctiletot |
118 |
& + tmpx*cos(yc(i,j,bi,bj)*deg2rad) |
& + tmpx* _rA(i,j,bi,bj)/rhoConstFresh |
119 |
sumtot = sumtot |
sumtot = sumtot |
120 |
& + cos(yc(i,j,bi,bj)*deg2rad) |
& + _rA(i,j,bi,bj) |
121 |
num_sfluxmm(bi,bj) = num_sfluxmm(bi,bj) + 1 |
num_sfluxmm(bi,bj) = num_sfluxmm(bi,bj) + 1 |
122 |
endif |
endif |
123 |
enddo |
enddo |
125 |
enddo |
enddo |
126 |
enddo |
enddo |
127 |
|
|
128 |
_GLOBAL_SUM_R8( sumtot , myThid ) |
_GLOBAL_SUM_RL( sumtot , myThid ) |
129 |
_GLOBAL_SUM_R8( fctiletot , myThid ) |
_GLOBAL_SUM_RL( fctiletot , myThid ) |
130 |
|
|
131 |
if (sumtot.eq.0.) sumtot = 1. |
if (sumtot.eq.0.) sumtot = 1. |
132 |
|
|
141 |
write(standardmessageunit,'(A,I5,2(X,D22.14))') |
write(standardmessageunit,'(A,I5,2(X,D22.14))') |
142 |
& ' --> bal_sfluxmm =', irec, |
& ' --> bal_sfluxmm =', irec, |
143 |
& fctiletot/sumtot, |
& fctiletot/sumtot, |
144 |
& ((fctiletot/sumtot)/wmean_sflux)**2 |
& objf_sfluxmm |
145 |
|
|
146 |
enddo |
enddo |
147 |
|
|