Line | |
---|
1 | |
---|
2 | setstate u 0 |
---|
3 | setstate v 0 |
---|
4 | setstate u_d 0 |
---|
5 | setstate v_d 0 |
---|
6 | setstate u_norm 1 |
---|
7 | setstate v_norm 1 |
---|
8 | setstate uc_d 0 |
---|
9 | setstate vc_d 0 |
---|
10 | |
---|
11 | outoebx uc_d 1 0 |
---|
12 | outoebx vc_d 1 0 |
---|
13 | |
---|
14 | set_bcoef uc_d 2 |
---|
15 | set_bcoef vc_d 2 |
---|
16 | |
---|
17 | set_scoef cout_div 8e7 |
---|
18 | #set_scoef cout_div 0 |
---|
19 | |
---|
20 | #wind coefficient |
---|
21 | #Poulain, JAOT, 2009 |
---|
22 | setstate au 0 |
---|
23 | setstate av 0 |
---|
24 | #setstate bu 0.0061 |
---|
25 | #setstate bv -0.0031 |
---|
26 | setstate bu 0 |
---|
27 | setstate bv 0 |
---|
28 | load_wind wind_interp.nc |
---|
29 | |
---|
30 | #Conditions initiales : |
---|
31 | setstate r_float_d 0 |
---|
32 | |
---|
33 | loadstate cfil 0 ij 0 A 1 ./CC.dat D |
---|
34 | loadstate dfil 0 ij 0 A 1 ./DD.dat D |
---|
35 | loadstate u_norm 0 ij 0 A 1 ./NU.dat D |
---|
36 | loadstate v_norm 0 ij 0 A 1 ./NV.dat D |
---|
37 | |
---|
38 | #Observations : |
---|
39 | read_obs_2 ./obs.dat |
---|
40 | read_bck ./uv_back_aviso.dat |
---|
41 | read_mask ./mask_aviso.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_geos_nodiv.dat |
---|
58 | save_output_uvtot ./uv_tot_nodiv.dat |
---|
59 | save_output_rfloat ./rfloat_nodiv.dat |
---|
60 | savestate Div 0 ij 0 A 3 ./Div_ana_nodiv.dat |
---|
61 | |
---|
62 | goto fin |
---|
63 | |
---|
64 | |
---|
65 | fin |
---|
Note: See
TracBrowser
for help on using the repository browser.