- Timestamp:
- 09/25/07 17:09:12 (17 years ago)
- Location:
- trunk/SRC
- Files:
-
- 38 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Calendar/julday.pro
r292 r297 11 11 ; Calendar 12 12 ; 13 ; @param MONTH {in}{ required} {type=scalar (integer or double) or array of scalars}13 ; @param MONTH {in}{optional} {type=scalar (integer or double) or array of scalars} 14 14 ; Number of the desired month (1 = January, ..., 12 = December). 15 15 ; 16 ; @param DAY {in}{ required} {type=scalar (integer or double) or array of scalars}16 ; @param DAY {in}{optional} {type=scalar (integer or double) or array of scalars} 17 17 ; Number of day of the month. 18 18 ; 19 ; @param YEARin {in}{ required} {type=scalar (integer or double) or array of scalars}19 ; @param YEARin {in}{optional} {type=scalar (integer or double) or array of scalars} 20 20 ; Number of the desired year.Year parameters must be valid 21 21 ; values from the civil calendar. Years B.C.E. are represented … … 93 93 IF n_elements(key_caltype) EQ 0 THEN key_caltype = 'greg' 94 94 if keyword_set(ndayspm) then key_caltype = '360d' 95 96 NP = n_params() 97 ; Process the input, if all are missing, use today's date. 98 IF (np EQ 0) THEN RETURN, SYSTIME(/JULIAN) 99 IF (np LT 3) THEN ras = report('Incorrect number of arguments.') 95 100 ; 96 101 YEAR = long(yearin) 97 102 zero = where(year EQ 0, cnt) 98 103 IF cnt NE 0 THEN YEAR[zero] = 654321L 99 100 NP = n_params()101 ; Process the input, if all are missing, use todays date.102 IF (np EQ 0) THEN RETURN, SYSTIME(/JULIAN)103 IF (np LT 3) THEN ras = report('Incorrect number of arguments.')104 104 ; 105 105 CASE key_caltype OF -
trunk/SRC/Grid/computegrid.pro
r287 r297 509 509 ; check y periodicity... Only according to ORCA grid 510 510 ;==================================================== 511 ; check the peri dicity if iyminmesh and iymaxmesh have the default definitions...511 ; check the periodicity if iyminmesh and iymaxmesh have the default definitions... 512 512 IF NOT keyword_set(plain) AND key_onearth EQ 1 $ 513 513 AND iyminmesh EQ 0l AND iymaxmesh eq jpjglo-1 AND jpj GE 3 AND jpi GE 2 THEN BEGIN … … 546 546 ;==================================================== 547 547 IF n_elements(periodic) NE 0 THEN forcenoperio = 1 - keyword_set(periodic) 548 ; check the peri dicity if ixminmesh and ixmaxmesh have the default definitions...548 ; check the periodicity if ixminmesh and ixmaxmesh have the default definitions... 549 549 IF NOT keyword_set(plain) AND NOT keyword_set(forcenoperio) AND key_onearth EQ 1 $ 550 550 AND ixminmesh EQ 0l AND ixmaxmesh eq jpiglo-1 AND jpi GE 3 THEN BEGIN -
trunk/SRC/Grid/micromeshmask.pro
r240 r297 101 101 ENDELSE 102 102 ENDELSE 103 ; ;get the variables list related to the partial steps103 ; get the variables list related to the partial steps 104 104 varlist_ps = ncdf_listvars(cdfid) 105 105 varlist_ps = strtrim(strlowcase(varlist_ps), 2) -
trunk/SRC/Interpolation/compute_fromirr_bilinear_weigaddr.pro
r295 r297 114 114 IF out[0] NE -1 THEN olon[out] = olon[out]+360 115 115 ; 116 ; we work only on the ou put grid water points116 ; we work only on the output grid water points 117 117 awater = where(amsk EQ 1) 118 118 nawater = n_elements(awater) -
trunk/SRC/Interpolation/inquad.pro
r242 r297 29 29 ; to specify that the quadrilateral are on a sphere and 30 30 ; that their coordinates are longitude-latitude coordinates. In this 31 ; case, e st-west periodicity, poles singularity and other pbs31 ; case, east-west periodicity, poles singularity and other pbs 32 32 ; related to longitude-latitude coordinates are managed 33 33 ; automatically. … … 345 345 !p = save.p 346 346 ENDIF 347 ; ;347 ; 348 348 RETURN, found 349 349 END -
trunk/SRC/Interpolation/lbcorca.pro
r262 r297 53 53 res = 0 54 54 55 ; e st-west periodicity checks55 ; east-west periodicity checks 56 56 57 57 IF array_equal(arr[0, *, *, *], arr[jpi-2, *, *, *]) NE 1 THEN BEGIN … … 156 156 ENDCASE 157 157 158 ; e st-west periodicity correction (again)158 ; east-west periodicity correction (again) 159 159 160 160 IF keyword_set(correction) THEN BEGIN -
trunk/SRC/Interpolation/spl_incr.pro
r242 r297 194 194 ; more than 2 pieces... 195 195 IF cntbad GT 1 THEN BEGIN 196 ; we take care of the piece located wetween bad[ib-1]+1 and bad[ib]196 ; we take care of the piece located between bad[ib-1]+1 and bad[ib] 197 197 FOR ib = 1, cntbad-1 DO BEGIN 198 198 ; if there is x2 values smaller that x[bad[ib]], then the x2 values … … 287 287 ; we force yp0 = 0.0d and ypn_1 = 0.0d 288 288 IF cntbad GT 1 THEN BEGIN 289 ; we take care of the piece located wetween bad[ib-1] and bad[ib]289 ; we take care of the piece located between bad[ib-1] and bad[ib] 290 290 FOR ib = 1, cntbad-1 DO BEGIN 291 291 ; if there is x2 values smaller that x[bad[ib]], then the x2 values -
trunk/SRC/Interpolation/spl_keep_mean.pro
r242 r297 7 7 ; of X2. The interpolation method is based on cubic spline, corrected 8 8 ; in a way that integral of the interpolated values is the same as the 9 ; integral of the input values. (-> for ex emple to build daily data9 ; integral of the input values. (-> for example to build daily data 10 10 ; from monthly mean and keep the monthly mean of the computed daily 11 ; data equa to the original values)11 ; data equal to the original values) 12 12 ; 13 13 ; @param x {in}{required} -
trunk/SRC/Obsolete/lec.pro
r254 r297 283 283 , iyminmesh-iymindta:iymaxmesh-iymindta, izminmesh-izmindta:izmaxmesh-izmindta] 284 284 ;--------------------------------------------------------------------- 285 ; on shift z si key_shift est defini nit285 ; on shift z si key_shift est defini 286 286 ;--------------------------------------------------------------------- 287 287 if n_elements(key_shift) NE 0 THEN BEGIN … … 290 290 endif 291 291 ;--------------------------------------------------------------------- 292 ; si /TOUT n''est pas active, on coupe z pour q ''il soit a la taille292 ; si /TOUT n''est pas active, on coupe z pour qu''il soit a la taille 293 293 ; du zoom: nx,ny nz 294 294 ;--------------------------------------------------------------------- -
trunk/SRC/Obsolete/lect.pro
r231 r297 98 98 ficname=ficname+'.'+box+'.'+direc+'.hovdat' 99 99 ;------------------------------------------------------------- 100 ; Est ce que le fichier de hovmo ller existe ?100 ; Est ce que le fichier de hovmoeller existe ? 101 101 ;------------------------------------------------------------- 102 102 ; -
trunk/SRC/Obsolete/nlec5j.pro
r254 r297 151 151 ;------------------------------------------------------------ 152 152 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 153 ; + on recupere la dim du no 9uveau domaine153 ; + on recupere la dim du nouveau domaine 154 154 ;------------------------------------------------------------ 155 155 if keyword_set(boite) then BEGIN -
trunk/SRC/Obsolete/nlecan.pro
r254 r297 142 142 ;------------------------------------------------------------ 143 143 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 144 ; + on recupere la dim du no 9uveau domaine144 ; + on recupere la dim du nouveau domaine 145 145 ;------------------------------------------------------------ 146 146 if keyword_set(boite) then BEGIN -
trunk/SRC/Obsolete/nlecmois.pro
r254 r297 99 99 ;------------------------------------------------------------ 100 100 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 101 ; + on recupere la dim du no 9uveau domaine101 ; + on recupere la dim du nouveau domaine 102 102 ;------------------------------------------------------------ 103 103 if keyword_set(boite) then BEGIN -
trunk/SRC/Obsolete/nlecsaison.pro
r254 r297 141 141 ;------------------------------------------------------------ 142 142 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 143 ; + on recupere la dim du no 9uveau domaine143 ; + on recupere la dim du nouveau domaine 144 144 ;------------------------------------------------------------ 145 145 if keyword_set(boite) then BEGIN -
trunk/SRC/Obsolete/nlecserie.pro
r254 r297 137 137 ;------------------------------------------------------------ 138 138 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 139 ; + on recupere la dim du no 9uveau domaine139 ; + on recupere la dim du nouveau domaine 140 140 ;------------------------------------------------------------ 141 141 if keyword_set(boite) then BEGIN -
trunk/SRC/ReadWrite/ncdf_gettime.pro
r272 r297 34 34 ; a double 1D array of IDL Julian days 35 35 ; In case of error return -1 if the time dimension was not found 36 ; or return -jpt if it as been found that the time dime sion size is jpt36 ; or return -jpt if it as been found that the time dimension size is jpt 37 37 ; 38 38 ; @restrictions -
trunk/SRC/Textoidl/textable.pro
r231 r297 166 166 ] 167 167 ;Special symbols -- 168 ; NOTES -- You must leave \infty before \in in the translat ation168 ; NOTES -- You must leave \infty before \in in the translation 169 169 ; table to avoid having the \in part of \infty translated 170 170 ; away. -
trunk/SRC/ToBeReviewed/GRILLE/cmpgrid.pro
r232 r297 40 40 if n_elements(ccmeshparameters) EQ 0 then return, 1 41 41 ; 42 ; we compare the structure which c aracterise the grid whith42 ; we compare the structure which characterise the grid with 43 43 ; ccmeshparameters 44 44 ; -
trunk/SRC/ToBeReviewed/GRILLE/decoupeterre.pro
r268 r297 108 108 END 109 109 ENDCASE 110 ; however for the vertical section we don''t w hant to extent the domain110 ; however for the vertical section we don''t want to extent the domain 111 111 ; in the direction perpendicular to the vertical section 112 112 if type EQ 'xz' then begin -
trunk/SRC/ToBeReviewed/INIT/initncdf.pro
r271 r297 109 109 zvarid = (where(namevar EQ 'nav_lev' or namevar EQ zaxisname OR namevar EQ 'level' OR namevar EQ 'lev' OR strmid(namevar, 0, 5) EQ 'depth'))[0] 110 110 if zvarid EQ -1 AND inside.ndims GT 3 then begin 111 ras = report( 'initncdf: the zaxis was not found..., check the use of ZAXISNAME keyword if you w hant to find one...')111 ras = report( 'initncdf: the zaxis was not found..., check the use of ZAXISNAME keyword if you want to find one...') 112 112 ; stop 113 113 endif -
trunk/SRC/ToBeReviewed/LECTURE/changeread.pro
r232 r297 34 34 if size(newread, /type) NE 8 then return, 0 35 35 ; 36 ; we compare the two structure which c aracterise the read36 ; we compare the two structure which characterise the read 37 37 ; 38 38 case 1 of -
trunk/SRC/ToBeReviewed/LECTURE/xncdf_lec.pro
r295 r297 75 75 COMMON motcle, mcatt, mccount, mcoffset, mciodir, mcshift, mcstride, mcvar 76 76 ;------------------------------------------------------------ 77 ; Trick for using keywords (we pass by variables declar ated in a common)77 ; Trick for using keywords (we pass by variables declared in a common) 78 78 ;------------------------------------------------------------ 79 79 res = -1 … … 207 207 ; 208 208 ; @param EVENT {in}{required} 209 ; A structure c aracterizing the type of event which arrive to a widget number1 2209 ; A structure characterizing the type of event which arrive to a widget number1 2 210 210 ; 211 211 ; @uses … … 378 378 ; 379 379 ; @file_comments 380 ; This procedure manage the second created w higet when we call xncdf_lec.380 ; This procedure manage the second created widget when we call xncdf_lec. 381 381 ; This widget concern the reading of the variable. 382 382 ; … … 385 385 ; created by xncdf_lec and which has allowed to select the variable to be read. 386 386 ; 387 ; OUTPUTS: indirectement res (le tableau ou la structure resultat) 387 ; @results 388 ; indirectement res (le tableau ou la structure resultat) 388 389 ; 389 390 ; @uses … … 449 450 widbase2 = widget_base(widbase, /column) 450 451 ; To each attribute, we created a widget (widbase21) containing in line a button 451 ; yes/no (widbase211), and two wi gdet text (widbase212, widbase213) comprising the452 ; yes/no (widbase211), and two widget text (widbase212, widbase213) comprising the 452 453 ; name and the value of the attribute. 453 454 widbase21 = lonarr(varcontient.natts) … … 459 460 widbase212 = lonarr(varcontient.natts) 460 461 widbase213 = lonarr(varcontient.natts) 461 for attid = 0, varcontient.natts-1 do BEGIN ;Lo p on the number of attribute.462 for attid = 0, varcontient.natts-1 do BEGIN ;Loop on the number of attribute. 462 463 widbase21[attid] = widget_base(widbase2, /row) 463 464 name=ncdf_attname(cdfid,varid,attid) … … 491 492 ; 492 493 ; @param EVENT 493 ; A structure c aracterising the event type which arrive at the widget number 1.494 ; A structure characterising the event type which arrive at the widget number 1. 494 495 ; 495 496 ; @uses -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/plt1d.pro
r295 r297 8 8 ; 9 9 ; @param TAB {in}{required} 10 ; The field whose we want to make the hovmo ller map can be 2 kind of thing:10 ; The field whose we want to make the hovmoeller map can be 2 kind of thing: 11 11 ; 1) An array which can be: 12 12 ; * 2d, 3d or 4d: array xy, xyz. xyt or xyzt. In this case, the array will pass … … 84 84 ; 85 85 ; @keyword TYPEIN 86 ; allows to specify the type of hovmo ller we want to do86 ; allows to specify the type of hovmoeller we want to do 87 87 ; 'xt','yt','zt','t' 88 88 ; with help of a keyword rather than the argument type. If the argument and the -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltsc.pro
r295 r297 68 68 ; Comment: We do not reinitialize when we call back pltsc 69 69 if NOT keyword_set(ov1d) then reinitplt 70 ; ;reduce data xyzt domain70 ; reduce data xyzt domain 71 71 72 72 if keyword_set(boxzoom) then BEGIN -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltt.pro
r295 r297 2 2 ; 3 3 ; @file_comments 4 ; Trace hovmo ller graphs: xt,yt,zt,t4 ; Trace hovmoeller graphs: xt,yt,zt,t 5 5 ; 6 6 ; @categories … … 8 8 ; 9 9 ; @param TAB {in}{required} 10 ; The field whose we want to make the hovmo ller map can be 2 kind of thing:10 ; The field whose we want to make the hovmoeller map can be 2 kind of thing: 11 11 ; 1) An array which can be: 12 12 ; * 3d or 4d: array 'xt','yt','zt','t'. The last component is the time. In this case, the array will … … 166 166 ; 167 167 ; @keyword TYPEIN 168 ; allows to specify the type of hovmo ller we want to do168 ; allows to specify the type of hovmoeller we want to do 169 169 ; 'xt','yt','zt','t' 170 170 ; with help of a keyword rather than the argument type. If the argument and the … … 181 181 ; Sebastien Masson (smasson\@lodyc.jussieu.fr) 182 182 ; 27/5/98 183 ; Jerome Vialard (adapting plt to hovmo ller drawing)183 ; Jerome Vialard (adapting plt to hovmoeller drawing) 184 184 ; 2/7/98 185 185 ; Sebastien Masson 14/8/98 (continents,barres) -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltz.pro
r295 r297 161 161 ; 162 162 ; @todo 163 ; seb definition of par maters L.215 à 221163 ; seb definition of paramaters L.215 à 221 164 164 ; 165 165 ;- … … 299 299 !y.range = [-1, 0] 300 300 ;-------------------------------------------------------------- 301 ; extrapolation of data son lands and specifying of the min/max value301 ; extrapolation of data on lands and specifying of the min/max value 302 302 ;-------------------------------------------------------------- 303 303 ; define masknan -
trunk/SRC/ToBeReviewed/PLOTS/DIVERS/checkfield.pro
r296 r297 288 288 ENDCASE 289 289 ;-------------------------------------------------------------- 290 ; is the size of the array compatible with t ehdomain?290 ; is the size of the array compatible with the domain? 291 291 ;-------------------------------------------------------------- 292 292 arr = fitintobox(temporary(arr), nx, ny, nz, firstx, firsty $ -
trunk/SRC/ToBeReviewed/PLOTS/DIVERS/placedessin.pro
r296 r297 22 22 ; 23 23 ; @param dtasize {in}{optional}{type=2 elements vector} 24 ; used by pltv to specify the size of the array that wil be plotted24 ; used by <pro>pltv</pro> to specify the size of the array that will be plotted 25 25 ; and therefore the default yxaspect to be used. 26 26 ; -
trunk/SRC/ToBeReviewed/PLOTS/VECTEUR/ajoutvect.pro
r296 r297 152 152 ; can not calculate the average) 153 153 ;----------------------------------------------------------- 154 ; exten tion of the mask154 ; extension of the mask 155 155 u = u*msku*shift(msku,1,0) 156 156 v = v*mskv*shift(mskv,0,1) -
trunk/SRC/ToBeReviewed/PLOTS/VECTEUR/vecteur.pro
r292 r297 451 451 ; 452 452 if NOT keyword_set(novectref) then BEGIN 453 dx = cmref*cv_cm2normal(0) ; Length of the vector of reference in normal zed coordinates.453 dx = cmref*cv_cm2normal(0) ; Length of the vector of reference in normalized coordinates. 454 454 if keyword_set(vectrefformat) then $ 455 455 normelegende = strtrim(string(normeref, format = vectrefformat), 1)+' ' $ -
trunk/SRC/ToBeReviewed/TRIANGULATION/definetri.pro
r262 r297 4 4 ; Define a triangulation array like <proidl>TRIANGULATE</proidl>. 5 5 ; But in a VERY SIMPLE CASE: 6 ; the points are regular y-gridded on nx*ny array.6 ; the points are regularly-gridded on nx*ny array. 7 7 ; Find a Delaunay triangulation for this set of points is easy: 8 8 ; Points define (nx-1)*(ny-1) rectangles which we can cut in 2 triangles. … … 124 124 ; 125 125 trinumber = 2*(upward-upward/nx) 126 ; ;we define the right triangles126 ; we define the right triangles 127 127 triangles[0, trinumber] = upward 128 128 triangles[1, trinumber] = upward+1 -
trunk/SRC/ToBeReviewed/TRIANGULATION/tracemask.pro
r232 r297 10 10 ; 2d array specifying the mask 11 11 ; 12 ; @param XIN {in}{required} ,12 ; @param XIN {in}{required} 13 13 ; 2d array specifying longitude coordinates. 14 14 ; 15 ; @param YIN {in}{required} ,15 ; @param YIN {in}{required} 16 16 ; 2d array specifying latitude coordinates. 17 17 ; … … 98 98 ; 99 99 IF testvar(var = key_performance) EQ 2 THEN $ 100 print, 'temps tracemask: determination du mask et de s ses coordonnes', systime(1)-tempdeux100 print, 'temps tracemask: determination du mask et de ses coordonnees', systime(1)-tempdeux 101 101 ; 102 102 ; We trace vertical segments: -
trunk/SRC/ToBeReviewed/TRIANGULATION/triangule_c.pro
r231 r297 58 58 ; 59 59 ; @restrictions 60 ; Data swhose we want to do the contour must be disposed in a matrix.60 ; Data whose we want to do the contour must be disposed in a matrix. 61 61 ; On the other hand, in the matrix, the points's arrangement can not be 62 62 ; irregular. If it is, use TRIANGULE. -
trunk/SRC/ToBeReviewed/UTILITAIRE/oups.pro
r238 r297 60 60 61 61 ; 62 ; we erase the last line to reexecute all the program mbut the last order.62 ; we erase the last line to reexecute all the program but the last order. 63 63 ; 64 64 IF strpos(vectjournal[n_elements(vectjournal)-1],';,pos=[' ) NE -1 then $ -
trunk/SRC/ToBeReviewed/WIDGET/AUTOUR_de_XXX/buildcmd.pro
r262 r297 352 352 6:boxzoom = boxextra 353 353 ENDCASE 354 ; Writ ting of this one as a string354 ; Writing of this one as a string 355 355 box = '['+strtrim(boxzoom[0], 1) 356 356 for i = 1, (n_elements(boxzoom)-1) < (3+2*(strpos(type, 'z') EQ -1)) do $ -
trunk/SRC/ToBeReviewed/WIDGET/AUTOUR_de_XXX/changefield.pro
r262 r297 49 49 widget_control, specifieid, set_value = exextra 50 50 ; Change the variable -> Do we need to change the vertical axis 51 ; according to the ty e of points (T or W)?51 ; according to the type of points (T or W)? 52 52 dthlv1id = widget_info(base, find_by_uname = 'dthlv1') 53 53 widget_control, dthlv1id, get_uvalue = dthlv1_uval -
trunk/SRC/ToBeReviewed/WIDGET/AUTOUR_de_XXX/xcreateanim.pro
r231 r297 94 94 widget_control, base, /realize 95 95 ; 96 ; Beginning of the gif file. Writ ting of an emptypicture.96 ; Beginning of the gif file. Writing of an empty picture. 97 97 IF keyword_set(fakecal) THEN date = index1 ELSE date = jul2date(calendar[index1]) 98 98 xxx2ps, /noerase, date1in = date, date2in = date … … 156 156 wset, current_window 157 157 158 ; If we are under x, we try to lauch xanim...158 ; If we are under X, we try to launch xanim... 159 159 160 160 if thisOS NE 'MAC' AND thisOS NE 'WIN' then begin -
trunk/SRC/ToBeReviewed/WIDGET/COMPOUND_WIDGET/cw_domain.pro
r254 r297 265 265 boxzoom[5] = boxzoom[4]+1 266 266 endif 267 ; Now, values and indexes are prop rely defined.267 ; Now, values and indexes are properly defined. 268 268 ; So we can apply them 269 269 widget_control, dthlv1id, set_value = {combobox_select:indice1}
Note: See TracChangeset
for help on using the changeset viewer.