Changeset 76

Timestamp:

03/05/12 10:39:54 (12 years ago)

Author:

aaclod

Message:

55 pas de temps pour sinobad.d, details

Location:

trunk

Files:

• data_out/plot_results.m (modified) (2 diffs)
• scripts/sinobad.i (modified) (3 diffs)
• src/sinobad.d (modified) (1 diff)

trunk/data_out/plot_results.m

@@ -1 @@
-f=netcdf('exp_T_3D_all_ebxhz.nc');
+f=netcdf('exp_T.nc');
 
 addpath('../scripts/matlab_toolbox');
@@ -11 +11 @@
 suff=cell(0,2);
 suff(end+1,:)={'a_c_init','first guess (FG)'};
-suff(end+1,:)={'b_forw0_49','forward du FG'};
+suff(end+1,:)={'b_forw0_52','forward du FG'};
 suff(end+1,:)={'a_c_fin','param control final (AP)'};
-suff(end+1,:)={'b_forwfin_99','forward du PA'};
+suff(end+1,:)={'b_forwfin_52','forward du PA'};
 suff(end+1,:)={'a_c_true','parametre vrai (TP)'};
-suff(end+1,:)={'b_obs_99','forward du TP'};
+suff(end+1,:)={'b_obs_52','forward du TP'};
 
 suff2=cell(0,2);
 suff2(end+1,:)={'n_c_init','first guess (FG)'};
-suff2(end+1,:)={'b_forw0_49','forward du FG'};
+suff2(end+1,:)={'b_forw0_52','forward du FG'};
 suff2(end+1,:)={'n_c_fin','param control final (AP)'};
-suff2(end+1,:)={'b_forwfin_99','forward du PA'};
+suff2(end+1,:)={'b_forwfin_52','forward du PA'};
 suff2(end+1,:)={'n_c_true','parametre vrai (TP)'};
-suff2(end+1,:)={'b_obs_99','forward du TP'};
+suff2(end+1,:)={'b_obs_52','forward du TP'};
 
 

trunk/scripts/sinobad.i

@@ -24 +24 @@
 
 ##OBSERVATION
-xwriteout 50 obs_48 ../data_out/exp_T.nc
+xwriteout 50 obs_52 ../data_out/exp_T.nc
 
 xrst_save ../data_in/file_rest/GYRE_00000420_restart_forwardYAO.nc
 
-goto FINRUN
+#goto FINRUN
 
-savestate sshb  1   ij   50    A       1       ../data_out/sshb_obs_50at   
-savestate tb  1   ijk   50    A       1       ../data_out/tb_obs_50.dat   
-savestate sb  1   ijk   50    A       1       ../data_out/sb_obs_50.dat   
-savestate ub  1   ijk   50    A       1       ../data_out/ub_obs_50.dat   
-savestate vb  1   ijk   50    A       1       ../data_out/vb_obs_50.dat   
+savestate sshb  1   ij   52    A       1       ../data_out/sshb_obs_52.dat   
+savestate tb  1   ijk   52    A       1       ../data_out/tb_obs_52.dat   
+savestate sb  1   ijk   52    A       1       ../data_out/sb_obs_52.dat   
+savestate ub  1   ijk   52    A       1       ../data_out/ub_obs_52.dat   
+savestate vb  1   ijk   52    A       1       ../data_out/vb_obs_52.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   
+loadobs sshb  1   ij   52    A       1       ../data_out/sshb_obs_52.dat D   
+loadobs tb  1   ijk   52    A       1       ../data_out/tb_obs_52.dat D   
+loadobs sb  1   ijk   52    A       1       ../data_out/sb_obs_52.dat D   
+loadobs ub  1   ijk   52    A       1       ../data_out/ub_obs_52.dat D   
+loadobs vb  1   ijk   52    A       1       ../data_out/vb_obs_52.dat D   
 
 
@@ -52 +52 @@
 
 xwriteout 0 init ../data_out/exp_T.nc
-xwriteout 50 forw0_48 ../data_out/exp_T.nc
+xwriteout 52 forw0_52 ../data_out/exp_T.nc
 
 ##TEST DE LA FONCTION OBJECTIVE
@@ -69 +69 @@
 setm_nsim    100
 setm_dxmin   1.0e-4
-setm_epsg    1.0e-10
+setm_epsg    1.0e-3
 setm_ddf1    1
 runm
 ##SAUVEGARDE
 xwriteout 0 fin ../data_out/exp_T.nc
-xwriteout 50 forwfin_48 ../data_out/exp_T.nc
+xwriteout 50 forwfin_52 ../data_out/exp_T.nc
 
 xrst_save ../data_in/file_rest/GYRE_00000400_restart_CONTROL_T.nc

trunk/src/sinobad.d

@@ -20 +20 @@
 defval NZ      31 
 defval NMAX    800 //    !  maximum of iterations for the SOR solver
-defval TA      121 |__DATE__ | number of time steps
+defval TA      55 |__DATE__ | number of time steps
 defval TU       2 | start time step
 #ifdef FILTER