Changeset 371
- Timestamp:
- 08/07/08 11:32:02 (16 years ago)
- Location:
- trunk/SRC
- Files:
-
- 137 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Calendar/caldat.pro
r327 r371 37 37 ; 38 38 ; @uses 39 ; cm_4cal39 ; <pro>cm_4cal</pro> 40 40 ; 41 41 ; @restrictions -
trunk/SRC/Calendar/date2jul.pro
r329 r371 56 56 ; @examples 57 57 ; 58 ; IDL> jday = juldate(19930124)59 ; IDL> print, date2jul(19931205) EQ julday(12,5,1993)58 ; IDL> jday = juldate(19930124) 59 ; IDL> print, date2jul(19931205) EQ julday(12,5,1993) 60 60 ; 1 61 ; IDL> print, date2jul(931205,/grads) EQ julday(12,5,1993)61 ; IDL> print, date2jul(931205,/grads) EQ julday(12,5,1993) 62 62 ; 1 63 ; IDL> print, date2jul(19931205.5d) EQ julday(12,5,1993,12,0,0)63 ; IDL> print, date2jul(19931205.5d) EQ julday(12,5,1993,12,0,0) 64 64 ; 1 65 ; IDL> print, date2jul(19931205.0d) EQ julday(12,5,1993,0,0,0)65 ; IDL> print, date2jul(19931205.0d) EQ julday(12,5,1993,0,0,0) 66 66 ; 1 67 67 ; -
trunk/SRC/Calendar/date2string.pro
r327 r371 19 19 ; @examples 20 20 ; 21 ; IDL> print, date2string(19900123)21 ; IDL> print, date2string(19900123) 22 22 ; Jan 23, 1990 23 ; IDL> print, date2string(19900123, format = '(C(CMOA))')23 ; IDL> print, date2string(19900123, format = '(C(CMOA))') 24 24 ; JAN 23, 1990 25 25 ; -
trunk/SRC/Calendar/daysinmonth.pro
r325 r371 19 19 ; 20 20 ; @uses 21 ; cm_4cal22 ; updatenew21 ; <pro>cm_4cal</pro> 22 ; <pro>updatenew</pro> 23 23 ; 24 24 ; @examples 25 ; IDL> ndays = daysinmonth(2, 2000) 25 ; 26 ; IDL> ndays = daysinmonth(2, 2000) 26 27 ; 27 28 ; @history -
trunk/SRC/Calendar/jul2date.pro
r327 r371 52 52 ; @examples 53 53 ; 54 ; IDL> print, jul2date(julday(12,23,1999))54 ; IDL> print, jul2date(julday(12,23,1999)) 55 55 ; 19991223 56 ; IDL> print, jul2date(julday(12,23,1999,12,0,0))56 ; IDL> print, jul2date(julday(12,23,1999,12,0,0)) 57 57 ; 19991224. 58 ; IDL> print, jul2date(julday(12,23,1999,12,0,0)), format='(f11.2)'58 ; IDL> print, jul2date(julday(12,23,1999,12,0,0)), format='(f11.2)' 59 59 ; 19991223.50 60 ; IDL> print, jul2date(julday(12,23,1999,0,0,0)), format='(f11.2)'60 ; IDL> print, jul2date(julday(12,23,1999,0,0,0)), format='(f11.2)' 61 61 ; 19991223.00 62 62 ; -
trunk/SRC/Calendar/julday.pro
r327 r371 51 51 ; 52 52 ; @uses 53 ; cm_4cal53 ; <pro>cm_4cal</pro> 54 54 ; 55 55 ; @restrictions -
trunk/SRC/Calendar/leapyr.pro
r327 r371 17 17 ; 18 18 ; @uses 19 ; cm_4cal19 ; <pro>cm_4cal</pro> 20 20 ; 21 21 ; @examples 22 ; IDL> result = leapyr(2000)22 ; IDL> result = leapyr(2000) 23 23 ; 24 24 ; @history -
trunk/SRC/Colors/color24.pro
r327 r371 25 25 ; or the decimal number 65535, type: 26 26 ; 27 ; IDL> color = COLOR24([255, 255, 0]) 28 ; 27 ; IDL> color = COLOR24([255, 255, 0]) 29 28 ; 30 29 ; @history -
trunk/SRC/Colors/colorbar.pro
r358 r371 96 96 ; To display a horizontal color bar above a contour plot, type: 97 97 ; 98 ; IDL> LOADCT, 5, NCOLORS=10099 ; IDL> CONTOUR, DIST(31,41), POSITION=[0.15, 0.15, 0.95, 0.75], $100 ; IDL> C_COLORS=INDGEN(25)*4, NLEVELS=25101 ; IDL> COLORBAR, NCOLORS=10098 ; IDL> LOADCT, 5, NCOLORS=100 99 ; IDL> CONTOUR, DIST(31,41), POSITION=[0.15, 0.15, 0.95, 0.75], $ 100 ; IDL> C_COLORS=INDGEN(25)*4, NLEVELS=25 101 ; IDL> COLORBAR, NCOLORS=100 102 102 ; 103 103 ; @history -
trunk/SRC/Colors/getcolor.pro
r327 r371 114 114 ; 115 115 ; @examples 116 ; 116 117 ; To load a yellow color in color index 100 and plot in yellow, type: 117 118 ; 118 ; 119 ; 119 ; IDL> yellow = GETCOLOR('yellow', 100) 120 ; IDL> PLOT, data, COLOR=yellow 120 121 ; 121 122 ; or, 122 123 ; 123 ; 124 ; IDL> PLOT, data, COLOR=GETCOLOR('yellow', 100) 124 125 ; 125 126 ; To do the same thing on a 24-bit color system with decomposed color on, type: 126 127 ; 127 ; 128 ; IDL> PLOT, data, COLOR=GETCOLOR('yellow', /TRUE) 128 129 ; 129 130 ; or in IDL 5.2 and higher, 130 131 ; 131 ; 132 ; 132 ; IDL> DEVICE, Decomposed=1 133 ; IDL> PLOT, data, COLOR=GETCOLOR('yellow') 133 134 ; 134 135 ; To load all 16 colors into the current color table, starting at 135 136 ; color index 200, type: 136 137 ; 137 ; 138 ; IDL> TVLCT, GETCOLOR(), 200 138 139 ; 139 140 ; To add the color names to a list widget: 140 141 ; 141 ; 142 ; IDL> listID = Widget_List(baseID, Value=GetColor(/Names), YSize=16) 142 143 ; 143 144 ; To load all 16 colors and have the color indices returned in a structure: 144 145 ; 145 ; 146 ; 147 ; 148 ; 146 ; IDL> DEVICE, Decomposed=0 147 ; IDL> colors = GetColor(/Load, Start=1) 148 ; IDL> HELP, colors, /Structure 149 ; IDL> PLOT, data, COLOR=colors.yellow 149 150 ; 150 151 ; To get the direct color values as 24-bit integers in color structure fields: 151 152 ; 152 ; 153 ; 154 ; 153 ; IDL> DEVICE, Decomposed=1 154 ; IDL> colors = GetColor(/Load) 155 ; IDL> PLOT, data, COLOR=colors.yellow 155 156 ; 156 157 ; Note that the START keyword value is ignored if on a 24-bit device, … … 158 159 ; writing code like this: 159 160 ; 160 ; 161 ; 162 ; 163 ; 164 ; 161 ; IDL> colors = GetColor(/Load) 162 ; IDL> PLOT, data, Color=colors.yellow; IDL> DEVICE, Decomposed=0 163 ; IDL> colors = GetColor(/Load, Start=1) 164 ; IDL> HELP, colors, /Structure 165 ; IDL> PLOT, data, COLOR=colors.yellow 165 166 ; 166 167 ; To get the direct color values as 24-bit integers in color structure fields: 167 168 ; 168 ; 169 ; 170 ; 169 ; IDL> DEVICE, Decomposed=1 170 ; IDL> colors = GetColor(/Load) 171 ; IDL> PLOT, data, COLOR=colors.yellow 171 172 ; 172 173 ; Note that the START keyword value is ignored if on a 24-bit device, … … 174 175 ; writing code like this: 175 176 ; 176 ; 177 ; 177 ; IDL> colors = GetColor(/Load) 178 ; IDL> PLOT, data, Color=colors.yellow 178 179 ; 179 180 ; @history -
trunk/SRC/Computation/curl.pro
r327 r371 36 36 ; 37 37 ; @uses 38 ; cm_4cal39 ; cm_4data40 ; cm_4mesh38 ; <pro>cm_4cal</pro> 39 ; <pro>cm_4data</pro> 40 ; <pro>cm_4mesh</pro> 41 41 ; 42 42 ; @restrictions … … 56 56 ; 57 57 ; @examples 58 ; IDL> @tst_initorca2 59 ; IDL> plt, curl(dist(jpi,jpj), dist(jpi,jpj)) 58 ; 59 ; IDL> @tst_initorca2 60 ; IDL> plt, curl(dist(jpi,jpj), dist(jpi,jpj)) 60 61 ; 61 62 ; @history -
trunk/SRC/Computation/div.pro
r327 r371 35 35 ; 36 36 ; @uses 37 ; cm_4cal38 ; cm_4data39 ; cm_4mmesh37 ; <pro>cm_4cal</pro> 38 ; <pro>cm_4data</pro> 39 ; <pro>cm_4mmesh</pro> 40 40 ; 41 41 ; @restrictions … … 55 55 ; 56 56 ; @examples 57 ; IDL> \@tst_initorca2 58 ; IDL> plt, div(dist(jpi,jpj), dist(jpi,jpj)) 57 ; 58 ; IDL> \@tst_initorca2 59 ; IDL> plt, div(dist(jpi,jpj), dist(jpi,jpj)) 59 60 ; 60 61 ; @history -
trunk/SRC/Computation/grad.pro
r327 r371 27 27 ; 28 28 ; @uses 29 ; cm_4cal30 ; cm_4data31 ; cm_4mesh29 ; <pro>cm_4cal</pro> 30 ; <pro>cm_4data</pro> 31 ; <pro>cm_4mesh</pro> 32 32 ; 33 33 ; @restrictions … … 41 41 ; 42 42 ; @examples 43 ; IDL> \@tst_initorca2 44 ; IDL> plt, grad({arr:gphit,g:'T'}, 'x') 45 ; IDL> plt, grad({arr:gphit,g:'T'}, 'y') 43 ; 44 ; IDL> \@tst_initorca2 45 ; IDL> plt, grad({arr:gphit,g:'T'}, 'x') 46 ; IDL> plt, grad({arr:gphit,g:'T'}, 'y') 46 47 ; 47 48 ; @history -
trunk/SRC/Computation/norm.pro
r327 r371 15 15 ; @param VV {in}{required} 16 16 ; Matrix representing the meridional coordinates (at V/U point) of a field of vectors 17 ; A 2D (xy), 3D (xyz or yt), 4D (xyzt) or a structure readable by17 ; A 2D (xy), 3D (xyz or yt), 4D (xyzt) array or a structure readable by 18 18 ; <pro>litchamp</pro> and containing a 2D (xy), 3D (xyz or yt), 4D (xyzt) array. 19 19 ; Note that the dimension of the array must suit the domain dimension. … … 27 27 ; 28 28 ; @returns 29 ; A 2D (xy), 3D (xyz or yt), 4D (xyzt) Array29 ; a 2D (xy), 3D (xyz or yt), 4D (xyzt) array 30 30 ; 31 31 ; @uses 32 ; cm_4mesh33 ; cm_4data34 ; cm_4cal32 ; <pro>cm_4mesh</pro> 33 ; <pro>cm_4data</pro> 34 ; <pro>cm_4cal</pro> 35 35 ; 36 36 ; @restrictions … … 51 51 ; 52 52 ; @examples 53 ; 53 54 ; To calculate the average of the norm of streams on all the domain 54 55 ; between 0 and 50: 55 ; 56 ; 56 ; IDL> domdef, 0, 50 57 ; IDL> res = norm(un, vn, dir = 'xyz') 57 58 ; 58 59 ; @history -
trunk/SRC/ForOldVersion/keep_compatibility.pro
r325 r371 17 17 ; 18 18 ; @uses 19 ; cm_general19 ; <pro>cm_general</pro> 20 20 ; 21 21 ; @examples -
trunk/SRC/ForOldVersion/updatekwd.pro
r325 r371 10 10 ; 11 11 ; @examples 12 ; IDL> @updatekwd13 ; (must be insert as an include in a procedure or function)14 12 ; 15 ; @restrictions the compatibility can be made only if the keyword 16 ; given through _extra are given with their complete name 13 ; must be insert as an include in a procedure or function 14 ; 15 ; IDL> @updatekwd 16 ; 17 ; @restrictions 18 ; the compatibility can be made only if the keywords 19 ; given through _extra are given with their complete name 17 20 ; 18 21 ; @history … … 23 26 ; 24 27 ; @todo 25 ; find the write way to comment this empty file for idldoc26 ; to avoid those error messages27 ; IDLDOC: unknown tag "file_comments" in file updatekwd.pro28 ; IDLDOC: unknown tag "categories" in file updatekwd.pro29 ; IDLDOC: unknown tag "restrictions" in file updatekwd.pro30 28 ; 31 29 ;- -
trunk/SRC/Grid/changemsk.pro
r327 r371 18 18 ; 19 19 ; @examples 20 ; IDL> a = changemsk(tmask[*,*,0]) 20 ; 21 ; IDL> a = changemsk(tmask[*,*,0]) 22 ; 21 23 ; to add ocean points 22 ; IDL> a = 1 - changemsk(1 - tmask[*,*,0])24 ; IDL> a = 1 - changemsk(1 - tmask[*,*,0]) 23 25 ; 24 26 ; @history -
trunk/SRC/Grid/computegrid.pro
r327 r371 222 222 ; 223 223 ; @uses 224 ; cm_4mesh225 ; cm_4data226 ; cm_4cal224 ; <pro>cm_4mesh</pro> 225 ; <pro>cm_4data</pro> 226 ; <pro>cm_4cal</pro> 227 227 ; 228 228 229 229 ; @restrictions 230 230 ; if the grid has x/y periodicity overlap and/or if 231 ; 232 ; 233 ; 234 ; 235 ; 236 ; 231 ; the mask has 0 everywhere at the border (like a closed sea) and 232 ; if (we did not activate /plain and xminmesh, xmaxmesh, yminmesh, 233 ; ymaxmesh keywords are defined to their default values), we redefine 234 ; xminmesh, xmaxmesh, yminmesh, ymaxmesh in order to reove the 235 ; overlapping part and/or to open the domain (avoid it be forced 236 ; to use cell_fill = 1). 237 237 ; 238 238 ; FUV points definition is not exact if the grid is irregular … … 241 241 ; Sebastien Masson (smasson\@lodyc.jussieu.fr) 242 242 ; 2000-04-20 243 ; 244 ; 243 ; Sept 2004, several bug fix to suit C grid type... 244 ; Aug 2005, rewritte almost everything... 245 245 ; 246 246 ; @version -
trunk/SRC/Grid/micromeshmask.pro
r327 r371 52 52 ; 53 53 ; @examples 54 ; IDL> meshdir='/d1fes2-raid2/smasson/DATA/ORCA05/' 55 ; IDL> micromeshmask, 'meshmask_ORCA_R05.nc',iodir=meshdir 54 ; 55 ; IDL> meshdir='/d1fes2-raid2/smasson/DATA/ORCA05/' 56 ; IDL> micromeshmask, 'meshmask_ORCA_R05.nc',iodir=meshdir 56 57 ; 57 58 ;- -
trunk/SRC/Grid/n128gaussian.pro
r325 r371 9 9 ; 10 10 ; @examples 11 ; IDL> lat=n128gaussian() 11 ; 12 ; IDL> lat=n128gaussian() 12 13 ; 13 14 ; @returns -
trunk/SRC/Grid/n160gaussian.pro
r325 r371 9 9 ; 10 10 ; @examples 11 ; IDL> lat=n160gaussian() 11 ; 12 ; IDL> lat=n160gaussian() 12 13 ; 13 14 ; @returns -
trunk/SRC/Grid/n256gaussian.pro
r325 r371 9 9 ; 10 10 ; @examples 11 ; IDL> lat=n256gaussian() 11 ; 12 ; IDL> lat=n256gaussian() 12 13 ; 13 14 ; @returns -
trunk/SRC/Grid/n48gaussian.pro
r325 r371 9 9 ; 10 10 ; @examples 11 ; IDL> lat=n48gaussian() 11 ; 12 ; IDL> lat=n48gaussian() 12 13 ; 13 14 ; @returns -
trunk/SRC/Grid/n80gaussian.pro
r325 r371 9 9 ; 10 10 ; @examples 11 ; IDL> lat=n80gaussian() 11 ; 12 ; IDL> lat=n80gaussian() 12 13 ; 13 14 ; @returns -
trunk/SRC/Grid/ncdf_meshread.pro
r359 r371 8 8 ; 9 9 ; @examples 10 ; IDL> ncdf_meshread [,' filename'] 10 ; 11 ; IDL> ncdf_meshread [,' filename'] 11 12 ; 12 13 ; @param filename {in}{optional}{default='meshmask.nc'}{type=scalar string} … … 59 60 ; 60 61 ; @uses 61 ; cm_4mesh62 ; cm_4data63 ; cm_4cal62 ; <pro>cm_4mesh</pro> 63 ; <pro>cm_4data</pro> 64 ; <pro>cm_4cal</pro> 64 65 ; 65 66 ; @restrictions -
trunk/SRC/Grid/ncdf_meshroms.pro
r327 r371 8 8 ; 9 9 ; @examples 10 ; IDL> ncdf_meshroms [,' filename'] 10 ; 11 ; IDL> ncdf_meshroms [,' filename'] 11 12 ; 12 13 ; @param filename {in}{optional}{default='roms_grd.nc'}{type=scalar string} … … 60 61 ; 61 62 ; @uses 62 ; cm_4mesh63 ; cm_4data64 ; cm_4cal63 ; <pro>cm_4mesh</pro> 64 ; <pro>cm_4data</pro> 65 ; <pro>cm_4cal</pro> 65 66 ; 66 67 ; @restrictions -
trunk/SRC/Grid/restoreboxparam.pro
r325 r371 9 9 ; 10 10 ; @uses 11 ; cm_4mesh (and cm_demomode_usedif we are in demo mode)11 ; <pro>cm_4mesh</pro> (and <pro>cm_demomode</pro> if we are in demo mode) 12 12 ; 13 13 ; @restrictions … … 16 16 ; 17 17 ; @examples 18 ; IDL> restoreboxparam, filename 18 ; 19 ; IDL> restoreboxparam, filename 19 20 ; 20 21 ; @history -
trunk/SRC/Grid/saveboxparam.pro
r325 r371 10 10 ; 11 11 ; @uses 12 ; cm_4mesh (and cm_demomode_usedif we are in demo mode)12 ; <pro>cm_4mesh</pro> (and <pro>cm_demomode</pro> if we are in demo mode) 13 13 ; 14 14 ; @restrictions … … 17 17 ; 18 18 ; @examples 19 ; IDL> saveboxparam, filename 19 ; 20 ; IDL> saveboxparam, filename 20 21 ; 21 22 ; @history -
trunk/SRC/Grid/smallmeshmask.pro
r327 r371 52 52 ; 53 53 ; @examples 54 ; IDL> meshdir='/d1fes2-raid2/smasson/DATA/ORCA05/' 55 ; IDL> smallmeshmask, 'meshmask_ORCA_R05.nc',iodir=meshdir 54 ; 55 ; IDL> meshdir='/d1fes2-raid2/smasson/DATA/ORCA05/' 56 ; IDL> smallmeshmask, 'meshmask_ORCA_R05.nc',iodir=meshdir 56 57 ; 57 58 ; @categories -
trunk/SRC/Interpolation/clickincell.pro
r370 r371 43 43 ; 44 44 ; @examples 45 ; IDL> res = clickincell()45 ; IDL> res = clickincell() 46 46 ; 47 47 ; Click with the left button to select a cell. Clicking one more … … 49 49 ; Click on the right button to quit. 50 50 ; 51 ; IDL> plt, findgen(jpi,jpj),/nodata,map=[90,0,0],/ortho52 ; IDL> print, clickincell(/draw,color=150,/xy)51 ; IDL> plt, findgen(jpi,jpj),/nodata,map=[90,0,0],/ortho 52 ; IDL> print, clickincell(/draw,color=150,/xy) 53 53 ; 54 54 ; @history -
trunk/SRC/Interpolation/cutpar.pro
r327 r371 42 42 ; @examples 43 43 ; 44 ; IDL> x0 = [2,6,2]45 ; IDL> y0 = [0,2,6]46 ; IDL> x1 = [3,8,4]47 ; IDL> y1 = [4,4,6]48 ; IDL> x2 = [1,6,4]49 ; IDL> y2 = [5,6,8]50 ; IDL> x3 = [0,4,2]51 ; IDL> y3 = [1,4,8]52 ; IDL> n = 453 ; IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata54 ; IDL> for i=0,2 do oplot, [x0[i],x1[i],x2[i],x3[i],x0[i]],[y0[i],y1[i],y2[i],y3[i],y0[i]]55 ; IDL> res=cutpar(x0, y0, x1, y1, x2, y2, x3, y3, n)56 ; IDL> for i=0,2 do oplot, [res[0,*,i]], [res[1,*,i]], color = 20+10*i, psym = 1, thick = 344 ; IDL> x0 = [2,6,2] 45 ; IDL> y0 = [0,2,6] 46 ; IDL> x1 = [3,8,4] 47 ; IDL> y1 = [4,4,6] 48 ; IDL> x2 = [1,6,4] 49 ; IDL> y2 = [5,6,8] 50 ; IDL> x3 = [0,4,2] 51 ; IDL> y3 = [1,4,8] 52 ; IDL> n = 4 53 ; IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata 54 ; IDL> for i=0,2 do oplot, [x0[i],x1[i],x2[i],x3[i],x0[i]],[y0[i],y1[i],y2[i],y3[i],y0[i]] 55 ; IDL> res=cutpar(x0, y0, x1, y1, x2, y2, x3, y3, n) 56 ; IDL> for i=0,2 do oplot, [res[0,*,i]], [res[1,*,i]], color = 20+10*i, psym = 1, thick = 3 57 57 ; 58 58 ; @history -
trunk/SRC/Interpolation/cutsegment.pro
r327 r371 30 30 ; 31 31 ; @examples 32 ; IDL> x0=[2,5] 33 ; IDL> y0=[5,1] 34 ; IDL> x1=[9,3] 35 ; IDL> y1=[1,8] 36 ; IDL> res=cutsegment(x0, y0, x1, y1, 10) 37 ; IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata 38 ; IDL> oplot, [x0[0], x1[0]], [y0[0], y1[0]] 39 ; IDL> oplot, [res[0,*,0]], [res[1,*,0]], color = 20, psym = 1, thick = 3 40 ; IDL> oplot, [x0[1], x1[1]], [y0[1], y1[1]] 41 ; IDL> oplot, [res[0,*,1]], [res[1,*,1]], color = 40, psym = 1, thick = 3 32 ; 33 ; IDL> x0=[2,5] 34 ; IDL> y0=[5,1] 35 ; IDL> x1=[9,3] 36 ; IDL> y1=[1,8] 37 ; IDL> res=cutsegment(x0, y0, x1, y1, 10) 38 ; IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata 39 ; IDL> oplot, [x0[0], x1[0]], [y0[0], y1[0]] 40 ; IDL> oplot, [res[0,*,0]], [res[1,*,0]], color = 20, psym = 1, thick = 3 41 ; IDL> oplot, [x0[1], x1[1]], [y0[1], y1[1]] 42 ; IDL> oplot, [res[0,*,1]], [res[1,*,1]], color = 40, psym = 1, thick = 3 42 43 ; 43 44 ; @history -
trunk/SRC/Interpolation/extrapolate.pro
r327 r371 39 39 ; 40 40 ; @examples 41 ; IDL> a=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 42 ; IDL> tvplus, a 43 ; IDL> tvplus, a*(1-tmask[*,*,0]) 41 ; 42 ; IDL> a=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 43 ; IDL> tvplus, a 44 ; IDL> tvplus, a*(1-tmask[*,*,0]) 45 ; 44 46 ; get the coastline: 45 ; IDL> a=extrapolate(tmask[*,*,0],tmask[*,*,0],1,/x_periodic) 46 ; IDL> tvplus, a-tmask[*,*,0] 47 ; 48 ; IDL> a=extrapolate(tmask[*,*,0],tmask[*,*,0],1,/x_periodic) 49 ; IDL> tvplus, a-tmask[*,*,0] 47 50 ; 48 51 ; @history -
trunk/SRC/Interpolation/extrapsmooth.pro
r327 r371 37 37 ; 38 38 ; @examples 39 ; IDL> a=extrapsmooth(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 40 ; IDL> tvplus, a 39 ; 40 ; IDL> a=extrapsmooth(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 41 ; IDL> tvplus, a 42 ; 41 43 ; compare to extrapolate result: 42 ; IDL> b=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 43 ; IDL> tvplus, b, window = 1 44 ; 45 ; IDL> b=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 46 ; IDL> tvplus, b, window = 1 44 47 ; 45 48 ; @history -
trunk/SRC/Interpolation/file_interp.pro
r362 r371 232 232 ; @examples 233 233 ; 234 ; IDL> file_interp, filein, fileout, gridout, inxaxisname = 'lo', inyaxisname = 'la', keep = ['lo', 'la', 'cond_sed'] 235 ; 236 ; IDL> file_interp, in, out, gdout, inuseasmask = 'sst', inmissing_value = -1.00000e+30, missing_value = -1000.00 234 ; IDL> file_interp, filein, fileout, gridout, inxaxisname = 'lo', inyaxisname = 'la', keep = ['lo', 'la', 'cond_sed'] 235 ; IDL> file_interp, in, out, gdout, inuseasmask = 'sst', inmissing_value = -1.00000e+30, missing_value = -1000.00 237 236 ; 238 237 ; @history -
trunk/SRC/Interpolation/fromirr.pro
r327 r371 65 65 ; To interpolate 1 field: 66 66 ; 67 ; IDL> tncep = fromirr('bilinear', topa, glamt, gphit, tmask[*,*,0], lonout, latout, mskout)67 ; IDL> tncep = fromirr('bilinear', topa, glamt, gphit, tmask[*,*,0], lonout, latout, mskout) 68 68 ; 69 69 ; or if you have several fields to interpolate from the same source and target grids … … 72 72 ; (that must be undefined or equal to 0 before calling fromirr) 73 73 ; 74 ; IDL> t1ncep = fromirr('bilinear', topa, glamt, gphit, tmask[*,*,0], lonout, latout, mskout $74 ; IDL> t1ncep = fromirr('bilinear', topa, glamt, gphit, tmask[*,*,0], lonout, latout, mskout $ 75 75 ; , WEIG = a, ADDR = b) 76 ; IDL> help, a, b76 ; IDL> help, a, b 77 77 ; 78 78 ; 2) use a and b that are now defined to bypass the computation of the weights 79 79 ; and addresses and speed-up the computation! 80 80 ; 81 ; IDL> t2ncep = fromirr('bilinear', topa, WEIG = a, ADDR = b)81 ; IDL> t2ncep = fromirr('bilinear', topa, WEIG = a, ADDR = b) 82 82 ; 83 83 ; @history -
trunk/SRC/Interpolation/fromreg.pro
r327 r371 67 67 ; To interpolate 1 field: 68 68 ; 69 ; IDL> topa = fromreg('bilinear', tncep, xncep, yncep, glamt, gphit)69 ; IDL> topa = fromreg('bilinear', tncep, xncep, yncep, glamt, gphit) 70 70 ; 71 71 ; or if you have several fields to interpolate from the same source and target grids … … 74 74 ; (that must be undefined or equal to 0 before calling fromreg 75 75 ; 76 ; IDL> t1opa = fromreg('bilinear', t1ncep, xncep, yncep, glamt, gphit, WEIG = a, ADDR = b)77 ; IDL> help, a, b76 ; IDL> t1opa = fromreg('bilinear', t1ncep, xncep, yncep, glamt, gphit, WEIG = a, ADDR = b) 77 ; IDL> help, a, b 78 78 ; 79 79 ; 2) use a and b that are now defined to bypass the computation of the weights and addresses 80 80 ; and speed-up the computation! 81 81 ; 82 ; IDL> t2opa = fromreg('bilinear', t2ncep, xncep, yncep, glamt, gphit, WEIG = a, ADDR = b)82 ; IDL> t2opa = fromreg('bilinear', t2ncep, xncep, yncep, glamt, gphit, WEIG = a, ADDR = b) 83 83 ; 84 84 ; @history -
trunk/SRC/Interpolation/get_gridparams.pro
r327 r371 14 14 ; 15 15 ; Case 1: 16 ; IDL> get_gridparams, file name/id, lonname, latname, lon, lat, jpi, jpj, n_dimensions16 ; IDL> get_gridparams, file name/id, lonname, latname, lon, lat, jpi, jpj, n_dimensions 17 17 ; 18 18 ; Case 2: 19 ; IDL> get_gridparams, lon, lat, jpi, jpj, n_dimensions19 ; IDL> get_gridparams, lon, lat, jpi, jpj, n_dimensions 20 20 ; 21 21 ; @param in1 {in}{required} -
trunk/SRC/Interpolation/inquad.pro
r327 r371 54 54 ; @examples 55 55 ; 56 ; IDL> x = 1.*[1, 2, 6, 7, 3]57 ; IDL> y = 1.*[1, 3, 3, 4, 7]58 ; IDL> x1 = 1.*[0,4,2]59 ; IDL> y1 = 1.*[1,4,8]60 ; IDL> x2 = 1.*[1,6,4]61 ; IDL> y2 = 1.*[5,6,8]62 ; IDL> x3 = 1.*[3,8,4]63 ; IDL> y3 = 1.*[4,4,6]64 ; IDL> x4 = 1.*[2,6,2]65 ; IDL> y4 = 1.*[0,2,6]66 ; IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata67 ; IDL> for i=0,2 do oplot, [x4[i],x1[i],x2[i],x3[i],x4[i]],[y4[i],y1[i],y2[i],y3[i],y4[i]]68 ; IDL> oplot, x, y, color = 20, psym = 1, thick = 269 ; IDL> print, inquad(x, y, x1, y1, x2, y2, x3, y3, x4, y4)56 ; IDL> x = 1.*[1, 2, 6, 7, 3] 57 ; IDL> y = 1.*[1, 3, 3, 4, 7] 58 ; IDL> x1 = 1.*[0,4,2] 59 ; IDL> y1 = 1.*[1,4,8] 60 ; IDL> x2 = 1.*[1,6,4] 61 ; IDL> y2 = 1.*[5,6,8] 62 ; IDL> x3 = 1.*[3,8,4] 63 ; IDL> y3 = 1.*[4,4,6] 64 ; IDL> x4 = 1.*[2,6,2] 65 ; IDL> y4 = 1.*[0,2,6] 66 ; IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata 67 ; IDL> for i=0,2 do oplot, [x4[i],x1[i],x2[i],x3[i],x4[i]],[y4[i],y1[i],y2[i],y3[i],y4[i]] 68 ; IDL> oplot, x, y, color = 20, psym = 1, thick = 2 69 ; IDL> print, inquad(x, y, x1, y1, x2, y2, x3, y3, x4, y4) 70 70 ; 71 71 ; On a sphere see -
trunk/SRC/Interpolation/inrecgrid.pro
r327 r371 35 35 ; @examples 36 36 ; 37 ; IDL> a=indgen(5)38 ; IDL> b=indgen(7)39 ; IDL> r=inrecgrid([0.25,3.25,2],[4.25,2.8,1.4],a,b)40 ; IDL> print, r37 ; IDL> a=indgen(5) 38 ; IDL> b=indgen(7) 39 ; IDL> r=inrecgrid([0.25,3.25,2],[4.25,2.8,1.4],a,b) 40 ; IDL> print, r 41 41 ; 20 13 7 42 ; IDL> r=inrecgrid([0.25,3.25,2],[4.25,2.8,1.4],a,a+1,b,b+1,/output2d)43 ; IDL> print, r42 ; IDL> r=inrecgrid([0.25,3.25,2],[4.25,2.8,1.4],a,a+1,b,b+1,/output2d) 43 ; IDL> print, r 44 44 ; 0.00000 4.00000 45 45 ; 3.00000 2.00000 -
trunk/SRC/Interpolation/lbcorca.pro
r331 r371 29 29 ; 30 30 ; @examples 31 ; IDL> \@tst_initorca2 32 ; IDL> a = dist(jpiglo, jpjglo) 33 ; IDL> print, lbcorca(a, 'T') 34 ; IDL> dummy = lbcorca(a, 'T', /correction) 35 ; IDL> print, lbcorca(a, 'T') 31 ; 32 ; IDL> \@tst_initorca2 33 ; IDL> a = dist(jpiglo, jpjglo) 34 ; IDL> print, lbcorca(a, 'T') 35 ; IDL> dummy = lbcorca(a, 'T', /correction) 36 ; IDL> print, lbcorca(a, 'T') 36 37 37 38 ; @history -
trunk/SRC/Interpolation/ll_narcs_distances.pro
r327 r371 40 40 ; 41 41 ; @examples 42 ; IDL> Lon_lat0 = [1.0, 2.0]; Initial point specified in radians 43 ; IDL> Arc_Dist = 2.0; Arc distance in radians 44 ; IDL> Az = 1.0; Azimuth in radians 45 ; IDL> Result = LL_ARC_DISTANCE(Lon_lat0, Arc_Dist, Az) 46 ; IDL> PRINT, Result 42 ; 43 ; IDL> Lon_lat0 = [1.0, 2.0]; Initial point specified in radians 44 ; IDL> Arc_Dist = 2.0; Arc distance in radians 45 ; IDL> Az = 1.0; Azimuth in radians 46 ; IDL> Result = LL_ARC_DISTANCE(Lon_lat0, Arc_Dist, Az) 47 ; IDL> PRINT, Result 47 48 ; 2.91415 -0.622234 48 49 ; 49 ; IDL> lon0 = [-10, 20, 100]50 ; IDL> lat0 = [0, -10, 45]51 ; IDL> lon1 = [10, 60, 280]52 ; IDL> lat1 = [0, 10, 45]53 ; IDL> dist = map_npoints(lon0, lat0, lon1, lat1, azimuth = azi, /two_by_two)54 ; IDL> earthradius = 6378206.4d055 ; IDL> res = ll_narcs_distances(lon0, lat0, dist/earthradius, azi, /degrees)56 ; IDL> print, reform(res[0, *])50 ; IDL> lon0 = [-10, 20, 100] 51 ; IDL> lat0 = [0, -10, 45] 52 ; IDL> lon1 = [10, 60, 280] 53 ; IDL> lat1 = [0, 10, 45] 54 ; IDL> dist = map_npoints(lon0, lat0, lon1, lat1, azimuth = azi, /two_by_two) 55 ; IDL> earthradius = 6378206.4d0 56 ; IDL> res = ll_narcs_distances(lon0, lat0, dist/earthradius, azi, /degrees) 57 ; IDL> print, reform(res[0, *]) 57 58 ; 10.000000 60.000000 280.00000 58 ; IDL> print, reform(res[1, *])59 ; IDL> print, reform(res[1, *]) 59 60 ; 1.1999280e-15 10.000000 45.000000 60 61 ; -
trunk/SRC/Interpolation/map_npoints.pro
r327 r371 48 48 ; and P1[i] (in that case, we have np0 = np1 = np) ; if /MIDDLE see this keyword. 49 49 ; @examples 50 ; IDL> print, $51 ; IDL> map_npoints([-105.15,1],[40.02,1],[-0.07,100,50],[51.30,20,0])50 ; IDL> print, $ 51 ; IDL> map_npoints([-105.15,1],[40.02,1],[-0.07,100,50],[51.30,20,0]) 52 52 ; 7551369.3 5600334.8 53 53 ; 12864354. 10921254. 54 54 ; 14919237. 5455558.8 55 55 ; 56 ; IDL> lon0 = [-10, 20, 100]57 ; IDL> lat0 = [0, -10, 45]58 ; IDL> lon1 = [10, 60, 280]59 ; IDL> lat1 = [0, 10, 45]60 ; IDL> dist = map_npoints(lon0, lat0, lon1, lat1, AZIMUTH = azi)61 ; IDL> help, dist, azi56 ; IDL> lon0 = [-10, 20, 100] 57 ; IDL> lat0 = [0, -10, 45] 58 ; IDL> lon1 = [10, 60, 280] 59 ; IDL> lat1 = [0, 10, 45] 60 ; IDL> dist = map_npoints(lon0, lat0, lon1, lat1, AZIMUTH = azi) 61 ; IDL> help, dist, azi 62 62 ; DIST DOUBLE = Array[3, 3] 63 63 ; AZI DOUBLE = Array[3, 3] 64 ; IDL> print, dist[4*lindgen(3)], azi[4*lindgen(3)]64 ; IDL> print, dist[4*lindgen(3)], azi[4*lindgen(3)] 65 65 ; 2226414.0 4957944.5 10018863. 66 66 ; 90.000000 64.494450 4.9615627e-15 67 ; IDL> dist = map_npoints(lon0, lat0, lon1, lat1, AZIMUTH = azi, /TWO_BY_TWO)68 ; IDL> help, dist, azi67 ; IDL> dist = map_npoints(lon0, lat0, lon1, lat1, AZIMUTH = azi, /TWO_BY_TWO) 68 ; IDL> help, dist, azi 69 69 ; DIST DOUBLE = Array[3] 70 70 ; AZI DOUBLE = Array[3] 71 ; IDL> print, dist, azi71 ; IDL> print, dist, azi 72 72 ; 2226414.0 4957944.5 10018863. 73 73 ; 90.000000 64.494450 4.9615627e-15 74 ; IDL> print, map_2points(lon0[0], lat0[0], lon1[0], lat1[0])74 ; IDL> print, map_2points(lon0[0], lat0[0], lon1[0], lat1[0]) 75 75 ; 20.000000 90.000000 76 ; IDL> print, map_npoints(lon0[0], lat0[0], lon1[0], lat1[0], AZIMUTH=azi)/6378206.4d0 / !dtor, azi76 ; IDL> print, map_npoints(lon0[0], lat0[0], lon1[0], lat1[0], AZIMUTH=azi)/6378206.4d0 / !dtor, azi 77 77 ; 20.000000 78 78 ; 90.000000 79 79 ; 80 ; IDL> lon0 = [-10, 20, 100]81 ; IDL> lat0 = [0, -10, 45]82 ; IDL> lon1 = [10, 60, 280]83 ; IDL> lat1 = [0, 10, 45]84 ; IDL> mid = map_npoints(lon0, lat0, lon1, lat1, /MIDDLE, /TWO_BY_TWO)85 ; IDL> print, reform(mid[0,*]), reform(mid[1,*])80 ; IDL> lon0 = [-10, 20, 100] 81 ; IDL> lat0 = [0, -10, 45] 82 ; IDL> lon1 = [10, 60, 280] 83 ; IDL> lat1 = [0, 10, 45] 84 ; IDL> mid = map_npoints(lon0, lat0, lon1, lat1, /MIDDLE, /TWO_BY_TWO) 85 ; IDL> print, reform(mid[0,*]), reform(mid[1,*]) 86 86 ; 0.0000000 40.000000 190.00000 87 87 ; 0.0000000 -1.5902773e-15 90.000000 88 ; IDL> print, (map_2points(lon0[0], lat0[0], lon1[0], lat1[0], npath = 3))[*, 1]88 ; IDL> print, (map_2points(lon0[0], lat0[0], lon1[0], lat1[0], npath = 3))[*, 1] 89 89 ; 0.0000000 0.0000000 90 ; IDL> print, (map_2points(lon0[1], lat0[1], lon1[1], lat1[1], npath = 3))[*, 1]90 ; IDL> print, (map_2points(lon0[1], lat0[1], lon1[1], lat1[1], npath = 3))[*, 1] 91 91 ; 40.000000 -1.5902773e-15 92 ; IDL> print, (map_2points(lon0[2], lat0[2], lon1[2], lat1[2], npath = 3))[*, 1]92 ; IDL> print, (map_2points(lon0[2], lat0[2], lon1[2], lat1[2], npath = 3))[*, 1] 93 93 ; 190.00000 90.000000 94 94 ; -
trunk/SRC/Interpolation/neighbor.pro
r327 r371 32 32 ; 33 33 ; @examples 34 ; IDL> print, neighbor(-105.15,40.02,[-0.07,100,50],[51.30,20,0], $ 35 ; IDL> distance=dis) 34 ; 35 ; IDL> print, neighbor(-105.15,40.02,[-0.07,100,50],[51.30,20,0], $ 36 ; distance=dis) 36 37 ; 0 37 ; IDL> print, dis38 ; IDL> print, dis 38 39 ; 105.684 206.125 160.228 39 40 ; -
trunk/SRC/Interpolation/quadrilateral2square.pro
r358 r371 56 56 ; @examples 57 57 ; 58 ; IDL> splot,[0,5],[0,3],/nodata,xstyle=1,ystyle=159 ; IDL> tracegrille, findgen(11)*.1, findgen(11)*.1,color=indgen(12)*2060 ; IDL> xin = (findgen(11)*.1)#replicate(1, 11)61 ; IDL> yin = replicate(1, 11)#(findgen(11)*.1)62 ; IDL> out = square2quadrilateral(2,1,3,0,5,1,2,3, xin, yin)63 ; IDL> tracegrille, reform(out[0,*],11,11), reform(out[1,*],11,11),color=indgen(12)*2058 ; IDL> splot,[0,5],[0,3],/nodata,xstyle=1,ystyle=1 59 ; IDL> tracegrille, findgen(11)*.1, findgen(11)*.1,color=indgen(12)*20 60 ; IDL> xin = (findgen(11)*.1)#replicate(1, 11) 61 ; IDL> yin = replicate(1, 11)#(findgen(11)*.1) 62 ; IDL> out = square2quadrilateral(2,1,3,0,5,1,2,3, xin, yin) 63 ; IDL> tracegrille, reform(out[0,*],11,11), reform(out[1,*],11,11),color=indgen(12)*20 64 64 ; 65 ; IDL> inorg=quadrilateral2square(2,1,3,0,5,1,2,3,out[0,*],out[1,*])66 ; IDL> tracegrille, reform(inorg[0,*],11,11), reform(inorg[1,*],11,11),color=indgen(12)*2065 ; IDL> inorg=quadrilateral2square(2,1,3,0,5,1,2,3,out[0,*],out[1,*]) 66 ; IDL> tracegrille, reform(inorg[0,*],11,11), reform(inorg[1,*],11,11),color=indgen(12)*20 67 67 ; 68 68 ; @history -
trunk/SRC/Interpolation/spl_fstdrv.pro
r325 r371 13 13 ; 14 14 ; @examples 15 ; IDL> y2 = spl_fstdrv(x, y, yscd, x2) 15 ; 16 ; IDL> y2 = spl_fstdrv(x, y, yscd, x2) 16 17 ; 17 18 ; @param x {in}{required} -
trunk/SRC/Interpolation/spl_incr.pro
r327 r371 34 34 ; 35 35 ; @examples 36 ; IDL> n = 100L 37 ; IDL> x = (dindgen(n))^2 38 ; IDL> y = abs(randomn(0, n)) 39 ; IDL> y[n/2:n/2+1] = 0. 40 ; IDL> y[n-n/3] = 0. 41 ; IDL> y[n-n/6:n-n/6+5] = 0. 42 ; IDL> y = total(y, /cumulative, /double) 43 ; IDL> x2 = dindgen((n-1)^2) 44 ; IDL> n2 = n_elements(x2) 45 ; IDL> print, min(y[1:n-1]-y[0:n-2]) LT 0 46 ; IDL> y2 = spl_incr( x, y, x2) 47 ; IDL> splot, x, y, xstyle = 1, ystyle = 1, ysurx=.25, petit = [1, 2, 1], /land 48 ; IDL> oplot, x2, y2, color = 100 49 ; IDL> c = y2[1:n2-1] - y2[0:n2-2] 50 ; IDL> print, min(c) LT 0 51 ; IDL> print, min(c, max = ma), ma 52 ; IDL> splot,c,xstyle=1,ystyle=1, yrange=[-.01,.05], ysurx=.25, petit = [1, 2, 2], /noerase 53 ; IDL> oplot,[0, n_elements(c)], [0, 0], linestyle = 1 36 ; 37 ; IDL> n = 100L 38 ; IDL> x = (dindgen(n))^2 39 ; IDL> y = abs(randomn(0, n)) 40 ; IDL> y[n/2:n/2+1] = 0. 41 ; IDL> y[n-n/3] = 0. 42 ; IDL> y[n-n/6:n-n/6+5] = 0. 43 ; IDL> y = total(y, /cumulative, /double) 44 ; IDL> x2 = dindgen((n-1)^2) 45 ; IDL> n2 = n_elements(x2) 46 ; IDL> print, min(y[1:n-1]-y[0:n-2]) LT 0 47 ; IDL> y2 = spl_incr( x, y, x2) 48 ; IDL> splot, x, y, xstyle = 1, ystyle = 1, ysurx=.25, petit = [1, 2, 1], /land 49 ; IDL> oplot, x2, y2, color = 100 50 ; IDL> c = y2[1:n2-1] - y2[0:n2-2] 51 ; IDL> print, min(c) LT 0 52 ; IDL> print, min(c, max = ma), ma 53 ; IDL> splot,c,xstyle=1,ystyle=1, yrange=[-.01,.05], ysurx=.25, petit = [1, 2, 2], /noerase 54 ; IDL> oplot,[0, n_elements(c)], [0, 0], linestyle = 1 54 55 ; 55 56 ; @history -
trunk/SRC/Interpolation/spl_keep_mean.pro
r327 r371 49 49 ; 12 monthly values of precipitations into daily values: 50 50 ; 51 ; IDL> yr1 = 199052 ; IDL> yr2 = 199253 ; IDL> nyr = yr2-yr1+154 ; IDL> n1 = 12*nyr+155 ; IDL> x = julday(1+findgen(n1), replicate(1, n1) $56 ; IDL> , replicate(yr1, n1), fltarr(n1))57 ; IDL> n2 = 365*nyr + total(leapyr(yr1+indgen(nyr))) + 158 ; IDL> x2 = julday(replicate(1, n2), 1+findgen(n2) $59 ; IDL> , replicate(yr1, n2), fltarr(n2))60 ; IDL> y = abs(randomn(0, n1-1))61 ; IDL> y2 = spl_keep_mean(x, y, x2, /ge0)62 63 ; IDL> print, min(x, max = ma), ma64 ; IDL> print, min(x2, max = ma), ma65 ; IDL> print, vairdate([min(x, max = ma), ma])66 ; IDL> print, total(y*(x[1:n1-1]-x[0:n1-2]))67 ; IDL> print, total(y2*(x2[1:n2-1]-x2[0:n2-2]))51 ; IDL> yr1 = 1990 52 ; IDL> yr2 = 1992 53 ; IDL> nyr = yr2-yr1+1 54 ; IDL> n1 = 12*nyr+1 55 ; IDL> x = julday(1+findgen(n1), replicate(1, n1) $ 56 ; IDL> , replicate(yr1, n1), fltarr(n1)) 57 ; IDL> n2 = 365*nyr + total(leapyr(yr1+indgen(nyr))) + 1 58 ; IDL> x2 = julday(replicate(1, n2), 1+findgen(n2) $ 59 ; IDL> , replicate(yr1, n2), fltarr(n2)) 60 ; IDL> y = abs(randomn(0, n1-1)) 61 ; IDL> y2 = spl_keep_mean(x, y, x2, /ge0) 62 ; 63 ; IDL> print, min(x, max = ma), ma 64 ; IDL> print, min(x2, max = ma), ma 65 ; IDL> print, vairdate([min(x, max = ma), ma]) 66 ; IDL> print, total(y*(x[1:n1-1]-x[0:n1-2])) 67 ; IDL> print, total(y2*(x2[1:n2-1]-x2[0:n2-2])) 68 68 ; 69 69 ; @history -
trunk/SRC/Interpolation/square2quadrilateral.pro
r358 r371 51 51 ; @examples 52 52 ; 53 ; IDL> splot,[0,5],[0,3],/nodata,xstyle=1,ystyle=154 ; IDL> tracegrille, findgen(11)*.1, findgen(11)*.1,color=indgen(12)*2055 ; IDL> xin = (findgen(11)*.1)#replicate(1, 11)56 ; IDL> yin = replicate(1, 11)#(findgen(11)*.1)57 ; IDL> out = square2quadrilateral(2,1,3,0,5,1,2,3, xin, yin)58 ; IDL> tracegrille, reform(out[0,*],11,11), reform(out[1,*],11,11),color=indgen(12)*2053 ; IDL> splot,[0,5],[0,3],/nodata,xstyle=1,ystyle=1 54 ; IDL> tracegrille, findgen(11)*.1, findgen(11)*.1,color=indgen(12)*20 55 ; IDL> xin = (findgen(11)*.1)#replicate(1, 11) 56 ; IDL> yin = replicate(1, 11)#(findgen(11)*.1) 57 ; IDL> out = square2quadrilateral(2,1,3,0,5,1,2,3, xin, yin) 58 ; IDL> tracegrille, reform(out[0,*],11,11), reform(out[1,*],11,11),color=indgen(12)*20 59 59 ; 60 60 ; @history -
trunk/SRC/Matrix/different.pro
r325 r371 30 30 ; @examples 31 31 ; 32 ; IDL> a = [2,4,6,8]33 ; IDL> b = [6,1,3,2]34 ; IDL> different(a,b) = [ 4, 8] ; Elements in A but not in B32 ; IDL> a = [2,4,6,8] 33 ; IDL> b = [6,1,3,2] 34 ; IDL> different(a,b) = [ 4, 8] ; Elements in A but not in B 35 35 ; 36 36 ; @history -
trunk/SRC/Matrix/extrac2.pro
r325 r371 6 6 ; By default, IDL can make extractions of subdomain: 7 7 ; 8 ; 9 ; 8 ; IDL> a=indgen(5,5) 9 ; IDL> print, a 10 10 ; 0 1 2 3 4 11 11 ; 5 6 7 8 9 … … 13 13 ; 15 16 17 18 19 14 14 ; 20 21 22 23 24 15 ; 15 ; IDL> print, a[[0,2],3] 16 16 ; 15 17 17 ; 17 ; IDL> print, a[[0,2],*] 18 18 ; 0 2 19 19 ; 5 7 … … 22 22 ; 20 22 23 23 ; but 24 ; 24 ; IDL> print, a[[0,2],[3,4]] 25 25 ; 15 22 26 26 ; while 27 ; 27 ; IDL> print, extrac2(a,[0,2],[3,4]) 28 28 ; 15 17 29 29 ; 20 22 … … 57 57 ; @examples 58 58 ; I have a dim 2 matrix named A. I want extract a small intersection 59 ; matrix 2d of the line 2, 3 and 7 and of the column 0 and 1:59 ; matrix 2d of the line 2, 3 and 7 and of the column 0 and 1: 60 60 ; 61 ; IDL> res=extrac2(A,[2,3,7],[0,1])61 ; IDL> res=extrac2(A,[2,3,7],[0,1]) 62 62 ; 63 63 ; other ex: 64 ; IDL> print, a64 ; IDL> print, a 65 65 ; a b c 66 66 ; d e f 67 67 ; g h i 68 ; IDL> print, extrac2(a,[0,2],[0,2])68 ; IDL> print, extrac2(a,[0,2],[0,2]) 69 69 ; a c 70 70 ; g i -
trunk/SRC/Matrix/inter.pro
r325 r371 30 30 ; 31 31 ; @examples 32 ; IDL> a = [2,4,6,8] 33 ; IDL> b = [6,1,3,2] 34 ; IDL> inter(a,b) = [ 2, 6] ; Common elements 32 ; 33 ; IDL> a = [2,4,6,8] 34 ; IDL> b = [6,1,3,2] 35 ; IDL> inter(a,b) = [ 2, 6] ; Common elements 35 36 ; 36 37 ; @history -
trunk/SRC/Matrix/union.pro
r325 r371 29 29 ; 30 30 ; @examples 31 ; IDL> a = [2,4,6,8] 32 ; IDL> b = [6,1,3,2] 33 ; IDL> union(a,b) = [ 1, 2, 3, 4, 6, 8] ; Elements in either set 31 ; 32 ; IDL> a = [2,4,6,8] 33 ; IDL> b = [6,1,3,2] 34 ; IDL> union(a,b) = [ 1, 2, 3, 4, 6, 8] ; Elements in either set 34 35 ; 35 36 ; @history -
trunk/SRC/Obsolete/congridseb.pro
r325 r371 6 6 ; Useless now... 7 7 ; 8 ; @examples 9 ; 8 10 ; old example: based on a old version of IDL (5.??) 9 ; IDL> print, congrid([[1,2,3,4],[5,6,7,8]],12,4)11 ; IDL> print, congrid([[1,2,3,4],[5,6,7,8]],12,4) 10 12 ; 1 1 1 2 2 2 3 3 3 3 4 4 11 13 ; 1 1 1 2 2 2 3 3 3 3 4 4 12 14 ; 5 5 5 6 6 6 7 7 7 7 8 8 13 15 ; 5 5 5 6 6 6 7 7 7 7 8 8 14 ; IDL> print, rebin([[1,2,3,4],[5,6,7,8]],12,4)16 ; IDL> print, rebin([[1,2,3,4],[5,6,7,8]],12,4) 15 17 ; 1 1 1 2 2 2 3 3 3 4 4 4 16 18 ; 3 3 3 4 4 4 5 5 5 6 6 6 17 19 ; 5 5 5 6 6 6 7 7 7 8 8 8 18 20 ; 5 5 5 6 6 6 7 7 7 8 8 8 19 ; IDL> print, congridseb([[1,2,3,4],[5,6,7,8]],12,4)21 ; IDL> print, congridseb([[1,2,3,4],[5,6,7,8]],12,4) 20 22 ; 1 1 1 2 2 2 3 3 3 4 4 4 21 23 ; 1 1 1 2 2 2 3 3 3 4 4 4 … … 26 28 ; 27 29 ; @param tableau {in}{required} 28 ; 1 or 2d array30 ; 1d or 2d array 29 31 ; 30 32 ; @param x {in}{required} -
trunk/SRC/Obsolete/extrait.pro
r325 r371 6 6 ; By default, IDL can make extractions of subdomain: 7 7 ; 8 ; 9 ; 8 ; IDL> a=indgen(5,5) 9 ; IDL> print, a 10 10 ; 0 1 2 3 4 11 11 ; 5 6 7 8 9 … … 13 13 ; 15 16 17 18 19 14 14 ; 20 21 22 23 24 15 ; 15 ; IDL> print, a[[0,2],3] 16 16 ; 15 17 17 ; 17 ; IDL> print, a[[0,2],*] 18 18 ; 0 2 19 19 ; 5 7 … … 22 22 ; 20 22 23 23 ; but 24 ; 24 ; IDL> print, a[[0,2],[3,4]] 25 25 ; 15 22 26 26 ; while 27 ; 27 ; IDL> print, extrait(a,[0,2],[3,4]) 28 28 ; 15 17 29 29 ; 20 22 … … 62 62 ; matrix 2d of the line 2,3 and 7 and of the column 0 and 1: 63 63 ; 64 ; IDL> res=extrait(A,[2,3,7],[0,1])64 ; IDL> res=extrait(A,[2,3,7],[0,1]) 65 65 ; 66 66 ; other ex: 67 ; IDL> print, a67 ; IDL> print, a 68 68 ; a b c 69 69 ; d e f 70 70 ; g h i 71 ; IDL> print, extrait(a,[0,2],[0,2])71 ; IDL> print, extrait(a,[0,2],[0,2]) 72 72 ; a c 73 73 ; g i -
trunk/SRC/Obsolete/ficdate.pro
r370 r371 23 23 ; 24 24 ; @examples 25 ; IDL> fic = ficdate(930124) 25 ; 26 ; IDL> fic = ficdate(930124) 26 27 ; 27 28 ; @history -
trunk/SRC/Obsolete/fictype.pro
r370 r371 15 15 ; @examples 16 16 ; 17 ; IDL> fictype = fictype(930124)17 ; IDL> fictype = fictype(930124) 18 18 ; 19 19 ; @history -
trunk/SRC/Obsolete/juldate.pro
r370 r371 32 32 ; 33 33 ; @examples 34 ; IDL> date = juldate(930124) 34 ; 35 ; IDL> date = juldate(930124) 35 36 ; 36 37 ; @history -
trunk/SRC/Obsolete/lec.pro
r370 r371 17 17 ; 18 18 ; @examples 19 ; IDL> resultat=lec('nom_Vairmer'[,date[,'nom_experience']]) 19 ; 20 ; IDL> resultat=lec('nom_Vairmer'[,date[,'nom_experience']]) 20 21 ; 21 22 ; @param nomchamp {in}{required} -
trunk/SRC/Obsolete/lect.pro
r370 r371 15 15 ; 16 16 ; @examples 17 ; IDL> resultat=lec('nom_Vairmer',date1,date2,'nomexp','direc',BOITE=boite) 17 ; 18 ; IDL> resultat=lec('nom_Vairmer',date1,date2,'nomexp','direc',BOITE=boite) 18 19 ; 19 20 ; @param nomchamp {in}{required} -
trunk/SRC/Obsolete/meshlec.pro
r370 r371 9 9 ; 10 10 ; @examples 11 ; IDL> meshmask[,' nomfich'] 11 ; 12 ; IDL> meshmask[,' nomfich'] 12 13 ; 13 14 ; @param nomfich {in}{required} -
trunk/SRC/Obsolete/nlec5j.pro
r370 r371 13 13 ; 14 14 ; @examples 15 ; IDL> res=nlec5j( nom, numsortie, nbretps) 15 ; 16 ; IDL> res=nlec5j( nom, numsortie, nbretps) 16 17 ; 17 18 ; @param nom {in}{required} -
trunk/SRC/Obsolete/nlecan.pro
r370 r371 6 6 ; 7 7 ; @obsolete 8 ; 8 9 ; @categories 9 10 ; Reading 10 11 ; 11 12 ; @examples 12 ; IDL> res=nlecan(nom, numsortie79, nbretps) 13 ; 14 ; IDL> res=nlecan(nom, numsortie79, nbretps) 13 15 ; 14 16 ; @param nom {in}{required} -
trunk/SRC/Obsolete/nlecsaison.pro
r370 r371 11 11 ; 12 12 ; @examples 13 ; IDL> res=nlecsaison(nom, numsortie79, nbretps) 13 ; 14 ; IDL> res=nlecsaison(nom, numsortie79, nbretps) 14 15 ; 15 16 ; @param nom {in}{required} -
trunk/SRC/Obsolete/nlecserie.pro
r370 r371 11 11 ; 12 12 ; @examples 13 ; IDL> res=nlecserie(nom, date1, date2) 13 ; 14 ; IDL> res=nlecserie(nom, date1, date2) 14 15 ; 15 16 ; @param nom {in}{required} -
trunk/SRC/Obsolete/norme.pro
r370 r371 51 51 ; To calculate the average of the norme of streams on all the domain 52 52 ; between 0 and 50: 53 ; 53 ; IDL> res=norme(un,vn,boxzoom=[0,50],dir='xyz') 54 54 ; 55 55 ; @history -
trunk/SRC/Obsolete/strrepl.pro
r325 r371 25 25 ; 26 26 ; @examples 27 ; ; Convert one letter into upper case28 27 ; 29 ; abc = 'abcdefghijklmnopqrstuvwxyz' 30 ; print,strrepl(abc,strpos(abc,'m'),'M') 28 ; Convert one letter into upper case 31 29 ; 32 ; ; prints "abcdefghijklMnopqrstuvwxyz" 30 ; IDL> abc = 'abcdefghijklmnopqrstuvwxyz' 31 ; IDL> print,strrepl(abc,strpos(abc,'m'),'M') 32 ; abcdefghijklMnopqrstuvwxyz 33 33 ; 34 34 ; 35 ; ; Use with strwhere function 36 ; a = 'abcabcabc' 37 ; print,strrepl(a,strwhere(a,'a'),'#') 35 ; Use with strwhere function 36 ; IDL> a = 'abcabcabc' 37 ; IDL> print,strrepl(a,strwhere(a,'a'),'#') 38 ; #bc#bc#bc#bc#bc 38 39 ; 39 ; ; prints "#bc#bc#bc#bc#bc" 40 ; 41 ; IDL> print, strrepl(a,'bc','!eeee!') 40 ; IDL> print, strrepl(a,'bc','!eeee!') 42 41 ; a!eeee!a!eeee!a!eeee! 43 ; 42 ; IDL> print, strrepl(a,'b','0000') 44 43 ; a0000ca0000ca0000 45 ; 44 ; IDL> print, strrepl(a,'toto','0000') 46 45 ; abcabcabc 47 46 ; -
trunk/SRC/Obsolete/vairdate.pro
r370 r371 10 10 ; 11 11 ; @examples 12 ; IDL> vdate = vairdate(1755087) 12 ; 13 ; IDL> vdate = vairdate(1755087) 13 14 ; 14 15 ; @param jdate {in}{required} -
trunk/SRC/Obsolete/vraidate.pro
r327 r371 21 21 ; 22 22 ; @examples 23 ; IDL> res=vraidate(date) 24 ; IDL> vraidate(980703) 23 ; 24 ; IDL> res=vraidate(date) 25 ; IDL> vraidate(980703) 25 26 ; donne 19980703 qui est un long 26 27 ; -
trunk/SRC/Picture/saveimage.pro
r327 r371 62 62 ; @examples 63 63 ; 64 ; IDL> openr, lun, filepath('hurric.dat', subdir='examples/data'), /get_lun65 ; IDL> image = bytarr(440, 330)66 ; IDL> readu, lun, image67 ; IDL> free_lun, lun68 ; IDL> loadct, 1369 ; IDL> tvscl, image70 ; IDL> saveimage, 'hurric.gif'64 ; IDL> openr, lun, filepath('hurric.dat', subdir='examples/data'), /get_lun 65 ; IDL> image = bytarr(440, 330) 66 ; IDL> readu, lun, image 67 ; IDL> free_lun, lun 68 ; IDL> loadct, 13 69 ; IDL> tvscl, image 70 ; IDL> saveimage, 'hurric.gif' 71 71 ; 72 72 ; @history -
trunk/SRC/Picture/showimage.pro
r325 r371 35 35 ; 36 36 ; @examples 37 ; IDL> showimage, filepath('rose.jpg', subdir='examples/data') 37 ; 38 ; IDL> showimage, filepath('rose.jpg', subdir='examples/data') 38 39 ; 39 40 ; @history -
trunk/SRC/Postscript/closeps.pro
r327 r371 14 14 ; 15 15 ; @uses 16 ; cm_4ps16 ; <pro>cm_4ps</pro> 17 17 ; 18 18 ; @history -
trunk/SRC/Postscript/openps.pro
r327 r371 39 39 ; 40 40 ; @uses 41 ; cm_4ps41 ; <pro>cm_4ps</pro> 42 42 ; 43 43 ; @history -
trunk/SRC/Postscript/printps.pro
r325 r371 197 197 ; 198 198 ; @uses 199 ; cm_4ps199 ; <pro>cm_4ps</pro> 200 200 ; 201 201 ; @examples 202 ; IDL> printps 202 ; 203 ; IDL> printps 203 204 ; 204 205 ; @history -
trunk/SRC/ReadWrite/idl-NetCDF/ncdf_listvars.pro
r308 r371 16 16 ; 17 17 ; @examples 18 ; IDL> varnames = ncdf_listvars('toto.nc') 18 ; 19 ; IDL> varnames = ncdf_listvars('toto.nc') 20 ; 19 21 ; or 20 ; IDL> i=ncdf_open('toto.nc') 21 ; IDL> varnames = ncdf_listvars(id) 22 ; 23 ; IDL> i=ncdf_open('toto.nc') 24 ; IDL> varnames = ncdf_listvars(id) 22 25 ; 23 26 ; @history -
trunk/SRC/ReadWrite/ncdf_getatt.pro
r344 r371 62 62 ; 63 63 ; @examples 64 ; IDL> ncdf_getatt, cdfid, 'sst', add_offset = add_offset, scale_factor = scale_factor, units = units 64 ; 65 ; IDL> ncdf_getatt, cdfid, 'sst', add_offset = add_offset, scale_factor = scale_factor, units = units 65 66 ; 66 67 ; @history -
trunk/SRC/ReadWrite/ncdf_getmask.pro
r361 r371 62 62 ; @examples 63 63 ; 64 ; IDL> mask = ncdf_getmask('HadISST1_1m_187001_200702_sst_reg1m.nc',useasmask = 'sst', missing_value = -1.00000e+30)64 ; IDL> mask = ncdf_getmask('HadISST1_1m_187001_200702_sst_reg1m.nc',useasmask = 'sst', missing_value = -1.00000e+30) 65 65 ; 66 ; IDL> mask = ncdf_getmask('meshmaskORCA2.nc', maskname = 'tmask')66 ; IDL> mask = ncdf_getmask('meshmaskORCA2.nc', maskname = 'tmask') 67 67 ; 68 ; IDL> mask = ncdf_getmask('t106.nc', useasmask = 'SLM', testop = 'le 0.5')68 ; IDL> mask = ncdf_getmask('t106.nc', useasmask = 'SLM', testop = 'le 0.5') 69 69 ; 70 70 ; @history -
trunk/SRC/ReadWrite/read_oasis.pro
r327 r371 28 28 ; 29 29 ; @examples 30 ; IDL> a=read_oasis('grids_orca_t106','a106.lon',320,160)31 ; IDL> m=read_oasis('masks_orca_t106','or1t.msk',182,149,/i4)32 30 ; 33 ; see also IDL> scanoasis,'grids_orca_t106' 31 ; IDL> a=read_oasis('grids_orca_t106','a106.lon',320,160) 32 ; IDL> m=read_oasis('masks_orca_t106','or1t.msk',182,149,/i4) 33 ; 34 ; see also 35 ; 36 ; IDL> scanoasis,'grids_orca_t106' 34 37 ; 35 38 ; @history -
trunk/SRC/ReadWrite/scanoasis.pro
r325 r371 13 13 ; 14 14 ; @examples 15 ; IDL> scanoasis,'grids_orca_t106' 15 ; 16 ; IDL> scanoasis,'grids_orca_t106' 16 17 ; 17 18 ; @history -
trunk/SRC/ReadWrite/write_ncdf.pro
r364 r371 574 574 ; 575 575 ; @uses 576 ; cm_general --> pour la var iodir576 ; <pro>cm_general</pro> for iodir variable 577 577 ; 578 578 ; @examples 579 579 ; 580 580 ; 1) 581 ; write_ncdf, vvsst, btoa581 ; IDL> write_ncdf, vvsst, btoa 582 582 ; 583 583 ; 2) 584 ; write_ncdf, a1, {var:a2,name:'rain'}, ['up','down','fix'], a4, a5, {var:a6,unlim:1,name:'sst',dname:['x','y','z','time']}584 ; IDL> write_ncdf, a1, {var:a2,name:'rain'}, ['up','down','fix'], a4, a5, {var:a6,unlim:1,name:'sst',dname:['x','y','z','time']} 585 585 ; 586 586 ; 3) 587 ; write_ncdf,v031tr,msf031dn,msf031up,v031,vmaskloc,btoa,runame,titplo1,vargrid,nytt $587 ; IDL> write_ncdf,v031tr,msf031dn,msf031up,v031,vmaskloc,btoa,runame,titplo1,vargrid,nytt $ 588 588 ; ,varname = ['v031tr','msf031dn','msf031up','v031','vmaskloc','btoa','runame','titplo1','vargrid','nytt'] $ 589 589 ; ,filename = iodir+'IDL_DATA/'+'waza3.nc' 590 590 ; 591 591 ; 4) 592 ; fileoutnc = iodir+'OBS/TRMM/'+'trmm_1d_'+iyystr+'0101_'+iyystr+'1231_reg0.25.nc'593 ; vv1 = {var:xaxis,name:'nav_lon',dname:'x',at1:{a:'units',b:'degrees_east'},at2:{a:'title',b:'longitude'}}594 ; vv2 = {var:yaxis,name:'nav_lat',dname:'y',at1:{a:'units',b:'degrees_north'},at2:{a:'title',b:'latitude'}}595 ; vv3 = {var:ttt,name:'time',dname:'time_counter',unlim:1,at1:{a:'units',b:timeunit},at2:{a:'title',b:'julian days'} $596 ; vv4 = {var:rain_i2, name:'rain'$597 ; ,at1:{a:'units',b:'mm/day'} 592 ; IDL> fileoutnc = iodir+'OBS/TRMM/'+'trmm_1d_'+iyystr+'0101_'+iyystr+'1231_reg0.25.nc' 593 ; IDL> vv1 = {var:xaxis,name:'nav_lon',dname:'x',at1:{a:'units',b:'degrees_east'},at2:{a:'title',b:'longitude'}} 594 ; IDL> vv2 = {var:yaxis,name:'nav_lat',dname:'y',at1:{a:'units',b:'degrees_north'},at2:{a:'title',b:'latitude'}} 595 ; IDL> vv3 = {var:ttt,name:'time',dname:'time_counter',unlim:1,at1:{a:'units',b:timeunit},at2:{a:'title',b:'julian days'} $ 596 ; IDL> vv4 = {var:rain_i2, name:'rain' $ 597 ; ,at1:{a:'units',b:'mm/day'} $ 598 598 ; ,at2:{a:'title',b:'trmm daily accumulated rainfall derived from the 3-hourly product (mm)'} $ 599 ; ,at3:{a:'scale_factor',b:scala_factor} 600 ; ,at4:{a:'add_offset',b:adda_offset} 601 ; ,at5:{a:'missing_value',b:missaval_i2} 599 ; ,at3:{a:'scale_factor',b:scala_factor} $ 600 ; ,at4:{a:'add_offset',b:adda_offset} $ 601 ; ,at5:{a:'missing_value',b:missaval_i2} $ 602 602 ; } 603 ; glbatt={ gb1:{a:'File_Name', b:'trmm_1d_'+iyystr+'0101_'+iyystr+'1231_reg0.25.nc'}$604 ; ,gb2:{a:'Model_Name', b:'TRMM 3B42_V6 derived product'} 603 ; IDL> glbatt={ gb1:{a:'File_Name', b:'trmm_1d_'+iyystr+'0101_'+iyystr+'1231_reg0.25.nc'} $ 604 ; ,gb2:{a:'Model_Name', b:'TRMM 3B42_V6 derived product'} $ 605 605 ; ,gb3:{a:'Source_File',b:'ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/Derived_Products/3B42_V6/Daily/'+iyystr+'/*.bin'} $ 606 606 ; ,gb4:{a:'IDL_Program_Name', b:'zz08_read_plot_row_trmm_precip.pro (clement@jamstec.go.jp)'} $ 607 ; ,gb5:{a:'Grid', b:'regular 0.25 degres resolution'} 607 ; ,gb5:{a:'Grid', b:'regular 0.25 degres resolution'}$ 608 608 ; } 609 ; write_ncdf, vv1, vv2, vv3, vv4, filename=fileoutnc, globattr=glbatt 610 ; 609 ; IDL> write_ncdf, vv1, vv2, vv3, vv4, filename=fileoutnc, globattr=glbatt 611 610 ; 612 611 ; @side effects … … 615 614 ; CBM 2007-09-10 616 615 ; 617 ; @todo clem 616 ; @todo 617 ; clem 618 618 ; 619 619 ;- -
trunk/SRC/Tests/tst_plt.pro
r368 r371 10 10 ; 11 11 ; To run test with Levitus temperature on a regular 1x1 grid : 12 ; IDL> @tst_initlev13 ; IDL> tst_plt12 ; IDL> @tst_initlev 13 ; IDL> tst_plt 14 14 ; 15 15 ; @history -
trunk/SRC/Tests/tst_pltt.pro
r368 r371 10 10 ; 11 11 ; To run test with Levitus temperature on a regular 1x1 grid : 12 ; IDL> @tst_initlev13 ; IDL> tst_pltt12 ; IDL> @tst_initlev 13 ; IDL> tst_pltt 14 14 ; 15 15 ; @history -
trunk/SRC/Tests/tst_pltz.pro
r368 r371 11 11 ; 12 12 ; To run test with Levitus temperature on a regular 1x1 grid : 13 ; IDL> @tst_initlev14 ; IDL> tst_pltz13 ; IDL> @tst_initlev 14 ; IDL> tst_pltz 15 15 ; 16 16 ; @history -
trunk/SRC/ToBeReviewed/CALCULS/depth2floatlevel.pro
r370 r371 26 26 ; 27 27 ; @examples 28 ; IDL> a=(jpk-1)/(1.*jpi*jpj)*findgen(jpi,jpj) 29 ; IDL> plt, 1e6*(a-floatlevel2depth(depth2floatlevel(a))),/nocontour 28 ; 29 ; IDL> a=(jpk-1)/(1.*jpi*jpj)*findgen(jpi,jpj) 30 ; IDL> plt, 1e6*(a-floatlevel2depth(depth2floatlevel(a))),/nocontour 30 31 ; 31 32 ; ->champ nul a 1e-6 pres -
trunk/SRC/ToBeReviewed/CALCULS/determ2.pro
r325 r371 27 27 ; 28 28 ; @examples 29 ; IDL> a=findgen(2,2,5) 30 ; IDL> print, determ2(a) 31 ; IDL> FOR i=0,4 DO print, determ(a[*,*,i]) 29 ; 30 ; IDL> a=findgen(2,2,5) 31 ; IDL> print, determ2(a) 32 ; IDL> FOR i=0,4 DO print, determ(a[*,*,i]) 32 33 ; 33 34 ; @history -
trunk/SRC/ToBeReviewed/CALCULS/determ3.pro
r231 r371 48 48 ; 49 49 ; @examples 50 ; a=findgen(3,3,5) 51 ; print, determ3(a^2) 52 ; FOR i=0,4 DO print, determ((a[*,*,i])^2) ; IDL solution 50 ; 51 ; IDL> a=findgen(3,3,5) 52 ; IDL> print, determ3(a^2) 53 ; IDL> FOR i=0,4 DO print, determ((a[*,*,i])^2) 53 54 ; 54 55 ; @history -
trunk/SRC/ToBeReviewed/CALCULS/floatlevel2depth.pro
r370 r371 21 21 ; 22 22 ; @examples 23 ; 24 ; 23 ; IDL> a=gdept[jpk-1]/(1.*jpi*jpj)*findgen(jpi,jpj) 24 ; IDL> plt, 1e6*(a-floatlevel2depth(depth2floatlevel(a))),/nocontour 25 25 ; 26 26 ; ->null field at 1e-6 pres -
trunk/SRC/ToBeReviewed/GRILLE/grille.pro
r327 r371 88 88 ; 89 89 ; @uses 90 ; cm_4mesh91 ; cm_4data90 ; <pro>cm_4mesh</pro> 91 ; <pro>cm_4data</pro> 92 92 ; 93 93 ; @restrictions -
trunk/SRC/ToBeReviewed/GRILLE/tracegrille.pro
r370 r371 39 39 ; @examples 40 40 ; 41 ; IDL> plt,indgen(jpi,jpj),/nocontour,/nofill42 ; IDL> vargrid='T'43 ; IDL> tracegrille,/ocean,color=2044 ; IDL> tracegrille,/earth,color=8041 ; IDL> plt,indgen(jpi,jpj),/nocontour,/nofill 42 ; IDL> vargrid='T' 43 ; IDL> tracegrille,/ocean,color=20 44 ; IDL> tracegrille,/earth,color=80 45 45 ; 46 46 ; @history -
trunk/SRC/ToBeReviewed/HOPE/read_hope.pro
r370 r371 286 286 ; @examples 287 287 ; 288 ; 289 ; 288 ; IDL> a=read_hope('xy','ocpt',filename='CLIM_CNT_1993-1998.nc') 289 ; IDL> help, a,/struct 290 290 ; ** Structure <82ec344>, 6 tags, length=1860176, refs=1: 291 291 ; ARRAY FLOAT Array[128, 242, 15] … … 295 295 ; GRID STRING 'T' 296 296 ; HOPEGRID STRUCT -> <Anonymous> Array[1] 297 ; 297 ; IDL> help, a.hopegrid,/struct 298 298 ; ** Structure <82eb9cc>, 8 tags, length=1588, refs=2: 299 299 ; XAXIS FLOAT Array[128] … … 305 305 ; LINETYPE STRING 'odd-even' 306 306 ; PTTYPE STRING 'T' 307 ; 307 ; IDL> help, litchamp(a) 308 308 ; <Expression> FLOAT = Array[128, 242, 15] 309 309 ; -
trunk/SRC/ToBeReviewed/IMAGE/animgif.pro
r370 r371 16 16 ; @restrictions 17 17 ; We use the Z buffer to go faster!. 18 ; 19 ; 20 ; 21 ; 18 ; So, if it failed in animgif, you have to do: 19 ; IDL> device, /close 20 ; IDL> set_plot, 'x' 21 ; IDL> retall 22 22 ; 23 23 ; If we want to delete the \@common, 24 ; 25 ; 24 ; we have to define manually the picture's size 25 ; (variables xsize and ysize) just as iodir. 26 26 ; 27 27 ; @history -
trunk/SRC/ToBeReviewed/INIT/initncdf.pro
r370 r371 43 43 ; 44 44 ; @examples 45 ; IDL> initncdf,'toto.nc',glam=[-180,180] 45 ; 46 ; IDL> initncdf,'toto.nc',glam=[-180,180] 46 47 ; 47 48 ; @history -
trunk/SRC/ToBeReviewed/LECTURE/GRIB/read_gribtable.pro
r325 r371 32 32 ; 33 33 ; @examples 34 ; IDL> tablename='ectab_128' 34 ; 35 ; IDL> tablename='ectab_128' 35 36 ; 36 37 ; @history -
trunk/SRC/ToBeReviewed/LECTURE/binary.pro
r325 r371 11 11 ; 12 12 ; @examples 13 ; 13 14 ; Binary representation of 11b: 14 ; 15 ; IDL> print, binary(11b) 15 16 ; 0 0 0 0 1 0 1 1 17 ; 16 18 ; Binary representation of pi (x86: Little-endian IEEE representation): 17 ; 19 ; IDL> print, format='(z9.8,5x,4(1x,8i1))', long(!pi,0), binary(!pi) 18 20 ; 40490fdb 01000000 01001001 00001111 11011011 (x86 Linux) 19 21 ; 0fdb4149 00001111 11011011 01000001 01001001 (Alpha OpenVMS) 20 ; 22 ; IDL> print, format='(8(1x,8i0))', binary(!dpi) 21 23 ; 01000000 00001001 00100001 11111011 01010100 01000100 00101101 00011000 24 ; 22 25 ; Some first tests before type double was added: 23 ; 26 ; IDL> print, format='(2a6,4x,2z9.8,4x,8z3.2)', $ 24 27 ; !version.arch, !version.os, long(!dpi,0,2), byte(!dpi,0,8) 25 28 ; x86 linux 54442d18 400921fb 18 2d 44 54 fb 21 09 40 -
trunk/SRC/ToBeReviewed/LECTURE/inverse_binary.pro
r325 r371 27 27 ; @examples 28 28 ; 29 ; 30 ; 31 ; 29 ; IDL> a=indgen(5) 30 ; IDL> b=binary(a) 31 ; IDL> help, b 32 32 ; B BYTE = Array[8, 2, 5] 33 ; 33 ; IDL> print, b 34 34 ; 0 0 0 0 0 0 0 0 35 35 ; 0 0 0 0 0 0 0 0 … … 46 46 ; 0 0 0 0 0 0 0 0 47 47 ; 0 0 0 0 0 1 0 0 48 ; 48 ; IDL> help, inverse_binary(b) 49 49 ; <Expression> INT = Array[5] 50 ; 50 ; IDL> print, inverse_binary(b) 51 51 ; 0 1 2 3 4 52 52 ; -
trunk/SRC/ToBeReviewed/LECTURE/litchamp.pro
r370 r371 79 79 ; @examples 80 80 ; 81 ; 81 ; IDL> print, vargrid,', ', varname,', ', varunit,', ', vardate,', ', varexp 82 82 ; T, , , 0, 83 ; 83 ; IDL> help, litchamp({a:indgen(5), u:'C', name:'toto'}) 84 84 ; <Expression> INT = Array[5] 85 ; 85 ; IDL> print, vargrid,', ', varname,', ', varunit,', ', vardate,', ', varexp 86 86 ; T, toto, C, 0, 87 ; 87 ; IDL> help, litchamp({a:indgen(5), da:'1999'}) 88 88 ; <Expression> INT = Array[5] 89 ; 89 ; IDL> print, vargrid,', ', varname,', ', varunit,', ', vardate,', ', varexp 90 90 ; T, toto, C, 1999, 91 91 ; -
trunk/SRC/ToBeReviewed/LECTURE/read_ftp.pro
r325 r371 143 143 ; in the current working directory: 144 144 ; 145 ; 146 ; 147 ; 148 ; 145 ; IDL> files = string(lindgen(10)+50,format='(%"MGB370.%3.3d.gz")') 146 ; IDL> READ_FTP, 'podaac.jpl.nasa.gov', files, $ 147 ; IDL> 'pub/sea_surface_height/topex_poseidon/mgdrb/data/MGB_370', /FILE 148 ; IDL> spawn,'dir MGB*',/log_output 149 149 ; Volume in drive C is Local Disk 150 150 ; Volume Serial Number is 34CE-24DF -
trunk/SRC/ToBeReviewed/LECTURE/xncdf_lec.pro
r370 r371 56 56 ; 57 57 ; @examples 58 ; IDL> help, xncdf_lec() 58 ; 59 ; IDL> help, xncdf_lec() 59 60 ; 60 61 ; @history -
trunk/SRC/ToBeReviewed/MATRICE/colle.pro
r327 r371 56 56 ; 57 57 ; @examples 58 ; IDL> print, colle(replicate(1,2,2,2),indgen(2,2,2),2) 58 ; 59 ; IDL> print, colle(replicate(1,2,2,2),indgen(2,2,2),2) 59 60 ; 1 1 60 61 ; 1 1 -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/plt.pro
r370 r371 106 106 ; map = [0, (lon1+lon2)/2., 0] 107 107 ; Comment: A good way to choose the type of the projection we want to do is to have a look at IDL demo: 108 ; IDL> demo108 ; IDL> demo 109 109 ; Then choose "Mapping" and "Global Mapping Demo" 110 110 ; Comment2: By default it is a cylindrical projection which is done (with or without the keyword map). 111 111 ; If we want an other projection, MAP must be activated and we have to add the keyword: /nom_projection. 112 112 ; For example, for a polar projection centered on the south pole: 113 ; IDL> domdef,-180,180,-90,-45114 ; IDL> plt, tab, /stereo,map=[-90,0,0]113 ; IDL> domdef,-180,180,-90,-45 114 ; IDL> plt, tab, /stereo,map=[-90,0,0] 115 115 ; 116 116 ; @keyword LABMAP {default=labmap=1} -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltv.pro
r327 r371 64 64 ; 65 65 ; @examples 66 ; IDL> tvplus, dist(100)66 ; IDL> tvplus, dist(100) 67 67 ; 68 68 ; @history -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/tvplus.pro
r327 r371 76 76 ; 77 77 ; @examples 78 ; IDL> tvplus, dist(100) 78 ; 79 ; IDL> tvplus, dist(100) 79 80 ; 80 81 ; @history -
trunk/SRC/ToBeReviewed/PLOTS/DIVERS/placedessin.pro
r370 r371 37 37 ; map = [0, (lon1+lon2)/2., 0] 38 38 ; Comment 1: A good way to choose the type of the projection we want to do is to have a look at IDL demo: 39 ; IDL> demo39 ; IDL> demo 40 40 ; Then choose "Mapping" and "Global Mapping Demo" 41 41 ; Comment 2: By default it is a cylindrical projection which is effectuated (with or without the keyword map). 42 42 ; If we want an other projection, MAP must be activated and we have to add the keyword: /nom_projection. 43 43 ; For example, for a polar projection centered on the south pole: 44 ; IDL> domdef,-180,180,-90,-4545 ; IDL> plt, tab, /stereo,map=[-90,0,0]44 ; IDL> domdef,-180,180,-90,-45 45 ; IDL> plt, tab, /stereo,map=[-90,0,0] 46 46 ; 47 47 ; @keyword NOCOLORBAR -
trunk/SRC/ToBeReviewed/PLOTS/LABEL/label_date.pro
r369 r371 53 53 ; For example, to plot from Jan 1, 1993, to July 12, 1994: 54 54 ; 55 ; IDL> Start_date = julday(1, 1, 1993)56 ; IDL> End_date = julday(7, 12, 1994)57 ; IDL> ; Simple mm/dd58 ; IDL> Dummy = LABEL_DATE(DATE_FORMAT='%N/%D')59 ; IDL> ;Time axis60 ; IDL> x = findgen(end_date+1 - start_date) + start_date61 ; IDL> ; Plot with X axis style set to exact62 ; IDL> PLOT, x, sqrt(x), XTICKFORMAT = 'LABEL_DATE', XSTYLE=155 ; IDL> Start_date = julday(1, 1, 1993) 56 ; IDL> End_date = julday(7, 12, 1994) 57 ; IDL> ; Simple mm/dd 58 ; IDL> Dummy = LABEL_DATE(DATE_FORMAT='%N/%D') 59 ; IDL> ;Time axis 60 ; IDL> x = findgen(end_date+1 - start_date) + start_date 61 ; IDL> ; Plot with X axis style set to exact 62 ; IDL> PLOT, x, sqrt(x), XTICKFORMAT = 'LABEL_DATE', XSTYLE=1 63 63 ; 64 64 ; Example with times: … … 66 66 ; For example, to plot from 3PM, Jan 1, 1993, to 5AM, Jan 3, 1993: 67 67 ; 68 ; IDL> ; Also starting offset69 ; IDL> Start_date = Julday(1,1,1993)70 ; IDL> ; Starting_time less offset71 ; IDL> Start_time = (3+12)/24.72 ; IDL> End_time = (Julday(1,3,1993) - Start_date) + 5./24.73 ; IDL> ; Ending date/time - offset, note that the order of operations is74 ; IDL> ; important to avoid loss of precision75 ; IDL> ; MMM NN <new line> HH:MM format76 ; IDL> Dummy = LABEL_DATE(DATE_FORMAT='%D %M!C%H:%I', offset=Start_date)77 ; IDL> ;Time axis78 ; IDL> x = findgen(20) * (End_time - Start_time) / 19 + start_time79 ; IDL> PLOT, x, sqrt(x), XTICKFORMAT = 'LABEL_DATE', XSTYLE=168 ; IDL> ; Also starting offset 69 ; IDL> Start_date = Julday(1,1,1993) 70 ; IDL> ; Starting_time less offset 71 ; IDL> Start_time = (3+12)/24. 72 ; IDL> End_time = (Julday(1,3,1993) - Start_date) + 5./24. 73 ; IDL> ; Ending date/time - offset, note that the order of operations is 74 ; IDL> ; important to avoid loss of precision 75 ; IDL> ; MMM NN <new line> HH:MM format 76 ; IDL> Dummy = LABEL_DATE(DATE_FORMAT='%D %M!C%H:%I', offset=Start_date) 77 ; IDL> ;Time axis 78 ; IDL> x = findgen(20) * (End_time - Start_time) / 19 + start_time 79 ; IDL> PLOT, x, sqrt(x), XTICKFORMAT = 'LABEL_DATE', XSTYLE=1 80 80 ; 81 81 ; @history -
trunk/SRC/ToBeReviewed/STATISTICS/a_timecorrelate.pro
r342 r371 99 99 ; 100 100 ; @examples 101 ; 101 102 ; Define an n-element sample population. 102 ; IDL> x = [3.73, 3.67, 3.77, 3.83, 4.67, 5.87, 6.70, 6.97, 6.40, 5.57]103 ; IDL> x = [3.73, 3.67, 3.77, 3.83, 4.67, 5.87, 6.70, 6.97, 6.40, 5.57] 103 104 ; 104 105 ; Compute the autocorrelation of X for LAG = -3, 0, 1, 3, 4, 8 105 ; IDL> lag = [-3, 0, 1, 3, 4, 8]106 ; IDL> result = a_correlate(x, lag)106 ; IDL> lag = [-3, 0, 1, 3, 4, 8] 107 ; IDL> result = a_correlate(x, lag) 107 108 ; 108 109 ; The result should be: -
trunk/SRC/ToBeReviewed/STATISTICS/c_timecorrelate.pro
r327 r371 89 89 ; 90 90 ; Define two n-elements sample populations. 91 ; IDL> x = [3.73, 3.67, 3.77, 3.83, 4.67, 5.87, 6.70, 6.97, 6.40, 5.57]92 ; IDL> y = [2.31, 2.76, 3.02, 3.13, 3.72, 3.88, 3.97, 4.39, 4.34, 3.95]91 ; IDL> x = [3.73, 3.67, 3.77, 3.83, 4.67, 5.87, 6.70, 6.97, 6.40, 5.57] 92 ; IDL> y = [2.31, 2.76, 3.02, 3.13, 3.72, 3.88, 3.97, 4.39, 4.34, 3.95] 93 93 ; 94 94 ; Compute the cross correlation of X and Y for LAG = -5, 0, 1, 5, 6, 7 95 ; IDL> lag = [-5, 0, 1, 5, 6, 7]96 ; IDL> result = c_timecorrelate(x, y, lag)95 ; IDL> lag = [-5, 0, 1, 5, 6, 7] 96 ; IDL> result = c_timecorrelate(x, y, lag) 97 97 ; 98 98 ; The result should be: -
trunk/SRC/ToBeReviewed/STRING/strkeywd.pro
r325 r371 26 26 ; 27 27 ; @examples 28 ; 28 29 ; We create a structure 29 ; 30 ; 30 ; IDL> b=get_extra(ok=111, year=[1997,1998,1999], age_capitaine=35) 31 ; IDL> help, b,/struct 31 32 ; ** Structure <817f4bc>, 3 tags, length=10, refs=1: 32 33 ; AGE_CAPITAINE INT 35 33 34 ; OK INT 111 34 35 ; YEAR INT Array[3] 36 ; 35 37 ; We put this structure as a string 36 ; 37 ; 38 ; IDL> a=strkeywd(b) 39 ; IDL> print, a 38 40 ; AGE_CAPITAINE=35, OK=111, YEAR=[1997,1998,1999] 41 ; 39 42 ; Now we can use the string a to pass keywords in a function thanks to execute!! 40 ; 41 ; 43 ; IDL> test=execute('c=get_extra('+a+')') 44 ; IDL> help, c,/struct 42 45 ; ** Structure <817f2dc>, 3 tags, length=10, refs=1: 43 46 ; AGE_CAPITAINE INT 35 -
trunk/SRC/ToBeReviewed/STRING/strsed.pro
r327 r371 25 25 ; @examples 26 26 ; 27 ; 28 ; 27 ; IDL> abc = 'abcdefghijklmnopqrstuvwxyz' 28 ; IDL> print, strsed(abc, 'm', 'M') 29 29 ; abcdefghijklMnopqrstuvwxyz 30 ; 30 ; IDL> print, strsed(abc, 'm.*t', 'M_T') 31 31 ; abcdefghijklM_Tuvwxyz 32 ; 33 ; 32 ; IDL> a = 'abcabcabc' 33 ; IDL> help, strsed([abc, a, 'eee'], 'abc', 'XXX_') 34 34 ; <Expression> STRING = Array[3] 35 ; 35 ; IDL> print, strsed([abc, a, 'eee'], 'abc', 'XXX_') 36 36 ; XXX_defghijklmnopqrstuvwxyz XXX_XXX_XXX_ eee 37 37 ; -
trunk/SRC/ToBeReviewed/STRING/strwhere.pro
r325 r371 16 16 ; The number of matches that were found 17 17 ; 18 ; 19 ; 18 ; The function returns an index array similar to the 19 ; result of the where function 20 20 ; 21 21 ; @examples 22 ; IDL> ind = strwhere('abcabcabc','a') 23 ; IDL> print, ind 22 ; 23 ; find position of string "a" in the string "abcabcabc" 24 ; 25 ; IDL> ind = strwhere('abcabcabc','a') 26 ; IDL> print, ind 24 27 ; [ 0, 3, 6 ] 25 28 ; -
trunk/SRC/ToBeReviewed/STRING/tostr.pro
r325 r371 26 26 ; 27 27 ; @examples 28 ; 28 ; IDL> help, tostr(1),tostr('a'),tostr(indgen(4)),tostr(['a','jkfjo']) 29 29 ; <Expression> STRING = '1' 30 30 ; <Expression> STRING = ''a'' 31 31 ; <Expression> STRING = '[0,1,2,3]' 32 32 ; <Expression> STRING = '['a','jkfjo']' 33 ; 33 ; IDL> print, tostr(['c''est bon','c''est bon']) 34 34 ; ['c''est bon','c''est bon'] 35 35 ; -
trunk/SRC/ToBeReviewed/STRUCTURE/chkstru.pro
r327 r371 37 37 ; 38 38 ; @examples 39 ; test = { a:1, b:2, c:3 } 40 ; required = ['a','c'] 41 ; if CHKSTRU(test,required) then print,'found a and c.' 42 ; IDL> print, CHKSTRU(test,'b') 39 ; 40 ; IDL> test = { a:1, b:2, c:3 } 41 ; IDL> required = ['a','c'] 42 ; IDL> if CHKSTRU(test,required) then print,'found a and c.' 43 ; IDL> print, CHKSTRU(test,'b') 43 44 ; 1 44 ; 45 ; IDL> print, CHKSTRU(test,'b',/extract) 45 46 ; 2 46 47 ; -
trunk/SRC/ToBeReviewed/STRUCTURE/extractstru.pro
r327 r371 25 25 ; @examples 26 26 ; 27 ; 28 ; 27 ; IDL> extra=get_extra(/ok, year=1999, age_capitaine=35 ) 28 ; IDL> help, extra,/struct 29 29 ; ** Structure <83e66bc>, 3 tags, length=6, refs=1: 30 30 ; AGE_CAPITAINE INT 35 31 31 ; OK INT 1 32 32 ; YEAR INT 1999 33 ; 33 ; IDL> help, extractstru(extra,['ok','hhuihi','YEAR']),/stru 34 34 ; ** Structure <831afac>, 1 tags, length=2, refs=1: 35 35 ; AGE_CAPITAINE INT 35 36 ; 36 ; IDL> help, extractstru(extra,['ok','hhuihi','YEAR'],/keep),/stru 37 37 ; ** Structure <834bbc4>, 2 tags, length=4, refs=1: 38 38 ; OK INT 1 -
trunk/SRC/ToBeReviewed/STRUCTURE/mixstru.pro
r325 r371 26 26 ; @examples 27 27 ; 28 ; 29 ; 30 ; 28 ; IDL> a=get_extra(/toto,ok=123) 29 ; IDL> b=get_extra(ok=111, year=1999, age_capitaine=35) 30 ; IDL> help, a,b,/struct 31 31 ; ** Structure <8334424>, 2 tags, length=4, refs=1: 32 32 ; OK INT 123 … … 36 36 ; OK INT 111 37 37 ; YEAR INT 1999 38 ; 38 ; IDL> help, mixstru(a,b),/struct 39 39 ; ** Structure <82f25ac>, 4 tags, length=8, refs=1: 40 40 ; AGE_CAPITAINE INT 35 … … 42 42 ; OK INT 123 43 43 ; TOTO INT 1 44 ; 44 ; IDL> help, mixstru(b,a),/struct 45 45 ; ** Structure <834604c>, 4 tags, length=8, refs=1: 46 46 ; TOTO INT 1 -
trunk/SRC/ToBeReviewed/STRUCTURE/struct2string.pro
r327 r371 25 25 ; @examples 26 26 ; 27 ; 27 ; IDL> print, struct2string(!d) 28 28 ; create_struct('NAME','X','X_SIZE',891,'Y_SIZE',630,'X_VSIZE' 29 29 ; ,891,'Y_VSIZE',630,'X_CH_SIZE',6,'Y_CH_SIZE',10,'X_PX_CM' -
trunk/SRC/ToBeReviewed/STRUCTURE/where_tag.pro
r262 r371 41 41 ; 42 42 ; @examples 43 ; Suppose STR is a structure with tags CAT_NO:indgen(10), and44 ; NAME:strarr(10). Find the indices where STR.CAT_NO is45 ; between 3 and 5.46 43 ; 47 ; IDL> print, WHERE_TAG( str, TAG_NAME = 'CAT_NO', VALUE = [3,4,5] ) ;or 48 ; IDL> print, WHERE_TAG( str, TAG_NUM = 0, RANGE = [3,5]) 44 ; Suppose STR is a structure with tags CAT_NO:indgen(10), and 45 ; NAME:strarr(10). Find the indices where STR.CAT_NO is between 3 and 5. 46 ; 47 ; IDL> print, WHERE_TAG( str, TAG_NAME = 'CAT_NO', VALUE = [3,4,5] ) 48 ; or 49 ; IDL> print, WHERE_TAG( str, TAG_NUM = 0, RANGE = [3,5]) 49 50 ; 50 51 ; @history -
trunk/SRC/ToBeReviewed/TRIANGULATION/definetri.pro
r325 r371 62 62 ; @examples 63 63 ; 64 ; IDL> triangles=definetri(3,3,[1,3])64 ; IDL> triangles=definetri(3,3,[1,3]) 65 65 ; 66 66 ; triangles will be this kind of triangulation: -
trunk/SRC/ToBeReviewed/WIDGET/COMPOUND_WIDGET/cw_calendar.pro
r327 r371 1 ; IDL> testwid,julday(1,1,1980)+lindgen(100)*51 ; IDL> testwid,julday(1,1,1980)+lindgen(100)*5 2 2 ; 3 3 ; PRO testwid_event, event -
trunk/SRC/ToBeReviewed/WIDGET/COMPOUND_WIDGET/cw_domain.pro
r327 r371 1 ; IDL> testwid1 ; IDL> testwid 2 2 ; PRO testwid_event, event 3 3 ; ; help, event, /struct -
trunk/SRC/ToBeReviewed/WIDGET/xnotice.pro
r327 r371 24 24 ; 25 25 ; @examples 26 ; IDL> id=xnotice('ca marche !C ou pas?') 27 ; IDL> widget_control, id, /destroy 26 ; 27 ; IDL> id=xnotice('ca marche !C ou pas?') 28 ; IDL> widget_control, id, /destroy 28 29 ; 29 30 ; @history -
trunk/SRC/ToBeReviewed/WIDGET/xquestion.pro
r327 r371 72 72 ; 73 73 ; @examples 74 ; IDL> help, xquestion('Postscript name') 74 ; 75 ; IDL> help, xquestion('Postscript name') 75 76 ; <Expression> STRING = 'toto.ps' 76 77 ; -
trunk/SRC/Utilities/createfunc.pro
r327 r371 26 26 ; 27 27 ; @examples 28 ; IDL> print, createfunc('3*2', FILENAMEIN='test') 29 ; IDL> print, createfunc('3*two', FILENAMEIN='test' $ 30 ; IDL> , KWDLIST='two = two', TWO=2) 28 ; 29 ; IDL> print, createfunc('3*2', FILENAMEIN='test') 30 ; IDL> print, createfunc('3*two', FILENAMEIN='test' $ 31 ; IDL> , KWDLIST='two = two', TWO=2) 31 32 ; 32 33 ; @history -
trunk/SRC/Utilities/createpro.pro
r327 r371 29 29 ; 30 30 ; @examples 31 ; IDL> createpro, ['print,''OK'''], FILENAMEIN='test' 32 ; IDL> createpro, ['if keyword_set(ok) then print,''OK'' else print, ''No'''] $ 33 ; IDL> , FILENAMEIN='test', KWDLIST =', ok = ok' 34 ; IDL> createpro, ['if keyword_set(ok) then print,''OK'' else print, ''No'''] $ 35 ; IDL> , FILENAMEIN='test', KWDLIST= ', ok = ok', /ok 31 ; 32 ; IDL> createpro, ['print,''OK'''], FILENAMEIN='test' 33 ; IDL> createpro, ['if keyword_set(ok) then print,''OK'' else print, ''No'''] $ 34 ; IDL> , FILENAMEIN='test', KWDLIST =', ok = ok' 35 ; IDL> createpro, ['if keyword_set(ok) then print,''OK'' else print, ''No'''] $ 36 ; IDL> , FILENAMEIN='test', KWDLIST= ', ok = ok', /ok 36 37 ; 37 38 ; @history -
trunk/SRC/Utilities/def_myuniquetmpdir.pro
r325 r371 9 9 ; 10 10 ; @examples 11 ; IDL> def_myuniquetmpdir 11 ; 12 ; IDL> def_myuniquetmpdir 12 13 ; 13 14 ; @uses 14 ; cm_general15 ; <pro>cm_general</pro> 15 16 ; 16 17 ; @history -
trunk/SRC/Utilities/demomode_compatibility.pro
r325 r371 5 5 ; 6 6 ; @uses 7 ; cm_general7 ; <pro>cm_general</pro> 8 8 ; 9 9 ; @version -
trunk/SRC/Utilities/find.pro
r327 r371 72 72 ; 73 73 ; @examples 74 ; IDL> print, find('*loadct') 74 ; 75 ; IDL> print, find('*loadct') 75 76 ; /usr/local/rsi/idl_6.0/lib/utilities/xloadct.pro 76 77 ; /usr/local/rsi/idl_6.0/lib/loadct.pro 77 ; IDL> print, find('*loadct', iodir=!dir,/recursive)78 ; IDL> print, find('*loadct', iodir=!dir,/recursive) 78 79 ; /usr/local/rsi/idl_6.0/lib/loadct.pro 79 80 ; /usr/local/rsi/idl_6.0/lib/utilities/xloadct.pro 80 ; IDL> print, find('*loadct.pro')81 ; IDL> print, find('*loadct.pro') 81 82 ; /usr/local/rsi/idl_6.0/lib/utilities/xloadct.pro 82 83 ; /usr/local/rsi/idl_6.0/lib/loadct.pro 83 ; IDL> print, find('*loadct',/nopro)84 ; IDL> print, find('*loadct',/nopro) 84 85 ; NOT FOUND 85 ; IDL> print, find('*loadct', iodir = '/usr/local/rsi/idl_6.0/lib')86 ; IDL> print, find('*loadct', iodir = '/usr/local/rsi/idl_6.0/lib') 86 87 ; /usr/local/rsi/idl_6.0/lib/loadct.pro 87 ; IDL> print, find('*loadct', iodir = '/usr/local/rsi/idl_6.0/lib', /test_write)88 ; IDL> print, find('*loadct', iodir = '/usr/local/rsi/idl_6.0/lib', /test_write) 88 89 ; NOT FOUND 89 ; IDL> print, find('*loadct', iodir = '/usr/local/rsi/idl_6.0/lib', /recursive)90 ; IDL> print, find('*loadct', iodir = '/usr/local/rsi/idl_6.0/lib', /recursive) 90 91 ; /usr/local/rsi/idl_6.0/lib/loadct.pro 91 92 ; /usr/local/rsi/idl_6.0/lib/utilities/xloadct.pro 92 ; IDL> print, find('mesh*', iodirectory = [iodir, !path])93 ; IDL> print, find('mesh*', iodirectory = [iodir, !path]) 93 94 ; /Users/sebastie/DATA/ORCA2/meshmaskORCA2closea.nc 94 95 ; /Users/sebastie/IDL/meshmaskclosesea.pro -
trunk/SRC/Utilities/fitintobox.pro
r327 r371 93 93 ; 94 94 ; @returns 95 ; an array with dimensions matching the domain 96 ; or -1 if there is an error... 95 ; an array with dimensions matching the domain or -1 if there is an error. 97 96 ; 98 97 ; @uses 99 ; cm_4mesh100 ; cm_4cal98 ; <pro>cm_4mesh</pro> 99 ; <pro>cm_4cal</pro> 101 100 ; 102 101 ; @examples 103 ; IDL> help, fitintobox(findgen(jpi,jpj)) 102 ; 103 ; IDL> help, fitintobox(findgen(jpi,jpj)) 104 104 ; <Expression> FLOAT = Array[41, 3] 105 ; IDL> help, fitintobox(findgen(jpi,jpj,78))105 ; IDL> help, fitintobox(findgen(jpi,jpj,78)) 106 106 ; Error: 107 107 ; the array dimensions [180,148,78] are incompatible -
trunk/SRC/Utilities/isadirectory.pro
r327 r371 27 27 ; @examples 28 28 ; 29 ; IDL> print, !dir29 ; IDL> print, !dir 30 30 ; /usr/local/rsi/idl_6.0 31 ; IDL> print, isadirectory(!dir)31 ; IDL> print, isadirectory(!dir) 32 32 ; /usr/local/rsi/idl_6.0/ 33 ; IDL> print, isadirectory(!dir+'notgood')33 ; IDL> print, isadirectory(!dir+'notgood') 34 34 ; 35 35 ; @history -
trunk/SRC/Utilities/isafile.pro
r327 r371 45 45 ; 46 46 ; @examples 47 ; IDL> print, isafile('/Users/sebastie/SAXO_RD/Commons/cm_4mesh.pro') 47 ; 48 ; IDL> print, isafile('/Users/sebastie/SAXO_RD/Commons/cm_4mesh.pro') 48 49 ; /Users/sebastie/SAXO_RD/Commons/cm_4mesh.pro 49 ; IDL> print, isafile('cm_4mesh.pro', iodir = '/Users/sebastie/SAXO_RD/Commons')50 ; IDL> print, isafile('cm_4mesh.pro', iodir = '/Users/sebastie/SAXO_RD/Commons') 50 51 ; /Users/sebastie/SAXO_RD/Commons/cm_4mesh.pro 51 ; IDL> print, isafile('cm_4mesh.pro', iodir = !path)52 ; IDL> print, isafile('cm_4mesh.pro', iodir = !path) 52 53 ; /Users/sebastie/SAXO_RD/Commons/cm_4mesh.pro 53 ; IDL> print, isafile('cm_4mesh.pro', iodir = '/Users/sebastie/SAXO_RD', /recursive)54 ; IDL> print, isafile('cm_4mesh.pro', iodir = '/Users/sebastie/SAXO_RD', /recursive) 54 55 ; /Users/sebastie/SAXO_RD/Commons/cm_4mesh.pro 55 ; IDL> print, isafile('cm_4mesh.pro', iodir = getenv('HOME'), /recursive)56 ; IDL> print, isafile('cm_4mesh.pro', iodir = getenv('HOME'), /recursive) 56 57 ; /Users/sebastie/SAXO_RD/Commons/cm_4mesh.pro 57 ; IDL> print, isafile('fake_file.pro')58 ; IDL> print, isafile('fake_file.pro') 58 59 ; 59 60 ; @history -
trunk/SRC/Utilities/linearequation.pro
r325 r371 28 28 ; 29 29 ; @examples 30 ; IDL> abc=linearequation(complex(1,2),[3,4]) 31 ; IDL> print, abc[0]*1+abc[1]*2+abc[2] 30 ; 31 ; IDL> abc=linearequation(complex(1,2),[3,4]) 32 ; IDL> print, abc[0]*1+abc[1]*2+abc[2] 32 33 ; 0.00000 33 34 ; -
trunk/SRC/Utilities/lineintersection.pro
r327 r371 39 39 ; 40 40 ; @examples 41 ; IDL> abc1=linearequation(complex(1,2),[3,4]) 42 ; IDL> abc2=linearequation(complex(1,2),[8,15]) 43 ; IDL> print, lineintersection(abc1, abc2) 41 ; 42 ; IDL> abc1=linearequation(complex(1,2),[3,4]) 43 ; IDL> abc2=linearequation(complex(1,2),[8,15]) 44 ; IDL> print, lineintersection(abc1, abc2) 44 45 ; ( 1.00000, 2.00000) 45 ; IDL> print, lineintersection(abc1, abc2,/float)46 ; IDL> print, lineintersection(abc1, abc2,/float) 46 47 ; 1.00000 2.00000 47 48 ; -
trunk/SRC/Utilities/protype.pro
r325 r371 18 18 ; @examples 19 19 ; 20 ; IDL> print, protype('protype')20 ; IDL> print, protype('protype') 21 21 ; func 22 ; IDL> print, protype('protype.pro')22 ; IDL> print, protype('protype.pro') 23 23 ; func 24 ; IDL> print, protype('init')24 ; IDL> print, protype('init') 25 25 ; batch 26 ; IDL> print, protype('plt')26 ; IDL> print, protype('plt') 27 27 ; proc 28 28 ; -
trunk/SRC/Utilities/report.pro
r327 r371 39 39 ; If there is not any widget activated: 40 40 ; 41 ; IDL> help, report('toto tata')41 ; IDL> help, report('toto tata') 42 42 ; % $MAIN$: toto tata 43 43 ; <Expression> INT = -1 44 ; IDL> help, report('does it works ?',/question)44 ; IDL> help, report('does it works ?',/question) 45 45 ; does it works ? y/n (default answer is y) 46 46 ; <Expression> BYTE = 1 47 ; IDL> help, report('question1: !C does it works ?',/question)47 ; IDL> help, report('question1: !C does it works ?',/question) 48 48 ; question1: 49 49 ; does it works ? y/n (default answer is y) -
trunk/SRC/Utilities/routine_name.pro
r325 r371 22 22 ; 23 23 ; @examples 24 ; IDL> print, routine_name() 24 ; 25 ; IDL> print, routine_name() 25 26 ; /usr1/com/smasson/IDL_RD/UTILITAIRE/report.pro 26 ; IDL> print, routine_name(1)27 ; IDL> print, routine_name(1) 27 28 ; /usr1/com/smasson/IDL_RD/PLOTS/DIVERS/determineminmax.pro 28 ; IDL> print, routine_name(2)29 ; IDL> print, routine_name(2) 29 30 ; /usr1/com/smasson/IDL_RD/PLOTS/DESSINE/plt.pro 30 ; IDL> print, routine_name(3)31 ; IDL> print, routine_name(3) 31 32 ; $MAIN$ 32 ; IDL> print, routine_name(4)33 ; IDL> print, routine_name(4) 33 34 ; $MAIN$ 34 35 ; -
trunk/SRC/Utilities/testvar.pro
r327 r371 15 15 ; 16 16 ; @examples 17 ; IDL> print, testvar(var=toto)17 ; IDL> print, testvar(var=toto) 18 18 ; 0 19 ; IDL> print, testvar(var='toto')19 ; IDL> print, testvar(var='toto') 20 20 ; toto 21 21 ; -
trunk/SRC/Utilities/undefine.pro
r325 r371 13 13 ; 14 14 ; @examples 15 ; IDL> a=116 ; IDL> undefine,a15 ; IDL> a=1 16 ; IDL> undefine,a 17 17 ; % Compiled module: UNDEFINE. 18 ; IDL> help, a18 ; IDL> help, a 19 19 ; A UNDEFINED = <Undefined> 20 20 ; -
trunk/SRC/Utilities/xfile.pro
r327 r371 23 23 ; 24 24 ; @examples 25 ; IDL> xfile,'plt' 25 ; 26 ; IDL> xfile,'plt' 26 27 ; 27 28 ; @history -
trunk/SRC/Utilities/xhelp.pro
r327 r371 20 20 ; 21 21 ; @examples 22 ; IDL> xhelp,'plt' 22 ; 23 ; IDL> xhelp,'plt' 23 24 ; 24 25 ; @history
Note: See TracChangeset
for help on using the changeset viewer.