138 |
|
|
139 |
if (sumtot.eq.0.) sumtot = 1. |
if (sumtot.eq.0.) sumtot = 1. |
140 |
|
|
141 |
do bj = jtlo,jthi |
if ( wmean_sflux .NE. 0. ) then |
142 |
do bi = itlo,ithi |
do bj = jtlo,jthi |
143 |
|
do bi = itlo,ithi |
144 |
fctilemm(bi,bj) = fctilemm(bi,bj) / sumtot |
fctilemm(bi,bj) = fctilemm(bi,bj) / sumtot |
145 |
objf_sfluxmm(bi,bj) = objf_sfluxmm(bi,bj) |
objf_sfluxmm(bi,bj) = objf_sfluxmm(bi,bj) |
146 |
& + (fctilemm(bi,bj)/wmean_sflux/nyearsrec)**2 |
& + (fctilemm(bi,bj)/wmean_sflux/nyearsrec)**2 |
147 |
enddo |
enddo |
148 |
enddo |
enddo |
149 |
|
endif |
150 |
|
|
151 |
c-- diagnostic: imbalance per year: |
c-- diagnostic: imbalance per year: |
152 |
tmpx2 = wsfluxmm(1,1) * fctiletot / sumtot |
tmpx2 = wsfluxmm(1,1) * fctiletot / sumtot |