[1] | 1 | #========================================================== |
---|
| 2 | # fichier de D'INSTRUCTIONS skimpy.i |
---|
| 3 | #========================================================== |
---|
| 4 | echo on |
---|
| 5 | |
---|
| 6 | #========================================================== |
---|
| 7 | !echo "INITIALISATION DU MODELE" |
---|
| 8 | xset_delta dt 1.50e3 |
---|
| 9 | xset_delta dx 5000 |
---|
| 10 | xset_delta dy 5000 |
---|
| 11 | xset_lamda 0.000832 |
---|
| 12 | setstate Alpha 3.14 |
---|
| 13 | |
---|
| 14 | #========================================================== |
---|
| 15 | !echo "TEST des DERIVEES" |
---|
| 16 | #testdf i j k t x pdx ptol [nmmod-ko-maxF] |
---|
| 17 | testdf 10 10 0 2 1 %0.01 0.01 |
---|
| 18 | |
---|
| 19 | #========================================================== |
---|
| 20 | # EXPERIENCE JUMELLE : |
---|
| 21 | # passe avant une fois sur tous les pas de temps |
---|
| 22 | # pour générer des observations avec les paramètres VRAIS |
---|
| 23 | print_time ON |
---|
| 24 | set_modeltime 0 |
---|
| 25 | FORWARD |
---|
| 26 | set_modeltime 0 |
---|
| 27 | print_time OFF |
---|
| 28 | |
---|
| 29 | #===> sauvegarde de la trajectoire calculee (pour voir) |
---|
| 30 | savestate Omega 1 ij 10% A 0 ./Otraj |
---|
| 31 | #y2dsurf Otraj 10 10 ij 1 1 |
---|
| 32 | |
---|
| 33 | #===> récupération des sorties obtenues comme observations |
---|
| 34 | OUTOOBS Omega 1 SZA |
---|
| 35 | |
---|
| 36 | #===> perturbation des paramètres à contrôler (first guess) |
---|
| 37 | setstate Alpha 2.718 |
---|
| 38 | |
---|
| 39 | #---------------------------------------------------------- |
---|
| 40 | # si on a des observations : |
---|
| 41 | # LOADOBS Omega 1 ij SZA A 0 ./Otraj |
---|
| 42 | |
---|
| 43 | #========================================================== |
---|
| 44 | # PREPARATION ET ASSIMILATION |
---|
| 45 | # choix de la fonction de cout |
---|
| 46 | cost lms 0.5 |
---|
| 47 | print_cost ON |
---|
| 48 | |
---|
| 49 | #---------------------------------------------------------- |
---|
| 50 | # soit le RUN standard de YAO : |
---|
| 51 | #=> choix du (ou des) pas de gradient |
---|
| 52 | setepsi Alpha 5000.5 |
---|
| 53 | setepsi Alpha 5.5 |
---|
| 54 | #=>et en route pour le run |
---|
| 55 | set_nbiter 20 |
---|
| 56 | RUN |
---|
| 57 | |
---|
| 58 | # soit le RUN sous M1QN3 .................................. |
---|
| 59 | setstate Alpha 2.718 |
---|
| 60 | setm_impres 5 |
---|
| 61 | setm_io 6 |
---|
| 62 | setm_mode 0 |
---|
| 63 | set_nbiter 5 |
---|
| 64 | setm_nsim 5 |
---|
| 65 | setm_dxmin 1.0e-20 |
---|
| 66 | setm_epsg 1.0e-20 |
---|
| 67 | setm_ddf1 1.0 |
---|
| 68 | runm |
---|
| 69 | #=========================== FIN ========================= |
---|
| 70 | |
---|
| 71 | |
---|