source: branches/branche-mb/scripts/sinobad_2D.i @ 68

Last change on this file since 68 was 68, checked in by berrada, 13 years ago

commit final avec filtre

File size: 2.2 KB
Line 
1#==============================================================================
2#               fichier de D'INSTRUCTIONS de référence
3#Julien Brajard
4# modified by mohamed Berrada (filtre August 2011)
5#==============================================================================
6#script faisant l'assimilation de référence
7#paramÚtres à régler : TA=50 (48 pdt effectifs)
8#NX=32,NY=22,NZ=31
9#paramÚres de contrÎle : ta_c
10#obs (2D a t=48 Yt=50) : tb
11#
12
13##INITIALISATION
14xistate_init 1 ../data_in/file_rest/GYRE_00000424_restart.nc
15load_eb ta_c   // pour tenir en compte ta_c sinon sera 0 cause filtre
16true_target_in_tab ta_c
17
18xset_fdt_fco th 1000.
19xset_fdt_fco tz 10.
20xset_fdt_fco kh 5.e+05
21xset_fdt_fco kz 1.2e-03
22
23xinitnc ../data_out/exp_T_2D.nc
24
25xwriteout 0 true ../data_out/exp_T_2D.nc
26
27print_time OFF
28set_modeltime 0
29FORWARD
30
31
32##OBSERVATION
33xwriteout 50 obs_48 ../data_out/exp_T_2D.nc
34
35#Pour ne sauvegarder que la surface
36xchangesavemode 1
37
38savestate tb  1   ijk   50    A       1       ../data_out/tb_obs_48_2D.dat   
39
40loadobs tb  1   ijk   50    A       1       ../data_out/tb_obs_48_2D.dat D   
41
42#goto FINRUN
43
44##INITIALISATION DE L'ASSIMILATION'
45xistate_init 1 ../data_in/file_rest/GYRE_00000400_restart.nc
46load_eb ta_c  // charger l'ebauche '
47setstate nu 0
48outoebx nu 1 0
49set_bcoef nu 1.e-8
50set_modeltime 0
51print_time OFF
52FORWARD
53
54xwriteout 0 init ../data_out/exp_T_2D.nc
55xwriteout 50 forw0_48 ../data_out/exp_T_2D.nc
56
57##TEST DE LA FONCTION OBJECTIVE
58goto FINTEST
59testof 1.e+00 0.5e-01 2. 10 1
60goto FINRUN
61FINTEST
62##DEBUT DE L'ASSIMILATION'
63print_cost ON
64set_modeltime 0
65!echo "RUN OPTIMIZATION WITH M1QN3 ......."
66setm_impres  5
67setm_io      6
68setm_mode    1
69set_nbiter   100
70setm_nsim    100
71setm_dxmin   1.0e-4
72setm_epsg    1.0e-10
73setm_ddf1    1
74runm
75
76goto nomulti
77load_eb ta_c  // charger l'ebauche '
78set_bcoef nu 5.e-05
79runm
80
81load_eb ta_c  // charger l'ebauche '
82set_bcoef nu 1.e-04
83runm
84
85##SAUVEGARDE
86xwriteout 0 fin ../data_out/exp_T_2D.nc
87xwriteout 50 forwfin_48 ../data_out/exp_T_2D.nc
88
89xrst_save ../data_in/file_rest/GYRE_00000400_restart_CONTROL_T_2D.nc
90
91
92nomulti
93
94##SAUVEGARDE
95xwriteout 0 fin ../data_out/exp_T_2D.nc
96xwriteout 50 forwfin_48 ../data_out/exp_T_2D.nc
97
98xrst_save ../data_in/file_rest/GYRE_00000400_restart_CONTROL_T_2D.nc
99
100
101
102goto FINRUN
103
104
105FINRUN
106
Note: See TracBrowser for help on using the repository browser.