68 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
69 |
CADJ STORE surf_el_streamice = comlev1, key = ikey_dynamics, |
CADJ STORE surf_el_streamice = comlev1, key = ikey_dynamics, |
70 |
CADJ & kind = isbyte |
CADJ & kind = isbyte |
71 |
|
CADJ STORE H_streamice_prev = comlev1, key = ikey_dynamics, |
72 |
|
CADJ & kind = isbyte |
73 |
#endif |
#endif |
74 |
|
|
75 |
C-- Calculate cost function on tile of this instance |
C-- Calculate cost function on tile of this instance |
84 |
cost_func1_streamice (bi,bj) = |
cost_func1_streamice (bi,bj) = |
85 |
& cost_func1_streamice (bi,bj) + |
& cost_func1_streamice (bi,bj) + |
86 |
& 0.5 * (S_obs(i,j,bi,bj)- |
& 0.5 * (S_obs(i,j,bi,bj)- |
87 |
& surf_el_streamice(i,j,bi,bj))**2 |
& surf_el_streamice(i,j,bi,bj))**2 + |
88 |
|
& 0.5 * (H_streamice(i,j,bi,bj)- |
89 |
|
& H_streamice_prev(i,j,bi,bj))**2 |
90 |
|
|
91 |
end do |
end do |
92 |
end do |
end do |