| 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 |