#==============================================================================
# 		fichier de D'INSTRUCTIONS sinobad1.i
#Mohamed Berrada 05/03/2009
#==============================================================================

xeos_init
xistate_init 2| init de ta_c et sa_c (GYRE) s'il n y a pas loadstate
xtraqsr_init

#loadstate wa   1 ijk 2 a 1 data_in/start_in_tb/GYRE_sortie_mohamed_pwn_198.dat D
#loadstate ua_c 1 ijk 0 a 1 data_in/start_in_tb/GYRE_sortie_mohamed_pub_198.dat D
#loadstate va_c 1 ijk 0 a 1 data_in/start_in_tb/GYRE_sortie_mohamed_pvb_198.dat D
#loadstate ta_c 1 ijk 0 a 1 data_in/start_in_tb/GYRE_sortie_mohamed_ptb_198.dat D
#loadstate sa_c 1 ijk 0 a 1 data_in/start_in_tb/GYRE_sortie_mohamed_psb_200.dat D

print_time OFF
set_modeltime 0
FORWARD
#savestate mod_name       ns  ijk   pdt  coding  format  file_name
savestate ta  1   ijk   50    A       3       ./data_out/test.pof
xcomparYF
end
print_cost ON
#GOTO TESTAD |TESTLT |TESTAD |TESTDF
outoobs sshb 1 22


GOTO TESTEND |TESTOF
FINTESTDF
#TEST DE LA PROGRAMMATION DES DERIVEES
TESTDF
setstate ua_c 0.00001
testdf 3  3  1  4    f 1.e-10  0.001  10
stop
FINTESTDF
#set_bcoef ua_c .001
#setstate ua_c 0.000001
#outoebx ua_c 1 0

GOTO TESTEND
GOTO TESTEND
#-----------------
# TEST DE LINEAIRE TANGENT
TESTLT
 testlt 1.e-05 1.e-07 1.1 20
stop
#-----------------
#TEST DE L'ADJOINT
TESTAD

testad 0.001
stop
testad 0.01 1.0e-3  1000   1.0e-16
stop
#------------------
#TEST DE LA FONCTION OBJECTIVE
TESTOF
setstate ta_c 100.e-06
testof 1.e-6 1.e-00 10. 10 1
TESTEND
#---------------------------------------------------------------------
#----------------------> choice of the run  <-------------------------
#---------------------------------------------------------------------
#-------> 1. RUNYAO
#-------> 2. RUNM2QN1
#-------> 3. RUNM1QN3


loadstate ta_c 1 ijk 0 a 1 data_in/start_in_tb/GYRE_sortie_mohamed_ptb_198.dat D
GOTO 3
#---------------------------------------------------------------------
3
!echo "RUN OPTIMIZATION WITH M1QN3 ......."
setm_impres  5
setm_io      6
setm_mode    1
set_nbiter   100
setm_nsim    100
setm_dxmin   1.0e-9
setm_epsg    1.0e-15
setm_ddf1    1.
runm

savestate ua_c1  1   ijk   0    A       3       ./uac1.pof
goto FINRUN


FINRUN





end
savestate ua  1   ijk   50    A       3       ./data_out/test_ua.dat
savestate va  1   ijk   50    A       3       ./data_out/test_va.dat
savestate wa  1   ijk   50    A       3       ./data_out/test_wa.dat
savestate ta  1   ijk   50    A       3       ./data_out/test_ta.dat
savestate sa  1   ijk   50    A       3       ./data_out/test_sa.dat
savestate sshn  1   ijk   50    A       3       ./data_out/test_sshn.dat
#savestate wa 1 ijk 0 A 3   ./data_out/testua.pof