Changeset 9 for trunk/scripts

Timestamp:

03/14/11 19:31:32 (13 years ago)

Author:

jbrlod

Message:

begin ref run

File:

• trunk/scripts/sinobad.i (modified) (2 diffs)

trunk/scripts/sinobad.i

@@ -6 +6 @@
 #paramÚtres à régler : TA=50 (48 pdt effectifs)
 #NX=32,NY=22,NZ=31
-#
+#paramÚres de contrÎle : sshn_c, ta_c, sa_c, ua_c, va_c
+#obs (3D a t=48 Yt=50) : sshb, tb, sb,ub,vb
 
 ##INITIALISATION
+xistate_init 1 ../data_in/file_rest/GYRE_00000424_restart.nc
+true_target_in_tab ta_c
+
+
+savestate sshn_c  1   ij   0    A       1       ../data_out/sshn_c_true.dat
+savestate ta_c  1   ijk   0    A       1       ../data_out/ta_c_true.dat
+savestate sa_c  1   ijk   0    A       1       ../data_out/sa_c_true.dat
+savestate ua_c  1   ijk   0    A       1       ../data_out/ua_c_true.dat
+savestate va_c  1   ijk   0    A       1       ../data_out/va_c_true.dat
+
+print_time OFF
+set_modeltime 0
+FORWARD
+
 
 ##OBSERVATION
-
-##INITIALISATION DE L'ASSIMILATION
-
-##DEBUT DE L'ASSIMILATION
-
+savestate sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat   
+savestate tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat   
+savestate sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat   
+savestate ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat   
+savestate vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat   
+
+loadobs sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat D   
+loadobs tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat D   
+loadobs sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat D   
+loadobs ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat D   
+loadobs vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat D   
+
+
+#goto FINRUN
+
+##INITIALISATION DE L'ASSIMILATION'
+xistate_init 1 ../data_in/file_rest/GYRE_00000400_restart.nc
+set_modeltime 0
+print_time OFF
+FORWARD
+savestate sshn_c  1   ij   0    A       1       ../data_out/sshn_c_init.dat
+savestate ta_c  1   ijk   0    A       1       ../data_out/ta_c_init.dat
+savestate sa_c  1   ijk   0    A       1       ../data_out/sa_c_init.dat
+savestate ua_c  1   ijk   0    A       1       ../data_out/ua_c_init.dat
+savestate va_c  1   ijk   0    A       1       ../data_out/va_c_init.dat
+
+##DEBUT DE L'ASSIMILATION'
+
+set_modeltime 0
+!echo "RUN OPTIMIZATION WITH M1QN3 ......."
+setm_impres  5
+setm_io      6
+setm_mode    1
+set_nbiter   100
+setm_nsim    100
+setm_dxmin   1.0e-4
+setm_epsg    1.0e-10
+setm_ddf1    1
+runm
 ##SAUVEGARDE
+savestate sshn_c  1   ij   0    A       1       ../data_out/sshn_c_fin.dat
+savestate ta_c  1   ijk   0    A       1       ../data_out/ta_c_fin.dat
+savestate sa_c  1   ijk   0    A       1       ../data_out/sa_c_fin.dat
+savestate ua_c  1   ijk   0    A       1       ../data_out/ua_c_fin.dat
+savestate va_c  1   ijk   0    A       1       ../data_out/va_c_fin.dat
+
+goto FINRUN
 
 ## analytiquement (2) 
@@ -44 +100 @@
 
 #setstate pca_ta 0
-GOTO NOF
+GOTO NOF 
 ## Faire une simulation direct 
 print_time OFF