54 |
#ifdef ALLOW_AUTODIFF |
#ifdef ALLOW_AUTODIFF |
55 |
|
|
56 |
adfc = 1.0 |
adfc = 1.0 |
57 |
|
#ifdef ALLOW_DIC |
58 |
|
adtotcost = 1.0 |
59 |
|
#endif |
60 |
|
|
61 |
do bj = jtlo,jthi |
do bj = jtlo,jthi |
62 |
do bi = itlo,ithi |
do bi = itlo,ithi |
63 |
#ifdef ALLOW_COST_VECTOR |
#ifdef ALLOW_COST_VECTOR |
67 |
end do |
end do |
68 |
#endif |
#endif |
69 |
#ifdef ALLOW_COST_STATE_FINAL |
#ifdef ALLOW_COST_STATE_FINAL |
|
c do k=1,nr |
|
70 |
do j=1,sNy |
do j=1,sNy |
71 |
do i=1,sNx |
do i=1,sNx |
72 |
objf_state_final(i,j,bi,bj) = 0. _d 0 |
do k=1,4*Nr |
73 |
|
objf_state_final(i,j,bi,bj,k) = 0. _d 0 |
74 |
|
enddo |
75 |
|
objf_state_final(i,j,bi,bj,4*Nr+1) = 0. _d 0 |
76 |
cph No init. of cost_state_final here, |
cph No init. of cost_state_final here, |
77 |
cph because we need it in ADM*TLM |
cph because we need it in ADM*TLM |
78 |
end do |
end do |
79 |
end do |
end do |
|
c end do |
|
80 |
#endif |
#endif |
81 |
enddo |
enddo |
82 |
enddo |
enddo |