54 |
|
|
55 |
adfc = 1.0 |
adfc = 1.0 |
56 |
do bj = jtlo,jthi |
do bj = jtlo,jthi |
57 |
do bi = itlo,ithi |
do bi = itlo,ithi |
58 |
#ifdef ALLOW_COST_VECTOR |
#ifdef ALLOW_COST_VECTOR |
59 |
|
do i=1,sNx |
60 |
|
objf_vector(i,bi,bj) = 0. _d 0 |
61 |
|
adobjf_vector(i,bi,bj) = 1. _d 0 |
62 |
|
end do |
63 |
|
#endif |
64 |
|
#ifdef ALLOW_COST_STATE_FINAL |
65 |
|
c do k=1,nr |
66 |
|
do j=1,sNy |
67 |
do i=1,sNx |
do i=1,sNx |
68 |
objf_vector(i,bi,bj) = 0. _d 0 |
objf_state_final(i,j,bi,bj) = 0. _d 0 |
69 |
adobjf_vector(i,bi,bj) = 1. _d 0 |
cph No init. of cost_state_final here, |
70 |
|
cph because we need it in ADM*TLM |
71 |
end do |
end do |
72 |
|
end do |
73 |
|
c end do |
74 |
#endif |
#endif |
75 |
enddo |
enddo |
76 |
enddo |
enddo |