Changeset 114 for trunk/mesoscale/PLOT


Ignore:
Timestamp:
Apr 22, 2011, 8:14:18 PM (14 years ago)
Author:
aslmd
Message:

LMD_MM_MARS: update graphic tools for GW, water cycle + generic graphic stuff (ps_start, map_latlon) + update notes

Location:
trunk/mesoscale/PLOT
Files:
1 added
8 edited
1 moved

Legend:

Unmodified
Added
Removed

  • trunk/mesoscale/PLOT/MINIMAL/fsc/ps_start.pro

    r85 r114  
    183183;            cmd = 'convert +antialias -density 300 ' + ps_struct.filename + ' -resize 25% ' + outfilename
    184184;            cmd = 'convert -density 600 ' + ps_struct.filename + ' -resize 12% ' + outfilename
    185 cmd = 'ps2epsi '+ps_struct.filename+' 2> /dev/null ; gs -sDEVICE=png16m -dTextAlphaBits=4 -sOutputFile='+outfilename+' -dNOPAUSE -dBATCH -dEPSCrop -r400x400 '+basename+'.epsi'
     185cmd = 'ps2epsi '+ps_struct.filename+' 2> /dev/null ; gs -sDEVICE=png16m -dTextAlphaBits=4 -sOutputFile='+outfilename+' -dNOPAUSE -dBATCH -dEPSCrop -r400x400 '+basename+'.epsi > /dev/null 2> /dev/null'
    186186            SPAWN, cmd
    187187        ENDIF
    … …  
    224224   
    225225   ; Change any parameters you feel like changing.
    226    !P.Charsize = 1.75
    227    !P.Thick = 2
    228    !X.Thick = 2
    229    !Y.Thick = 2
     226;;;;;;; AYMERIC
     227;   !P.Charsize = 1.75
     228;   !P.Thick = 2
     229;   !X.Thick = 2
     230;   !Y.Thick = 2
    230231   !Z.Thick = 2
    231232   !P.Symsize = 1.25
    232233   !P.Font = 1
     234
     235;;;;;;; AYMERIC
     236  !P.Charsize = 1.2
     237  !p.charthick = 2.0
     238  !p.thick = 2.0
     239  !x.thick = 2.0
     240  !y.thick = 2.0
     241
     242
    233243   Device, Set_Font='Helvetica', /TT_FONT
    234244

  • trunk/mesoscale/PLOT/MINIMAL/map_latlon.pro

    r91 r114  
    219219loadct,pal,/silent
    220220;reverse_ct
    221 ;if (pal eq 0) then reverse_ct
     221if (pal eq 0) then reverse_ct
    222222contour, what_I_plot, $
    223223        lon,lat, $

  • trunk/mesoscale/PLOT/MINIMAL/profile.pro

    r85 r114  
    5555
    5656if (n_elements(column) eq 0) then column=findgen(n_elements(what_I_plot))
    57 if (n_elements(mention) eq 0) then mention='' else mention='('+mention+')'
     57if (n_elements(mention) eq 0) then mention=''
     58if (mention ne '') then mention='('+mention+')'
    5859
    5960;------------------
    … …  
    124125        ytitle=ytitle,$
    125126        title=title_user,$
    126 ;/ylog
     127;/ylog, $
    127128        subtitle=mention
    128129

  • trunk/mesoscale/PLOT/SPEC/GW/gravitwave2.pro

    r85 r114  
    22;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    33;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    4 file='../EN_COURS/gw/ok_hires_wrfout_d01_9999-09-09_09:00:00_z'
    5 file='../EN_COURS/gw/_wrfout_d01_9999-09-09_09:00:00_z'
    6 file='../LMD_MM_MARS/TESTGW/_wrfout_d01_9999-09-09_09:00:00_z'
    7 ;file='../LMD_MM_MARS/TESTGW/wind15_wrfout_d01_9999-09-09_09:00:00_z'
    8 file='../TESTGW/wind20c_wrfout_d01_9999-09-09_09:00:00_z'
     4        ;;;; KEEP REFERENCES of FIRST PLOTS
     5        ;file='../EN_COURS/gw/ok_hires_wrfout_d01_9999-09-09_09:00:00_z'
     6        ;file='../EN_COURS/gw/_wrfout_d01_9999-09-09_09:00:00_z'
     7        ;file='../LMD_MM_MARS/TESTGW/_wrfout_d01_9999-09-09_09:00:00_z'
     8        ;;file='../LMD_MM_MARS/TESTGW/wind15_wrfout_d01_9999-09-09_09:00:00_z'
     9        ;file='../TESTGW/wind20c_wrfout_d01_9999-09-09_09:00:00_z'
     10folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_3D.157077/'
     11folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_narrowmountain_2D.156878/'
     12;folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_narrowmountain_3D.157025/'
     13folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_narrowmountain_2D.morepoints.160818/'
     14folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_LT15.161544/'
     15file=folder+'wrfout_d01_9999-09-09_09:00:00_z'
    916charvar='W' & charvarc='W'
    10 ;charvar='tk' & charvarc='W'
    11 charvar='tk' & charvarc='W' & cond=1
     17charvar='tk' & charvarc='W'
     18;charvar='tk' & charvarc='W' & cond=1
    1219charvar='tk' & charvarc='tk' & cond=1
    13 ;charvar='tpot' & charvarc='tpot'
     20;;charvar='tpot' & charvarc='tpot'
    1421;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    1522;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    … …  
    4148;
    4249;
    43 for nt=1,15 do begin
     50;for nt=4,10 do begin
     51;for nt=27,27 do begin
     52for nt=1,40 do begin
    4453zefile='gravitwave_'+charvar+'_'+charvarc+'_'+string(nt+100,'(I0)')
    4554if (n_elements(cond) ne 0) then zefile = 'cond_' + zefile
    … …  
    5766;
    5867;
    59 what_I_plot = reform(invar(*,1,*,nt))
    60 overcontour = reform(invarc(*,1,*,nt))
     68s = size(invar) & middle = floor(s[2]/2) & print, 'PLOT at y subs ', middle
     69what_I_plot = reform(invar(*,middle,*,nt))
     70overcontour = reform(invarc(*,middle,*,nt))
    6171   if (n_elements(cond) ne 0) then begin
    6272     column = columnp
    6373     ;column = reform(columnp(*,1,*,1))
    6474     @tempcond.inc
    65      yeye = reform(overplot(*,1,*,nt))
     75     yeye = reform(overplot(*,middle,*,nt))
    6676     w = where(yeye le 0.) & yeye[w] = 0.
    6777     if (charvar eq 'tk') then what_I_plot = what_I_plot - yeye
    68      if (charvarc eq 'tk') then overcontour = overcontour - yeye
     78     if (charvarc eq 'tk') then overcontour = overcontour - yeye 
    6979   endif
    7080xx = findgen(n_elements(what_I_plot(*,0)))
    7181zz = vert / 1000.
    7282
    73 w = where(what_I_plot lt 0.)
    74 w = where(what_I_plot lt -1.)
    75 w = where(what_I_plot lt -3.)
    76 if (w(0) ne -1) then begin
    77 for i=0,n_elements(w)-1 do begin
    78  ind = ARRAY_INDICES(what_I_plot, w(i))
    79  print, 'x '+string(ind(0),'(I0)')+' z '+string(zz(ind(1)),'(I0)')+' t '+string(nt*925./3700.,'(F4.1)'), what_I_plot(ind[0],ind[1])
    80 endfor
    81 endif
     83        ;w = where(what_I_plot lt 0.)
     84        w = where(what_I_plot lt -1.)
     85        ;w = where(what_I_plot lt -3.)
     86        if (w(0) ne -1) then begin
     87        for i=0,n_elements(w)-1 do begin
     88        ind = ARRAY_INDICES(what_I_plot, w(i))
     89        print, 'x '+string(ind(0),'(I0)')+' z '+string(zz(ind(1)),'(I0)')+' t '+string(nt*925./3700.,'(F4.1)'), what_I_plot(ind[0],ind[1])
     90        endfor
     91        endif
    8292;
    8393;
    … …  
    100110;
    101111PS_End, /PNG
     112SPAWN, 'mv '+zefile+'.png '+folder+' ; \rm '+zefile+'.*ps*'
    102113endfor
    103114end

  • trunk/mesoscale/PLOT/SPEC/GW/gravitwave_inc.pro

    r85 r114  
    99title_user    = 'Vertical velocity (m s!U-1!N)'
    1010;;;;;;;; tk
    11 ;minfield_init =  110.
    12 ;maxfield_init =  210.
    13 ;format        = '(I0)'
    14 ;pal           = 33
    15 ;title_user    = 'Temperature (K)'
    16 ;;;;;;;; tk - tsat
    17 minfield_init =  -20.;-4.
    18 maxfield_init =  -10.
     11minfield_init =  110.
     12maxfield_init =  210.
     13format        = '(I0)'
     14pal            = 33
     15title_user    = 'Temperature (K)'
     16;;;;;;; tk - tsat
     17minfield_init =  -4. ;-8.
     18maxfield_init =  4. ;+8.
     19ndiv          =  8     
    1920format        = '(F4.1)'
    2021pal           =  0
    … …  
    3132intervalx = 5.
    3233intervaly = 5.
     34windowy = [0.,130.]
     35windowx = [10., 110.]
     36intervalx = 10.
     37intervaly = 10.
    3338title_axis = ['horizontal coordinate (x5km)','altitude (km)']
    3439poscb=0.55

  • trunk/mesoscale/PLOT/SPEC/GW/gravitwaveprof.pro

    r85 r114  
    77
    88;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    9 path='../TESTGW/'
    10 path='/tmp7/aslmd/'
    11 experiment='mcd_lt16_rcpcst'
    12 ;experiment='mcd_lt16'
    13 ;experiment='short_mcd_lt16_rcpcst'
    14 ;experiment='short_mcd_lt16'
    15 ;experiment='short_mcd_lt16_calcprho'
    16 ;experiment='mcd_lt16_calcprho'
    17 ;experiment='short_mcd_lt16_calcprho_tot'
    18 ;experiment='colder_gcm_add_short'
    19 ;experiment='colder_gcm_add'
    20 experiment=''
    21 ;experiment='wind15'
    22 ;experiment='wind20c'
    23 ;experiment='colder_gcm_add_3.9_fine_rad_mountain'
    24 experiment='test'
    25 file=path+experiment+'_wrfout_d01_9999-09-09_09:00:00_z'
     9        path='../TESTGW/'
     10        path='/tmp7/aslmd/'
     11path='/d5/aslmd/GRAVITWAVE/'
     12        experiment='mcd_lt16_rcpcst'
     13        ;experiment='mcd_lt16'
     14        ;experiment='short_mcd_lt16_rcpcst'
     15        ;experiment='short_mcd_lt16'
     16        ;experiment='short_mcd_lt16_calcprho'
     17        ;experiment='mcd_lt16_calcprho'
     18        ;experiment='short_mcd_lt16_calcprho_tot'
     19        ;experiment='colder_gcm_add_short'
     20        ;experiment='colder_gcm_add'
     21        experiment=''
     22        ;experiment='wind15'
     23        ;experiment='wind20c'
     24        ;experiment='colder_gcm_add_3.9_fine_rad_mountain'
     25        experiment='test'
     26experiment='GW_MARS_highwind_3D.157077'
     27experiment='GW_MARS_highwind_3D_lower_opacity.161436'
     28experiment='GW_MARS_highwind_3D_loweropacity_lotspoints.161468'
     29experiment='GW_MARS_highwind_3D_loweropacity_morepoints.161542'
     30experiment='GW_MARS_highwind_3D_loweropacity_morepoints_LT15.161544'
     31        file=path+experiment+'_wrfout_d01_9999-09-09_09:00:00_z'
     32file=path+experiment+'/wrfout_d01_9999-09-09_09:00:00_z'
     33;file=path+experiment+'/wrfout_d01_9999-09-09_10:00:00_z'
     34;file=path+experiment+'/wrfout_d01_9999-09-09_11:00:00_z'
    2635charvar='W' & charvarc='W'
    2736charvar='tk' & charvarc='W'
    … …  
    4352nts=1  & nte=9    & nxs=59  & nxe=59
    4453;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     54nts=1  & nte=10  & nxs=50   & nxe=50
     55;nts=1  & nte=1  & nxs=6   & nxe=6
     56nts=1  & nte=100  & nxs=50   & nxe=50
     57;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    4558;
    4659;
    … …  
    5568        charvar=charvar, $
    5669        invar=invar
    57 getcdf, $
    58         file=file, $
    59         charvar=charvarc, $
    60         invar=invarc
     70;getcdf, $
     71;        file=file, $
     72;        charvar=charvarc, $
     73;        invar=invarc
    6174getcdf, $
    6275        file=file, $
    … …  
    6780        charvar='PTOT', $
    6881        invar=columnp
    69 ;
    70 ;
    71 ;
     82;;
     83;;
     84;;
     85        s = size(invar) & middle = floor(s[2]/2) & print, 'PLOT at y subs ', middle
     86        nte = min(s[4]-1,nte) & print, 'nte ', nte
     87        mean_invar = TOTAL(invar(*,*,*,nts:nte),4) / float(nte-nts+1)
     88        ;mean_columnp = TOTAL(columnp(*,*,*,nts:nte),4) / float(nte-nts+1)
     89        anomal_invar = invar & for i=0,s[4]-1 do anomal_invar(*,*,*,i) = anomal_invar(*,*,*,i) - mean_invar
     90;
     91;
     92;
     93        ;goto, perturb
    7294for nt=nts,nte do begin
    7395for nx=nxs,nxe do begin
    … …  
    87109;
    88110;
    89 ;what_I_plot = reform(invar(*,1,*,nt))
    90 what_I_plot = reform(invar(*,59,*,nt))
     111;;what_I_plot = reform(invar(*,1,*,nt))
     112;what_I_plot = reform(invar(*,59,*,nt))
     113what_I_plot = reform(invar(*,middle,*,nt))
    91114
    92115   ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    … …  
    112135   endif else begin
    113136     what_I_plot = TOTAL(what_I_plot,1) / n_elements(reform(what_I_plot(*,0)))
    114    endelse
     137   endelse
     138
    115139     ;;;;;;;;;;;
    116140     ;;;;;;;;;;;
    … …  
    121145     alt = [40.,100.]
    122146     alt = [40.,80.]
    123         ;hache = 8.
    124         ;column = -hache*alog(reform(columnp(nx,1,*,nt))/610.)
    125         ;overplot_column = -hache*alog(pressure/610.)
    126 ;column = reform(columnp(nx,1,*,nt))
     147     alt = [0.,120.]
     148     ;;;;;;;;;;;
     149     ;;;;;;;;;;;
     150     ;mention = ''
     151     SPAWN, '\rm param_plot.idl'
     152     minfield_init = 85.  ;100.
     153     maxfield_init = 185. ;250.
     154     title_axis = ['Temperature (K)','Altitude (km)']
     155     overplot = mcd_tpot
     156     ;;;;;;;;;;;
     157     SPAWN, "echo 'intervaly=10.' >> param_plot.idl"
     158     SPAWN, "echo 'intervalx=5.' >> param_plot.idl"
     159     ;;;;;;;;;;;
     160
     161        ;overplot = reform(mean_invar(nx,middle,*))
     162        ;overplot_column = column
     163
     164;;;;;;;;;;;;;;;;;;;;;;;;;;;; PLOT PRESSION mettre profile en log
     165;column = reform(columnp(nx,middle,*,nt))
    127166;overplot_column = pressure
     167;       ;;overplot = mean_invar
     168;       ;;overplot_column = mean_columnp
    128169;alt = [1.e2,1.e-3]
    129 ;;;; mettre profile en log
    130      ;;;;;;;;;;;
    131      ;;;;;;;;;;;
    132      mention = ''
    133      SPAWN, '\rm param_plot.idl'
    134      minfield_init = 80.  ;100.
    135      maxfield_init = 190. ;250.
    136      title_axis = ['Temperature (K)','Altitude (km)']
    137170;title_axis = ['Temperature (K)','Pressure (Pa)']
    138      overplot = mcd_tpot
    139      ;;;;;;;;;;;
    140 
     171;;;;;;;;;;;;;;;;;;;;;;;;;;;;; PLOT PRESSION
     172
     173;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    141174;; calculer la CAPE ? Tsat - Tenv / Tenv
    142175
    … …  
    169202        mention=mention                         ; add text precision within the plot window (default is nothing or '')
    170203
    171 ;;; sponge layer
    172 ;w=where(FINITE(overplot) ne 0)
    173 ;oplot, [minfield_init,maxfield_init], [max(overplot_column[w])-50.,max(overplot_column[w])-50.], linestyle=1
    174 
    175 ;column = columnp
    176 ;@tempcond.inc
    177 ;yeye = reform(overplot(nx,1,*,nt))
    178 ;;w = where(yeye le 0.) & yeye[w] = !VALUES.F_NAN
    179 ;;oplot, yeye, vert/1000., linestyle=1
     204        ;;; sponge layer
     205        ;w=where(FINITE(overplot) ne 0)
     206        ;oplot, [minfield_init,maxfield_init], [max(overplot_column[w])-50.,max(overplot_column[w])-50.], linestyle=1
     207
     208        ;column = columnp
     209        ;@tempcond.inc
     210        ;yeye = reform(overplot(nx,1,*,nt))
     211        ;;w = where(yeye le 0.) & yeye[w] = !VALUES.F_NAN
     212        ;;oplot, yeye, vert/1000., linestyle=1
    180213
    181214@tempcond.inc
    182 overplot = reform(overplot(nx,1,*,nt))
     215overplot = reform(overplot(nx,middle,*,nt))
    183216oplot, overplot, overplot_column, linestyle=1
    184217
    185 ;calc_tk = calc_tpot * (pressure/610.)^(1/4.4)
    186 ;oplot, calc_tk, column, linestyle=3
     218        ;calc_tk = calc_tpot * (pressure/610.)^(1/4.4)
     219        ;oplot, calc_tk, column, linestyle=3
    187220
    188221PS_End, /PNG
    189222endfor
    190223endfor
     224
     225perturb:
     226PS_Start, filename=experiment+'_'+charvar+'_perturb_'+string(nxs+100,'(I0)')+'.ps'
     227
     228     what_I_plot_loop = reform(anomal_invar(nxs,middle,*,*))
     229
     230     what_I_plot = reform(what_I_plot_loop(*,nts))
     231     column = vert/1000.
     232     ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     233     alt = [0.,120.]
     234     minfield_init = -30.
     235     maxfield_init = 30.
     236     title_axis = ['Temperature anomaly (K)','Altitude (km)']
     237                ;; peut pas changer     
     238                ;SPAWN, "echo 'intervaly=10.' >> param_plot.idl"
     239                ;SPAWN, "echo 'intervalx=5.' >> param_plot.idl"
     240     ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     241     profile, $
     242        what_I_plot, $                          ; 1D vertical profile
     243        column, $                               ; altitudes     
     244        alt=alt, $                              ; altitude range [altmin, altmax]
     245        minfield=minfield_init, $               ; minimum value of plotted field (=0: calculate)
     246        maxfield=maxfield_init, $               ; maximum value of plotted field (=0: calculate)
     247        inprofile=overplot, $                   ; another vertical profile to overplot
     248        incolumn=overplot_column, $             ; altitudes of the other vertical profile (in case /= column)
     249;        discrete=discrete, $                    ; show the profile points (= type of points in !psym)
     250;        title_plot=title_user, $                ; title of the plot ('Profile' is default)
     251        title_axis=title_axis, $                ; title of the [x,y] axis (['Field','Altitude'] is default)
     252        mention=mention                         ; add text precision within the plot window (default is nothing or '')
     253
     254     for iii=nts+1, nte do begin
     255        overplot = reform(what_I_plot_loop(*,iii))
     256        overplot_column = column
     257        oplot, overplot, overplot_column
     258     endfor
     259PS_End, /PNG
     260
    191261end

  • trunk/mesoscale/PLOT/SPEC/SERIES/mawd.pro

    r112 r114  
    22
    33;;;;;;;
    4 ;f_user='/donnees/aslmd/SVN/trunk/mesoscale/TMPDIR/GCMINI/mawd.nc'
    5 ;f_user='/tmp7/aslmd/mawd.nc'
    6 f_user='/d5/aslmd/mawd.nc'
     4experiment='a035'
     5;experiment='a040'
     6;experiment='a045'
     7;;f_user='/donnees/aslmd/SVN/trunk/mesoscale/TMPDIR/GCMINI/mawd.nc'
     8;;f_user='/tmp7/aslmd/mawd.nc'
     9f_user='/d5/aslmd/WATERCYCLE/mawd_'+experiment+'.nc'
    710v_user = 'mtot'
    811x_user = 'Time'
    … …  
    2629SPAWN, '\rm param_plot.idl ; cp mawd_inc.pro param_plot.idl'
    2730
    28 PS_START, file='mawd.ps'
    29   !P.Charsize = 1.2
    30   !p.charthick = 2.0
    31   !p.thick = 2.0
    32   !x.thick = 2.0
    33   !y.thick = 2.0
     31PS_START, file='mawd_'+experiment+'.ps'
     32;  !P.Charsize = 1.2
     33;  !p.charthick = 2.0
     34;  !p.thick = 2.0
     35;  !x.thick = 2.0
     36;  !y.thick = 2.0
    3437
    3538cdfid = ncdf_open(f_user)
    … …  
    4144ncdf_varget, cdfid, varid, yy
    4245
     46varid=ncdf_varid(cdfid,x_user)
     47ncdf_varget, cdfid, varid, xx
     48
    4349if (tes eq 'no') then begin
    4450  ;;; en precip-microns
    4551  champ = champ * 1.e6 / 917.
    46   ;;; prendre Time depuis un autre fichier
    47   cdfid = ncdf_open('/d5/aslmd/gcm_LT14_a035.nc')
    48   varid=ncdf_varid(cdfid,x_user)
    49   ncdf_varget, cdfid, varid, xx
    5052  ;;;; entre 0 et 360
    5153  ;xx = xx MOD 360
     54  xx = xx - 360.*3.
    5255endif else begin
    53   varid=ncdf_varid(cdfid,x_user)
    54   ncdf_varget, cdfid, varid, xx
    5556  xx = xx - 360.*2.
    5657endelse
    … …  
    7677endif
    7778
     79what_I_plot = smooth(what_I_plot,[10,1])
     80
    7881map_latlon, $
    7982        what_I_plot, $                          ; 2D field

  • trunk/mesoscale/PLOT/SPEC/SERIES/mawd_inc.pro

    r112 r114  
    88        format          =       '(I0)'
    99        colors          =       32
    10         pal             =       22;4;33              ;; GOOD: 4, 18, 22, 16, 37, 33, 39, 6, 11
     10        pal             =       33;22;4;33              ;; GOOD: 4, 18, 22, 16, 37, 33, 39, 6, 11
    1111
    1212;;;;;*************************************** FILL LIMITS
    1313        minfield_init   =       0.
    14         maxfield_init   =       65.
    15         ndiv            =       13
     14        maxfield_init   =       100. ;120. ;65.
     15        ndiv            =       10 ;12 ;13
    1616
    1717        ;;;;;*************************************** LIMIT TRICKS
    … …  
    2020                lim_max = maxfield_init & w=where((what_I_plot ge lim_max) and (what_I_plot le 1e9)) & if (w[0] ne -1) then what_I_plot[w]=lim_max
    2121                lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min
    22                 ;lim_blank = 0.1 & w=where(abs(what_I_plot) le lim_blank) & if (w[0] ne -1) then what_I_plot[w]=missing_value
     22                lim_blank = 1.0 & w=where(abs(what_I_plot) le lim_blank) & if (w[0] ne -1) then what_I_plot[w]=missing_value
    2323
    2424;;;;;*************************************** CONTOUR
    … …  
    3838;;;;; MULTIYEAR
    3939        isotropic       =       'false'
    40         windowx         =       [1500.,4000.]
    41         intervalx       =       500.0
     40        windowx         =       [0.,2880.] ;[0.,3600.] ;[1500.,4000.]
     41        intervalx       =       360.0
    4242        windowy         =       [-90.,90.]
    4343        intervaly       =       30.0
    44         lev             =       [100.]
     44        lev             =       [600.] ;[100.]
Note: See TracChangeset for help on using the changeset viewer.