Philippe 23/02/2016: lfi2cdf: modif for JPHEXT/=1
[MNH-git_open_source-lfs.git] / readme / obs2mesonh.LISEZMOI
1 #
2 # obs2mesonh
3 #
4 # Intégration des observations 1D,2D,3D dans la grille Mesonh et écriture d'un fichier
5 # diachronique pouvant être utilisé via diaprog ou extractdia (readvar/writevar)
6 #
7 #  1) préparer un(des) fichier(s) ascii contenant toutes les observations
8 #  au format suivant: une obs par ligne, la valeur indéfinie prévue étant 999.0
9 #       lon lat altitude_en_mètres valeur   
10 #    ou lat lon altitude_en_mètres valeur   
11 #
12 #  2) disposer d un fichier diachronique (fic_diachro_avec_zs): les observations seront integrees sur la grille definie dans ce fichier
13 #
14 #  3) initialiser (facultatif si les fichiers sont dans le repertoire courant)
15     export DIROBS=directory_fichiers_obs
16     export DIRLFI=directory_fichier_grille_mesonh
17 #     puis executer
18     obs2mesonh #et répondre aux questions en interactif
19                # Un fichier "dirobs2mnh" consignera toutes vos réponses 
20                # rentrées au clavier
21 #    ou
22     obs2mesonh < dirobs2mnh_execution_precedente
23 #    ou
24     obs2mesonh << eof
25 fic_diachro_avec_zs     # pour initialiser la grille Mesonh et les dates/heures
26 0/1/2/3                 # + ou - de prints de contrôle
27 LL                      # format du fichier d obs (LL=lon lat alt valeur, 
28                         #                          ll=lat lon alt valeur)
29 fichier1_obs            # format lon lat altitude valeur (indef=999.0)
30 nom_nouveau_champ       # nom des obs du fichier1 suivant la logique diaprog
31 unité_nouveau_champ     # chaine de caractères libre
32 1D/2D/3D                # champ créé de type 1D ou 3D ou 2D 
33                         # pour le cas 2D, seul K=1 sera initialisé
34 LL                      # format du fichier d obs (LL=lon lat alt valeur, ll= lat lon alt valeur)
35 fichier2_obs
36 nom_nouveau_champ2
37 unité_nouveau_champ2
38 1D/2D/3D
39 END                     # fin donc fermeture du fichier diachro résultat
40 eof
41 #
42 #  4) le fichier résultat est nommé fic_diachro_avec_zs+'obs':
43 #    sa grille et ses dates/heures sont celles de fic_diachro_avec_zs,
44 #    le champ ZS est celui de fic_diachro_avec_zs
45 #    il contient autant de champs que de fichiers d'obs traités lors de 
46 #   l'exécution de obs2mesonh.
47 #   + un champ ALT_nom_champ si le type=2D : altitude des obs
48 #   + un champ N_nom_champ                 : nombre d obs par point de grille
49 #
50 #     diaprog ou extractdia ou ... sont utilisables
51 #
52 #     Pour les tracés, utiliser la version postérieure à diaprog_LX_020204
53 #                      activer l'option LSPOT=T pour obtenir le tracé de
54 #     toutes les valeurs même si celles-ci sont trés peu denses.
55 #
56
57 #Méthode utilisée:
58 #=================
59 #
60 #  pour chaque obs lue, 
61 #   - recherche du point de grille Mesonh  I,J contenant cette obs.
62 #   - puis recherche du niveau vertical en tenant compte de la 
63 #     grille verticale au point I,J ( grille verticale W si le nom du champ
64 #     commence par la lettre W, grille de masse pour tous les autres noms
65 #   - stockage de l'obs au point de grille (I,J,K)
66 #
67 # Pour les composantes du vent, on suppose que les observations sont
68 # lues en Uzonal et Vméridien et on transforme ces composantes pour les
69 # utiliser dans la grille Mesonh.
70 #
71 #   Moyenne arithmétique des obs contenues dans la même maille du modèle
72 #   Valeur indefinie si aucune obs.
73 #   Pour les obs dont l unité est "dBz" : passage dbZe à Ze pour effectuer
74 # la moyenne puis retour à dbZe pour l'écriture
75 #   Les observations dont l'altitude est inférieure à l'altitude du premier
76 # niveau de modèle sont stockées au niveau k=1 et un message est affiché lors de
77 # l'exécution.
78 #  
79 #  Mise à jour le 02/04/2004
80 #  Mise à jour le 01/03/2005