132 |
|
|
133 |
if ( wmean_sflux .NE. 0. ) then |
if ( wmean_sflux .NE. 0. ) then |
134 |
objf_sfluxmm = objf_sfluxmm |
objf_sfluxmm = objf_sfluxmm |
135 |
& + ( (fctiletot/sumtot-1.D0)/wmean_sflux )**2 |
& + ( (fctiletot/sumtot)/wmean_sflux )**2 |
136 |
else |
else |
137 |
objf_sfluxmm = 0. _d 0 |
objf_sfluxmm = 0. _d 0 |
138 |
endif |
endif |
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 |
& ((fctiletot/sumtot)/wmean_sflux)**2 |
145 |
|
|
146 |
enddo |
enddo |
147 |
|
|