source: trunk/POSTSCRIPT/plein2dessin.pro @ 2

Last change on this file since 2 was 2, checked in by opalod, 22 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 6.0 KB
Line 
1;------------------------------------------------------------
2;------------------------------------------------------------
3;------------------------------------------------------------
4;+
5; NAME:plein2dessin
6;
7; PURPOSE:permet de faire plusieurs dessin sur un Postscript
8;
9; CATEGORY:sorties graphiques
10;
11; CALLING SEQUENCE:plein2dessin
12;
13; INPUTS:
14;
15; KEYWORD PARAMETERS:
16;
17;       /EPS: pour selectionner le format encapsulated postscript
18;
19;       /PORTRAIT: oblige la feuille a etre en position "debout".
20;
21;       /LANDSCAPE: oblige la feuille a etre en position "allongee".
22;       
23;       /TRANS: fait un postscript et l'imprime -si on le desire- sur un
24; transparant
25;
26; OUTPUTS:
27;
28; COMMON BLOCKS:
29;       common.pro
30;
31; SIDE EFFECTS:
32;
33; RESTRICTIONS:
34;
35; EXAMPLE:
36;
37; MODIFICATION HISTORY: Sebastien Masson (smasson@lodyc.jussieu.fr)
38;                       5/7/98
39;-
40;------------------------------------------------------------
41;------------------------------------------------------------
42;------------------------------------------------------------
43pro plein2dessin,PORTRAIT=portrait,LANDSCAPE=landscape, FILENAME = filename
44@common
45;------------------------------------------------------------
46reinitplt, /z,/invert
47if keyword_set(portrait) then key_portrait = 1
48if keyword_set(landscape) then key_portrait = 0
49openps, FILENAME = filename
50;------------------------------------------------------------
51; partie a changer
52;------------------------------------------------------------
53;plt,lec('VOTEMPER'),petit=[2,1,19,10]
54;plt,lec('VOTEMPER'),petit=[2,1,19,10], /rempli
55;pltz,lec('VOTEMPER'),petit=[5,10,16,28],boite=[100,300,-1,1,500],/rempli
56;plt,lec('VOTEMPER') ,petit=[1,2,14,19]
57;plt,lec('VOSALINE') ,petit=[14,2,27,19]
58;plt,1000000.* lec('VORESTO '),0,TITRE='10E6*VORESTO',petit=[1,1,19,10],/rem
59;plt,1000000.* lec('VORESTO '),0,TITRE='10E6*VORESTO',petit=[1,10,19,19],/rem
60;plt,1000000.* lec('VORESTO '),0,TITRE='10E6*VORESTO',petit=[1,19,19,28],/rem
61;placecolor,[1,1,19,10]
62
63
64; p1=pression(nlec('sn',820803,'AT1',dat=820831,dir='t'),nlec('tn',820803,'AT1',dat=820831,dir='t'),grad='x')
65; p2=pression(nlec('sn',820803,'RE2',dat=820831,dir='t'),nlec('tn',820803,'RE2',dat=820831,dir='t'),grad='x')
66;vardate = 'AUG-82'
67;varname = 'gradiant zonal de pression *1000'
68;pltz,p2*1e3 ,boite=[310,370,-.5,.5,300],-1.,1,int=.2 ,zoom=2000,/nocolor,petit=[1,2,20,16], /post 
69;varexp = 'RE1'
70;pltz,p1*1e3 ,boite=[310,370,-.5,.5,300],-1.,1,int=.2,zoom=2000,/nocolor,petit=[1,15,20,29],/ post
71
72;bte = [305, 315, 3, 10]
73;;a = nlec('htoth',810101,'AT1',dat=821231,boite=bte,dir='xy')-nlec('hpycn',810101,'AT1',dat=821231,boite=bte,dir='xy')
74;a = nlec('sss',810101,'AT1',dat=821231,boite=bte,dir='xy')
75;varexp = 'REF, RE1 et RE3'
76;pltt, nlec('sss',810101,'RE3',dat=821231,boite=bte,dir='xy'),sty1d=2 
77;pltt, a, /ov1d
78;pltt, nlec('sss',810101,'ATF',dat=821231,boite=bte,dir='xy'),/ov1d,sty1d=1
79;pltt, nlec('htoth',810101,'RE3',dat=820631,boite=bte,dir='xy')- nlec('hpycn',810101,'RE3',dat=820631,boite=bte,dir='xy'),/ov1d,sty1d=2 
80;pltt, nlec('htoth',810101,'ATF',dat=821231,boite=bte,dir='xy')-nlec('hpycn',810101,'ATF',dat=821231,boite=bte,dir='xy'),/ov1d,sty1d=1
81
82;a = nlecs('utrend6',810115,'RE5',821215,dir='yt',boite=[305,375,-0.5,0.5])*1e6
83;a = nlecs('un',810115,'RE5',821215,dir='yt',boite=[305,375,-0.5,0.5])
84;vardate = '81-82'
85;varname = varname+'*1e6'
86;pltz, a, -1, 1,petit=[1,2,1],/nocolor, nite = 5
87;a = nlecs('utrend6',810115,'RE6',821215,dir='yt',boite=[305,375,-0.5,0.5])*1e6
88;a = nlecs('un',810115,'RE6',821215,dir='yt',boite=[305,375,-0.5,0.5])
89;varname = varname+'*1e6'
90;pltz,a,-1,1,/por,petit=[1,2,2],/nocolor,/noerase , nite = 5
91;placecolor,[4,1.75,17,2.25]
92;plt,(nlecs('vtrend1',810815,'RE5')-reform((nlecs('spgv',810815,'RE5'))[*]#replicate(1,jpk),nxv,nyv,jpk))*1e6,-5,5,/nocolor ,petit=[1,2,1]
93; plt,(nlecs('vtrend1',810815,'RE6')-reform((nlecs('spgv',810815,'RE6'))[*]#replicate(1,jpk),nxv,nyv,jpk))*1e6,-5,5,/nocolor ,petit=[1,2,2]
94
95;plt, nlecs('sst',810715,'RE5'),20,28,int=.5,petit=[1,3,1],/nocolo,/por,/rempli         
96;plt, nlecs('sst',810715,'RE7'),20,28,int=.5,petit=[1,3,2],/nocolo,/por,/rempli,/noerase         
97
98;a = nlecs('sst',810715,'RE5')-nlecs('sst',810715,'RE7')
99;varexp = 'RE5-RE7'
100;set_plot, 'x'
101;lct,  64
102;set_plot, 'ps'
103;plt, a,-1,1,petit=[1,3,3],/nocolo,/por,/rempli,/noerase           
104
105; plt,hdyn(nlecs('sn',810615,'RE5')<1e6,nlecs('tn',810615,'RE5')<1e6)-hdyn(nlecs('sn',810615,'RE6')<1e6,nlecs('tn',810615,'RE6')<1e6) , 0,10,/inv ,/por,petit=[2,2,1],/nocolor,ysurx=2
106; plt,hdyn(nlecs('sn',810715,'RE5')<1e6,nlecs('tn',810715,'RE5')<1e6)-hdyn(nlecs('sn',810715,'RE6')<1e6,nlecs('tn',810715,'RE6')<1e6) , 0,10,/inv ,/por,petit=[2,2,2],/nocolor,ysurx=2,/noerase
107; plt,hdyn(nlecs('sn',810815,'RE5')<1e6,nlecs('tn',810815,'RE5')<1e6)-hdyn(nlecs('sn',810815,'RE6')<1e6,nlecs('tn',810815,'RE6')<1e6) , 0,10,/inv ,/por,petit=[2,2,3],/nocolor,ysurx=2,/noerase
108; plt,hdyn(nlecs('sn',810915,'RE5')<1e6,nlecs('tn',810915,'RE5')<1e6)-hdyn(nlecs('sn',810915,'RE6')<1e6,nlecs('tn',810915,'RE6')<1e6) , 0,10,/inv ,/por,petit=[2,2,4],/nocolor,ysurx=2,/noerase
109;placecolor,[4,15,17,15.5]
110
111
112;a = nlecs('utrend1',810715,'RE5',boite=[305,375,-1, 1],dir='y')*1e6
113;;b = nlecs('spgu',810815,'RE5',boite=[300,330,7,8],dir='y')*1e6
114;;b = b#replicate(1, jpk)
115;;a = a-b
116;;varname = 'gradient, zonal de pression (totale)'
117;pltz, a, -1, 1,/nocolor, /rempli,petit=[1,3,1] 
118
119;c = nlecs('utrend1',810715,'RE7',boite=[305,375,-1, 1],dir='y')*1e6
120;;d = nlecs('spgu',810815,'RE7',boite=[300,330,7,8],dir='y')*1e6
121;;d = d#replicate(1, jpk)
122;;c = c-d
123;;varname = 'gradient, zonal de pression (totale)'
124;pltz,c,-1, 1,/nocolor, /rempli,petit=[1,3,2], /noerase 
125;varexp = 'RE5-RE7'
126;pltz, a-c, -1, 1,/nocolor, /rempli,petit=[1,3,3],/noerase 
127
128p = nlec('gps',810813,'REF')*1000
129plt, p-moyenne(p,'xy') ,/land,petit=[1,2,1],-8,8,/rempli
130 h=hdyn(nlec('sn'),nlec('tn'))
131plt, h-moyenne(h, 'xy'), petit=[1,2,2], /rempli
132;------------------------------------------------------------
133; fin de la partie a changer
134;------------------------------------------------------------
135closeps
136imprime,TRANS=trans
137return
138end
Note: See TracBrowser for help on using the repository browser.