source: altifloat/exp_twin_aviso/assim.i @ 147

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 
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#Observations :
28read_obs_2 ./obs_jum.dat
29read_bck ./uv_bck.dat
30
31read_mask ./mask.dat
32set_modeltime 0
33print_cost OFF
34
35setm_impres 1
36setm_io 6
37setm_mode 0
38#set_nbiter 100
39set_nbiter 20
40setm_nsim 200
41setm_dxmin 1.0e-12
42setm_epsg 1.0e-12
43setm_ddf1 1.0
44
45multirun
46
47save_output_uv ./uv_total_jum1.dat
48save_output_rfloat ./rfloat_total_jum1.dat
49savestate Div 0 ij 0 A 3 ./Div_ana.dat
50
51goto fin
52
53
54fin
Note: See TracBrowser for help on using the repository browser.