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 | #wind coefficient |
---|
28 | #Poulain, JAOT, 2009 |
---|
29 | setstate au 0 |
---|
30 | setstate av 0 |
---|
31 | #setstate bu 0.0061 |
---|
32 | #setstate bv -0.0031 |
---|
33 | setstate bu 0 |
---|
34 | setstate bv 0 |
---|
35 | load_wind wind_interp.nc |
---|
36 | |
---|
37 | #Observations : |
---|
38 | read_obs_2 ./obs_jum.dat |
---|
39 | read_bck ./uv_bck.dat |
---|
40 | |
---|
41 | read_mask ./mask.dat |
---|
42 | set_modeltime 0 |
---|
43 | print_cost OFF |
---|
44 | |
---|
45 | setm_impres 1 |
---|
46 | setm_io 6 |
---|
47 | setm_mode 0 |
---|
48 | #set_nbiter 100 |
---|
49 | set_nbiter 20 |
---|
50 | setm_nsim 200 |
---|
51 | setm_dxmin 1.0e-12 |
---|
52 | setm_epsg 1.0e-12 |
---|
53 | setm_ddf1 1.0 |
---|
54 | |
---|
55 | multirun |
---|
56 | |
---|
57 | save_output_uv ./uv_total_jum1.dat |
---|
58 | save_output_rfloat ./rfloat_total_jum1.dat |
---|
59 | savestate Div 0 ij 0 A 3 ./Div_ana.dat |
---|
60 | |
---|
61 | goto fin |
---|
62 | |
---|
63 | |
---|
64 | fin |
---|
Note: See
TracBrowser
for help on using the repository browser.