| 1 | #-*- coding:UTF-8 -*- |
|---|
| 2 | Case='BOMEX' |
|---|
| 3 | from datetime import datetime |
|---|
| 4 | |
|---|
| 5 | DIRDATA='./' |
|---|
| 6 | |
|---|
| 7 | nomvar='theta' |
|---|
| 8 | |
|---|
| 9 | # couleurs des courbes |
|---|
| 10 | #--------------------- |
|---|
| 11 | # cycle de couleurs |
|---|
| 12 | listcoul=['black','red','green','fuchsia','blue','lime','darkviolet','cyan','darkorange','slateblue','brown','gold'] |
|---|
| 13 | |
|---|
| 14 | # exemple dictionnaire de couleurs en fonction du nom de fichier (sans le prefixe time_ ou prof_) |
|---|
| 15 | #black=reference; blue=bb; violet=turb, orange=pas de temps, Diffusion=vert, resolution=grey, advection=red, resv=brown,domain= microphysique=turquoise |
|---|
| 16 | dicocoul={\ |
|---|
| 17 | 'ensmin_SCM.nc' : 'violet',\ |
|---|
| 18 | 'ensmax_SCM.nc' : 'red',\ |
|---|
| 19 | 'ensavg_SCM.nc' : 'pink',\ |
|---|
| 20 | 'Out_klevel.nc' : 'black',\ |
|---|
| 21 | 'LES0.nc' : 'blue',\ |
|---|
| 22 | 'LES1.nc' : 'blue',\ |
|---|
| 23 | 'LES2.nc' : 'blue',\ |
|---|
| 24 | 'LES3.nc' : 'blue',\ |
|---|
| 25 | 'LES4.nc' : 'blue',\ |
|---|
| 26 | 'LES5.nc' : 'blue',\ |
|---|
| 27 | 'LES6.nc' : 'blue',\ |
|---|
| 28 | 'LES7.nc' : 'blue',\ |
|---|
| 29 | 'LES8.nc' : 'blue',\ |
|---|
| 30 | } |
|---|
| 31 | |
|---|
| 32 | dicostyl={\ |
|---|
| 33 | 'ensmin_SCM.nc' : '-',\ |
|---|
| 34 | 'ensmax_SCM.nc' : '-',\ |
|---|
| 35 | 'ensavg_SCM.nc' : '-',\ |
|---|
| 36 | 'Out_klevel.nc' : '-',\ |
|---|
| 37 | 'LES0.nc' : '-.',\ |
|---|
| 38 | 'LES1.nc' : '-.',\ |
|---|
| 39 | 'LES2.nc' : '-.',\ |
|---|
| 40 | 'LES3.nc' : '-.',\ |
|---|
| 41 | 'LES4.nc' : '-.',\ |
|---|
| 42 | 'LES5.nc' : '-.',\ |
|---|
| 43 | 'LES6.nc' : '-.',\ |
|---|
| 44 | 'LES7.nc' : '-.',\ |
|---|
| 45 | 'LES8.nc' : '-.',\ |
|---|
| 46 | } |
|---|
| 47 | |
|---|
| 48 | vmintab={\ |
|---|
| 49 | 'theta' : 299,\ |
|---|
| 50 | 'wu_shcon' : 0,\ |
|---|
| 51 | 'alphau_shcon':0,\ |
|---|
| 52 | 'mu_shcon':0,\ |
|---|
| 53 | 'u':4,\ |
|---|
| 54 | 'v':0,\ |
|---|
| 55 | 'wth':-0.1,\ |
|---|
| 56 | 'tke':0,\ |
|---|
| 57 | 'th2':0,\ |
|---|
| 58 | 'thu_shcon':299,\ |
|---|
| 59 | 'ql':0,\ |
|---|
| 60 | 'qlu_shcon':0,\ |
|---|
| 61 | 'qt':0,\ |
|---|
| 62 | 'qv':0,\ |
|---|
| 63 | 'qtu_shcon':0,\ |
|---|
| 64 | 'rneb':0,\ |
|---|
| 65 | 'thvu_shcon':299,\ |
|---|
| 66 | } |
|---|
| 67 | |
|---|
| 68 | vmaxtab={\ |
|---|
| 69 | 'theta' : 319,\ |
|---|
| 70 | 'wu_shcon' : 3,\ |
|---|
| 71 | 'alphau_shcon':0.5,\ |
|---|
| 72 | 'mu_shcon':1.,\ |
|---|
| 73 | 'u':16,\ |
|---|
| 74 | 'v':7,\ |
|---|
| 75 | 'wth':0.3,\ |
|---|
| 76 | 'tke':5,\ |
|---|
| 77 | 'th2':1.5,\ |
|---|
| 78 | 'thu_shcon':319,\ |
|---|
| 79 | 'ql':0.00005,\ |
|---|
| 80 | 'qlu_shcon':0.005,\ |
|---|
| 81 | 'qt':0.020,\ |
|---|
| 82 | 'qv':0.020,\ |
|---|
| 83 | 'qtu_shcon':0.020,\ |
|---|
| 84 | 'rneb':1.,\ |
|---|
| 85 | 'thvu_shcon':319,\ |
|---|
| 86 | } |
|---|
| 87 | # listcoul ou dicocoul ou dicohightune ou ... |
|---|
| 88 | #coul1d=listcoul |
|---|
| 89 | styl1d=dicostyl |
|---|
| 90 | coul1d=dicocoul |
|---|
| 91 | |
|---|
| 92 | # liste des fichiers sélectionnés |
|---|
| 93 | listfic=[\ |
|---|
| 94 | 'LES/BOMEX/REF/LES0.nc',\ |
|---|
| 95 | 'LES/BOMEX/REF/LES1.nc',\ |
|---|
| 96 | 'LES/BOMEX/REF/LES2.nc',\ |
|---|
| 97 | 'LES/BOMEX/REF/LES3.nc',\ |
|---|
| 98 | 'LES/BOMEX/REF/LES4.nc',\ |
|---|
| 99 | 'LES/BOMEX/REF/LES5.nc',\ |
|---|
| 100 | 'LES/BOMEX/REF/LES6.nc',\ |
|---|
| 101 | 'LES/BOMEX/REF/LES7.nc',\ |
|---|
| 102 | 'LES/BOMEX/REF/LES8.nc',\ |
|---|
| 103 | 'WAVE1/BOMEX/REF/ensmin_SCM.nc',\ |
|---|
| 104 | 'WAVE1/BOMEX/REF/ensavg_SCM.nc',\ |
|---|
| 105 | 'WAVE1/BOMEX/REF/ensmax_SCM.nc',\ |
|---|
| 106 | '../../../../MUSC/simulations/arp631/CMIP6/L91_300s/BOMEX/REF/Output/netcdf/Out_klevel.nc',\ |
|---|
| 107 | ] |
|---|
| 108 | #'ensmin_A24SC.nc',\ |
|---|
| 109 | #'ensmax_A24SC.nc' |
|---|
| 110 | |
|---|
| 111 | |
|---|
| 112 | # date de la forme : datetime(YYYY,MM,DD,H,M,S) sans 0 à gauche dans les nombres |
|---|
| 113 | dateprof=datetime(1969,06,24,05,30,0) |
|---|
| 114 | datedeb=datetime(1969,06,24,00,0,0) |
|---|
| 115 | datefin=datetime(1969,06,24,06,0,0) |
|---|
| 116 | |
|---|
| 117 | # axe vertical altitude ou pression : 'z' ou 'p' |
|---|
| 118 | axevert='z' |
|---|
| 119 | zmin=0 |
|---|
| 120 | zmax=3000 |
|---|
| 121 | |
|---|
| 122 | niv=500. |
|---|
| 123 | |
|---|
| 124 | # limites valeurs |
|---|
| 125 | #pour Theta |
|---|
| 126 | vmin=299 |
|---|
| 127 | vmax=319 |
|---|
| 128 | #pour wu |
|---|
| 129 | |
|---|
| 130 | #pour alpha_u |
|---|
| 131 | |
|---|
| 132 | #pour mu |
|---|
| 133 | |
|---|
| 134 | #pour thu |
|---|
| 135 | |
|---|
| 136 | #pour u |
|---|
| 137 | |
|---|
| 138 | # pour v |
|---|