Last change
on this file since 147 was
147,
checked in by jbrlod, 9 years ago
|
add divergence term in cost function
|
File size:
832 bytes
|
Line | |
---|
1 | setstate u 0 |
---|
2 | setstate v 0 |
---|
3 | setstate u_d 0 |
---|
4 | setstate v_d 0 |
---|
5 | setstate u_norm 1 |
---|
6 | setstate v_norm 1 |
---|
7 | setstate uc_d 0 |
---|
8 | setstate vc_d 0 |
---|
9 | |
---|
10 | outoebx uc_d 1 0 |
---|
11 | outoebx vc_d 1 0 |
---|
12 | |
---|
13 | set_bcoef uc_d 1 |
---|
14 | set_bcoef vc_d 1 |
---|
15 | |
---|
16 | #set_scoef cout_div 1e7 |
---|
17 | set_scoef cout_div 0 |
---|
18 | |
---|
19 | #Conditions initiales : |
---|
20 | setstate r_float_d 0 |
---|
21 | |
---|
22 | loadstate cfil 0 ij 0 A 1 ./CC.dat D |
---|
23 | loadstate dfil 0 ij 0 A 1 ./DD.dat D |
---|
24 | loadstate u_norm 0 ij 0 A 1 ./NU.dat D |
---|
25 | loadstate v_norm 0 ij 0 A 1 ./NV.dat D |
---|
26 | |
---|
27 | #Observations : |
---|
28 | read_obs_2 ./obs_jum.dat |
---|
29 | read_bck ./uv_bck.dat |
---|
30 | |
---|
31 | read_mask ./mask.dat |
---|
32 | set_modeltime 0 |
---|
33 | print_cost OFF |
---|
34 | |
---|
35 | setm_impres 1 |
---|
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 | multirun |
---|
46 | |
---|
47 | save_output_uv ./uv_total_jum1.dat |
---|
48 | save_output_rfloat ./rfloat_total_jum1.dat |
---|
49 | savestate Div 0 ij 0 A 3 ./Div_ana.dat |
---|
50 | |
---|
51 | goto fin |
---|
52 | |
---|
53 | |
---|
54 | fin |
---|
Note: See
TracBrowser
for help on using the repository browser.