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 |
|
CADJ STORE H_streamice = comlev1, key = ikey_dynamics, |
74 |
|
CADJ & kind = isbyte |
75 |
#endif |
#endif |
76 |
|
|
77 |
C-- Calculate cost function on tile of this instance |
C-- Calculate cost function on tile of this instance |
86 |
cost_func1_streamice (bi,bj) = |
cost_func1_streamice (bi,bj) = |
87 |
& cost_func1_streamice (bi,bj) + |
& cost_func1_streamice (bi,bj) + |
88 |
& 0.5 * (S_obs(i,j,bi,bj)- |
& 0.5 * (S_obs(i,j,bi,bj)- |
89 |
& surf_el_streamice(i,j,bi,bj))**2 |
& surf_el_streamice(i,j,bi,bj))**2 + |
90 |
|
& 0.5 * (H_streamice(i,j,bi,bj)- |
91 |
|
& H_streamice_prev(i,j,bi,bj))**2 |
92 |
|
|
93 |
end do |
end do |
94 |
end do |
end do |