[6] | 1 | #!/usr/bin/env python |
---|
| 2 | # -*- coding: utf-8 -*- |
---|
| 3 | import string |
---|
| 4 | import numpy as np |
---|
| 5 | import matplotlib.pyplot as plt |
---|
| 6 | from pylab import * |
---|
| 7 | from mpl_toolkits.basemap import Basemap |
---|
| 8 | from mpl_toolkits.basemap import shiftgrid, cm |
---|
| 9 | import netCDF4 |
---|
| 10 | |
---|
| 11 | |
---|
| 12 | f11='/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUB_CH' |
---|
| 13 | f2='_ANTARC_JUNE2010.DAT' |
---|
| 14 | amsubch=np.array([1,2,3,4,5]) |
---|
| 15 | amsub_numlines=np.zeros([5],float) |
---|
| 16 | |
---|
| 17 | #for i in range(0,5): |
---|
| 18 | #f=f1+str(amsuach[i])+f2 |
---|
| 19 | #fichier=open(f,'r') |
---|
| 20 | #amsub_numlines[i] = 0 |
---|
| 21 | #for line in fichier: amsub_numlines[i] += 1 |
---|
| 22 | |
---|
| 23 | #fichier.close |
---|
| 24 | |
---|
| 25 | amsub_numlines=np.array([ 24061450., 24061453., 19295174., 19295173., 19295172.]) |
---|
| 26 | |
---|
| 27 | ambch1=np.zeros([20,amsub_numlines[0]],float) |
---|
| 28 | ambch2=np.zeros([20,amsub_numlines[1]],float) |
---|
| 29 | ambch3=np.zeros([20,amsub_numlines[2]],float) |
---|
| 30 | ambch4=np.zeros([20,amsub_numlines[3]],float) |
---|
| 31 | ambch5=np.zeros([20,amsub_numlines[4]],float) |
---|
| 32 | |
---|
| 33 | |
---|
| 34 | #ch1 |
---|
| 35 | i=0 |
---|
| 36 | f=f11+str(amsubch[i])+f2 |
---|
| 37 | fichier=open(f,'r') |
---|
| 38 | a=amsub_numlines[i]-1 |
---|
| 39 | |
---|
| 40 | for iligne in range(0,a): |
---|
| 41 | line=fichier.readline() |
---|
| 42 | liste = line.split() |
---|
| 43 | for j in range(0,19): |
---|
| 44 | ambch1[j,iligne] = float(liste[j]) |
---|
| 45 | |
---|
| 46 | |
---|
| 47 | fichier.close |
---|
| 48 | |
---|
| 49 | |
---|
| 50 | #ch2 |
---|
| 51 | i=1 |
---|
| 52 | f=f11+str(amsubch[i])+f2 |
---|
| 53 | fichier=open(f,'r') |
---|
| 54 | a=amsub_numlines[i]-1 |
---|
| 55 | |
---|
| 56 | for iligne in range(0,a): |
---|
| 57 | line=fichier.readline() |
---|
| 58 | liste = line.split() |
---|
| 59 | for j in range(0,19): |
---|
| 60 | ambch2[j,iligne] = float(liste[j]) |
---|
| 61 | |
---|
| 62 | |
---|
| 63 | fichier.close |
---|
| 64 | |
---|
| 65 | #ch3 |
---|
| 66 | i=2 |
---|
| 67 | f=f11+str(amsubch[i])+f2 |
---|
| 68 | fichier=open(f,'r') |
---|
| 69 | a=amsub_numlines[i]-1 |
---|
| 70 | |
---|
| 71 | for iligne in range(0,a): |
---|
| 72 | line=fichier.readline() |
---|
| 73 | liste = line.split() |
---|
| 74 | for j in range(0,19): |
---|
| 75 | ambch3[j,iligne] = float(liste[j]) |
---|
| 76 | |
---|
| 77 | |
---|
| 78 | fichier.close |
---|
| 79 | |
---|
| 80 | #ch4 |
---|
| 81 | i=3 |
---|
| 82 | f=f11+str(amsubch[i])+f2 |
---|
| 83 | fichier=open(f,'r') |
---|
| 84 | a=amsub_numlines[i]-1 |
---|
| 85 | |
---|
| 86 | for iligne in range(0,a): |
---|
| 87 | line=fichier.readline() |
---|
| 88 | liste = line.split() |
---|
| 89 | for j in range(0,19): |
---|
| 90 | ambch4[j,iligne] = float(liste[j]) |
---|
| 91 | |
---|
| 92 | |
---|
| 93 | fichier.close |
---|
| 94 | |
---|
| 95 | #ch5 |
---|
| 96 | i=4 |
---|
| 97 | f=f11+str(amsubch[i])+f2 |
---|
| 98 | fichier=open(f,'r') |
---|
| 99 | a=amsub_numlines[i]-1 |
---|
| 100 | |
---|
| 101 | for iligne in range(0,a): |
---|
| 102 | line=fichier.readline() |
---|
| 103 | liste = line.split() |
---|
| 104 | for j in range(0,19): |
---|
| 105 | ambch5[j,iligne] = float(liste[j]) |
---|
| 106 | |
---|
| 107 | |
---|
| 108 | fichier.close |
---|
| 109 | |
---|
| 110 | #####################################" |
---|
| 111 | #altitude |
---|
| 112 | x=ambch1[0,:] |
---|
| 113 | y=ambch1[1,:] |
---|
| 114 | z=ambch1[14,:] |
---|
| 115 | m=ambch1[14,:] |
---|
| 116 | z0=0 |
---|
| 117 | z1=1 |
---|
| 118 | y11 = -50 |
---|
| 119 | |
---|
| 120 | import newgrid_zml |
---|
| 121 | mgrid_tab1, mppgrid_tab1, mngrid_tab1, mnngrid_tab1, msigma_grid_tab1, xvec, yvec, m_tab1, mm_tupab1 = newgrid_zml.newgrid(x, y, z, m, vnx1,z0, z1, y11) |
---|
| 122 | |
---|
| 123 | coef=0.6 |
---|
| 124 | cax1=0 |
---|
| 125 | cax2=1 |
---|
| 126 | cax3=0.01 |
---|
| 127 | my_cmap=cm.s3pcpn_l |
---|
| 128 | |
---|
| 129 | mask=mppgrid_tab1.round() |
---|
| 130 | |
---|
| 131 | zz=mask |
---|
| 132 | txt='mask amsub' |
---|
| 133 | outfile='mask_amsub.png' |
---|
| 134 | drawmapf.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 135 | |
---|
| 136 | #################### |
---|
| 137 | #tup |
---|
| 138 | tupamsub_max=np.array([ambch1[18,:].max(),ambch2[18,:].max(),ambch3[18,:].max(),ambch4[18,:].max(),ambch5[18,:].max()]) |
---|
| 139 | |
---|
| 140 | #ch1 |
---|
| 141 | |
---|
| 142 | x=ambch1[0,:] |
---|
| 143 | y=ambch1[1,:] |
---|
| 144 | z=ambch1[18,:] |
---|
| 145 | m=ambch1[14,:] |
---|
| 146 | z0=0 |
---|
| 147 | z1=300 |
---|
| 148 | y11 = -50 |
---|
| 149 | |
---|
| 150 | import newgrid_zml |
---|
| 151 | tupgrid_tab1, tuppgrid_tab1, tupngrid_tab1, tupnngrid_tab1, tupsigma_grid_tab1, xvec, yvec, tup_tab1, mm_tupab1 = newgrid_zml.newgrid(x, y, z, m, vnx1,z0, z1, y11) |
---|
| 152 | |
---|
| 153 | #ch2 |
---|
| 154 | |
---|
| 155 | x=ambch2[0,:] |
---|
| 156 | y=ambch2[1,:] |
---|
| 157 | z=ambch2[18,:] |
---|
| 158 | m=ambch2[14,:] |
---|
| 159 | z0=0 |
---|
| 160 | z1=300 |
---|
| 161 | y11 = -50 |
---|
| 162 | |
---|
| 163 | import newgrid_zml |
---|
| 164 | tupgrid_tab2, tuppgrid_tab2, tupngrid_tab2, tupnngrid_tab2, tupsigma_grid_tab2, xvec, yvec, tup_tab2, mm_tupab2 = newgrid_zml.newgrid(x, y, z, m, vnx1,z0, z1, y11) |
---|
| 165 | |
---|
| 166 | #ch3 |
---|
| 167 | |
---|
| 168 | x=ambch3[0,:] |
---|
| 169 | y=ambch3[1,:] |
---|
| 170 | z=ambch3[18,:] |
---|
| 171 | m=ambch3[14,:] |
---|
| 172 | z0=0 |
---|
| 173 | z1=300 |
---|
| 174 | y11 = -50 |
---|
| 175 | |
---|
| 176 | import newgrid_zml |
---|
| 177 | tupgrid_tab3, tuppgrid_tab3, tupngrid_tab3, tupnngrid_tab3, tupsigma_grid_tab3, xvec, yvec, tup_tab3, mm_tupab3 = newgrid_zml.newgrid(x, y, z, m, vnx1,z0, z1, y11) |
---|
| 178 | |
---|
| 179 | #ch4 |
---|
| 180 | |
---|
| 181 | x=ambch4[0,:] |
---|
| 182 | y=ambch4[1,:] |
---|
| 183 | z=ambch4[18,:] |
---|
| 184 | m=ambch4[14,:] |
---|
| 185 | z0=0 |
---|
| 186 | z1=300 |
---|
| 187 | y11 = -50 |
---|
| 188 | |
---|
| 189 | import newgrid_zml |
---|
| 190 | tupgrid_tab4, tuppgrid_tab4, tupngrid_tab4, tupnngrid_tab4, tupsigma_grid_tab4, xvec, yvec, tup_tab4, mm_tupab4 = newgrid_zml.newgrid(x, y, z, m, vnx1,z0, z1, y11) |
---|
| 191 | |
---|
| 192 | #ch5 |
---|
| 193 | |
---|
| 194 | x=ambch5[0,:] |
---|
| 195 | y=ambch5[1,:] |
---|
| 196 | z=ambch5[18,:] |
---|
| 197 | m=ambch5[14,:] |
---|
| 198 | z0=0 |
---|
| 199 | z1=300 |
---|
| 200 | y11 = -50 |
---|
| 201 | |
---|
| 202 | import newgrid_zml |
---|
| 203 | tupgrid_tab5, tuppgrid_tab5, tupngrid_tab5, tupnngrid_tab5, tupsigma_grid_tab5, xvec, yvec, tup_tab5, mm_tupab5 = newgrid_zml.newgrid(x, y, z, m, vnx1,z0, z1, y11) |
---|