| 86 |
|
|
| 87 |
|
|
| 88 |
! S_obs(i,j,bi,bj) = 0.0 |
! S_obs(i,j,bi,bj) = 0.0 |
| 89 |
|
if (streamice_cost_mask(i,j,bi,bj).eq.1.0) THEN |
| 90 |
cost_func1_streamice (bi,bj) = |
cost_func1_streamice (bi,bj) = |
| 91 |
& cost_func1_streamice (bi,bj) + |
& cost_func1_streamice (bi,bj) + |
| 92 |
& streamice_wgt_surf * (S_obs(i,j,bi,bj)- |
& streamice_wgt_surf * (S_obs(i,j,bi,bj)- |
| 93 |
& surf_el_streamice(i,j,bi,bj))**2 / Nx / Ny + |
& surf_el_streamice(i,j,bi,bj))**2 / Nx / Ny + |
| 94 |
& streamice_wgt_drift * (H_streamice(i,j,bi,bj)- |
& streamice_wgt_drift * (H_streamice(i,j,bi,bj)- |
| 95 |
& H_streamice_prev(i,j,bi,bj))**2 / Nx / Ny |
& H_streamice_prev(i,j,bi,bj))**2 / Nx / Ny |
| 96 |
|
endif |
| 97 |
|
|
| 98 |
|
|
| 99 |
end do |
end do |