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 | |
---|