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

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

commit final avec filtre

  • Property svn:executable set to *
File size: 2.9 KB
Line 
1#==============================================================================
2#               fichier de D'INSTRUCTIONS de référence
3#Julien Brajard
4#==============================================================================
5#script faisant l'assimilation de référence
6#paramÚtres à régler : TA=50 (48 pdt effectifs)
7#NX=32,NY=22,NZ=31
8#paramÚres de contrÎle : ta_c
9#obs (3D a t=48 Yt=50) : sshb, tb, sb,ub,vb
10#
11
12##INITIALISATION
13xistate_init 1 ../data_in/file_rest/GYRE_00000424_restart.nc
14load_eb ta_c   // pour tenir en compte ta_c sinon sera 0 cause filtre
15true_target_in_tab ta_c
16xinitnc ../data_out/exp_T.nc
17
18xwriteout 0 true ../data_out/exp_T.nc
19#setstate nu 1
20init_nu
21xset_fdt_fco th 100.
22xset_fdt_fco tz 100.
23xset_fdt_fco kh 5.e+02
24xset_fdt_fco kz 1.2e-05
25
26print_time OFF
27set_modeltime 0
28FORWARD
29#xwriteout 50 obs_48 ../data_out/exp_T.nc
30#savestate gcx_dynspg_flt 1 ijk 0 A 3 ../data_out/gcx_dynspg_flt_jb.dat
31#savestate gcb_dynspg_flt 1 ijk 0 A 3 ../data_out/gcb_dynspg_flt_jb.dat
32#savestate gcx2 1 ijk 0 A 3 ../data_out/gcx2_jb.dat
33
34#goto FINRUN
35
36##OBSERVATION
37xwriteout 50 obs_48 ../data_out/exp_T.nc
38savestate sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat   
39savestate tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat   
40savestate sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat   
41savestate ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat   
42savestate vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat   
43GOTO TESTAD
44
45loadobs sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat D   
46loadobs tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat D   
47loadobs sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat D   
48loadobs ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat D   
49loadobs vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat D   
50
51GOTO FINTEST
52#TEST DE L'ADJOINT'
53TESTAD
54testad  1.e-01
55stop
56#-----------------
57#TEST DE LA PROGRAMMATION DES DERIVEES
58TESTDF
59testdf 3  3  1  4    f 1.e-10  0.001  10
60stop
61#-----------------
62# TEST DE LINEAIRE TANGENT
63TESTLT
64 testlt 1.e+00 0.5e-01 2. 15
65stop
66FINTEST
67#goto FINRUN
68
69##INITIALISATION DE L'ASSIMILATION'
70xistate_init 1 ../data_in/file_rest/GYRE_00000400_restart.nc
71load_eb ta_c  // charger l'ebauche '
72setstate nu 0
73outoebx nu 1 0
74set_bcoef nu 1.e-16
75set_modeltime 0
76print_time OFF
77FORWARD
78
79xwriteout 0 init ../data_out/exp_T.nc
80xwriteout 50 forw0_48 ../data_out/exp_T.nc
81
82##TEST DE LA FONCTION OBJECTIVE
83goto FINTEST
84testof 1.e+00 0.5e-01 2. 10 1
85goto FINRUN
86FINTEST
87##DEBUT DE L'ASSIMILATION'
88print_cost ON
89set_modeltime 0
90!echo "RUN OPTIMIZATION WITH M1QN3 ......."
91setm_impres  5
92setm_io      6
93setm_mode    1
94set_nbiter   100
95setm_nsim    100
96setm_dxmin   1.0e-4
97setm_epsg    1.0e-10
98setm_ddf1    1
99runm
100##SAUVEGARDE
101xwriteout 0 fin ../data_out/exp_T.nc
102xwriteout 50 forwfin_48 ../data_out/exp_T.nc
103
104xrst_save ../data_in/file_rest/GYRE_00000400_restart_CONTROL_T.nc
105
106
107
108goto FINRUN
109
110
111FINRUN
112
Note: See TracBrowser for help on using the repository browser.