1 | setstate u 0 |
---|
2 | setstate v 0 |
---|
3 | setstate u_d 0 |
---|
4 | setstate v_d 0 |
---|
5 | |
---|
6 | #loadstate u 0 ij 0 A 1 ../obs_float/uzero.dat D |
---|
7 | #loadstate v 0 ij 0 A 1 ../obs_float/vzero.dat D |
---|
8 | |
---|
9 | |
---|
10 | #Conditions initiales : |
---|
11 | loadstate r_float 0 i 0 A 1 ../obs_float/init.dat D |
---|
12 | setstate r_float_d 0 |
---|
13 | |
---|
14 | #Observations : |
---|
15 | loadobs r_cout_d 0 i 0 A 1 ../obs_float/obs.dat D |
---|
16 | |
---|
17 | set_modeltime 0 |
---|
18 | print_cost ON |
---|
19 | |
---|
20 | setm_impres 3 |
---|
21 | setm_io 6 |
---|
22 | setm_mode 0 |
---|
23 | #set_nbiter 100 |
---|
24 | set_nbiter 20 |
---|
25 | setm_nsim 200 |
---|
26 | setm_dxmin 1.0e-12 |
---|
27 | setm_epsg 1.0e-12 |
---|
28 | setm_ddf1 1.0 |
---|
29 | |
---|
30 | |
---|
31 | run_inc |
---|
32 | goto fin |
---|
33 | |
---|
34 | #uniquement activer les trajectoires de type M model |
---|
35 | activ M only |
---|
36 | |
---|
37 | forward |
---|
38 | |
---|
39 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_forward.dat |
---|
40 | |
---|
41 | activ D only |
---|
42 | |
---|
43 | #testdf 1 29 21 4 r 0.0001 0.0001 |
---|
44 | |
---|
45 | |
---|
46 | |
---|
47 | |
---|
48 | print_cost ON |
---|
49 | |
---|
50 | setm_impres 3 |
---|
51 | setm_io 6 |
---|
52 | setm_mode 0 |
---|
53 | #set_nbiter 100 |
---|
54 | set_nbiter 20 |
---|
55 | setm_nsim 200 |
---|
56 | setm_dxmin 1.0e-12 |
---|
57 | setm_epsg 1.0e-12 |
---|
58 | setm_ddf1 1.0 |
---|
59 | |
---|
60 | RUNM |
---|
61 | savestate r_cout_d 0 i 0 A 3 ../obs_float/r_cout_d_optim_1.dat |
---|
62 | |
---|
63 | |
---|
64 | |
---|
65 | update_uv |
---|
66 | setstate u_d 0 |
---|
67 | setstate v_d 0 |
---|
68 | |
---|
69 | |
---|
70 | |
---|
71 | savestate u 0 ij 0 A 3 ../obs_float/uoptim1.dat |
---|
72 | savestate v 0 ij 0 A 3 ../obs_float/voptim1.dat |
---|
73 | |
---|
74 | #FIN 1ere OPTIM |
---|
75 | |
---|
76 | set_modeltime 0 |
---|
77 | activ M only |
---|
78 | forward |
---|
79 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim1.dat |
---|
80 | |
---|
81 | |
---|
82 | activ D only |
---|
83 | RUNM |
---|
84 | |
---|
85 | update_uv |
---|
86 | |
---|
87 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim.dat |
---|
88 | savestate r_cout_d 0 i 0 A 3 ../obs_float/r_cout_d_optim_2.dat |
---|
89 | savestate u 0 ij 0 A 3 ../obs_float/uoptim2.dat |
---|
90 | savestate v 0 ij 0 A 3 ../obs_float/voptim2.dat |
---|
91 | savestate u_d 0 ij 0 A 3 ../obs_float/udoptim.dat |
---|
92 | savestate v_d 0 ij 0 A 3 ../obs_float/vdoptim.dat |
---|
93 | |
---|
94 | |
---|
95 | #FIN 2eme OPTIM |
---|
96 | |
---|
97 | |
---|
98 | set_modeltime 0 |
---|
99 | activ M only |
---|
100 | forward |
---|
101 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim2.dat |
---|
102 | |
---|
103 | |
---|
104 | activ D only |
---|
105 | RUNM |
---|
106 | |
---|
107 | update_uv |
---|
108 | savestate u 0 ij 0 A 3 ../obs_float/uoptim3.dat |
---|
109 | savestate v 0 ij 0 A 3 ../obs_float/voptim3.dat |
---|
110 | |
---|
111 | |
---|
112 | #FIN 3eme OPTIM |
---|
113 | |
---|
114 | |
---|
115 | set_modeltime 0 |
---|
116 | activ M only |
---|
117 | forward |
---|
118 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim3.dat |
---|
119 | |
---|
120 | |
---|
121 | activ D only |
---|
122 | RUNM |
---|
123 | |
---|
124 | update_uv |
---|
125 | savestate u 0 ij 0 A 3 ../obs_float/uoptim4.dat |
---|
126 | savestate v 0 ij 0 A 3 ../obs_float/voptim4.dat |
---|
127 | |
---|
128 | |
---|
129 | #FIN 4eme OPTIM |
---|
130 | |
---|
131 | |
---|
132 | set_modeltime 0 |
---|
133 | activ M only |
---|
134 | forward |
---|
135 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim4.dat |
---|
136 | |
---|
137 | activ D only |
---|
138 | RUNM |
---|
139 | |
---|
140 | update_uv |
---|
141 | savestate u 0 ij 0 A 3 ../obs_float/uoptim5.dat |
---|
142 | savestate v 0 ij 0 A 3 ../obs_float/voptim5.dat |
---|
143 | |
---|
144 | |
---|
145 | #FIN 5eme OPTIM |
---|
146 | |
---|
147 | |
---|
148 | set_modeltime 0 |
---|
149 | activ M only |
---|
150 | forward |
---|
151 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim5.dat |
---|
152 | |
---|
153 | |
---|
154 | activ D only |
---|
155 | RUNM |
---|
156 | |
---|
157 | update_uv |
---|
158 | savestate u 0 ij 0 A 3 ../obs_float/uoptim6.dat |
---|
159 | savestate v 0 ij 0 A 3 ../obs_float/voptim6.dat |
---|
160 | |
---|
161 | |
---|
162 | #FIN 6eme OPTIM |
---|
163 | |
---|
164 | |
---|
165 | set_modeltime 0 |
---|
166 | activ M only |
---|
167 | forward |
---|
168 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim6.dat |
---|
169 | |
---|
170 | activ D only |
---|
171 | RUNM |
---|
172 | |
---|
173 | update_uv |
---|
174 | savestate u 0 ij 0 A 3 ../obs_float/uoptim7.dat |
---|
175 | savestate v 0 ij 0 A 3 ../obs_float/voptim7.dat |
---|
176 | |
---|
177 | |
---|
178 | #FIN 7eme OPTIM |
---|
179 | |
---|
180 | |
---|
181 | set_modeltime 0 |
---|
182 | activ M only |
---|
183 | forward |
---|
184 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim7.dat |
---|
185 | |
---|
186 | activ D only |
---|
187 | RUNM |
---|
188 | |
---|
189 | update_uv |
---|
190 | savestate u 0 ij 0 A 3 ../obs_float/uoptim8.dat |
---|
191 | savestate v 0 ij 0 A 3 ../obs_float/voptim8.dat |
---|
192 | |
---|
193 | |
---|
194 | #FIN 8eme OPTIM |
---|
195 | |
---|
196 | |
---|
197 | set_modeltime 0 |
---|
198 | activ M only |
---|
199 | forward |
---|
200 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim8.dat |
---|
201 | |
---|
202 | |
---|
203 | activ D only |
---|
204 | RUNM |
---|
205 | |
---|
206 | update_uv |
---|
207 | savestate u 0 ij 0 A 3 ../obs_float/uoptim9.dat |
---|
208 | savestate v 0 ij 0 A 3 ../obs_float/voptim9.dat |
---|
209 | |
---|
210 | |
---|
211 | #FIN 9eme OPTIM |
---|
212 | |
---|
213 | |
---|
214 | set_modeltime 0 |
---|
215 | activ M only |
---|
216 | forward |
---|
217 | savestate r_float 0 i 0 A 3 ../obs_float/r_float_optim9.dat |
---|
218 | |
---|
219 | |
---|
220 | |
---|
221 | print_cost OFF |
---|
222 | #testof 0.1 1 10 10 1 |
---|
223 | #testing the gradients of each module |
---|
224 | |
---|
225 | fin |
---|