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 |