#==========================================================
#           fichier de D'INSTRUCTIONS skimpy.i
#==========================================================
echo on

#==========================================================
!echo "INITIALISATION DU MODELE"
xset_delta dt 1.50e3
xset_delta dx 5000
xset_delta dy 5000
xset_lamda 0.000832
setstate Alpha  3.14

#==========================================================
!echo "TEST des DERIVEES"
#testdf  i  j  k  t  x   pdx   ptol [nmmod-ko-maxF]
testdf  10 10  0  2  1  %0.01  0.01

#==========================================================
# EXPERIENCE JUMELLE :
# passe avant une fois sur tous les pas de temps
# pour générer des observations avec les paramètres VRAIS
print_time ON
set_modeltime 0
FORWARD
set_modeltime 0
print_time OFF

#===> sauvegarde de la trajectoire calculee (pour voir)
savestate Omega 1 ij 10% A 0 ./Otraj
#y2dsurf Otraj 10 10 ij 1 1

#===> récupération des sorties obtenues comme observations
OUTOOBS Omega 1 SZA

#===> perturbation des paramètres à contrôler (first guess)
setstate Alpha  2.718

#----------------------------------------------------------
# si on a des observations :
# LOADOBS Omega 1 ij SZA A 0 ./Otraj

#==========================================================
# PREPARATION ET ASSIMILATION
# choix de la fonction de cout
cost lms 0.5
print_cost ON

#----------------------------------------------------------
# soit le RUN standard de YAO :
#=> choix du (ou des) pas de gradient
setepsi Alpha    5000.5
setepsi Alpha    5.5
#=>et en route pour le run
set_nbiter 20
RUN

# soit le RUN sous M1QN3 ..................................
setstate Alpha  2.718
setm_impres  5
setm_io      6
setm_mode    0
set_nbiter   5
setm_nsim    5
setm_dxmin   1.0e-20
setm_epsg    1.0e-20
setm_ddf1    1.0
runm
#=========================== FIN =========================