source: altifloat/exp_twin_aviso/assim.i @ 199

Last change on this file since 199 was 160, checked in by jbrlod, 9 years ago

finalize wind effect

File size: 994 bytes
Line 
1setstate u 0
2setstate v 0
3setstate u_d 0
4setstate v_d 0
5setstate u_norm 1
6setstate v_norm 1
7setstate uc_d 0
8setstate vc_d 0
9
10outoebx uc_d 1 0
11outoebx vc_d 1 0
12
13set_bcoef uc_d 1
14set_bcoef vc_d 1
15
16#set_scoef cout_div 1e7
17set_scoef cout_div 0
18
19#Conditions initiales :
20setstate r_float_d 0
21
22loadstate cfil 0 ij 0 A 1 ./CC.dat D
23loadstate dfil 0 ij 0 A 1 ./DD.dat D
24loadstate u_norm 0 ij 0 A 1 ./NU.dat D
25loadstate v_norm 0 ij 0 A 1 ./NV.dat D
26
27#wind coefficient
28#Poulain, JAOT, 2009
29setstate au 0
30setstate av 0
31#setstate bu 0.0061
32#setstate bv -0.0031
33setstate bu 0
34setstate bv 0
35load_wind wind_interp.nc
36
37#Observations :
38read_obs_2 ./obs_jum.dat
39read_bck ./uv_bck.dat
40
41read_mask ./mask.dat
42set_modeltime 0
43print_cost OFF
44
45setm_impres 1
46setm_io 6
47setm_mode 0
48#set_nbiter 100
49set_nbiter 20
50setm_nsim 200
51setm_dxmin 1.0e-12
52setm_epsg 1.0e-12
53setm_ddf1 1.0
54
55multirun
56
57save_output_uv ./uv_total_jum1.dat
58save_output_rfloat ./rfloat_total_jum1.dat
59savestate Div 0 ij 0 A 3 ./Div_ana.dat
60
61goto fin
62
63
64fin
Note: See TracBrowser for help on using the repository browser.