Index: /trunk/scripts/run_batch_sinobad.sh
===================================================================
--- /trunk/scripts/run_batch_sinobad.sh	(revision 40)
+++ /trunk/scripts/run_batch_sinobad.sh	(revision 41)
@@ -1,3 +1,3 @@
-#!/bin/tcsh
+#!/bin/sh
 #PBS -e ../scripts/script_outputs/
 #PBS -o ../scripts/script_outputs/
@@ -8,7 +8,15 @@
 #Command to submit the job:
 #qsub run_batch_sinobad.sh
-cd /usr/home/jbrlod/these/postdoc/collaborations/NEMO/nemoyao/svn/trunk/src/
+cd $PBS_O_WORKDIR
+cd ../src/
+
+echo "****************************************" 
+echo "*       run reference experiment       *"
+echo "*             (sequential)             *"
+echo "*              26/04/2011              *"
+echo "*  julien.brajard@locean-ipsl.upmc.fr  *"
+echo "****************************************"
 
 $YAODIR/etc/YaoI -p +O3 sinobad
-./sinobad -i ../scripts/sinobad.i
+time ./sinobad -i ../scripts/sinobad.i
 
Index: /trunk/scripts/run_batch_sinobad_seq.sh
===================================================================
--- /trunk/scripts/run_batch_sinobad_seq.sh	(revision 41)
+++ /trunk/scripts/run_batch_sinobad_seq.sh	(revision 41)
@@ -0,0 +1,23 @@
+#!/bin/sh
+#PBS -e ../scripts/script_outputs/
+#PBS -o ../scripts/script_outputs/
+#PBS -q short
+#PBS -l ncpus=16,nodes=1:x64
+
+#run reference experience
+#Command to submit the job:
+#qsub run_batch_sinobad.sh
+cd $PBS_O_WORKDIR
+cd ../src/
+
+
+echo "****************************************" 
+echo "*       run reference experiment       *"
+echo "*             (sequential)             *"
+echo "*              26/04/2011              *"
+echo "*  julien.brajard@locean-ipsl.upmc.fr  *"
+echo "****************************************"
+
+$YAODIR/etc/YaoI +O3 sinobad
+time ./sinobad -i ../scripts/sinobad.i
+
Index: unk/scripts/script_outputs/run_batch_sinobad.sh.o321
===================================================================
--- /trunk/scripts/script_outputs/run_batch_sinobad.sh.o321	(revision 40)
+++ 	(revision )
@@ -1,130 +1,0 @@
-Warning: no access to tty (Bad file descriptor).
-Thus no job control in this shell.
-appending ferret path
-appending ferret path
- CSHRC: BINTYPE linux86_64 ...
- --> test d'imprimante :  hpl5.
-
-***>> DISPLAY n'est pas défini. Pour le définir tapez d<NOMDUTERMINAL>
-        <NOMDUTERMINAL> correspond à l'écran où vous voulez faire afficher
-        vos fenêtres. Par exemple :
-
-  dboreas pour la Station Solaris boreas.lodyc.jussieu.fr:0 134.157.9.52:0
-  ddedale pour la Station Linux dedale.lodyc.jussieu.fr:0   134.157.9.123:0
-  calypso pour le Terminal X calypso.lodyc.jussieu.fr:0     134.157.9.3:0
-
-                 Bienvenue au LODYC
-    Actuellement vous êtes logué sur zelos.locean-ipsl.upmc.fr
-
- 09:47:25 up 196 days, 21:47,  0 users,  load average: 0,00, 0,01, 0,00
-USER     TTY        LOGIN@   IDLE   JCPU   PCPU WHAT
-
-TERM ... vt100
-TTY .... n'est pas un tty
- > MACHINE TYPE determined now ... <linux>.
- > MQNDIR directory variable not declared...standard directory.
-
-       ====================================================================
-	            YaoI V.9 : PROJECT CREATION PROCEDURE                   
-	            01.11.2008 (release 01.11.2008)                         
-       ====================================================================
-        LOCAL MACHINE BINARY TYPE FOUND ..... <linux>                      
-       ====================================================================
-
- > Run of YaoI for project : sinobad
- > HEADERS generation not required: 
- >   the description file has not changed from the last generation of Y1sinobad.h and Y2sinobad.h.
- >   Pay attention all the warning showed in the generation process are not been displayed this way because of the non regeneration.
- >   If you want to force the regeneration remove the Yworkdir/ directory.
- > enter proc construction 
- Extern object : 
-                 ../object/ncutil.o 
- > end proc construction 
- > start compilation 
- > end of compilation, start linking 
-
-       =======================  PROJECT IS NOW READY  =====================
-                           YaoI V.9 : END OF PROCEDURE                     
-                           01.11.2008 (release 01.11.2008)                 
-       ====================================================================
-                           TOTAL COMPILATION TIME : 14 second(s)            
-       ====================================================================
-
-
-Starting Yao : YaoI V.9 01.11.2008 (release 01.11.2008)  -O3  
-Project sinobad. 
-(You are in double mode)
-This is a parallel version of the application, the max number of threads that the system supports are: 24.
-RUN OPTIMIZATION WITH M1QN3 .......
-////////////////////////////////////////////////////////////////////////
-//                       NEMO/YAO PROJECT                             //
-//                     M. Berrada 02-2009                            //
-//                  LOCEAN-IPSL.UPMC (Paris 6)                        //
-//====================================================================//
-norm_sup_diff_ta_c :  0.0000000000000000e+00
-norm_2_diff_ta_c :  0.0000000000000000e+00
-Yao#>
-Yao#>xistate_init 1 ../data_in/file_rest/GYRE_00000424_restart.nc
-start from rest_file
- Yao#>true_target_in_tab ta_c
-Yao#>xinitnc ../data_out/exp_T.nc
-Yao#>
-Yao#>xwriteout 0 true ../data_out/exp_T.nc
-Yao#>
-Yao#>print_time OFF
-Yao#>set_modeltime 0
-Yao#>FORWARD
-Yao#>
-Yao#>
-Yao#>xwriteout 50 obs_48 ../data_out/exp_T.nc
-Yao#>savestate sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat   
- 704 state(s) selected saving module sshb
-Yao#>savestate tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat   
- 21824 state(s) selected saving module tb
-Yao#>savestate sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat   
- 21824 state(s) selected saving module sb
-Yao#>savestate ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat   
- 21824 state(s) selected saving module ub
-Yao#>savestate vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat   
- 21824 state(s) selected saving module vb
-Yao#>
-Yao#>loadobs sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat D   
-modul sshb: 704 data read ; 704 state(s) selected as obs (704 inserted + 0 overrided)
-Yao#>loadobs tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat D   
-modul tb: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
-Yao#>loadobs sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat D   
-modul sb: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
-Yao#>loadobs ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat D   
-modul ub: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
-Yao#>loadobs vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat D   
-modul vb: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
-Yao#>
-Yao#>
-Yao#>
-Yao#>xistate_init 1 ../data_in/file_rest/GYRE_00000400_restart.nc
-start from rest_file
- Yao#>set_modeltime 0
-Yao#>print_time OFF
-Yao#>FORWARD
-Yao#>
-Yao#>xwriteout 0 init ../data_out/exp_T.nc
-Yao#>xwriteout 50 forw0_48 ../data_out/exp_T.nc
-Yao#>
-Yao#>goto FINTEST
-Yao#>
-Yao#>set_modeltime 0
-Yao#>setm_impres  5
-Yao#>setm_io      6
-Yao#>setm_mode    1
-Yao#>set_nbiter   100
-Yao#>setm_nsim    100
-Yao#>setm_dxmin   1.0e-4
-Yao#>setm_epsg    1.0e-10
-Yao#>setm_ddf1    1
-Yao#>runm
-
-start run time: Tue Apr 26 09:47:47 2011
-
-norm_sup_diff_ta_c :  5.6939927493938196e+00
-norm_2_diff_ta_c :  2.9456334868054643e+01
-Terminated
Index: unk/scripts/script_outputs/run_batch_sinobad.sh.o322
===================================================================
--- /trunk/scripts/script_outputs/run_batch_sinobad.sh.o322	(revision 40)
+++ 	(revision )
@@ -1,389 +1,0 @@
-Warning: no access to tty (Bad file descriptor).
-Thus no job control in this shell.
-appending ferret path
-appending ferret path
- CSHRC: BINTYPE linux86_64 ...
- --> test d'imprimante :  hpl5.
-
-***>> DISPLAY n'est pas défini. Pour le définir tapez d<NOMDUTERMINAL>
-        <NOMDUTERMINAL> correspond à l'écran où vous voulez faire afficher
-        vos fenêtres. Par exemple :
-
-  dboreas pour la Station Solaris boreas.lodyc.jussieu.fr:0 134.157.9.52:0
-  ddedale pour la Station Linux dedale.lodyc.jussieu.fr:0   134.157.9.123:0
-  calypso pour le Terminal X calypso.lodyc.jussieu.fr:0     134.157.9.3:0
-
-                 Bienvenue au LODYC
-    Actuellement vous êtes logué sur zelos.locean-ipsl.upmc.fr
-
- 09:50:09 up 196 days, 21:50,  0 users,  load average: 1,83, 1,71, 0,72
-USER     TTY        LOGIN@   IDLE   JCPU   PCPU WHAT
-
-TERM ... vt100
-TTY .... n'est pas un tty
- > MACHINE TYPE determined now ... <linux>.
- > MQNDIR directory variable not declared...standard directory.
-
-       ====================================================================
-	            YaoI V.9 : PROJECT CREATION PROCEDURE                   
-	            01.11.2008 (release 01.11.2008)                         
-       ====================================================================
-        LOCAL MACHINE BINARY TYPE FOUND ..... <linux>                      
-       ====================================================================
-
- > Run of YaoI for project : sinobad
- > HEADERS generation not required: 
- >   the description file has not changed from the last generation of Y1sinobad.h and Y2sinobad.h.
- >   Pay attention all the warning showed in the generation process are not been displayed this way because of the non regeneration.
- >   If you want to force the regeneration remove the Yworkdir/ directory.
- > enter proc construction 
- Extern object : 
-                 ../object/ncutil.o 
- > end proc construction 
- > start compilation 
- > end of compilation, start linking 
-
-       =======================  PROJECT IS NOW READY  =====================
-                           YaoI V.9 : END OF PROCEDURE                     
-                           01.11.2008 (release 01.11.2008)                 
-       ====================================================================
-                           TOTAL COMPILATION TIME : 13 second(s)            
-       ====================================================================
-
-
-Starting Yao : YaoI V.9 01.11.2008 (release 01.11.2008)  -O3  
-Project sinobad. 
-(You are in double mode)
-This is a parallel version of the application, the max number of threads that the system supports are: 24.
-RUN OPTIMIZATION WITH M1QN3 .......
-////////////////////////////////////////////////////////////////////////
-//                       NEMO/YAO PROJECT                             //
-//                     M. Berrada 02-2009                            //
-//                  LOCEAN-IPSL.UPMC (Paris 6)                        //
-//====================================================================//
-norm_sup_diff_ta_c :  0.0000000000000000e+00
-norm_2_diff_ta_c :  0.0000000000000000e+00
-Yao#>
-Yao#>xistate_init 1 ../data_in/file_rest/GYRE_00000424_restart.nc
-start from rest_file
- Yao#>true_target_in_tab ta_c
-Yao#>xinitnc ../data_out/exp_T.nc
-Yao#>
-Yao#>xwriteout 0 true ../data_out/exp_T.nc
-Yao#>
-Yao#>print_time OFF
-Yao#>set_modeltime 0
-Yao#>FORWARD
-Yao#>
-Yao#>
-Yao#>xwriteout 50 obs_48 ../data_out/exp_T.nc
-Yao#>savestate sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat   
- 704 state(s) selected saving module sshb
-Yao#>savestate tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat   
- 21824 state(s) selected saving module tb
-Yao#>savestate sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat   
- 21824 state(s) selected saving module sb
-Yao#>savestate ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat   
- 21824 state(s) selected saving module ub
-Yao#>savestate vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat   
- 21824 state(s) selected saving module vb
-Yao#>
-Yao#>loadobs sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat D   
-modul sshb: 704 data read ; 704 state(s) selected as obs (704 inserted + 0 overrided)
-Yao#>loadobs tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat D   
-modul tb: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
-Yao#>loadobs sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat D   
-modul sb: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
-Yao#>loadobs ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat D   
-modul ub: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
-Yao#>loadobs vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat D   
-modul vb: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
-Yao#>
-Yao#>
-Yao#>
-Yao#>xistate_init 1 ../data_in/file_rest/GYRE_00000400_restart.nc
-start from rest_file
- Yao#>set_modeltime 0
-Yao#>print_time OFF
-Yao#>FORWARD
-Yao#>
-Yao#>xwriteout 0 init ../data_out/exp_T.nc
-Yao#>xwriteout 50 forw0_48 ../data_out/exp_T.nc
-Yao#>
-Yao#>goto FINTEST
-Yao#>
-Yao#>set_modeltime 0
-Yao#>setm_impres  5
-Yao#>setm_io      6
-Yao#>setm_mode    1
-Yao#>set_nbiter   100
-Yao#>setm_nsim    100
-Yao#>setm_dxmin   1.0e-4
-Yao#>setm_epsg    1.0e-10
-Yao#>setm_ddf1    1
-Yao#>runm
-
-start run time: Tue Apr 26 09:50:30 2011
-
-norm_sup_diff_ta_c :  5.6928683673747766e+00
-norm_2_diff_ta_c :  2.9456930296697998e+01
-norm_sup_diff_ta_c :  6.4456128245485527e+00
-norm_2_diff_ta_c :  2.5358057811903969e+01
-norm_sup_diff_ta_c :  6.5684672748454798e+00
-norm_2_diff_ta_c :  2.5910609236589693e+01
-norm_sup_diff_ta_c :  6.4682904636272180e+00
-norm_2_diff_ta_c :  2.5272654891909998e+01
-norm_sup_diff_ta_c :  6.4241515552409876e+00
-norm_2_diff_ta_c :  2.5273839573783327e+01
-norm_sup_diff_ta_c :  6.4976221477397669e+00
-norm_2_diff_ta_c :  2.6027933104251819e+01
-norm_sup_diff_ta_c :  6.4420472537883509e+00
-norm_2_diff_ta_c :  2.5352006728036727e+01
-norm_sup_diff_ta_c :  6.4418684398540247e+00
-norm_2_diff_ta_c :  2.5350883241422451e+01
-norm_sup_diff_ta_c :  6.4409824764140833e+00
-norm_2_diff_ta_c :  2.5345420593959393e+01
-norm_sup_diff_ta_c :  6.6385418331035364e+00
-norm_2_diff_ta_c :  2.5933307557810213e+01
-norm_sup_diff_ta_c :  6.4823833858379114e+00
-norm_2_diff_ta_c :  2.5383526478692996e+01
-norm_sup_diff_ta_c :  6.4819694911845911e+00
-norm_2_diff_ta_c :  2.5382919579971695e+01
-norm_sup_diff_ta_c :  6.4799200450786341e+00
-norm_2_diff_ta_c :  2.5379978851884822e+01
-norm_sup_diff_ta_c :  6.4701858913310755e+00
-norm_2_diff_ta_c :  2.5367544567196965e+01
-norm_sup_diff_ta_c :  6.7977506076701868e+00
-norm_2_diff_ta_c :  2.6472100059374544e+01
-norm_sup_diff_ta_c :  6.5451647197612512e+00
-norm_2_diff_ta_c :  2.5472124089082534e+01
-norm_sup_diff_ta_c :  6.5444146549113977e+00
-norm_2_diff_ta_c :  2.5470631271124056e+01
-norm_sup_diff_ta_c :  6.5407034313091028e+00
-norm_2_diff_ta_c :  2.5463373364657119e+01
-norm_sup_diff_ta_c :  6.5230738078960169e+00
-norm_2_diff_ta_c :  2.5431910965081048e+01
-norm_sup_diff_ta_c :  6.5072075282939661e+00
-norm_2_diff_ta_c :  2.5407864869996192e+01
-norm_sup_diff_ta_c :  6.4961010372050989e+00
-norm_2_diff_ta_c :  2.5393444667264557e+01
-norm_sup_diff_ta_c :  6.4883266841777552e+00
-norm_2_diff_ta_c :  2.5384534932549929e+01
-norm_sup_diff_ta_c :  6.4828845416911829e+00
-norm_2_diff_ta_c :  2.5378879096119331e+01
-norm_sup_diff_ta_c :  6.4845167552837122e+00
-norm_2_diff_ta_c :  2.5380524855728378e+01
-norm_sup_diff_ta_c :  6.4871837055095423e+00
-norm_2_diff_ta_c :  2.5383307339528361e+01
-norm_sup_diff_ta_c :  6.4863835727580774e+00
-norm_2_diff_ta_c :  2.5382460911736413e+01
-norm_sup_diff_ta_c :  6.4858237659343470e+00
-norm_2_diff_ta_c :  2.5381874029423425e+01
-norm_sup_diff_ta_c :  6.4854318057903040e+00
-norm_2_diff_ta_c :  2.5381466493927661e+01
-norm_sup_diff_ta_c :  6.4851571475871790e+00
-norm_2_diff_ta_c :  2.5381183594325236e+01
-norm_sup_diff_ta_c :  6.4849649822124231e+00
-norm_2_diff_ta_c :  2.5380984144114944e+01
-norm_sup_diff_ta_c :  6.4850226795085657e+00
-norm_2_diff_ta_c :  2.5381043802479155e+01
-norm_sup_diff_ta_c :  6.4851170932658899e+00
-norm_2_diff_ta_c :  2.5381141000079062e+01
-norm_sup_diff_ta_c :  6.4850884830363977e+00
-norm_2_diff_ta_c :  2.5381112228820555e+01
-norm_sup_diff_ta_c :  6.4850970661052454e+00
-norm_2_diff_ta_c :  2.5381120227344230e+01
-Yao#>(6) Arret sur dxmin pendant le line-search 
-
-start run time ....: Tue Apr 26 09:50:30 2011
-  end run time ....: Tue Apr 26 10:01:31 2011
-  run duration ....: 661s
-
-Yao#>xwriteout 0 fin ../data_out/exp_T.nc
-Yao#>xwriteout 50 forwfin_48 ../data_out/exp_T.nc
-Yao#>
-Yao#>xrst_save ../data_in/file_rest/GYRE_00000400_restart_CONTROL_T.nc
-Yao#>
-Yao#>
-Yao#>
-Yao#>goto FINRUN
-Yao#>
-
- m1qn3 (version 2.0b, december 1993): entry point
-     dimension of the problem (n):    21824
-     absolute precision on x (dxmin): 0.10E-03
-     expected decrease for f (df1): 0.30E+04
-     relative precision on g (epsg): 0.10E-09
-     maximal number of iterations (niter):   100
-     maximal number of simulations (nsim):   100
-     printing level (impres):   5
-
- m1qn3: scalar initial scaling mode
-
-     allocated memory (nrz) :   436488
-     used memory :              414664
-     number of updates :             8
-     (y,s) pairs are stored in core memory
-
- m1qn3: cold start
-
-     f         =  0.29733594E+04
-     norm of g =  0.11778759E+03
-
- m1qn3a: descent direction -g: precon =  0.429E+00
-
- -------------------------------------------------------------------------------
-
- m1qn3: iter   1, simul   1, f= 0.29733594E+04, h'(0)=-0.59467D+04
-
- m1qn3: line search
-
-     mlis0       fpn=-0.595E+04 d2= 0.25E+04  tmin= 0.96E-05 tmax= 0.10E+21
-     mlis3     0.100E+01 -0.243E+04 -0.421E+03
-
- m1qn3: matrix update:
-     oren-spedicato factor =  0.411E+00
-
- m1qn3: stopping criterion on g:  0.33231E+00
-
- m1qn3: descent direction d: angle(-g,d) =  17.7 degrees
-
- -------------------------------------------------------------------------------
-
- m1qn3: iter   2, simul   2, f= 0.54576740E+03, h'(0)=-0.72852D+03
-
- m1qn3: line search
-
-     mlis0       fpn=-0.729E+03 d2= 0.38E+03  tmin= 0.37E-04 tmax= 0.10E+21
-     mlis3     0.100E+01 -0.301E+03  0.723E+02
-
- m1qn3: convergence rate, s(k)/s(k-1) =  0.38691E+00
-
- m1qn3: matrix update:
-     oren-spedicato factor =  0.346E+00
-
- m1qn3: stopping criterion on g:  0.20463E+00
-
- m1qn3: descent direction d: angle(-g,d) =  25.3 degrees
-
- -------------------------------------------------------------------------------
-
- m1qn3: iter   3, simul   3, f= 0.24504416E+03, h'(0)=-0.17172D+03
-
- m1qn3: line search
-
-     mlis0       fpn=-0.172E+03 d2= 0.62E+02  tmin= 0.10E-03 tmax= 0.10E+21
-     mlis2                                      0.100E+01  0.178E+02  0.774E+02
-     mlis3     0.339E+00 -0.301E+02 -0.810E+02
-
- m1qn3: convergence rate, s(k)/s(k-1) =  0.13676E+00
-
- m1qn3: matrix update:
-     oren-spedicato factor =  0.126E+00
-
- m1qn3: stopping criterion on g:  0.12375E+00
-
- m1qn3: descent direction d: angle(-g,d) =  21.1 degrees
-
- -------------------------------------------------------------------------------
-
- m1qn3: iter   4, simul   5, f= 0.21490956E+03, h'(0)=-0.64216D+02
-
- m1qn3: line search
-
-     mlis0       fpn=-0.642E+02 d2= 0.22E+02  tmin= 0.14E-03 tmax= 0.10E+21
-     mlis3     0.100E+01 -0.126E+02 -0.176E+02
-
- m1qn3: convergence rate, s(k)/s(k-1) =  0.17682E+01
-
- m1qn3: matrix update:
-     oren-spedicato factor =  0.216E+00
-
- m1qn3: stopping criterion on g:  0.10386E+00
-
- m1qn3: descent direction d: angle(-g,d) =  42.7 degrees
-
- -------------------------------------------------------------------------------
-
- m1qn3: iter   5, simul   6, f= 0.20227124E+03, h'(0)=-0.48986D+02
-
- m1qn3: line search
-
-     mlis0       fpn=-0.490E+02 d2= 0.30E+02  tmin= 0.95E-04 tmax= 0.10E+21
-     mlis2                                      0.100E+01  0.150E+02  0.215E+02
-     mlis2                                      0.244E+00  0.773E-01 -0.314E+02
-     mlis2                                      0.241E+00  0.151E+00 -0.316E+02
-     mlis3     0.229E+00 -0.436E+01 -0.328E+02
-
- m1qn3: convergence rate, s(k)/s(k-1) =  0.26415E+00
-
- m1qn3: matrix update:
-     oren-spedicato factor =  0.180E+00
-
- m1qn3: stopping criterion on g:  0.92139E-01
-
- m1qn3: descent direction d: angle(-g,d) =  38.2 degrees
-
- -------------------------------------------------------------------------------
-
- m1qn3: iter   6, simul  10, f= 0.19791341E+03, h'(0)=-0.43879D+02
-
- m1qn3: line search
-
-     mlis0       fpn=-0.439E+02 d2= 0.26E+02  tmin= 0.74E-04 tmax= 0.10E+21
-     mlis2                                      0.100E+01  0.163E+02  0.116E+02
-     mlis2                                      0.210E+00  0.136E+00 -0.335E+02
-     mlis2                                      0.207E+00  0.210E+00 -0.336E+02
-     mlis2                                      0.197E+00  0.607E+00 -0.342E+02
-     mlis3     0.148E+00 -0.515E+01 -0.368E+02
-
- m1qn3: convergence rate, s(k)/s(k-1) =  0.60951E+00
-
- m1qn3: matrix update:
-     oren-spedicato factor =  0.236E+00
-
- m1qn3: stopping criterion on g:  0.82234E-01
-
- m1qn3: descent direction d: angle(-g,d) =  35.6 degrees
-
- -------------------------------------------------------------------------------
-
- m1qn3: iter   7, simul  15, f= 0.19276703E+03, h'(0)=-0.52331D+02
-
- m1qn3: line search
-
-     mlis0       fpn=-0.523E+02 d2= 0.44E+02  tmin= 0.59E-04 tmax= 0.10E+21
-     mlis2                                      0.100E+01  0.171E+02  0.181E+02
-     mlis2                                      0.229E+00  0.216E+00 -0.350E+02
-     mlis2                                      0.227E+00  0.297E+00 -0.352E+02
-     mlis2                                      0.215E+00  0.970E+00 -0.373E+02
-     mlis2                                      0.161E+00  0.344E+00 -0.410E+02
-     mlis2                                      0.113E+00  0.259E+01 -0.443E+02
-     mlis2                                      0.791E-01  0.615E+01 -0.460E+02
-     mlis2                                      0.554E-01  0.587E+01 -0.481E+02
-     mlis3     0.388E-01 -0.199E+01 -0.497E+02
-     mlis3     0.438E-01 -0.223E+01 -0.493E+02
-     mlis2                                      0.519E-01  0.487E+01 -0.488E+02
-     mlis2                                      0.494E-01  0.498E+01 -0.490E+02
-     mlis2                                      0.477E-01  0.507E+01 -0.491E+02
-     mlis2                                      0.465E-01  0.513E+01 -0.491E+02
-     mlis2                                      0.457E-01  0.518E+01 -0.492E+02
-     mlis3     0.451E-01 -0.230E+01 -0.493E+02
-     mlis3     0.453E-01 -0.231E+01 -0.492E+02
-     mlis2                                      0.456E-01  0.518E+01 -0.492E+02
-     mlis3     0.455E-01 -0.232E+01 -0.492E+02
-     mlis2                                      0.455E-01  0.518E+01 -0.492E+02
-
-     mlis1   STOP on tmin        step           functions     derivatives
-     mlis5                  0.45495603E-01    0.19044658E+03 -0.492E+02
-     mlis5                  0.45521464E-01    0.19795140E+03 -0.492E+02
-
- -------------------------------------------------------------------------------
-
- m1qn3: output mode is  6
-     number of iterations:    7
-     number of simulations:     35
-     realized relative precision on g:  0.82E-01
-     norm of x =  0.19350844E+04
-     f         =  0.19044658E+03
-     norm of g =  0.94057922E+01
Index: /trunk/scripts/script_outputs/run_batch_sinobad_seq.sh.e329
===================================================================
--- /trunk/scripts/script_outputs/run_batch_sinobad_seq.sh.e329	(revision 41)
+++ /trunk/scripts/script_outputs/run_batch_sinobad_seq.sh.e329	(revision 41)
@@ -0,0 +1,160 @@
+In file included from Ysinobad.cpp:561:
+Y2sinobad.h:3919: attention : ignoring #pragma omp parallel
+Y2sinobad.h:3939: attention : ignoring #pragma omp parallel
+In file included from Ysinobad.cpp:561:
+Y2sinobad.h:4010: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4082: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4094: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4106: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4171: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4194: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4217: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4270: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4354: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4373: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4387: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4445: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4474: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4497: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4631: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4732: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4747: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4765: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4782: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4854: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4952: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4962: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4980: attention : ignoring #pragma omp parallel
+Y2sinobad.h:4997: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5015: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5034: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5103: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5154: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5168: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5180: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5187: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5200: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5212: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5219: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5231: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5241: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5248: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5267: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5269: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5273: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5278: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5287: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5423: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5452: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5456: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5461: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5599: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5626: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5643: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5647: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5651: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5656: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5667: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5675: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5809: attention : ignoring #pragma omp parallel
+Y2sinobad.h:5841: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5845: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5855: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5859: attention : ignoring #pragma omp atomic
+Y2sinobad.h:5938: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6083: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6087: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6089: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6093: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6100: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6115: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6119: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6131: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6135: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6141: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6169: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6173: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6177: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6179: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6183: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6185: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6189: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6194: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6269: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6273: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6277: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6279: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6283: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6287: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6294: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6304: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6316: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6329: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6331: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6333: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6343: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6345: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6347: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6464: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6476: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6484: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6494: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6498: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6503: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6567: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6571: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6575: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6589: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6593: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6597: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6602: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6618: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6622: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6634: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6638: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6643: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6658: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6660: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6664: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6675: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6677: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6681: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6686: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6758: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6760: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6762: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6766: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6770: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6786: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6788: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6790: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6794: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6798: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6805: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6821: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6834: attention : ignoring #pragma omp atomic
+Y2sinobad.h:6839: attention : ignoring #pragma omp parallel
+Y2sinobad.h:6977: attention : ignoring #pragma omp parallel
+Y2sinobad.h:7056: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7060: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7062: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7066: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7070: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7074: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7076: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7080: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7084: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7086: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7090: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7094: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7098: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7100: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7104: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7108: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7113: attention : ignoring #pragma omp parallel
+Y2sinobad.h:7132: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7136: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7138: attention : ignoring #pragma omp atomic
+Y2sinobad.h:7142: attention : ignoring #pragma omp atomic
+995.76user 0.69system 16:47.55elapsed 98%CPU (0avgtext+0avgdata 6011792maxresident)k
+0inputs+39632outputs (0major+749524minor)pagefaults 0swaps
Index: /trunk/scripts/script_outputs/run_batch_sinobad_seq.sh.o329
===================================================================
--- /trunk/scripts/script_outputs/run_batch_sinobad_seq.sh.o329	(revision 41)
+++ /trunk/scripts/script_outputs/run_batch_sinobad_seq.sh.o329	(revision 41)
@@ -0,0 +1,582 @@
+Warning: no access to tty (Bad file descriptor).
+Thus no job control in this shell.
+appending ferret path
+appending ferret path
+ CSHRC: BINTYPE linux86_64 ...
+ --> test d'imprimante :  hpl5.
+
+***>> DISPLAY n'est pas défini. Pour le définir tapez d<NOMDUTERMINAL>
+        <NOMDUTERMINAL> correspond à l'écran où vous voulez faire afficher
+        vos fenêtres. Par exemple :
+
+  dboreas pour la Station Solaris boreas.lodyc.jussieu.fr:0 134.157.9.52:0
+  ddedale pour la Station Linux dedale.lodyc.jussieu.fr:0   134.157.9.123:0
+  calypso pour le Terminal X calypso.lodyc.jussieu.fr:0     134.157.9.3:0
+
+                 Bienvenue au LODYC
+    Actuellement vous êtes logué sur zelos.locean-ipsl.upmc.fr
+
+ 15:25:01 up 197 days,  3:25,  0 users,  load average: 0,65, 0,24, 0,08
+USER     TTY        LOGIN@   IDLE   JCPU   PCPU WHAT
+
+TERM ... vt100
+TTY .... n'est pas un tty
+****************************************
+*       run reference experiment       *
+*             (sequential)             *
+*              26/04/2011              *
+*  julien.brajard@locean-ipsl.upmc.fr  *
+****************************************
+ > MACHINE TYPE determined now ... <linux>.
+ > MQNDIR directory variable not declared...standard directory.
+
+       ====================================================================
+	            YaoI V.9 : PROJECT CREATION PROCEDURE                   
+	            01.11.2008 (release 01.11.2008)                         
+       ====================================================================
+        LOCAL MACHINE BINARY TYPE FOUND ..... <linux>                      
+       ====================================================================
+
+ > Run of YaoI for project : sinobad
+ > HEADERS generation not required: 
+ >   the description file has not changed from the last generation of Y1sinobad.h and Y2sinobad.h.
+ >   Pay attention all the warning showed in the generation process are not been displayed this way because of the non regeneration.
+ >   If you want to force the regeneration remove the Yworkdir/ directory.
+ > enter proc construction 
+ Extern object : 
+                 ../object/ncutil.o 
+ > end proc construction 
+ > start compilation 
+ > end of compilation, start linking 
+
+       =======================  PROJECT IS NOW READY  =====================
+                           YaoI V.9 : END OF PROCEDURE                     
+                           01.11.2008 (release 01.11.2008)                 
+       ====================================================================
+                           TOTAL COMPILATION TIME : 15 second(s)            
+       ====================================================================
+
+
+Starting Yao : YaoI V.9 01.11.2008 (release 01.11.2008)  -O3  
+Project sinobad. 
+(You are in double mode)
+This is a parallel version of the application, the max number of threads that the system supports are: 1.
+RUN OPTIMIZATION WITH M1QN3 .......
+////////////////////////////////////////////////////////////////////////
+//                       NEMO/YAO PROJECT                             //
+//                     M. Berrada 02-2009                            //
+//                  LOCEAN-IPSL.UPMC (Paris 6)                        //
+//====================================================================//
+norm_sup_diff_ta_c :  0.0000000000000000e+00
+norm_2_diff_ta_c :  0.0000000000000000e+00
+Yao#>
+Yao#>xistate_init 1 ../data_in/file_rest/GYRE_00000424_restart.nc
+start from rest_file
+ Yao#>true_target_in_tab ta_c
+Yao#>xinitnc ../data_out/exp_T.nc
+Yao#>
+Yao#>xwriteout 0 true ../data_out/exp_T.nc
+Yao#>
+Yao#>print_time OFF
+Yao#>set_modeltime 0
+Yao#>FORWARD
+Yao#>
+Yao#>
+Yao#>xwriteout 50 obs_48 ../data_out/exp_T.nc
+Yao#>savestate sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat   
+ 704 state(s) selected saving module sshb
+Yao#>savestate tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat   
+ 21824 state(s) selected saving module tb
+Yao#>savestate sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat   
+ 21824 state(s) selected saving module sb
+Yao#>savestate ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat   
+ 21824 state(s) selected saving module ub
+Yao#>savestate vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat   
+ 21824 state(s) selected saving module vb
+Yao#>
+Yao#>loadobs sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat D   
+modul sshb: 704 data read ; 704 state(s) selected as obs (704 inserted + 0 overrided)
+Yao#>loadobs tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat D   
+modul tb: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
+Yao#>loadobs sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat D   
+modul sb: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
+Yao#>loadobs ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat D   
+modul ub: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
+Yao#>loadobs vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat D   
+modul vb: 21824 data read ; 21824 state(s) selected as obs (21824 inserted + 0 overrided)
+Yao#>
+Yao#>
+Yao#>
+Yao#>xistate_init 1 ../data_in/file_rest/GYRE_00000400_restart.nc
+start from rest_file
+ Yao#>set_modeltime 0
+Yao#>print_time OFF
+Yao#>FORWARD
+Yao#>
+Yao#>xwriteout 0 init ../data_out/exp_T.nc
+Yao#>xwriteout 50 forw0_48 ../data_out/exp_T.nc
+Yao#>
+Yao#>goto FINTEST
+Yao#>
+Yao#>set_modeltime 0
+Yao#>setm_impres  5
+Yao#>setm_io      6
+Yao#>setm_mode    1
+Yao#>set_nbiter   100
+Yao#>setm_nsim    100
+Yao#>setm_dxmin   1.0e-4
+Yao#>setm_epsg    1.0e-10
+Yao#>setm_ddf1    1
+Yao#>runm
+
+start run time: Tue Apr 26 15:25:29 2011
+
+norm_sup_diff_ta_c :  5.6957761203654993e+00
+norm_2_diff_ta_c :  2.9458927804041817e+01
+norm_sup_diff_ta_c :  6.4525837069643242e+00
+norm_2_diff_ta_c :  2.5364756069289729e+01
+norm_sup_diff_ta_c :  6.5694796001323450e+00
+norm_2_diff_ta_c :  2.5910201309997589e+01
+norm_sup_diff_ta_c :  6.4730238085878593e+00
+norm_2_diff_ta_c :  2.5277352029028023e+01
+norm_sup_diff_ta_c :  6.4239994441875208e+00
+norm_2_diff_ta_c :  2.5277655022396658e+01
+norm_sup_diff_ta_c :  6.4981752788432825e+00
+norm_2_diff_ta_c :  2.6011383203801639e+01
+norm_sup_diff_ta_c :  6.4431740199931848e+00
+norm_2_diff_ta_c :  2.5361084473246660e+01
+norm_sup_diff_ta_c :  6.4429823314555872e+00
+norm_2_diff_ta_c :  2.5359877359486145e+01
+norm_sup_diff_ta_c :  6.4420334255107630e+00
+norm_2_diff_ta_c :  2.5354011151211974e+01
+norm_sup_diff_ta_c :  6.4375249301799524e+00
+norm_2_diff_ta_c :  2.5328662139641139e+01
+norm_sup_diff_ta_c :  6.6143342410930384e+00
+norm_2_diff_ta_c :  2.6066496113131752e+01
+norm_sup_diff_ta_c :  6.4752818500407923e+00
+norm_2_diff_ta_c :  2.5372534736558389e+01
+norm_sup_diff_ta_c :  6.7013689433940637e+00
+norm_2_diff_ta_c :  2.6234486389933608e+01
+norm_sup_diff_ta_c :  6.5553413783916223e+00
+norm_2_diff_ta_c :  2.5539378625795980e+01
+norm_sup_diff_ta_c :  6.5545407688029993e+00
+norm_2_diff_ta_c :  2.5536944541846530e+01
+norm_sup_diff_ta_c :  6.5505777751811731e+00
+norm_2_diff_ta_c :  2.5525120713987512e+01
+norm_sup_diff_ta_c :  6.9657641803630579e+00
+norm_2_diff_ta_c :  2.7331841099141510e+01
+norm_sup_diff_ta_c :  6.6724058544048059e+00
+norm_2_diff_ta_c :  2.5872523591600963e+01
+norm_sup_diff_ta_c :  6.9938589488872278e+00
+norm_2_diff_ta_c :  2.6992441893623258e+01
+norm_sup_diff_ta_c :  7.1218420421489466e+00
+norm_2_diff_ta_c :  2.7509855675717919e+01
+norm_sup_diff_ta_c :  7.0236097728618372e+00
+norm_2_diff_ta_c :  2.7094694365769563e+01
+norm_sup_diff_ta_c :  7.3259350215800989e+00
+norm_2_diff_ta_c :  2.9205030909929718e+01
+norm_sup_diff_ta_c :  7.0799638187297571e+00
+norm_2_diff_ta_c :  2.7362841880350924e+01
+norm_sup_diff_ta_c :  7.0794001972087610e+00
+norm_2_diff_ta_c :  2.7359860482350378e+01
+norm_sup_diff_ta_c :  7.0766106998332727e+00
+norm_2_diff_ta_c :  2.7345192094444286e+01
+norm_sup_diff_ta_c :  7.0633603488811243e+00
+norm_2_diff_ta_c :  2.7277530261132259e+01
+norm_sup_diff_ta_c :  7.0514351283916223e+00
+norm_2_diff_ta_c :  2.7219490206551317e+01
+norm_sup_diff_ta_c :  7.1741358196147669e+00
+norm_2_diff_ta_c :  2.8160627646491843e+01
+norm_sup_diff_ta_c :  7.0728222286113489e+00
+norm_2_diff_ta_c :  2.7324005716924308e+01
+
+ m1qn3 (version 2.0b, december 1993): entry point
+     dimension of the problem (n):    21824
+     absolute precision on x (dxmin): 0.10E-03
+     expected decrease for f (df1): 0.30E+04
+     relative precision on g (epsg): 0.10E-09
+     maximal number of iterations (niter):   100
+     maximal number of simulations (nsim):   100
+     printing level (impres):   5
+
+ m1qn3: scalar initial scaling mode
+
+     allocated memory (nrz) :   436488
+     used memory :              414664
+     number of updates :             8
+     (y,s) pairs are stored in core memory
+
+ m1qn3: cold start
+
+     f         =  0.29733594E+04
+     norm of g =  0.11776870E+03
+
+ m1qn3a: descent direction -g: precon =  0.429E+00
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter   1, simul   1, f= 0.29733594E+04, h'(0)=-0.59467D+04
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.595E+04 d2= 0.25E+04  tmin= 0.96E-05 tmax= 0.10E+21
+     mlis3     0.100E+01 -0.243E+04 -0.420E+03
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.411E+00
+
+ m1qn3: stopping criterion on g:  0.33221E+00
+
+ m1qn3: descent direction d: angle(-g,d) =  17.7 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter   2, simul   2, f= 0.54561414E+03, h'(0)=-0.72774D+03
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.728E+03 d2= 0.38E+03  tmin= 0.37E-04 tmax= 0.10E+21
+     mlis3     0.100E+01 -0.301E+03  0.724E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.38662E+00
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.346E+00
+
+ m1qn3: stopping criterion on g:  0.20451E+00
+
+ m1qn3: descent direction d: angle(-g,d) =  25.3 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter   3, simul   3, f= 0.24475352E+03, h'(0)=-0.17134D+03
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.171E+03 d2= 0.62E+02  tmin= 0.10E-03 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.188E+02  0.777E+02
+     mlis3     0.335E+00 -0.295E+02 -0.818E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.13521E+00
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.125E+00
+
+ m1qn3: stopping criterion on g:  0.12434E+00
+
+ m1qn3: descent direction d: angle(-g,d) =  21.1 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter   4, simul   5, f= 0.21527711E+03, h'(0)=-0.65060D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.651E+02 d2= 0.23E+02  tmin= 0.14E-03 tmax= 0.10E+21
+     mlis3     0.100E+01 -0.132E+02 -0.174E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.18044E+01
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.215E+00
+
+ m1qn3: stopping criterion on g:  0.10449E+00
+
+ m1qn3: descent direction d: angle(-g,d) =  42.9 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter   5, simul   6, f= 0.20204126E+03, h'(0)=-0.48310D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.483E+02 d2= 0.29E+02  tmin= 0.97E-04 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.121E+02  0.187E+02
+     mlis2                                      0.259E+00  0.308E-01 -0.304E+02
+     mlis2                                      0.256E+00  0.692E-01 -0.305E+02
+     mlis2                                      0.243E+00  0.316E+00 -0.313E+02
+     mlis3     0.182E+00 -0.184E+01 -0.362E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.20514E+00
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.139E+00
+
+ m1qn3: stopping criterion on g:  0.95869E-01
+
+ m1qn3: descent direction d: angle(-g,d) =  38.0 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter   6, simul  11, f= 0.20019969E+03, h'(0)=-0.52585D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.526E+02 d2= 0.35E+02  tmin= 0.66E-04 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.227E+02  0.275E+02
+     mlis3     0.214E+00 -0.142E+01 -0.360E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.12925E+01
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.150E+00
+
+ m1qn3: stopping criterion on g:  0.82567E-01
+
+ m1qn3: descent direction d: angle(-g,d) =  35.8 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter   7, simul  13, f= 0.19878349E+03, h'(0)=-0.44408D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.444E+02 d2= 0.32E+02  tmin= 0.69E-04 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.192E+01  0.133E+02
+     mlis2                                      0.354E+00  0.136E+01 -0.273E+02
+     mlis2                                      0.351E+00  0.945E+00 -0.276E+02
+     mlis3     0.333E+00 -0.585E+01 -0.275E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.14855E+01
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.133E+00
+
+ m1qn3: stopping criterion on g:  0.78750E-01
+
+ m1qn3: descent direction d: angle(-g,d) =  43.1 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter   8, simul  17, f= 0.19293538E+03, h'(0)=-0.47398D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.474E+02 d2= 0.49E+02  tmin= 0.67E-04 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.152E+02  0.454E+02
+     mlis3     0.293E+00 -0.108E+02 -0.231E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.10953E+01
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.242E+00
+
+ m1qn3: stopping criterion on g:  0.68289E-01
+
+ m1qn3: descent direction d: angle(-g,d) =  47.4 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter   9, simul  19, f= 0.18213156E+03, h'(0)=-0.23784D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.238E+02 d2= 0.19E+02  tmin= 0.12E-03 tmax= 0.10E+21
+     mlis3     0.100E+01 -0.262E+01  0.157E+01
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.21273E+01
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.252E+00
+
+ m1qn3: stopping criterion on g:  0.77526E-01
+
+ m1qn3: descent direction d: angle(-g,d) =  50.4 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter  10, simul  20, f= 0.17950673E+03, h'(0)=-0.13960D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.140E+02 d2= 0.57E+01  tmin= 0.21E-03 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.130E+01 -0.613E+01
+     mlis3     0.232E+00 -0.746E+00 -0.119E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.12746E+00
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.341E-01
+
+ m1qn3: stopping criterion on g:  0.77492E-01
+
+ m1qn3: descent direction d: angle(-g,d) =  47.6 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter  11, simul  22, f= 0.17876028E+03, h'(0)=-0.43376D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.434E+02 d2= 0.50E+02  tmin= 0.69E-04 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.238E+02  0.228E+02
+     mlis2                                      0.186E+00  0.508E-01 -0.281E+02
+     mlis2                                      0.185E+00  0.110E+00 -0.282E+02
+     mlis2                                      0.175E+00  0.315E+00 -0.289E+02
+     mlis2                                      0.131E+00  0.169E+01 -0.321E+02
+     mlis3     0.920E-01 -0.307E+01 -0.353E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.11652E+01
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.151E+00
+
+ m1qn3: stopping criterion on g:  0.68888E-01
+
+ m1qn3: descent direction d: angle(-g,d) =  37.4 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter  12, simul  28, f= 0.17569373E+03, h'(0)=-0.31251D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.313E+02 d2= 0.24E+02  tmin= 0.11E-03 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.154E+02  0.398E+01
+     mlis3     0.174E+00 -0.524E+00 -0.234E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.13020E+01
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.759E-01
+
+ m1qn3: stopping criterion on g:  0.63335E-01
+
+norm_sup_diff_ta_c :  7.3367048656352747e+00
+norm_2_diff_ta_c :  2.9583294862566490e+01
+norm_sup_diff_ta_c :  7.1169577991374720e+00
+norm_2_diff_ta_c :  2.7602124338220236e+01
+norm_sup_diff_ta_c :  7.3071843540080774e+00
+norm_2_diff_ta_c :  3.0043100486932495e+01
+norm_sup_diff_ta_c :  7.1400934612163294e+00
+norm_2_diff_ta_c :  2.7814698337905291e+01
+norm_sup_diff_ta_c :  7.1398621951946009e+00
+norm_2_diff_ta_c :  2.7812449424164779e+01
+norm_sup_diff_ta_c :  7.1387168323405970e+00
+norm_2_diff_ta_c :  2.7801352960701074e+01
+norm_sup_diff_ta_c :  7.1332770740398157e+00
+norm_2_diff_ta_c :  2.7749472997647597e+01
+norm_sup_diff_ta_c :  7.1283813869365442e+00
+norm_2_diff_ta_c :  2.7703958029693364e+01
+norm_sup_diff_ta_c :  7.2231961643108118e+00
+norm_2_diff_ta_c :  2.8743583930850225e+01
+norm_sup_diff_ta_c :  7.1424633418925989e+00
+norm_2_diff_ta_c :  2.7834458045384906e+01
+norm_sup_diff_ta_c :  7.1423226749309290e+00
+norm_2_diff_ta_c :  2.7833110842932065e+01
+norm_sup_diff_ta_c :  7.1416255390056360e+00
+norm_2_diff_ta_c :  2.7826453189497322e+01
+norm_sup_diff_ta_c :  7.1383143817790735e+00
+norm_2_diff_ta_c :  2.7795115690263486e+01
+norm_sup_diff_ta_c :  7.1353346263774622e+00
+norm_2_diff_ta_c :  2.7767318265964022e+01
+norm_sup_diff_ta_c :  7.1332484638103235e+00
+norm_2_diff_ta_c :  2.7748088832739079e+01
+norm_sup_diff_ta_c :  7.1338745509990442e+00
+norm_2_diff_ta_c :  2.7753837924537869e+01
+norm_sup_diff_ta_c :  7.1343122875102747e+00
+norm_2_diff_ta_c :  2.7757872199172485e+01
+norm_sup_diff_ta_c :  7.1341811572917688e+00
+norm_2_diff_ta_c :  2.7756660616271883e+01
+norm_sup_diff_ta_c :  7.1342202579387415e+00
+norm_2_diff_ta_c :  2.7757024531087740e+01
+norm_sup_diff_ta_c :  7.1342088138469446e+00
+norm_2_diff_ta_c :  2.7756914985620227e+01
+norm_sup_diff_ta_c :  7.1342002307780970e+00
+norm_2_diff_ta_c :  2.7756839169203612e+01
+norm_sup_diff_ta_c :  7.1341945087321985e+00
+norm_2_diff_ta_c :  2.7756785474136279e+01
+Yao#>(6) Arret sur dxmin pendant le line-search 
+
+start run time ....: Tue Apr 26 15:25:29 2011
+  end run time ....: Tue Apr 26 15:42:03 2011
+  run duration ....: 994s
+
+Yao#>xwriteout 0 fin ../data_out/exp_T.nc
+Yao#>xwriteout 50 forwfin_48 ../data_out/exp_T.nc
+Yao#>
+Yao#>xrst_save ../data_in/file_rest/GYRE_00000400_restart_CONTROL_T.nc
+Yao#>
+Yao#>
+Yao#>
+Yao#>goto FINRUN
+Yao#>
+ m1qn3: descent direction d: angle(-g,d) =  49.6 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter  13, simul  30, f= 0.17516956E+03, h'(0)=-0.32339D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.323E+02 d2= 0.45E+02  tmin= 0.77E-04 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.203E+02  0.143E+02
+     mlis3     0.167E+00 -0.616E+00 -0.249E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.13233E+01
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.178E+00
+
+ m1qn3: stopping criterion on g:  0.60300E-01
+
+ m1qn3: descent direction d: angle(-g,d) =  49.2 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter  14, simul  32, f= 0.17455342E+03, h'(0)=-0.32767D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.328E+02 d2= 0.50E+02  tmin= 0.81E-04 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.362E+02  0.248E+02
+     mlis2                                      0.122E+00  0.184E+00 -0.235E+02
+     mlis2                                      0.120E+00  0.213E+00 -0.236E+02
+     mlis2                                      0.114E+00  0.363E+00 -0.240E+02
+     mlis2                                      0.858E-01  0.792E+00 -0.253E+02
+     mlis3     0.601E-01 -0.308E+01 -0.258E+02
+
+ m1qn3: convergence rate, s(k)/s(k-1) =  0.37937E+00
+
+ m1qn3: matrix update:
+     oren-spedicato factor =  0.383E-01
+
+ m1qn3: stopping criterion on g:  0.54590E-01
+
+ m1qn3: descent direction d: angle(-g,d) =  50.5 degrees
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: iter  15, simul  38, f= 0.17146970E+03, h'(0)=-0.13935D+02
+
+ m1qn3: line search
+
+     mlis0       fpn=-0.139E+02 d2= 0.12E+02  tmin= 0.16E-03 tmax= 0.10E+21
+     mlis2                                      0.100E+01  0.772E+01 -0.313E+01
+     mlis2                                      0.149E+00  0.165E+01 -0.131E+02
+     mlis2                                      0.147E+00  0.174E+01 -0.131E+02
+     mlis2                                      0.140E+00  0.175E+01 -0.132E+02
+     mlis2                                      0.105E+00  0.181E+01 -0.136E+02
+     mlis2                                      0.733E-01  0.244E+01 -0.138E+02
+     mlis3     0.513E-01 -0.551E+00 -0.133E+02
+     mlis3     0.579E-01 -0.727E+00 -0.132E+02
+     mlis2                                      0.626E-01  0.277E+01 -0.140E+02
+     mlis3     0.612E-01 -0.779E+00 -0.131E+02
+     mlis2                                      0.616E-01  0.278E+01 -0.140E+02
+     mlis2                                      0.615E-01  0.279E+01 -0.140E+02
+     mlis2                                      0.614E-01  0.279E+01 -0.140E+02
+     mlis2                                      0.613E-01  0.279E+01 -0.140E+02
+
+     mlis1   STOP on tmin        step           functions     derivatives
+     mlis5                  0.61167818E-01    0.17069031E+03 -0.131E+02
+     mlis5                  0.61310437E-01    0.17425839E+03 -0.140E+02
+
+ -------------------------------------------------------------------------------
+
+ m1qn3: output mode is  6
+     number of iterations:   15
+     number of simulations:     52
+     realized relative precision on g:  0.55E-01
+     norm of x =  0.19348639E+04
+     f         =  0.17069031E+03
+     norm of g =  0.70167274E+01
Index: /trunk/src/Yworkdir/Y2sinobad.h
===================================================================
--- /trunk/src/Yworkdir/Y2sinobad.h	(revision 40)
+++ /trunk/src/Yworkdir/Y2sinobad.h	(revision 41)
@@ -3873,5 +3873,4 @@
 	 Yting[0]=YS1_tb( Yi, Yj, 0, YTemps-1);
 	 Ysbcgyre_qsr_qns(Yi, Yj)->forward( Yting[0] );
-	 Ydbg_nanf ("F: ", YNBI_sbcgyre_qsr_qns, "sbcgyre_qsr_qns");
 	 Ysbcgyre_emp1(Yi, Yj)->forward();
 	 if (1==0 || Yi-1<0)
@@ -3893,5 +3892,4 @@
 	 Yting[4]=YS1_sbcgyre_emp1( Yi, Yj, YTemps);
 	 Ysbcgyre_zsumemp_zsurf(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_sbcgyre_zsumemp_zsurf, "sbcgyre_zsumemp_zsurf");
 
 }
@@ -3903,8 +3901,6 @@
 	 Yting[2]=YS2_sbcgyre_zsumemp_zsurf( 31, 21, YTemps);
 	 Ysbcgyre_emp(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_sbcgyre_emp, "sbcgyre_emp");
 	 Yting[0]=YS1_sbcgyre_emp( Yi, Yj, YTemps);
 	 Ysbcgyre_emps(Yi, Yj)->forward( Yting[0] );
-	 Ydbg_nanf ("F: ", YNBI_sbcgyre_emps, "sbcgyre_emps");
 	 Ysbcgyre_utau_vtau(Yi, Yj)->forward();
 	 if (1==0 || Yi-1<0)
@@ -3919,5 +3915,4 @@
 	 Yting[3]=YS2_sbcgyre_utau_vtau( Yi, Yj, YTemps);
 	 Ysbcgyre_wndm(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_sbcgyre_wndm, "sbcgyre_wndm");
 
 }
@@ -3939,5 +3934,4 @@
 	 Yting[3]=YS1_sb( Yi, Yj, Yk-1, YTemps-1);
 	 Ybn2_rn2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_bn2_rn2, "bn2_rn2");
 
 }
@@ -4011,5 +4005,4 @@
 	 Yting[18]=YS1_va( Yi+1, Yj, 1, YTemps-1);
 	 Yzdf_avt_avmu_avmv(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17], Yting[18] );
-	 Ydbg_nanf ("F: ", YNBI_zdf_avt_avmu_avmv, "zdf_avt_avmu_avmv");
 
 }
@@ -4052,5 +4045,4 @@
 	 Yting[30]=YS1_zdf_avt_avmu_avmv( Yi, Yj, 30, YTemps);
 	 Yzdfmxl_hmld(Yi, Yj)->forward();
-	 Ydbg_nanf ("F: ", YNBI_zdfmxl_hmld, "zdfmxl_hmld");
 	 Yting[0]=YS1_rhop( Yi, Yj, 0, YTemps-1);
 	 Yting[1]=YS1_rhop( Yi, Yj, 1, YTemps-1);
@@ -4085,5 +4077,4 @@
 	 Yting[30]=YS1_rhop( Yi, Yj, 30, YTemps-1);
 	 Yzdfmxl_nmln_hmlp_hmlpt(Yi, Yj)->forward();
-	 Ydbg_nanf ("F: ", YNBI_zdfmxl_nmln_hmlp_hmlpt, "zdfmxl_nmln_hmlp_hmlpt");
 
 }
@@ -4098,5 +4089,4 @@
 	 Yting[2]=YS1_sa( Yi, Yj, 0, YTemps-1);
 	 Yta_sa_trasbc(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_ta_sa_trasbc, "ta_sa_trasbc");
 
 }
@@ -4111,5 +4101,4 @@
 	 Yting[1]=YS1_ta_sa_trasbc( Yi, Yj, YTemps);
 	 Yta_traqsr(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_ta_traqsr, "ta_traqsr");
 
 }
@@ -4134,5 +4123,4 @@
 	 Yting[4]=YS1_ta( Yi+1, Yj, Yk, YTemps-1);
 	 Yzcent_hor_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_zcent_hor_traadv_cen2, "zcent_hor_traadv_cen2");
 	 Yting[0]=YS1_ua( Yi, Yj, Yk, YTemps-1);
 	 Yting[1]=YS1_va( Yi, Yj, Yk, YTemps-1);
@@ -4147,5 +4135,4 @@
 	 Yting[4]=YS1_sa( Yi+1, Yj, Yk, YTemps-1);
 	 Yzcens_hor_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_zcens_hor_traadv_cen2, "zcens_hor_traadv_cen2");
 
 }
@@ -4165,5 +4152,4 @@
 	 Yting[4]=YS2_zcent_hor_traadv_cen2( Yi, Yj, Yk, YTemps);
 	 Yta_hor_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_ta_hor_traadv_cen2, "ta_hor_traadv_cen2");
 	 Yting[0]=YS2_ta_sa_trasbc( Yi, Yj, YTemps);
 	 if (1==0 || Yi-1<0)
@@ -4178,5 +4164,4 @@
 	 Yting[4]=YS2_zcens_hor_traadv_cen2( Yi, Yj, Yk, YTemps);
 	 Ysa_hor_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_sa_hor_traadv_cen2, "sa_hor_traadv_cen2");
 
 }
@@ -4197,5 +4182,4 @@
 	 Yting[2]=YS1_ta( Yi, Yj, Yk-1, YTemps-1);
 	 Yzcent_ver_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_zcent_ver_traadv_cen2, "zcent_ver_traadv_cen2");
 	 Yting[0]=YS1_wa( Yi, Yj, Yk, YTemps-1);
 	 Yting[1]=YS1_sa( Yi, Yj, Yk, YTemps-1);
@@ -4205,5 +4189,4 @@
 	 Yting[2]=YS1_sa( Yi, Yj, Yk-1, YTemps-1);
 	 Yzcens_ver_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_zcens_ver_traadv_cen2, "zcens_ver_traadv_cen2");
 
 }
@@ -4222,5 +4205,4 @@
 	 Yting[2]=YS1_zcent_ver_traadv_cen2( Yi, Yj, Yk+1, YTemps);
 	 Yta_ver_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_ta_ver_traadv_cen2, "ta_ver_traadv_cen2");
 	 Yting[0]=YS1_sa_hor_traadv_cen2( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_zcens_ver_traadv_cen2( Yi, Yj, Yk, YTemps);
@@ -4230,5 +4212,4 @@
 	 Yting[2]=YS1_zcens_ver_traadv_cen2( Yi, Yj, Yk+1, YTemps);
 	 Ysa_ver_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_sa_ver_traadv_cen2, "sa_ver_traadv_cen2");
 
 }
@@ -4250,5 +4231,4 @@
 	 Yting[2]=YS1_tb( Yi, Yj+1, Yk, YTemps-1);
 	 Yztuv_traldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_ztuv_traldf_lap, "ztuv_traldf_lap");
 	 Yting[0]=YS1_sb( Yi, Yj, Yk, YTemps-1);
 	 if (1==0 || Yi+1>YA1_S3dt-1)
@@ -4261,5 +4241,4 @@
 	 Yting[2]=YS1_sb( Yi, Yj+1, Yk, YTemps-1);
 	 Yzsuv_traldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_zsuv_traldf_lap, "zsuv_traldf_lap");
 	 Yting[0]=YS1_ta_ver_traadv_cen2( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_ztuv_traldf_lap( Yi, Yj, Yk, YTemps);
@@ -4274,5 +4253,4 @@
 	 Yting[4]=YS2_ztuv_traldf_lap( Yi, Yj-1, Yk, YTemps);
 	 Yta_traldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_ta_traldf_lap, "ta_traldf_lap");
 	 Yting[0]=YS1_sa_ver_traadv_cen2( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_zsuv_traldf_lap( Yi, Yj, Yk, YTemps);
@@ -4287,5 +4265,4 @@
 	 Yting[4]=YS2_zsuv_traldf_lap( Yi, Yj-1, Yk, YTemps);
 	 Ysa_traldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_sa_traldf_lap, "sa_traldf_lap");
 
 }
@@ -4303,13 +4280,10 @@
 	 Yting[1]=YS1_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps);
 	 Yzwisd_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_zwisd_trazdf_imp, "zwisd_trazdf_imp");
 	 Yting[0]=YS1_ta_traldf_lap( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_tb( Yi, Yj, Yk, YTemps-1);
 	 Yzrhs_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_zrhs_trazdf_imp, "zrhs_trazdf_imp");
 	 Yting[0]=YS1_sa_traldf_lap( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_sb( Yi, Yj, Yk, YTemps-1);
 	 Yzrhs_sa_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_zrhs_sa_trazdf_imp, "zrhs_sa_trazdf_imp");
 
 }
@@ -4330,5 +4304,4 @@
 	 Yting[3]=YS1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps);
 	 Yzwt_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_zwt_trazdf_imp, "zwt_trazdf_imp");
 	 Yting[0]=YS1_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -4342,5 +4315,4 @@
 	 Yting[3]=YS1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps);
 	 Yta_de_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_ta_de_trazdf_imp, "ta_de_trazdf_imp");
 	 Yting[0]=YS1_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -4354,5 +4326,4 @@
 	 Yting[3]=YS1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps);
 	 Ysa_de_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_sa_de_trazdf_imp, "sa_de_trazdf_imp");
 
 }
@@ -4370,5 +4341,4 @@
 	 Yting[4]=YS1_ta_c( Yi, Yj, Yk);
 	 Yta(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_ta, "ta");
 	 Yting[0]=YS2_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -4380,5 +4350,4 @@
 	 Yting[4]=YS1_sa_c( Yi, Yj, Yk);
 	 Ysa(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_sa, "sa");
 
 }
@@ -4394,5 +4363,4 @@
 	 Yting[3]=YS1_ta_c( Yi, Yj, Yk);
 	 Ytb(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_tb, "tb");
 	 Yting[0]=YS1_sb( Yi, Yj, Yk, YTemps-1);
 	 Yting[1]=YS1_sa( Yi, Yj, Yk, YTemps-1);
@@ -4400,5 +4368,4 @@
 	 Yting[3]=YS1_sa_c( Yi, Yj, Yk);
 	 Ysb(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_sb, "sb");
 
 }
@@ -4413,8 +4380,6 @@
 	 Yting[1]=YS1_sb( Yi, Yj, Yk, YTemps);
 	 Yrhd(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_rhd, "rhd");
 	 Yting[0]=YS1_rhd( Yi, Yj, Yk, YTemps);
 	 Yrhop(Yi, Yj, Yk)->forward( Yting[0] );
-	 Ydbg_nanf ("F: ", YNBI_rhop, "rhop");
 
 }
@@ -4446,5 +4411,4 @@
 	 Yting[5]=YS1_va( Yi, Yj+1, Yk, YTemps-1);
 	 Yhmf_TF_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5] );
-	 Ydbg_nanf ("F: ", YNBI_hmf_TF_dynadv_cen2, "hmf_TF_dynadv_cen2");
 
 }
@@ -4463,5 +4427,4 @@
 	 Yting[3]=YS2_hmf_TF_dynadv_cen2( Yi, Yj-1, Yk, YTemps);
 	 Yzua_dhmf_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_zua_dhmf_dynadv_cen2, "zua_dhmf_dynadv_cen2");
 	 Yting[0]=YS3_hmf_TF_dynadv_cen2( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yi-1<0)
@@ -4475,5 +4438,4 @@
 	 Yting[3]=YS4_hmf_TF_dynadv_cen2( Yi, Yj-1, Yk, YTemps);
 	 Yzva_dhmf_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_zva_dhmf_dynadv_cen2, "zva_dhmf_dynadv_cen2");
 
 }
@@ -4507,5 +4469,4 @@
 	 Yting[6]=YS1_va( Yi, Yj, Yk-1, YTemps-1);
 	 Yzfuv_UW_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6] );
-	 Ydbg_nanf ("F: ", YNBI_zfuv_UW_dynadv_cen2, "zfuv_UW_dynadv_cen2");
 
 }
@@ -4524,5 +4485,4 @@
 	 Yting[2]=YS1_zua_dhmf_dynadv_cen2( Yi, Yj, Yk, YTemps);
 	 Yua_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_ua_dynadv_cen2, "ua_dynadv_cen2");
 	 Yting[0]=YS2_zfuv_UW_dynadv_cen2( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -4532,5 +4492,4 @@
 	 Yting[2]=YS1_zva_dhmf_dynadv_cen2( Yi, Yj, Yk, YTemps);
 	 Yva_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_va_dynadv_cen2, "va_dynadv_cen2");
 
 }
@@ -4554,5 +4513,4 @@
 	 Yting[3]=YS1_va( Yi+1, Yj, Yk, YTemps-1);
 	 Yzwxyz_vor_ene(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_zwxyz_vor_ene, "zwxyz_vor_ene");
 
 }
@@ -4580,5 +4538,4 @@
 	 Yting[6]=YS2_zwxyz_vor_ene( Yi+1, Yj, Yk, YTemps);
 	 Yua_vor_ene(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6] );
-	 Ydbg_nanf ("F: ", YNBI_ua_vor_ene, "ua_vor_ene");
 	 Yting[0]=YS1_va_dynadv_cen2( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yi-1<0)
@@ -4601,5 +4558,4 @@
 	 Yting[6]=YS1_zwxyz_vor_ene( Yi, Yj+1, Yk, YTemps);
 	 Yva_vor_ene(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6] );
-	 Ydbg_nanf ("F: ", YNBI_va_vor_ene, "va_vor_ene");
 	 Yting[0]=YS1_ua_vor_ene( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yj-1<0)
@@ -4614,5 +4570,4 @@
 	 Yting[4]=YS1_hdivn( Yi+1, Yj, Yk, YTemps-2);
 	 Yua_dyn_ldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_ua_dyn_ldf_lap, "ua_dyn_ldf_lap");
 	 Yting[0]=YS1_va_vor_ene( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yi-1<0)
@@ -4627,5 +4582,4 @@
 	 Yting[4]=YS1_hdivn( Yi, Yj+1, Yk, YTemps-2);
 	 Yva_dyn_ldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
-	 Ydbg_nanf ("F: ", YNBI_va_dyn_ldf_lap, "va_dyn_ldf_lap");
 
 }
@@ -4667,13 +4621,10 @@
 	 Yting[7]=YS2_zhpij_dyn_hpg_zco( Yi, Yj, Yk-1, YTemps);
 	 Yzhpij_dyn_hpg_zco(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7] );
-	 Ydbg_nanf ("F: ", YNBI_zhpij_dyn_hpg_zco, "zhpij_dyn_hpg_zco");
 	 Yting[0]=YS1_ua_dyn_ldf_lap( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps);
 	 Yua_dyn_hpg_zco(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_ua_dyn_hpg_zco, "ua_dyn_hpg_zco");
 	 Yting[0]=YS1_va_dyn_ldf_lap( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS2_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps);
 	 Yva_dyn_hpg_zco(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_va_dyn_hpg_zco, "va_dyn_hpg_zco");
 
 }
@@ -4690,5 +4641,4 @@
 	 Yting[1]=YS2_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps);
 	 Yzwisd_ua_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_zwisd_ua_dynzdf_imp, "zwisd_ua_dynzdf_imp");
 	 Yting[0]=YS3_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -4697,15 +4647,12 @@
 	 Yting[1]=YS3_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps);
 	 Yzwisd_va_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_zwisd_va_dynzdf_imp, "zwisd_va_dynzdf_imp");
 	 Yting[0]=YS1_ua_dyn_hpg_zco( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1);
 	 Yting[2]=YS1_sbcgyre_utau_vtau( Yi, Yj, YTemps);
 	 Yzrhs_ua_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_zrhs_ua_dynzdf_imp, "zrhs_ua_dynzdf_imp");
 	 Yting[0]=YS1_va_dyn_hpg_zco( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_vb( Yi, Yj, Yk, YTemps-1);
 	 Yting[2]=YS2_sbcgyre_utau_vtau( Yi, Yj, YTemps);
 	 Yzrhs_va_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_zrhs_va_dynzdf_imp, "zrhs_va_dynzdf_imp");
 
 }
@@ -4726,5 +4673,4 @@
 	 Yting[3]=YS1_zwt_ua_dynzdf_imp( Yi, Yj, Yk-1, YTemps);
 	 Yzwt_ua_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_zwt_ua_dynzdf_imp, "zwt_ua_dynzdf_imp");
 	 Yting[0]=YS1_zwisd_ua_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -4738,5 +4684,4 @@
 	 Yting[3]=YS1_zwt_ua_dynzdf_imp( Yi, Yj, Yk-1, YTemps);
 	 Yua_de_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_ua_de_dynzdf_imp, "ua_de_dynzdf_imp");
 	 Yting[0]=YS1_zwisd_va_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -4750,5 +4695,4 @@
 	 Yting[3]=YS1_zwt_va_dynzdf_imp( Yi, Yj, Yk-1, YTemps);
 	 Yzwt_va_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_zwt_va_dynzdf_imp, "zwt_va_dynzdf_imp");
 	 Yting[0]=YS1_zwisd_va_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -4762,5 +4706,4 @@
 	 Yting[3]=YS1_zwt_va_dynzdf_imp( Yi, Yj, Yk-1, YTemps);
 	 Yva_de_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_va_de_dynzdf_imp, "va_de_dynzdf_imp");
 
 }
@@ -4777,5 +4720,4 @@
 	 Yting[3]=YS1_zwt_ua_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 Yua_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_ua_dynzdf_imp, "ua_dynzdf_imp");
 	 Yting[0]=YS2_zwisd_va_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -4786,5 +4728,4 @@
 	 Yting[3]=YS1_zwt_va_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 Yva_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_va_dynzdf_imp, "va_dynzdf_imp");
 
 }
@@ -4798,9 +4739,7 @@
 	 Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1);
 	 Yua_norm_dynzdf(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_ua_norm_dynzdf, "ua_norm_dynzdf");
 	 Yting[0]=YS1_va_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_vb( Yi, Yj, Yk, YTemps-1);
 	 Yva_norm_dynzdf(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_va_norm_dynzdf, "va_norm_dynzdf");
 
 }
@@ -4821,5 +4760,4 @@
 	 Yting[2]=YS1_sshn( Yi+1, Yj, YTemps-1);
 	 Yspguv_now_dynspg(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_spguv_now_dynspg, "spguv_now_dynspg");
 
 }
@@ -4835,10 +4773,8 @@
 	 Yting[2]=YS1_spguv_now_dynspg( Yi, Yj, YTemps);
 	 Yua_1_dynspg_flt(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_ua_1_dynspg_flt, "ua_1_dynspg_flt");
 	 Yting[0]=YS1_va_norm_dynzdf( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_vb( Yi, Yj, Yk, YTemps-1);
 	 Yting[2]=YS2_spguv_now_dynspg( Yi, Yj, YTemps);
 	 Yva_1_dynspg_flt(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_va_1_dynspg_flt, "va_1_dynspg_flt");
 
 }
@@ -4881,5 +4817,4 @@
 	 Yting[30]=YS1_ua_1_dynspg_flt( Yi, Yj, 30, YTemps);
 	 Ysum_spgu_dynspg_flt(Yi, Yj)->forward();
-	 Ydbg_nanf ("F: ", YNBI_sum_spgu_dynspg_flt, "sum_spgu_dynspg_flt");
 	 Yting[0]=YS1_va_1_dynspg_flt( Yi, Yj, 0, YTemps);
 	 Yting[1]=YS1_va_1_dynspg_flt( Yi, Yj, 1, YTemps);
@@ -4914,5 +4849,4 @@
 	 Yting[30]=YS1_va_1_dynspg_flt( Yi, Yj, 30, YTemps);
 	 Ysum_spgv_dynspg_flt(Yi, Yj)->forward();
-	 Ydbg_nanf ("F: ", YNBI_sum_spgv_dynspg_flt, "sum_spgv_dynspg_flt");
 
 }
@@ -4934,9 +4868,7 @@
 	 Yting[3]=YS1_sum_spgv_dynspg_flt( Yi, Yj-1, YTemps);
 	 Ygcb_dynspg_flt(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_gcb_dynspg_flt, "gcb_dynspg_flt");
 	 Yting[0]=YS1_gcx2( Yi, Yj, YTemps-1);
 	 Yting[1]=YS1_gcx2( Yi, Yj, YTemps-2);
 	 Ygcx_dynspg_flt(Yi, Yj)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_gcx_dynspg_flt, "gcx_dynspg_flt");
 
 }
@@ -5016,5 +4948,4 @@
 	 Yting[18]=YS1_solsor_dynspg_flt( Yi-1, Yj-1, YTemps);
 	 Ysolsor_dynspg_flt(Yi, Yj)->forward();
-	 Ydbg_nanf ("F: ", YNBI_solsor_dynspg_flt, "solsor_dynspg_flt");
 
 }
@@ -5026,5 +4957,4 @@
 	 Yting[0]=YS1_solsor_dynspg_flt( Yi, Yj, YTemps);
 	 Ygcx2(Yi, Yj)->forward();
-	 Ydbg_nanf ("F: ", YNBI_gcx2, "gcx2");
 
 }
@@ -5045,5 +4975,4 @@
 	 Yting[2]=YS1_gcx2( Yi, Yj+1, YTemps);
 	 Yspguv_dynspg_flt(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_spguv_dynspg_flt, "spguv_dynspg_flt");
 
 }
@@ -5059,10 +4988,8 @@
 	 Yting[2]=YS1_ua_c( Yi, Yj, Yk);
 	 Yua(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_ua, "ua");
 	 Yting[0]=YS1_va_1_dynspg_flt( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS2_spguv_dynspg_flt( Yi, Yj, YTemps);
 	 Yting[2]=YS1_va_c( Yi, Yj, Yk);
 	 Yva(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] );
-	 Ydbg_nanf ("F: ", YNBI_va, "va");
 
 }
@@ -5078,5 +5005,4 @@
 	 Yting[3]=YS1_sshn_c( Yi, Yj);
 	 Ysshn(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_sshn, "sshn");
 	 Yting[0]=YS1_sshb( Yi, Yj, YTemps-1);
 	 Yting[1]=YS1_sshn( Yi, Yj, YTemps-1);
@@ -5084,5 +5010,4 @@
 	 Yting[3]=YS1_sshn_c( Yi, Yj);
 	 Ysshb(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_sshb, "sshb");
 
 }
@@ -5099,5 +5024,4 @@
 	 Yting[3]=YS1_ua_c( Yi, Yj, Yk);
 	 Yub(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_ub, "ub");
 	 Yting[0]=YS1_vb( Yi, Yj, Yk, YTemps-1);
 	 Yting[1]=YS1_va( Yi, Yj, Yk, YTemps-1);
@@ -5105,5 +5029,4 @@
 	 Yting[3]=YS1_va_c( Yi, Yj, Yk);
 	 Yvb(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_vb, "vb");
 
 }
@@ -5126,5 +5049,4 @@
 	 Yting[3]=YS1_va( Yi, Yj-1, Yk, YTemps);
 	 Yhdivn(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_hdivn, "hdivn");
 	 if (1==0 || Yi+1>YA1_S3dt-1)
 	 Yting[0]=0; 
@@ -5138,5 +5060,4 @@
 	 Yting[3]=YS1_ua( Yi, Yj, Yk, YTemps);
 	 Yrotn(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
-	 Ydbg_nanf ("F: ", YNBI_rotn, "rotn");
 
 }
@@ -5152,5 +5073,4 @@
 	 Yting[1]=YS1_hdivn( Yi, Yj, Yk, YTemps);
 	 Ywa(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
-	 Ydbg_nanf ("F: ", YNBI_wa, "wa");
 
 }
@@ -5179,5 +5099,4 @@
 	 YG1_wa( Yi, Yj, Yk+1, YTemps) += Ytbeta[0];
 	 YG1_hdivn( Yi, Yj, Yk, YTemps) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_wa, "wa");
 
 }
@@ -5209,5 +5128,4 @@
 	 YG1_ua( Yi, Yj+1, Yk, YTemps) += Ytbeta[2];
 	 YG1_ua( Yi, Yj, Yk, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_rotn, "rotn");
 	 Yting[0]=YS1_ua( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yi-1<0)
@@ -5231,5 +5149,4 @@
 	 else 
 	 YG1_va( Yi, Yj-1, Yk, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_hdivn, "hdivn");
 
 }
@@ -5253,5 +5170,4 @@
 	 YG1_va( Yi, Yj, Yk, YTemps) += Ytbeta[2];
 	 YG1_va_c( Yi, Yj, Yk) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_vb, "vb");
 	 Yting[0]=YS1_ub( Yi, Yj, Yk, YTemps-1);
 	 Yting[1]=YS1_ua( Yi, Yj, Yk, YTemps-1);
@@ -5266,5 +5182,4 @@
 	 YG1_ua( Yi, Yj, Yk, YTemps) += Ytbeta[2];
 	 YG1_ua_c( Yi, Yj, Yk) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_ub, "ub");
 
 }
@@ -5287,5 +5202,4 @@
 	 YG1_sshn( Yi, Yj, YTemps) += Ytbeta[2];
 	 YG1_sshn_c( Yi, Yj) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_sshb, "sshb");
 	 Yting[0]=YS1_sshb( Yi, Yj, YTemps-1);
 	 Yting[1]=YS1_wa( Yi, Yj, 0, YTemps-1);
@@ -5300,5 +5214,4 @@
 	 YG1_sbcgyre_emp( Yi, Yj, YTemps) += Ytbeta[2];
 	 YG1_sshn_c( Yi, Yj) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_sshn, "sshn");
 
 }
@@ -5320,5 +5233,4 @@
 	 YG2_spguv_dynspg_flt( Yi, Yj, YTemps) += Ytbeta[1];
 	 YG1_va_c( Yi, Yj, Yk) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_va, "va");
 	 Yting[0]=YS1_ua_1_dynspg_flt( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_spguv_dynspg_flt( Yi, Yj, YTemps);
@@ -5331,5 +5243,4 @@
 	 YG1_spguv_dynspg_flt( Yi, Yj, YTemps) += Ytbeta[1];
 	 YG1_ua_c( Yi, Yj, Yk) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_ua, "ua");
 
 }
@@ -5362,5 +5273,4 @@
 	 #pragma omp atomic
 	 YG1_gcx2( Yi, Yj+1, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_spguv_dynspg_flt, "spguv_dynspg_flt");
 
 }
@@ -5377,5 +5287,4 @@
 	 #pragma omp atomic
 	 YG1_solsor_dynspg_flt( Yi, Yj, YTemps) += Ytbeta[0];
-	 Ydbg_nanf ("B: ", YNBI_gcx2, "gcx2");
 
 }
@@ -5510,5 +5419,4 @@
 	 else 
 	 YG1_solsor_dynspg_flt( Yi-1, Yj-1, YTemps) += Ytbeta[18];
-	 Ydbg_nanf ("B: ", YNBI_solsor_dynspg_flt, "solsor_dynspg_flt");
 
 }
@@ -5525,5 +5433,4 @@
 	 YG1_gcx2( Yi, Yj, YTemps-1) += Ytbeta[0];
 	 YG1_gcx2( Yi, Yj, YTemps-2) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_gcx_dynspg_flt, "gcx_dynspg_flt");
 	 Yting[0]=YS1_sum_spgu_dynspg_flt( Yi, Yj, YTemps);
 	 if (1==0 || Yi-1<0)
@@ -5549,5 +5456,4 @@
 	 #pragma omp atomic
 	 YG1_sum_spgv_dynspg_flt( Yi, Yj-1, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_gcb_dynspg_flt, "gcb_dynspg_flt");
 
 }
@@ -5623,5 +5529,4 @@
 	 YG1_va_1_dynspg_flt( Yi, Yj, 29, YTemps) += Ytbeta[29];
 	 YG1_va_1_dynspg_flt( Yi, Yj, 30, YTemps) += Ytbeta[30];
-	 Ydbg_nanf ("B: ", YNBI_sum_spgv_dynspg_flt, "sum_spgv_dynspg_flt");
 	 Yting[0]=YS1_ua_1_dynspg_flt( Yi, Yj, 0, YTemps);
 	 Yting[1]=YS1_ua_1_dynspg_flt( Yi, Yj, 1, YTemps);
@@ -5689,9 +5594,8 @@
 	 YG1_ua_1_dynspg_flt( Yi, Yj, 29, YTemps) += Ytbeta[29];
 	 YG1_ua_1_dynspg_flt( Yi, Yj, 30, YTemps) += Ytbeta[30];
-	 Ydbg_nanf ("B: ", YNBI_sum_spgu_dynspg_flt, "sum_spgu_dynspg_flt");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -5708,5 +5612,4 @@
 	 YG1_vb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 YG2_spguv_now_dynspg( Yi, Yj, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_va_1_dynspg_flt, "va_1_dynspg_flt");
 	 Yting[0]=YS1_ua_norm_dynzdf( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1);
@@ -5718,9 +5621,8 @@
 	 YG1_ub( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 YG1_spguv_now_dynspg( Yi, Yj, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_ua_1_dynspg_flt, "ua_1_dynspg_flt");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA2_S3dt-1; YY>=0; --YY){
   Yj=YY;
@@ -5739,19 +5641,18 @@
 	 Yspguv_now_dynspg(Yi, Yj)->backward( Yting[0], Yting[1], Yting[2] );
 	 Yvsmatt (YNBS_spguv_now_dynspg, YNBI_spguv_now_dynspg, YMAX_JAC_NBI, &YG1Y_spguv_now_dynspg, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sshn( Yi, Yj, YTemps-1) += Ytbeta[0];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sshn( Yi, Yj+1, YTemps-1) += Ytbeta[1];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sshn( Yi+1, Yj, YTemps-1) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_spguv_now_dynspg, "spguv_now_dynspg");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -5764,8 +5665,7 @@
 	 Yva_norm_dynzdf(Yi, Yj, Yk)->backward( Yting[0], Yting[1] );
 	 Yvsmatt (YNBS_va_norm_dynzdf, YNBI_va_norm_dynzdf, YMAX_JAC_NBI, &YG1Y_va_norm_dynzdf, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 YG1_vb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_va_norm_dynzdf, "va_norm_dynzdf");
 	 Yting[0]=YS1_ua_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1);
@@ -5773,8 +5673,7 @@
 	 Yua_norm_dynzdf(Yi, Yj, Yk)->backward( Yting[0], Yting[1] );
 	 Yvsmatt (YNBS_ua_norm_dynzdf, YNBI_ua_norm_dynzdf, YMAX_JAC_NBI, &YG1Y_ua_norm_dynzdf, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 YG1_ub( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_ua_norm_dynzdf, "ua_norm_dynzdf");
 
 }
@@ -5800,5 +5699,4 @@
 	 YG1_va_de_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[2];
 	 YG1_zwt_va_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_va_dynzdf_imp, "va_dynzdf_imp");
 	 Yting[0]=YS2_zwisd_ua_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -5817,5 +5715,4 @@
 	 YG1_ua_de_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[2];
 	 YG1_zwt_ua_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_ua_dynzdf_imp, "ua_dynzdf_imp");
 
 }
@@ -5845,5 +5742,4 @@
 	 else 
 	 YG1_zwt_va_dynzdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_va_de_dynzdf_imp, "va_de_dynzdf_imp");
 	 Yting[0]=YS1_zwisd_va_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -5867,5 +5763,4 @@
 	 else 
 	 YG1_zwt_va_dynzdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_zwt_va_dynzdf_imp, "zwt_va_dynzdf_imp");
 	 Yting[0]=YS1_zwisd_ua_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -5889,5 +5784,4 @@
 	 else 
 	 YG1_zwt_ua_dynzdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_ua_de_dynzdf_imp, "ua_de_dynzdf_imp");
 	 Yting[0]=YS1_zwisd_ua_dynzdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -5911,8 +5805,7 @@
 	 else 
 	 YG1_zwt_ua_dynzdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_zwt_ua_dynzdf_imp, "zwt_ua_dynzdf_imp");
-
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -5929,5 +5822,4 @@
 	 YG1_vb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 YG2_sbcgyre_utau_vtau( Yi, Yj, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_zrhs_va_dynzdf_imp, "zrhs_va_dynzdf_imp");
 	 Yting[0]=YS1_ua_dyn_hpg_zco( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1);
@@ -5939,5 +5831,4 @@
 	 YG1_ub( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 YG1_sbcgyre_utau_vtau( Yi, Yj, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_zrhs_ua_dynzdf_imp, "zrhs_ua_dynzdf_imp");
 	 Yting[0]=YS3_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -5948,11 +5839,10 @@
 	 Yzwisd_va_dynzdf_imp(Yi, Yj, Yk)->backward( Yting[0], Yting[1] );
 	 Yvsmatt (YNBS_zwisd_va_dynzdf_imp, YNBI_zwisd_va_dynzdf_imp, YMAX_JAC_NBI, &YG1Y_zwisd_va_dynzdf_imp, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG3_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 if (1==0 || Yk+1>YA3_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG3_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_zwisd_va_dynzdf_imp, "zwisd_va_dynzdf_imp");
 	 Yting[0]=YS2_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -5963,11 +5853,10 @@
 	 Yzwisd_ua_dynzdf_imp(Yi, Yj, Yk)->backward( Yting[0], Yting[1] );
 	 Yvsmatt (YNBS_zwisd_ua_dynzdf_imp, YNBI_zwisd_ua_dynzdf_imp, YMAX_JAC_NBI, &YG1Y_zwisd_ua_dynzdf_imp, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG2_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 if (1==0 || Yk+1>YA3_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG2_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_zwisd_ua_dynzdf_imp, "zwisd_ua_dynzdf_imp");
 
 }
@@ -5984,5 +5873,4 @@
 	 YG1_va_dyn_ldf_lap( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 YG2_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_va_dyn_hpg_zco, "va_dyn_hpg_zco");
 	 Yting[0]=YS1_ua_dyn_ldf_lap( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps);
@@ -5992,5 +5880,4 @@
 	 YG1_ua_dyn_ldf_lap( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 YG1_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_ua_dyn_hpg_zco, "ua_dyn_hpg_zco");
 	 if (1==0 || Yi+1>YA1_S3dt-1)
 	 Yting[0]=0; 
@@ -6047,8 +5934,7 @@
 	 else 
 	 YG2_zhpij_dyn_hpg_zco( Yi, Yj, Yk-1, YTemps) += Ytbeta[7];
-	 Ydbg_nanf ("B: ", YNBI_zhpij_dyn_hpg_zco, "zhpij_dyn_hpg_zco");
-
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6080,5 +5966,4 @@
 	 else 
 	 YG1_hdivn( Yi, Yj+1, Yk, YTemps-2) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_va_dyn_ldf_lap, "va_dyn_ldf_lap");
 	 Yting[0]=YS1_ua_vor_ene( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yj-1<0)
@@ -6104,5 +5989,4 @@
 	 else 
 	 YG1_hdivn( Yi+1, Yj, Yk, YTemps-2) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_ua_dyn_ldf_lap, "ua_dyn_ldf_lap");
 	 Yting[0]=YS1_va_dynadv_cen2( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yi-1<0)
@@ -6142,5 +6026,4 @@
 	 else 
 	 YG1_zwxyz_vor_ene( Yi, Yj+1, Yk, YTemps) += Ytbeta[6];
-	 Ydbg_nanf ("B: ", YNBI_va_vor_ene, "va_vor_ene");
 	 Yting[0]=YS1_ua_dynadv_cen2( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yj-1<0)
@@ -6180,5 +6063,4 @@
 	 else 
 	 YG2_zwxyz_vor_ene( Yi+1, Yj, Yk, YTemps) += Ytbeta[6];
-	 Ydbg_nanf ("B: ", YNBI_ua_vor_ene, "ua_vor_ene");
 
 }
@@ -6199,23 +6081,22 @@
 	 Yzwxyz_vor_ene(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
 	 Yvsmatt (YNBS_zwxyz_vor_ene, YNBI_zwxyz_vor_ene, YMAX_JAC_NBI, &YG1Y_zwxyz_vor_ene, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[1];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[2];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_zwxyz_vor_ene, "zwxyz_vor_ene");
-
-}
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6232,12 +6113,11 @@
 	 Yva_dynadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] );
 	 Yvsmatt (YNBS_va_dynadv_cen2, YNBI_va_dynadv_cen2, YMAX_JAC_NBI, &YG1Y_va_dynadv_cen2, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG2_zfuv_UW_dynadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 if (1==0 || Yk+1>YA3_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG2_zfuv_UW_dynadv_cen2( Yi, Yj, Yk+1, YTemps) += Ytbeta[1];
 	 YG1_zva_dhmf_dynadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_va_dynadv_cen2, "va_dynadv_cen2");
 	 Yting[0]=YS1_zfuv_UW_dynadv_cen2( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -6249,16 +6129,15 @@
 	 Yua_dynadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] );
 	 Yvsmatt (YNBS_ua_dynadv_cen2, YNBI_ua_dynadv_cen2, YMAX_JAC_NBI, &YG1Y_ua_dynadv_cen2, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_zfuv_UW_dynadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 if (1==0 || Yk+1>YA3_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_zfuv_UW_dynadv_cen2( Yi, Yj, Yk+1, YTemps) += Ytbeta[1];
 	 YG1_zua_dhmf_dynadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_ua_dynadv_cen2, "ua_dynadv_cen2");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6288,31 +6167,30 @@
 	 Yzfuv_UW_dynadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6] );
 	 Yvsmatt (YNBS_zfuv_UW_dynadv_cen2, YNBI_zfuv_UW_dynadv_cen2, YMAX_JAC_NBI, &YG1Y_zfuv_UW_dynadv_cen2, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_wa( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_wa( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_wa( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[2];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[3];
 	 if (1==0 || Yk-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[4];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[5];
 	 if (1==0 || Yk-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[6];
-	 Ydbg_nanf ("B: ", YNBI_zfuv_UW_dynadv_cen2, "zfuv_UW_dynadv_cen2");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6342,5 +6220,4 @@
 	 else 
 	 YG4_hmf_TF_dynadv_cen2( Yi, Yj-1, Yk, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_zva_dhmf_dynadv_cen2, "zva_dhmf_dynadv_cen2");
 	 Yting[0]=YS1_hmf_TF_dynadv_cen2( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yi-1<0)
@@ -6364,5 +6241,4 @@
 	 else 
 	 YG2_hmf_TF_dynadv_cen2( Yi, Yj-1, Yk, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_zua_dhmf_dynadv_cen2, "zua_dhmf_dynadv_cen2");
 
 }
@@ -6391,31 +6267,30 @@
 	 Yhmf_TF_dynadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5] );
 	 Yvsmatt (YNBS_hmf_TF_dynadv_cen2, YNBI_hmf_TF_dynadv_cen2, YMAX_JAC_NBI, &YG1Y_hmf_TF_dynadv_cen2, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[2];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[3];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[4];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[5];
-	 Ydbg_nanf ("B: ", YNBI_hmf_TF_dynadv_cen2, "hmf_TF_dynadv_cen2");
-
-}
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6427,7 +6302,6 @@
 	 Yrhop(Yi, Yj, Yk)->backward( Yting[0] );
 	 Yvsmatt (YNBS_rhop, YNBI_rhop, YMAX_JAC_NBI, &YG1Y_rhop, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_rhd( Yi, Yj, Yk, YTemps) += Ytbeta[0];
-	 Ydbg_nanf ("B: ", YNBI_rhop, "rhop");
 	 Yting[0]=YS1_tb( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_sb( Yi, Yj, Yk, YTemps);
@@ -6437,9 +6311,8 @@
 	 YG1_tb( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 YG1_sb( Yi, Yj, Yk, YTemps) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_rhd, "rhd");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6454,12 +6327,11 @@
 	 Ysb(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
 	 Yvsmatt (YNBS_sb, YNBI_sb, YMAX_JAC_NBI, &YG1Y_sb, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sa( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sa( Yi, Yj, Yk, YTemps) += Ytbeta[2];
 	 YG1_sa_c( Yi, Yj, Yk) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_sb, "sb");
 	 Yting[0]=YS1_tb( Yi, Yj, Yk, YTemps-1);
 	 Yting[1]=YS1_ta( Yi, Yj, Yk, YTemps-1);
@@ -6469,12 +6341,11 @@
 	 Ytb(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
 	 Yvsmatt (YNBS_tb, YNBI_tb, YMAX_JAC_NBI, &YG1Y_tb, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_tb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ta( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ta( Yi, Yj, Yk, YTemps) += Ytbeta[2];
 	 YG1_ta_c( Yi, Yj, Yk) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_tb, "tb");
 
 }
@@ -6502,5 +6373,4 @@
 	 YG1_zwt_trazdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[3];
 	 YG1_sa_c( Yi, Yj, Yk) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_sa, "sa");
 	 Yting[0]=YS2_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -6521,5 +6391,4 @@
 	 YG1_zwt_trazdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[3];
 	 YG1_ta_c( Yi, Yj, Yk) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_ta, "ta");
 
 }
@@ -6549,5 +6418,4 @@
 	 else 
 	 YG1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_sa_de_trazdf_imp, "sa_de_trazdf_imp");
 	 Yting[0]=YS1_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -6571,5 +6439,4 @@
 	 else 
 	 YG1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_ta_de_trazdf_imp, "ta_de_trazdf_imp");
 	 Yting[0]=YS1_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk-1<0)
@@ -6593,8 +6460,7 @@
 	 else 
 	 YG1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_zwt_trazdf_imp, "zwt_trazdf_imp");
-
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6608,7 +6474,6 @@
 	 Yvsmatt (YNBS_zrhs_sa_trazdf_imp, YNBI_zrhs_sa_trazdf_imp, YMAX_JAC_NBI, &YG1Y_zrhs_sa_trazdf_imp, Yjac[0], Ytbeta);
 	 YG1_sa_traldf_lap( Yi, Yj, Yk, YTemps) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_zrhs_sa_trazdf_imp, "zrhs_sa_trazdf_imp");
 	 Yting[0]=YS1_ta_traldf_lap( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_tb( Yi, Yj, Yk, YTemps-1);
@@ -6617,7 +6482,6 @@
 	 Yvsmatt (YNBS_zrhs_trazdf_imp, YNBI_zrhs_trazdf_imp, YMAX_JAC_NBI, &YG1Y_zrhs_trazdf_imp, Yjac[0], Ytbeta);
 	 YG1_ta_traldf_lap( Yi, Yj, Yk, YTemps) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_tb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_zrhs_trazdf_imp, "zrhs_trazdf_imp");
 	 Yting[0]=YS1_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yk+1>YA3_S3dt-1)
@@ -6628,15 +6492,14 @@
 	 Yzwisd_trazdf_imp(Yi, Yj, Yk)->backward( Yting[0], Yting[1] );
 	 Yvsmatt (YNBS_zwisd_trazdf_imp, YNBI_zwisd_trazdf_imp, YMAX_JAC_NBI, &YG1Y_zwisd_trazdf_imp, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps) += Ytbeta[0];
 	 if (1==0 || Yk+1>YA3_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_zwisd_trazdf_imp, "zwisd_trazdf_imp");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6667,5 +6530,4 @@
 	 else 
 	 YG2_zsuv_traldf_lap( Yi, Yj-1, Yk, YTemps) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_sa_traldf_lap, "sa_traldf_lap");
 	 Yting[0]=YS1_ta_ver_traadv_cen2( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_ztuv_traldf_lap( Yi, Yj, Yk, YTemps);
@@ -6691,5 +6553,4 @@
 	 else 
 	 YG2_ztuv_traldf_lap( Yi, Yj-1, Yk, YTemps) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_ta_traldf_lap, "ta_traldf_lap");
 	 Yting[0]=YS1_sb( Yi, Yj, Yk, YTemps-1);
 	 if (1==0 || Yi+1>YA1_S3dt-1)
@@ -6704,15 +6565,14 @@
 	 Yzsuv_traldf_lap(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] );
 	 Yvsmatt (YNBS_zsuv_traldf_lap, YNBI_zsuv_traldf_lap, YMAX_JAC_NBI, &YG1Y_zsuv_traldf_lap, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sb( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sb( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_zsuv_traldf_lap, "zsuv_traldf_lap");
 	 Yting[0]=YS1_tb( Yi, Yj, Yk, YTemps-1);
 	 if (1==0 || Yi+1>YA1_S3dt-1)
@@ -6727,19 +6587,18 @@
 	 Yztuv_traldf_lap(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] );
 	 Yvsmatt (YNBS_ztuv_traldf_lap, YNBI_ztuv_traldf_lap, YMAX_JAC_NBI, &YG1Y_ztuv_traldf_lap, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_tb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_tb( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_tb( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_ztuv_traldf_lap, "ztuv_traldf_lap");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6757,11 +6616,10 @@
 	 Yvsmatt (YNBS_sa_ver_traadv_cen2, YNBI_sa_ver_traadv_cen2, YMAX_JAC_NBI, &YG1Y_sa_ver_traadv_cen2, Yjac[0], Ytbeta);
 	 YG1_sa_hor_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_zcens_ver_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[1];
 	 if (1==0 || Yk+1>YA3_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_zcens_ver_traadv_cen2( Yi, Yj, Yk+1, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_sa_ver_traadv_cen2, "sa_ver_traadv_cen2");
 	 Yting[0]=YS1_ta_hor_traadv_cen2( Yi, Yj, Yk, YTemps);
 	 Yting[1]=YS1_zcent_ver_traadv_cen2( Yi, Yj, Yk, YTemps);
@@ -6774,15 +6632,14 @@
 	 Yvsmatt (YNBS_ta_ver_traadv_cen2, YNBI_ta_ver_traadv_cen2, YMAX_JAC_NBI, &YG1Y_ta_ver_traadv_cen2, Yjac[0], Ytbeta);
 	 YG1_ta_hor_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_zcent_ver_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[1];
 	 if (1==0 || Yk+1>YA3_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_zcent_ver_traadv_cen2( Yi, Yj, Yk+1, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_ta_ver_traadv_cen2, "ta_ver_traadv_cen2");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6799,13 +6656,12 @@
 	 Yzcens_ver_traadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] );
 	 Yvsmatt (YNBS_zcens_ver_traadv_cen2, YNBI_zcens_ver_traadv_cen2, YMAX_JAC_NBI, &YG1Y_zcens_ver_traadv_cen2, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_wa( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sa( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 if (1==0 || Yk-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sa( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_zcens_ver_traadv_cen2, "zcens_ver_traadv_cen2");
 	 Yting[0]=YS1_wa( Yi, Yj, Yk, YTemps-1);
 	 Yting[1]=YS1_ta( Yi, Yj, Yk, YTemps-1);
@@ -6817,17 +6673,16 @@
 	 Yzcent_ver_traadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] );
 	 Yvsmatt (YNBS_zcent_ver_traadv_cen2, YNBI_zcent_ver_traadv_cen2, YMAX_JAC_NBI, &YG1Y_zcent_ver_traadv_cen2, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_wa( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ta( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
 	 if (1==0 || Yk-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ta( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_zcent_ver_traadv_cen2, "zcent_ver_traadv_cen2");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -6859,5 +6714,4 @@
 	 YG2_zcens_hor_traadv_cen2( Yi, Yj-1, Yk, YTemps) += Ytbeta[3];
 	 YG2_zcens_hor_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_sa_hor_traadv_cen2, "sa_hor_traadv_cen2");
 	 Yting[0]=YS1_ta_traqsr( Yi, Yj, Yk, YTemps);
 	 if (1==0 || Yi-1<0)
@@ -6883,5 +6737,4 @@
 	 YG2_zcent_hor_traadv_cen2( Yi, Yj-1, Yk, YTemps) += Ytbeta[3];
 	 YG2_zcent_hor_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_ta_hor_traadv_cen2, "ta_hor_traadv_cen2");
 
 }
@@ -6903,19 +6756,18 @@
 	 Yzcens_hor_traadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
 	 Yvsmatt (YNBS_zcens_hor_traadv_cen2, YNBI_zcens_hor_traadv_cen2, YMAX_JAC_NBI, &YG1Y_zcens_hor_traadv_cen2, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sa( Yi, Yj, Yk, YTemps-1) += Ytbeta[2];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sa( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[3];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sa( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_zcens_hor_traadv_cen2, "zcens_hor_traadv_cen2");
 	 Yting[0]=YS1_ua( Yi, Yj, Yk, YTemps-1);
 	 Yting[1]=YS1_va( Yi, Yj, Yk, YTemps-1);
@@ -6932,19 +6784,18 @@
 	 Yzcent_hor_traadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] );
 	 Yvsmatt (YNBS_zcent_hor_traadv_cen2, YNBI_zcent_hor_traadv_cen2, YMAX_JAC_NBI, &YG1Y_zcent_hor_traadv_cen2, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[1];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ta( Yi, Yj, Yk, YTemps-1) += Ytbeta[2];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ta( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[3];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ta( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_zcent_hor_traadv_cen2, "zcent_hor_traadv_cen2");
 
 }
@@ -6965,5 +6816,4 @@
 	 YG1_sbcgyre_qsr_qns( Yi, Yj, YTemps) += Ytbeta[0];
 	 YG1_ta_sa_trasbc( Yi, Yj, YTemps) += Ytbeta[1];
-	 Ydbg_nanf ("B: ", YNBI_ta_traqsr, "ta_traqsr");
 
 }
@@ -6984,5 +6834,4 @@
 	 #pragma omp atomic
 	 YG1_sa( Yi, Yj, 0, YTemps-1) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_ta_sa_trasbc, "ta_sa_trasbc");
 
 }
@@ -7058,5 +6907,4 @@
 	 YG1_rhop( Yi, Yj, 29, YTemps-1) += Ytbeta[29];
 	 YG1_rhop( Yi, Yj, 30, YTemps-1) += Ytbeta[30];
-	 Ydbg_nanf ("B: ", YNBI_zdfmxl_nmln_hmlp_hmlpt, "zdfmxl_nmln_hmlp_hmlpt");
 	 Yting[0]=YS1_zdf_avt_avmu_avmv( Yi, Yj, 0, YTemps);
 	 Yting[1]=YS1_zdf_avt_avmu_avmv( Yi, Yj, 1, YTemps);
@@ -7124,9 +6972,8 @@
 	 YG1_zdf_avt_avmu_avmv( Yi, Yj, 29, YTemps) += Ytbeta[29];
 	 YG1_zdf_avt_avmu_avmv( Yi, Yj, 30, YTemps) += Ytbeta[30];
-	 Ydbg_nanf ("B: ", YNBI_zdfmxl_hmld, "zdfmxl_hmld");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
+
+}
+}
+#pragma omp parallel for num_threads(16)
 for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
@@ -7207,64 +7054,63 @@
 	 if (1==0 || Yi-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi-1, Yj, 29, YTemps-1) += Ytbeta[3];
 	 if (1==0 || Yi-1<0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi-1, Yj+1, 29, YTemps-1) += Ytbeta[4];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj, 29, YTemps-1) += Ytbeta[5];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj+1, 29, YTemps-1) += Ytbeta[6];
 	 if (1==0 || Yi-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi-1, Yj, 1, YTemps-1) += Ytbeta[7];
 	 if (1==0 || Yi-1<0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi-1, Yj+1, 1, YTemps-1) += Ytbeta[8];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj, 1, YTemps-1) += Ytbeta[9];
 	 if (1==0 || Yj+1>YA2_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_ua( Yi, Yj+1, 1, YTemps-1) += Ytbeta[10];
 	 if (1==0 || Yj-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj-1, 29, YTemps-1) += Ytbeta[11];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj, 29, YTemps-1) += Ytbeta[12];
 	 if (1==0 || Yi+1>YA1_S3dt-1 || Yj-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi+1, Yj-1, 29, YTemps-1) += Ytbeta[13];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi+1, Yj, 29, YTemps-1) += Ytbeta[14];
 	 if (1==0 || Yj-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj-1, 1, YTemps-1) += Ytbeta[15];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi, Yj, 1, YTemps-1) += Ytbeta[16];
 	 if (1==0 || Yi+1>YA1_S3dt-1 || Yj-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi+1, Yj-1, 1, YTemps-1) += Ytbeta[17];
 	 if (1==0 || Yi+1>YA1_S3dt-1){} 
 	 else 
-	   // #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_va( Yi+1, Yj, 1, YTemps-1) += Ytbeta[18];
-	 Ydbg_nanf ("B: ", YNBI_zdf_avt_avmu_avmv, "zdf_avt_avmu_avmv");
-
-}
-}
-//#pragma omp parallel for num_threads(16)
-for(YY=YA3_S3dt-1; YY>=0; --YY){ 
+
+}
+}
+#pragma omp parallel for num_threads(16)
+for(YY=YA3_S3dt-1; YY>=0; --YY){
   Yk=YY;
 for(Yj=YA2_S3dt-1; Yj>=0; --Yj)
@@ -7284,17 +7130,16 @@
 	 Ybn2_rn2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
 	 Yvsmatt (YNBS_bn2_rn2, YNBI_bn2_rn2, YMAX_JAC_NBI, &YG1Y_bn2_rn2, Yjac[0], Ytbeta);
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_tb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0];
 	 if (1==0 || Yk-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_tb( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[1];
-	 //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sb( Yi, Yj, Yk, YTemps-1) += Ytbeta[2];
 	 if (1==0 || Yk-1<0){} 
 	 else 
-	   //	 #pragma omp atomic
+	 #pragma omp atomic
 	 YG1_sb( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_bn2_rn2, "bn2_rn2");
 
 }
@@ -7324,5 +7169,4 @@
 	 YG1_sbcgyre_utau_vtau( Yi, Yj, YTemps) += Ytbeta[2];
 	 YG2_sbcgyre_utau_vtau( Yi, Yj, YTemps) += Ytbeta[3];
-	 Ydbg_nanf ("B: ", YNBI_sbcgyre_wndm, "sbcgyre_wndm");
 	 Ysbcgyre_utau_vtau(Yi, Yj)->backward();
 	 Yting[0]=YS1_sbcgyre_emp( Yi, Yj, YTemps);
@@ -7331,5 +7175,4 @@
 	 Yvsmatt (YNBS_sbcgyre_emps, YNBI_sbcgyre_emps, YMAX_JAC_NBI, &YG1Y_sbcgyre_emps, Yjac[0], Ytbeta);
 	 YG1_sbcgyre_emp( Yi, Yj, YTemps) += Ytbeta[0];
-	 Ydbg_nanf ("B: ", YNBI_sbcgyre_emps, "sbcgyre_emps");
 	 Yting[0]=YS1_sbcgyre_emp1( Yi, Yj, YTemps);
 	 Yting[1]=YS1_sbcgyre_zsumemp_zsurf( 31, 21, YTemps);
@@ -7341,5 +7184,4 @@
 	 YG1_sbcgyre_zsumemp_zsurf( 31, 21, YTemps) += Ytbeta[1];
 	 YG2_sbcgyre_zsumemp_zsurf( 31, 21, YTemps) += Ytbeta[2];
-	 Ydbg_nanf ("B: ", YNBI_sbcgyre_emp, "sbcgyre_emp");
 
 }
@@ -7380,5 +7222,4 @@
 	 YG2_sbcgyre_zsumemp_zsurf( 31, Yj-1, YTemps) += Ytbeta[3];
 	 YG1_sbcgyre_emp1( Yi, Yj, YTemps) += Ytbeta[4];
-	 Ydbg_nanf ("B: ", YNBI_sbcgyre_zsumemp_zsurf, "sbcgyre_zsumemp_zsurf");
 	 Ysbcgyre_emp1(Yi, Yj)->backward();
 	 Yting[0]=YS1_tb( Yi, Yj, 0, YTemps-1);
@@ -7387,5 +7228,4 @@
 	 Yvsmatt (YNBS_sbcgyre_qsr_qns, YNBI_sbcgyre_qsr_qns, YMAX_JAC_NBI, &YG1Y_sbcgyre_qsr_qns, Yjac[0], Ytbeta);
 	 YG1_tb( Yi, Yj, 0, YTemps-1) += Ytbeta[0];
-	 Ydbg_nanf ("B: ", YNBI_sbcgyre_qsr_qns, "sbcgyre_qsr_qns");
 
 }
Index: /trunk/src/cost.dat
===================================================================
--- /trunk/src/cost.dat	(revision 40)
+++ /trunk/src/cost.dat	(revision 41)
@@ -2,2 +2,53 @@
  0.0000000000000000e+00
  0.0000000000000000e+00
+ 2.9733593674719423e+03
+ 5.4561410776679122e+02
+ 2.4475352252509134e+02
+ 2.6358555318068841e+02
+ 2.1527711430855038e+02
+ 2.0204126463990215e+02
+ 2.1416662011590833e+02
+ 2.0207208737645496e+02
+ 2.0211049147089304e+02
+ 2.0235703823123197e+02
+ 2.0019969413511828e+02
+ 2.2286623273509434e+02
+ 1.9878349444274735e+02
+ 2.0070292696001030e+02
+ 2.0014258993979595e+02
+ 1.9972801663599338e+02
+ 1.9293537153893590e+02
+ 2.0814944891652880e+02
+ 1.8213155756443817e+02
+ 1.7950672225933866e+02
+ 1.8080233872884580e+02
+ 1.7876027971923725e+02
+ 2.0259497236690703e+02
+ 1.7881103040469617e+02
+ 1.7886996041817454e+02
+ 1.7907554124113310e+02
+ 1.8045371934564810e+02
+ 1.7569372841627006e+02
+ 1.9113648527038990e+02
+ 1.7516955251891332e+02
+ 1.9550146370923514e+02
+ 1.7455342251319595e+02
+ 2.1077987977826763e+02
+ 1.7473743190883036e+02
+ 1.7476659805882767e+02
+ 1.7491607437730161e+02
+ 1.7534570209262407e+02
+ 1.7146968928574532e+02
+ 1.7918664114219453e+02
+ 1.7311809756377912e+02
+ 1.7320818032615955e+02
+ 1.7322190117805096e+02
+ 1.7328346667692077e+02
+ 1.7390681350015083e+02
+ 1.7091883954920155e+02
+ 1.7074257389516592e+02
+ 1.7423942813610816e+02
+ 1.7069031205380722e+02
+ 1.7425456057653693e+02
+ 1.7425631257707835e+02
+ 1.7425753342598696e+02
Index: /trunk/src/solsor_flag.h
===================================================================
--- /trunk/src/solsor_flag.h	(revision 40)
+++ /trunk/src/solsor_flag.h	(revision 41)
@@ -12,5 +12,9 @@
 	3	from	solsor_flag		1  i-1 	j    k  t
 	4	from	solsor_flag		1  i 	j-1  k  t
-} */
+} 
+output 1 : resmax
+oupput 2 : flag (0 -> continue; 1 -> stop)
+
+*/
 
 //jn (current iteration) = Yk+1
@@ -18,7 +22,9 @@
 if (Yk<nmin)
 {
-YS1=0;
+  YS1=0;
+ YS2=0;
 }
-
-if (Yi==0) {
-
+ else {
+   if (Yi==0 & Yj==0) YS1=0
+     
+ } //Yk>=nmin
