Changeset 118 for trunk/SRC/Obsolete
- Timestamp:
- 06/27/06 17:47:06 (18 years ago)
- Location:
- trunk/SRC/Obsolete
- Files:
-
- 19 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Obsolete/common.pro
r97 r118 1 ;+ 1 2 ; @obsolete 3 ; 4 ; @version $Id$ 5 ; 6 ;- 7 2 8 @all_cm -
trunk/SRC/Obsolete/cp.pro
r114 r118 4 4 ;+ 5 5 ; @file_comments copy files 6 ; 6 7 ; @obsolete 8 ; 7 9 ; @file_comments file_copy should be used instead... 8 10 ; 11 ; @param filenamein {in}{required} 12 ; @param filenameout {in}{required} 13 ; 14 ; @keyword _EXTRA used to pass your keywords 15 ; 9 16 ; @history 10 ; June 2005: Sebastien Masson => obsolete routine 17 ; June 2005: Sebastien Masson => obsolete routine 18 ; 19 ; @version $Id$ 11 20 ; 12 21 ;- -
trunk/SRC/Obsolete/ficdate.pro
r114 r118 21 21 ; @history Jerome Vialard (jv\@lodyc.jussieu.fr) 22 22 ;1/7/98 23 ; 24 ; @version $Id$ 25 ; 23 26 ;- 24 27 ;------------------------------------------------------------ … … 46 49 rien = juldate(vdate) 47 50 ; constitution de la date yymmdd 48 case 1 of 51 case 1 of 49 52 year lt 10: s_year ='0'+string(format='(i1)',year) 50 53 year ge 10 and year lt 100:s_year= string(format='(i2)',year) 51 year ge 100:BEGIN 54 year ge 100:BEGIN 52 55 year = year-1900 53 56 if year LT 10 then s_year ='0'+string(format='(i1)',year) $ -
trunk/SRC/Obsolete/fictype.pro
r114 r118 5 5 ; @file_comments gives fictype ('DA' 'MO' 'AN' 'SE' 'EX') corresponding 6 6 ; to the given date in vairmer format ( yymmdd ) 7 ; 8 ; @param vdate {in}{required} 9 ; @param dim unused 10 ; 7 11 ; @uses common.pro 8 12 ; … … 12 16 ; @history Jerome Vialard (jv\@lodyc.jussieu.fr) 13 17 ; 2/7/98 18 ; 19 ; @version $Id$ 20 ; 14 21 ;- 15 22 ;------------------------------------------------------------ -
trunk/SRC/Obsolete/imprime.pro
r114 r118 3 3 ;------------------------------------------------------------ 4 4 ;+ 5 ; @file_comments 5 ; @file_comments use printps instead... 6 6 ; @obsolete 7 ; 8 ; @param filename {in}{required} 9 ; 10 ; @keyword TRANS 11 ; 12 ; @keyword NB 13 ; 14 ; @uses printps 7 15 ; 8 16 ; @history 9 17 ; June 2005: Sebastien Masson, english version 18 ; 19 ; @version $Id$ 20 ; 10 21 ;- 11 22 ;------------------------------------------------------------ … … 19 30 ; 20 31 thisOS = strupcase(strmid(!version.os_family, 0, 3)) 21 CASE thisOS OF 32 CASE thisOS OF 22 33 'MAC':return 23 34 'WIN':return … … 30 41 0:printps 31 42 1:printps, filename 32 ELSE: BEGIN 43 ELSE: BEGIN 33 44 ras = report('imprime accept only one element: psfilename') 34 45 return 35 END 36 ENDCASE 46 END 47 ENDCASE 37 48 ;------------------------------------------------------------ 38 49 return -
trunk/SRC/Obsolete/jourdsmois.pro
r114 r118 6 6 ; @obsolete 7 7 ; 8 ; @param mois {in}{optional} 9 ; @param annee {in}{optional} 10 ; 8 11 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 9 12 ; June 2005: Sebastien Masson, english version 13 ; 14 ; @version $Id$ 15 ; 10 16 ;- 11 17 ;------------------------------------------------------------ -
trunk/SRC/Obsolete/juldate.pro
r114 r118 6 6 ; @obsolete 7 7 ; 8 ; @file_comments gives julian date equivalent of a date in vairmer 9 ; ( yymmdd or yyyymmdd ) format 8 ; @file_comments gives julian date equivalent of a date in vairmer 9 ; ( yymmdd or yyyymmdd ) format 10 10 ; sets month, day and year to the corresp values 11 11 ; … … 13 13 ; 14 14 ; @param vvdate {in} date de la forme yymmdd ou yyyymmdd 15 ; 16 ; @keyword _EXTRA used to pass your keywords 15 17 ; 16 18 ; @uses common.pro vraidate … … 20 22 ; jour julien de l'annee -1. 21 23 ; COMPATIBLE AVEC L'AN 2000 !!!: une date de la forme yymmdd est 22 ; donvertit sous la forme yyyymmdd a l''aide de vraidate.24 ; convertit sous la forme yyyymmdd a l''aide de vraidate. 23 25 ; Attention, les variables globales year, month day sont attribuees 24 26 ; 25 ; @examples 27 ; @examples 26 28 ; IDL> date = juldate(930124) 27 29 ; 28 30 ; @history Jerome Vialard (jv\@lodyc.jussieu.fr) 29 31 ; 2/7/98 32 ; 33 ; @version $Id$ 34 ; 30 35 ;- 31 36 ;------------------------------------------------------------ … … 47 52 mm = month & dd = day & yy = year 48 53 ;------------------------------------------------------------ 49 ndate = n_elements(vdate) 54 ndate = n_elements(vdate) 50 55 if total(mm EQ 0) EQ ndate then mm[*] = 6 51 56 if total(dd EQ 0) EQ ndate then dd[*] = 15 -
trunk/SRC/Obsolete/lec.pro
r114 r118 20 20 ; 21 21 ; @param nomchamp {in}{required} 2 choix possibles: 22 ; 1) nom de champ Vairmer (chaine de 8 caracteres en majuscule ou 22 ; 1) nom de champ Vairmer (chaine de 8 caracteres en majuscule ou 23 23 ; minuscule commencant par vo ou so). Dans cette methode on saute directement 24 24 ; d'en-tete en en-tete jusqu'a trouver le bon fichier. … … 31 31 ; 32 32 ; 33 ; @keyword /ANOM {in} type du fichier vairmer par rapport auquel on doit calculer 33 ; @keyword /ANOM {in} type du fichier vairmer par rapport auquel on doit calculer 34 34 ; l'anomalie ('EX','AN','SE','MO','') 35 35 ; … … 37 37 ; ds ce cas en input on met seulement 'vo' ou 'so' la fonction retourne le 38 38 ; nombre de fichiers lus. 39 ; 40 ; @keyword BOITE 39 41 ; 40 42 ; @keyword /EXPANOM {in} si on calcule l'anom par rapport a une exper … … 60 62 ; celle du sous domaine definit par domdef (nx,ny,nz). avec /TOUT le 61 63 ; champ a la taille du domaine qui a etait selectionne pour la 62 ; cession en cours (jpi,jpj,jpk). 63 ; pour les sous domaines cf: 64 ; cession en cours (jpi,jpj,jpk). 65 ; pour les sous domaines cf: 64 66 ; http://www.ipsl.jussieu.fr/~smlod/sousdomaine.html 65 67 ; Retourne -1 en cas d'erreur. 66 68 ; 67 ; @uses common.pro isnumber.pro fivardate.pro 69 ; @uses common.pro isnumber.pro fivardate.pro 68 70 ; 69 71 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) … … 76 78 ; Sebastien Masson (decoupe pour les sous domaines...) 77 79 ; 2/99 80 ; 81 ; @version $Id$ 82 ; 78 83 ;;- 79 84 ;-------------------------------------------------------------- … … 141 146 ; definition de la taille du fichier a aller chercher: jpidta,jpjdta,jpkdta... 142 147 ;--------------------------------------------------------------------- 143 if n_elements(jpidta) EQ 0 THEN BEGIN 148 if n_elements(jpidta) EQ 0 THEN BEGIN 144 149 if n_elements(ixmindta) EQ 0 OR n_elements(ixmaxdta) EQ 0 then $ 145 150 jpidta = jpiglo else jpidta = ixmaxdta-ixmindta+1 146 151 endif 147 if n_elements(jpjdta) EQ 0 THEN BEGIN 152 if n_elements(jpjdta) EQ 0 THEN BEGIN 148 153 if n_elements(iymindta) EQ 0 OR n_elements(iymaxdta) EQ 0 then $ 149 154 jpjdta = jpjglo else jpjdta = iymaxdta-iymindta+1 150 155 endif 151 if n_elements(jpkdta) EQ 0 THEN BEGIN 156 if n_elements(jpkdta) EQ 0 THEN BEGIN 152 157 if n_elements(izmindta) EQ 0 OR n_elements(izmaxdta) EQ 0 then $ 153 158 jpkdta = jpkglo else jpkdta = izmaxdta-izmindta+1 … … 267 272 ; on shift z si key_shift est defininit 268 273 ;--------------------------------------------------------------------- 269 if n_elements(key_shift) NE 0 THEN BEGIN 274 if n_elements(key_shift) NE 0 THEN BEGIN 270 275 if dim EQ 'SO' then z = shift(z,key_shift, 0) $ 271 276 ELSE z = shift(z,key_shift, 0, 0) … … 275 280 ; du zoom: nx,ny nz 276 281 ;--------------------------------------------------------------------- 277 if NOT keyword_set(tout) then BEGIN 282 if NOT keyword_set(tout) then BEGIN 278 283 ;------------------------------------------------------------- 279 284 ; changement de domaine 280 285 ;------------------------------------------------------------- 281 if keyword_set(boite) then BEGIN 286 if keyword_set(boite) then BEGIN 282 287 Case 1 Of 283 288 N_Elements(Boite) Eq 1:bte=[lon1, lon2, lat1, lat2, 0.,boite[0]] … … 290 295 oldboite = [lon1, lon2, lat1, lat2, prof1, prof2] 291 296 domdef, bte,GRILLE=vargrid 292 ENDIF 297 ENDIF 293 298 ;------------------------------------------------------------- 294 299 grille,mask,glam,gphi,gdep,nx,ny,nz,premierx,premiery,premierz,dernierx,derniery,dernierz … … 296 301 if dim EQ 'SO' then z = z[premierx:dernierx, premiery:derniery] $ 297 302 ELSE z = z[premierx:dernierx, premiery:derniery, premierz:dernierz] 298 ENDIF ELSE BEGIN 303 ENDIF ELSE BEGIN 299 304 case vargrid OF ; on recupere le mask en entier ds le cas ou /TOUT 300 305 'U':mask = umask() ; n''est pas active et on le choisit en fontion 301 306 'T':mask = tmask ; de la valeur de vargrid 302 'W':mask = tmask 307 'W':mask = tmask 303 308 'V':mask = vmask() 304 309 'F':mask = fmask() 305 310 ENDCASE 306 ENDELSE 311 ENDELSE 307 312 ;--------------------------------------------------------------------- 308 313 ; calcul d'une anomalie si le keyword anom est active … … 311 316 case anom of 312 317 'EX' : adate = 0 313 'AN' : adate = floor(date/10000)*10000 318 'AN' : adate = floor(date/10000)*10000 314 319 'SE' : adate = floor(date - floor(date/10000)*10000)/100 * 100 315 320 'MO' : adate = floor(date/100)*100 … … 327 332 ;--------------------------------------------------------------------- 328 333 IF n_elements(valmask) EQ 0 THEN valmask = 1e20 329 if dim EQ 'SO' then BEGIN 334 if dim EQ 'SO' then BEGIN 330 335 terre = where(mask[*,*,0] EQ 0) 331 if terre[0] NE -1 then z[terre] = valmask 332 ENDIF ELSE BEGIN 336 if terre[0] NE -1 then z[terre] = valmask 337 ENDIF ELSE BEGIN 333 338 terre = where(mask[*,*,0] EQ 0) 334 339 if terre[0] NE -1 then z[where(mask EQ 0)] = valmask 335 ENDELSE 340 ENDELSE 336 341 ;--------------------------------------------------------------------- 337 342 free_lun,numlec -
trunk/SRC/Obsolete/lect.pro
r114 r118 3 3 ;------------------------------------------------------------ 4 4 ;+ 5 ; 5 ; 6 6 ; @file_comments 7 7 ; lit les fichiers Vairmer de date1 a date2 et en sort un tableau … … 19 19 ; IDL> resultat=lec('nom_Vairmer',date1,date2,'nomexp','direc',BOITE=boite) 20 20 ; 21 ; @param nomchamp {in}{required} nom de champ Vairmer (chaine de 8 caracteres commencant22 ; 23 ; 24 ; 21 ; @param nomchamp {in}{required} nom de champ Vairmer (chaine de 8 caracteres 22 ; commencant par VO ou SO 23 ; @param date1 {in}{required} date de depart de la serie temporelle a 24 ; @param date2 {in}{required} date de fin (date2) de la serie temporelle a 25 25 ; extraire 26 ; @param nomexp {in}{required} nom de l' experience a lire (prefix pardefaut) 27 ; 28 ; @param direc {in}{required} 'x' 'y' 'z' 'xy' 'xz' 'yz' 'xyz' 'xt' 'yt' 'zt' 'xyt' 'xzt' 29 ; 'yzt' 'xyzt' directions selon lesquelles effectuer les moyennes 30 ; (si rien n'est donne on n'effectue pas de moyenne) 31 ; 32 ; 33 ; 34 ; @keyword boite {in} boite sur laquelle integrer (par defaut tt le domaine) 35 ; @keyword anom {in} type de fichiers ('SE','AN','') a relire pour calc une 36 ; anomalie 37 ; @keyword expanom {in} experience pour laquelle on veut calculer une anomalie 26 ; @param nomexp {in}{required} nom de l' experience a lire (prefix pardefaut) 27 ; 28 ; @param direc {in}{required} 29 ; 'x' 'y' 'z' 'xy' 'xz' 'yz' 'xyz' 'xt' 'yt' 'zt' 'xyt' 'xzt' 30 ; 'yzt' 'xyzt' directions selon lesquelles effectuer les moyennes 31 ; (si rien n'est donne on n'effectue pas de moyenne) 32 ; 33 ; @keyword boite {in} boite sur laquelle integrer (par defaut tt le domaine) 34 ; @keyword anom {in} type de fichiers ('SE','AN','') a relire pour calc une 35 ; anomalie 36 ; @keyword expanom {in} experience pour laquelle on veut calculer une anomalie 38 37 ; (par defaut la meme que nomexp) 39 38 ; … … 42 41 ; @history Jerome Vialard (jv\@lodyc.jussieu.fr) 43 42 ; 2/7/98 43 ; 44 ; @version $Id$ 45 ; 44 46 ;- 45 47 ;-------------------------------------------------------------- 46 48 ;-------------------------------------------------------------- 47 49 ;-------------------------------------------------------------- 48 ; 1 2 3 4 5 50 ; 1 2 3 4 5 49 51 function lect, nomchamp,date1,date2,nomexp,direc,BOITE=boite, $ 50 52 ANOM=anom,EXPANOM=expanom … … 68 70 if fictype(date1) ne fictype(date2) then $ 69 71 return, report('Les deux dates doivent correspondre au meme type de fic vairmer') 70 fictyp = fictype(date1) 72 fictyp = fictype(date1) 71 73 ;------------------------------------------------------------- 72 74 ; creation du nom du fichier … … 77 79 if (keyword_set(anom)) then ficname=ficname +'.'+anom 78 80 if (keyword_set(expanom)) then ficname=ficname +'-'+expanom 79 case n_elements(boite) of 81 case n_elements(boite) of 80 82 4 : box = strcompress(string(format='(i4,"_",i4,"_",i4,"_",i4)',boite) $ 81 83 ,/remove_all) … … 90 92 ; 91 93 ; structure du fichier : 92 ; jpt (valeur de la dim temporelle), dimtableau (dimension du tableau) 94 ; jpt (valeur de la dim temporelle), dimtableau (dimension du tableau) 93 95 ; dimttab[0], dimttab[1], ... (valeur des dim ) 94 96 ; time (axe des tps), ttab (tableau a lire) … … 101 103 dimtableau = long(1) 102 104 readu, numlec, jpt,dimtableau 103 case dimtableau of 105 case dimtableau of 104 106 1 : begin 105 107 n1 = long(1) … … 145 147 ; changement de domaine 146 148 ;------------------------------------------------------------- 147 if keyword_set(boite) then BEGIN 149 if keyword_set(boite) then BEGIN 148 150 Case 1 Of 149 151 N_Elements(Boite) Eq 1:bte=[lon1, lon2, lat1, lat2, 0.,boite[0]] … … 156 158 oldboite = [lon1, lon2, lat1, lat2, prof1, prof2] 157 159 domdef, bte 158 ENDIF 159 ;------------------------------------------------------------- 160 ; Boucle de lecture des fichiers 160 ENDIF 161 ;------------------------------------------------------------- 162 ; Boucle de lecture des fichiers 161 163 ;------------------------------------------------------------- 162 164 case fictyp of … … 190 192 ; Moyenne du champs tab 191 193 ;-------------------------------------------------------- 192 IF n_params() EQ 5 THEN if direc NE '' then BEGIN 193 if nx EQ 1 OR ny EQ 1 OR nz EQ 1 THEN BEGIN 194 IF n_params() EQ 5 THEN if direc NE '' then BEGIN 195 if nx EQ 1 OR ny EQ 1 OR nz EQ 1 THEN BEGIN 194 196 if string(format='(a2)',nomchamp) EQ 'SO' then tab = reform(tab, nx, ny, /over) $ 195 197 ELSE tab = reform(tab, nx, ny, nz, /over) 196 ENDIF 198 ENDIF 197 199 tab = moyenne(tab,direc) 198 200 endif 199 201 if (jpt eq 1) then begin 200 202 ttab = tab 201 endif else BEGIN 203 endif else BEGIN 202 204 ttab = colle(ttab, tab, (size(tab))[0]+1 ) 203 205 endelse … … 217 219 endcase 218 220 ;------------------------------------------------------------ 219 ; Fin de boucle de lecture des fichiers 221 ; Fin de boucle de lecture des fichiers 220 222 ;------------------------------------------------------------ 221 223 vdat=long(10000)*year+long(100)*month+day … … 229 231 taille = size(ttab) 230 232 writeu, numlec, long(jpt),long(taille[0]) 231 case taille[0] of 233 case taille[0] of 232 234 1 : writeu, numlec,long(taille[1]) 233 235 2 : writeu, numlec,long(taille[1]),long(taille[2]) … … 258 260 close, /all 259 261 ; 260 IF keyword_set(key_performance) THEN print, 'temps lect', systime(1)-tempsun 262 IF keyword_set(key_performance) THEN print, 'temps lect', systime(1)-tempsun 261 263 ; 262 264 return, ttab -
trunk/SRC/Obsolete/meshlec.pro
r114 r118 4 4 ;+ 5 5 ; 6 ; @file_comments lecture du mask de des sorties d'OPA. les sources se trouvent ds les6 ; @file_comments lecture du mask des sorties d'OPA. les sources se trouvent ds les 7 7 ; repertoires sur maia du type: 8 8 ; /nom_exp/RESTARTS 9 9 ; 10 10 ; @obsolete 11 ; 11 12 ; @examples 12 13 ; IDL> meshmask[,' nomfich'] 13 14 ; 14 ; @param nomfich {in}{required} string, c''est le nom du fichier a lire. Par defaut, c''est meshmask.15 ; 16 ; 17 ; 18 ; 19 ; 15 ; @param nomfich {in}{required} string, c''est le nom du fichier a lire. Par defaut, c''est meshmask. 16 ; 17 ; 18 ; @keyword GLAMBOUNDARY {in} un vecteur de 2 elements specifaint le min et le 19 ; max qui doivent etre imposes en longitude (obligatoire si le 20 ; tableau depasse 360 degres) 20 21 ; 21 22 ; @keyword /pasblabla {in} pour suprimer les blablas … … 24 25 ; precision. ce Mot clef est maintenant active 25 26 ; automatiquement. 27 ; 28 ; @keyword GETDIMENSIONS 26 29 ; 27 30 ; @uses common.pro … … 35 38 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 36 39 ; Marina Levy : lecture en double precision (cas calcul sur shine) 40 ; 41 ; @version $Id$ 42 ; 37 43 ;- 38 44 ;------------------------------------------------------------ … … 109 115 ; sizefile8:double = 1 110 116 ; sizefile4:double = 0 111 ; ELSE:BEGIN 117 ; ELSE:BEGIN 112 118 ; nothing = report('The OPA Mesh file as not the good size!') 113 119 ; free_lun,numlec 114 120 ; close, numlec 115 ; return 121 ; return 116 122 ; END 117 123 ; endcase … … 150 156 ;------------------------------------------------------- 151 157 152 IF NOT keyword_set(double) THEN BEGIN 158 IF NOT keyword_set(double) THEN BEGIN 153 159 z3d = fltarr(jpiglo, jpjglo, jpkglo) 154 160 z2d = fltarr(jpiglo, jpjglo) 155 161 z1d = fltarr(jpkglo) 156 ENDIF ELSE BEGIN 162 ENDIF ELSE BEGIN 157 163 z3d = dblarr(jpiglo, jpjglo, jpkglo) 158 164 z2d = dblarr(jpiglo, jpjglo) 159 165 z1d = dblarr(jpkglo) 160 ENDELSE 166 ENDELSE 161 167 162 168 if not keyword_set(pasblabla) then print, ' ' … … 175 181 176 182 if not keyword_set(pasblabla) then print, ' ' 177 readu, numlec, tab, z2d 183 readu, numlec, tab, z2d 178 184 GPHIT=float(z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh]) 179 185 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur GPHIT[25,31]: ',GPHIT[25,31] 180 readu, numlec, tab, z2d 186 readu, numlec, tab, z2d 181 187 GPHIU=float(z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh]) 182 188 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur GPHIU[25,31]: ',GPHIU[25,31] 183 readu, numlec, tab, z2d 189 readu, numlec, tab, z2d 184 190 GPHIV=float(z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh]) 185 191 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur GPHIV[25,31]: ',GPHIV[25,31] 186 readu, numlec, tab, z2d 192 readu, numlec, tab, z2d 187 193 GPHIF=float(z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh]) 188 194 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur GPHIF[25,31]: ',z2d[25,31] … … 217 223 218 224 if not keyword_set(pasblabla) then print, ' ' 219 readu, numlec, tab, z3d 225 readu, numlec, tab, z3d 220 226 TMASK=byte(z3d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh,izminmesh:izmaxmesh]) 221 227 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur TMASK[25,5,0]: ',TMASK[25,5,0] 222 readu, numlec, tab, z3d 228 readu, numlec, tab, z3d 223 229 UMASKred=byte(z3d[ixmaxmesh,iyminmesh:iymaxmesh,izminmesh:izmaxmesh]) 224 230 umaskred = reform(umaskred) 225 231 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur UMASK[25,5,0]: ',z3d[25,5,0] 226 readu, numlec, tab, z3d 232 readu, numlec, tab, z3d 227 233 VMASKred=byte(z3d[ixminmesh:ixmaxmesh,iymaxmesh,izminmesh:izmaxmesh]) 228 234 vmaskred = reform(vmaskred) 229 235 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur VMASK[25,5,0]: ',z3d[25,5,0] 230 readu, numlec, tab, z3d 236 readu, numlec, tab, z3d 231 237 fmaskredy=byte(z3d[ixmaxmesh,iyminmesh:iymaxmesh,izminmesh:izmaxmesh]) 232 238 coast = where(fmaskredy NE 0 and fmaskredy NE 1) … … 240 246 ; 241 247 if not keyword_set(pasblabla) then print, ' ' 242 readu, numlec, tab, z2d 248 readu, numlec, tab, z2d 243 249 ;FF=z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh] 244 250 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur FF[25,5]: ',z2d[25,5] … … 246 252 GDEPT=float(z1d[izminmesh:izmaxmesh]) 247 253 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur GDEPT[1]: ',GDEPT[1] 248 readu, numlec, tab, z1d 254 readu, numlec, tab, z1d 249 255 GDEPW=float(z1d[izminmesh:izmaxmesh]) 250 256 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur GDEPW[1]: ',GDEPW[1] 251 readu, numlec, tab, z1d 257 readu, numlec, tab, z1d 252 258 E3T=float(z1d[izminmesh:izmaxmesh]) 253 259 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur E3T[3]: ',E3T[3] 254 readu, numlec, tab, z1d 260 readu, numlec, tab, z1d 255 261 E3W=float(z1d[izminmesh:izmaxmesh]) 256 262 if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur E3W[3]: ',E3W[3] … … 260 266 ;------------------------------------------------------- 261 267 ;------------------------------------------------------- 262 ; bornes de glam qui ne doivent pas depasser 360 degres.. .268 ; bornes de glam qui ne doivent pas depasser 360 degres.. 263 269 ;------------------------------------------------------- 264 270 ; minglam = min(glamt, max = maxglam) … … 329 335 if not keyword_set(pasblabla) then print,'lecture '+nomfich+' finie' 330 336 widget_control, noticebase, bad_id = toto, /destroy 331 if keyword_set(key_performance) THEN print, 'temps meshlec', systime(1)-tempsun 337 if keyword_set(key_performance) THEN print, 'temps meshlec', systime(1)-tempsun 332 338 333 339 return -
trunk/SRC/Obsolete/ncdf_meshlec.pro
r114 r118 6 6 ; @obsolete 7 7 ; 8 ; @param filename 9 ; 10 ; @keyword _EXTRA used to pass your keywords 11 ; 8 12 ; @history 9 13 ; Aug. 2005, Sebastien Masson: switch to ncdf_meshread 14 ; 15 ; @version $Id$ 16 ; 10 17 ;- 11 18 ;------------------------------------------------------------ -
trunk/SRC/Obsolete/nlec.pro
r114 r118 15 15 ; 16 16 ; 17 ; 18 ; 19 ; @keyword DIREC:'x' 'y' 'z' 'xy' 'xz' 'yz' 'xyz' directions selon lesquelles17 ; @keyword BOITE boite sur laquelle integrer (par defaut tt le domaine) 18 ; 19 ; @keyword DIREC 'x' 'y' 'z' 'xy' 'xz' 'yz' 'xyz' directions selon lesquelles 20 20 ; effectuer les moyennes (si rien n'est donne on n'effectue pas de moyenne) 21 21 ; 22 22 ; @keyword GRILLE impose la grille a laquelle est rapporte le champ. rq permet 23 ; d'aller plus vite ds la lecture 23 ; d'aller plus vite ds la lecture 24 24 ; 25 25 ; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 26 ; 27 ; @uses common.pro vraidate juldate nlec5j nlecserie 26 ; @keyword STRUCTURE 27 ; @keyword SEUILMIN 28 ; @keyword SEUILMAX 29 ; @keyword NAN 30 ; 31 ; @keyword _EXTRA used to pass your keywords 32 ; 33 ; @uses common.pro vraidate juldate nlec5j nlecserie 28 34 ; 29 35 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 30 36 ; 14/8/98 37 ; 38 ; @version $Id$ 39 ; 31 40 ;- 32 41 ;------------------------------------------------------------ 33 42 ;REF.07.790101.grid.T.nc 34 ;REF.07.790101.grid.U.nc 35 ;REF.07.790101.grid.V.nc 43 ;REF.07.790101.grid.U.nc 44 ;REF.07.790101.grid.V.nc 36 45 ;REF.07.790101.grid.W.nc 37 46 ;------------------------------------------------------------ … … 51 60 ;------------------------------------------------------------ 52 61 case n_params() of 53 1:BEGIN 62 1:BEGIN 54 63 year=year+1900*(year ne 0 and year ne -1 and year lt 100) 55 64 date = day+100*month+10000*year 56 65 end 57 2:BEGIN 66 2:BEGIN 58 67 if size(debut, /tname) EQ 'STRING' then begin 59 68 prefix = strupcase(debut) … … 150 159 datejul = 5*numsortie-1+julday(1, 1, year) 151 160 if n_elements(date2) ne 0 then date2jul = 5*numsortie2-1+julday(1, 1, year2) $ 152 ELSE date2jul = datejul 161 ELSE date2jul = datejul 153 162 res = nlecserie(nom, datejul, date2jul,BOITE=boite,GRILLE=grille,TOUT=tout, _EXTRA = ex) 154 IF res[0] EQ -1 THEN BEGIN 163 IF res[0] EQ -1 THEN BEGIN 155 164 dejaserie = 1 156 165 GOTO, passerie 157 ENDIF 166 ENDIF 158 167 ENDIF ELSE BEGIN 159 168 passerie: … … 163 172 dejapasserie = 1 164 173 GOTO, serie 165 ENDIF 174 ENDIF 166 175 endelse 167 176 ENDIF ELSE BEGIN … … 233 242 ; moyenne eventuelle 234 243 ;------------------------------------------------------------ 235 IF keyword_set(direc) THEN BEGIN 244 IF keyword_set(direc) THEN BEGIN 236 245 IF jpt EQ 1 THEN res = moyenne(res,direc,BOITE=boite, nan = notan) $ 237 246 ELSE res = grossemoyenne(res,direc,BOITE=boite, nan = notan) 238 ENDIF 247 ENDIF 239 248 ;------------------------------------------------------------ 240 249 ; mise en placer des parametres pour le trace... … … 248 257 if keyword_set(structure) then res = {tab:res, grille:vargrid, unite:varunit, experience:varexp, nom:varname} 249 258 ;------------------------------------------------------------ 250 if keyword_set(key_performance) THEN print, 'temps nlec', systime(1)-tempsun 259 if keyword_set(key_performance) THEN print, 'temps nlec', systime(1)-tempsun 251 260 time = time-2l 252 261 return,res … … 254 263 ;------------------------------------------------------------ 255 264 end 256 257 258 259 260 261 262 263 264 265 266 267 268 269 -
trunk/SRC/Obsolete/nlec5j.pro
r114 r118 6 6 ; contiennent les sorties a 5j regroupees par type de grille par 7 7 ; paquets de 6 mois. 8 ; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS 8 ; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS 9 9 ; @obsolete 10 10 ; … … 16 16 ; @param nom {in}{required} nom du tableau contenant le champ que l'on veut (le 17 17 ; trouver avec ncdflec,fichier,/var 18 ; @param numsortie {in}{required}le numero du pas de temps que l'on veut sortir du 18 ; 19 ; @param numsortie {in}{required} 20 ; le numero du pas de temps que l'on veut sortir du 19 21 ; fichier (compte a partir de 1) a partir de year 22 ; 20 23 ; @param nbretps {in}{required} nombre de pas de temps a extraire 21 24 ; 22 ; 23 ; @keyword BOITE boite sur laquelle integrer (par defaut tt le domaine) 25 ; @keyword BOITE boite sur laquelle integrer (par defaut tt le domaine) 24 26 ; 25 27 ; @keyword GRILLE impose la grille a laquelle est rapporte le champ. rq permet 26 ; d'aller plus vite ds la lecture 28 ; d'aller plus vite ds la lecture 27 29 ; 28 30 ; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 29 31 ; 30 ; @returns tableau 2d (qd on ne demande pas de serie) ou 3d ou 4d ds le cas dune serie 31 ; 32 ; @keyword _EXTRA used to pass your keywords 33 34 ; @returns tableau 2d (qd on ne demande pas de serie) ou 3d ou 4d ds le cas dune serie 35 ; 32 36 ; @uses common.pro 33 37 ; … … 35 39 ; 36 40 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 41 ; 42 ; 43 ; @version $Id$ 37 44 ; 38 45 ;- … … 51 58 ; gestion du nom du fichier 52 59 ;------------------------------------------------------------ 53 if numsortie le 36 then BEGIN 60 if numsortie le 36 then BEGIN 54 61 mmdd='0101' 55 62 numsort=numsortie … … 84 91 ras = report('LES FICHIERS: '+iodir+prefix+'.'+numfich+'.'+s_date+'.grid... n''existe pas.') 85 92 return, -1 86 ENDIF 93 ENDIF 87 94 for i=0,n_elements(quelsfichiers)-1 do begin 88 95 IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+quelsfichiers[i]+' > /dev/null' … … 94 101 vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille 95 102 goto, grilletrouvee 96 ENDIF 103 ENDIF 97 104 endfor 98 105 ncdf_close,cdfid … … 132 139 premiery = 0 133 140 premierz = 0 134 endif else BEGIN 141 endif else BEGIN 135 142 ;------------------------------------------------------------ 136 143 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 137 144 ; + on recupere la dim du no9uveau domaine 138 145 ;------------------------------------------------------------ 139 if keyword_set(boite) then BEGIN 146 if keyword_set(boite) then BEGIN 140 147 Case N_Elements(Boite) Of 141 148 1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex … … 148 155 ENDIF 149 156 grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery, premierz, dernierx, derniery, dernierz 150 ENDELSE 157 ENDELSE 151 158 ;------------------------------------------------------------ 152 159 ; determination du nombre de pas de tps a extraire ds ce fichier … … 169 176 res = [[[res]] $ 170 177 ,[[nlec5j(nom,1+36*(mmdd eq '0101'),nbretps-nt,tout=tout,GRILLE=vargrid,BOITE=boite)]]] $ 171 ELSE BEGIN 178 ELSE BEGIN 172 179 res=[res[*] $ 173 180 ,(nlec5j(nom,1+36*(mmdd eq '0101'),nbretps-nt,tout=tout,GRILLE=vargrid, BOITE=boite))[*]] 174 res=reform(res,nx,ny,nz,nbretps, /over) 175 ENDELSE 181 res=reform(res,nx,ny,nz,nbretps, /over) 182 ENDELSE 176 183 ncdf_varget,cdfid,'time_counter',temps, offset = [numsort-1],count=[nt] 177 184 time = [long(temps)+julday(1, 5, 1979), time] -
trunk/SRC/Obsolete/nlecan.pro
r114 r118 4 4 ;+ 5 5 ; @file_comments lit les moyennes annuelles 6 ; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS 6 ; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS 7 7 ; 8 8 ; @obsolete 9 9 ; @categories lecture de NETCDF 10 10 ; 11 ; @examples 11 ; @examples 12 12 ; IDL> res=nlecan(nom, numsortie79, nbretps) 13 13 ; … … 21 21 ; 22 22 ; @keyword GRILLE impose la grille a laquelle est rapporte le champ. rq permet 23 ; d'aller plus vite ds la lecture 23 ; d'aller plus vite ds la lecture 24 24 ; 25 25 ; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 26 26 ; 27 ; @keyword _EXTRA used to pass your keywords 28 ; 27 29 ; @returns tableau 2d (qd la serie ne fait que 1 pas de temps) ou 28 30 ; 3d (valable ds 1 premier tps que pour les tableaux 2d) 29 ; 31 ; 30 32 ; @uses common.pro 31 33 ; … … 33 35 ; 34 36 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 37 ; 38 ; 39 ; @version $Id$ 35 40 ; 36 41 ;- … … 55 60 vargrid=grille 56 61 nomfich=prefix+'.moyenne.annuelle.grid.'+vargrid+'.nc' 57 IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+iodir+nomfich+' > /dev/null' 62 IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+iodir+nomfich+' > /dev/null' 58 63 cdfid=ncdf_open(iodir+nomfich) 59 64 varcontient=ncdf_varinq(cdfid,nom) … … 72 77 vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille 73 78 goto, grilletrouvee 74 ENDIF 79 ENDIF 75 80 endfor 76 81 ncdf_close,cdfid … … 121 126 premiery = 0 122 127 premierz = 0 123 endif else BEGIN 128 endif else BEGIN 124 129 ;------------------------------------------------------------ 125 130 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 126 131 ; + on recupere la dim du no9uveau domaine 127 132 ;------------------------------------------------------------ 128 if keyword_set(boite) then BEGIN 133 if keyword_set(boite) then BEGIN 129 134 Case N_Elements(Boite) Of 130 135 1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex … … 137 142 ENDIF 138 143 grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery, premierz 139 ENDELSE 144 ENDELSE 140 145 ;------------------------------------------------------------ 141 146 if varcontient.ndims eq 3 then $ -
trunk/SRC/Obsolete/nlecmois.pro
r114 r118 12 12 ; 13 13 ; @keyword GRILLE impose la grille a laquelle est rapporte le champ. rq permet 14 ; d'aller plus vite ds la lecture 14 ; d'aller plus vite ds la lecture 15 15 ; 16 16 ; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 17 ; @keyword _EXTRA used to pass your keywords 18 17 19 ; @uses common.pro 18 20 ; 19 21 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 22 ; 23 ; @version $Id$ 20 24 ; 21 25 ;- … … 59 63 vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille 60 64 goto, grilletrouvee 61 ENDIF 65 ENDIF 62 66 endfor 63 67 ncdf_close,cdfid … … 85 89 premiery=0 86 90 premierz=0 87 endif else BEGIN 91 endif else BEGIN 88 92 ;------------------------------------------------------------ 89 93 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 90 94 ; + on recupere la dim du no9uveau domaine 91 95 ;------------------------------------------------------------ 92 if keyword_set(boite) then BEGIN 96 if keyword_set(boite) then BEGIN 93 97 Case N_Elements(Boite) Of 94 98 1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex … … 101 105 ENDIF 102 106 grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery, premierz 103 ENDELSE 107 ENDELSE 104 108 ;------------------------------------------------------------ 105 109 ; determination du nombre de pas de tps a extraire ds ce fichier … … 121 125 res = [[[res]] $ 122 126 ,[[nlecmois(nom,numsortie79+nt,nbretps-nt,tout=tout,GRILLE=vargrid,BOITE=boite)]]] $ 123 else BEGIN 127 else BEGIN 124 128 res=[res[*] $ 125 ,(nlecmois(nom,numsortie79+nt,nbretps-nt,tout=tout,GRILLE=vargrid,BOITE=boite))[*]] 126 127 res=reform(res,nx,ny,nz,nbretps, /over) 128 ENDELSE 129 ,(nlecmois(nom,numsortie79+nt,nbretps-nt,tout=tout,GRILLE=vargrid,BOITE=boite))[*]] 130 131 res=reform(res,nx,ny,nz,nbretps, /over) 132 ENDELSE 129 133 ncdf_varget,cdfid,'time_counter',temps, offset = [numsortie-1],count=[nt] 130 134 time = [long(temps)+julday(12, 31, 1978), time] -
trunk/SRC/Obsolete/nlecsaison.pro
r114 r118 4 4 ;+ 5 5 ; @file_comments lit les moyennes annuelles 6 ; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS 6 ; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS 7 7 ; @obsolete 8 8 ; … … 21 21 ; 22 22 ; @keyword GRILLE impose la grille a laquelle est rapporte le champ. rq permet 23 ; d'aller plus vite ds la lecture 23 ; d'aller plus vite ds la lecture 24 24 ; 25 25 ; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 26 26 ; 27 ; @keyword _EXTRA used to pass your keywords 28 27 29 ; @returns tableau 2d (qd la serie ne fait que 1 pas de temps) ou 28 30 ; 3d (valable ds 1 premier tps que pour les tableaux 2d) 29 ; 31 ; 30 32 ; @uses common.pro 31 33 ; … … 33 35 ; 34 36 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 37 ; 38 ; 39 ; @version $Id$ 35 40 ; 36 41 ;- … … 72 77 vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille 73 78 goto, grilletrouvee 74 ENDIF 79 ENDIF 75 80 endfor 76 81 ncdf_close,cdfid … … 118 123 premiery = 0 119 124 premierz = 0 120 endif else BEGIN 125 endif else BEGIN 121 126 ;------------------------------------------------------------ 122 127 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 123 128 ; + on recupere la dim du no9uveau domaine 124 129 ;------------------------------------------------------------ 125 if keyword_set(boite) then BEGIN 130 if keyword_set(boite) then BEGIN 126 131 CASE N_Elements(Boite) Of 127 132 1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex … … 134 139 ENDIF 135 140 grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery,premierz 136 ENDELSE 141 ENDELSE 137 142 ;------------------------------------------------------------ 138 143 if varcontient.ndims eq 3 then $ -
trunk/SRC/Obsolete/nlecserie.pro
r114 r118 4 4 ;+ 5 5 ; @file_comments lit les series temporelles se rapportant a une variable 6 ; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS 6 ; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS 7 7 ; @obsolete 8 8 ; 9 9 ; @categories lecture de NETCDF 10 10 ; 11 ; @examples 11 ; @examples 12 12 ; IDL> res=nlecserie(nom, date1, date2) 13 13 ; … … 20 20 ; 21 21 ; @keyword GRILLE impose la grille a laquelle est rapporte le champ. rq permet 22 ; d'aller plus vite ds la lecture 22 ; d'aller plus vite ds la lecture 23 23 ; 24 24 ; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 25 25 ; 26 ; @keyword _EXTRA used to pass your keywords 27 ; 26 28 ; @returns tableau 2d (qd la serie ne fait que 1 pas de temps) ou 27 29 ; 3d (valable ds 1 premier tps que pour les tableaux 2d) 28 ; 30 ; 29 31 ; @uses common.pro 30 32 ; … … 32 34 ; 33 35 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 36 ; 37 ; @version $Id$ 34 38 ; 35 39 ;- … … 56 60 ; print, 'LES FICHIERS: '+iodir+prefix+'.'+nom+'.grid... n''existe pas.' 57 61 return, -1 58 ENDIF ELSE BEGIN 62 ENDIF ELSE BEGIN 59 63 vargrid=strmid(quelsfichiers[0],strpos(quelsfichiers[0],'grid')+5,1) ;nom de grille 60 64 IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+quelsfichiers[0]+' > /dev/null' 61 65 cdfid=ncdf_open(quelsfichiers[0]) 62 ENDELSE 66 ENDELSE 63 67 ;------------------------------------------------------------ 64 68 ; lecture de certains attributs et de l'axe des temps … … 116 120 premiery = 0 117 121 premierz = 0 118 endif else BEGIN 122 endif else BEGIN 119 123 ;------------------------------------------------------------ 120 124 ; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 121 125 ; + on recupere la dim du no9uveau domaine 122 126 ;------------------------------------------------------------ 123 if keyword_set(boite) then BEGIN 127 if keyword_set(boite) then BEGIN 124 128 Case N_Elements(Boite) Of 125 129 1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex -
trunk/SRC/Obsolete/vairdate.pro
r114 r118 9 9 ; @file_comments you better use jul2date 10 10 ; 11 ; @examples 11 ; @examples 12 12 ; IDL> vdate = vairdate(1755087) 13 13 ; 14 14 ; @param jdate {in}{required} date en jours juliens 15 ;16 15 ; 17 16 ; @keyword MENSUEL a activer si on veut pour que les dates dont le jour est … … 23 22 ; 0. par ex: 24 23 ; 19990601 -> 19990000 24 ; 25 ; @keyword _EXTRA used to pass your keywords 25 26 ; 26 27 ; @returns vdate date vairmer … … 36 37 ; d''utiliser vairdate avec des scalaire ou des 37 38 ; tableaux 39 ; 40 ; @version $Id$ 41 ; 38 42 ;- 39 43 ;------------------------------------------------------------ … … 49 53 index = where(year eq -1) 50 54 if index[0] NE -1 then year[index] = 0 51 if keyword_set(mensuel) THEN BEGIN 55 if keyword_set(mensuel) THEN BEGIN 52 56 index = where(day EQ 15) 53 57 if index[0] NE -1 then day[index] = 0 54 58 endif 55 if keyword_set(annuel) THEN BEGIN 59 if keyword_set(annuel) THEN BEGIN 56 60 index = where(day EQ 1 AND month EQ 6) 57 if index[0] NE -1 then BEGIN 61 if index[0] NE -1 then BEGIN 58 62 day[index] = 0 59 63 month[index] = 0 -
trunk/SRC/Obsolete/vraidate.pro
r114 r118 6 6 ; @obsolete 7 7 ; 8 ;9 8 ; @param date {in}{required} une date du type yyyymmdd 10 9 ; 11 ; 10 ; @keyword /GRADS if 1 le year le 49 then year = 2000+year 12 11 ; if 50 le year le 99 then year = 1900+year 12 ; 13 ; @keyword _EXTRA used to pass your keywords 13 14 ; 14 15 ; @returns une date vairmer du type yyyymmdd 15 16 ; si year est nulle ou egale a -1 ne change rien 16 17 ; 17 ; @examples 18 ; @examples 18 19 ; IDL> res=vraidate(date) 19 ; IDL> vraidate(980703) 20 ; IDL> vraidate(980703) 20 21 ; donne 19980703 qui est un long 21 22 ; … … 23 24 ; 3/7/98 24 25 ; remove automatic change from year 1 to 1901... Aug 2004 26 ; 27 ; @version $Id$ 28 ; 25 29 ;- 26 30 ;------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.