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