Changeset 118 for trunk/SRC/Obsolete


Ignore:
Timestamp:
06/27/06 17:47:06 (18 years ago)
Author:
pinsard
Message:

add $ in Calendar, Grid, Interpolation, Obsolete and Postscript *.pro files, add svn:keywords Id to all these files, some improvements in header

Location:
trunk/SRC/Obsolete
Files:
19 edited

Legend:

Unmodified
Added
Removed
  • trunk/SRC/Obsolete/common.pro

    r97 r118  
     1;+ 
    12; @obsolete 
     3; 
     4; @version $Id$ 
     5; 
     6;- 
     7 
    28@all_cm 
  • trunk/SRC/Obsolete/cp.pro

    r114 r118  
    44;+ 
    55; @file_comments copy files 
     6; 
    67; @obsolete 
     8; 
    79; @file_comments file_copy should be used instead... 
    810; 
     11; @param filenamein {in}{required} 
     12; @param filenameout {in}{required} 
     13; 
     14; @keyword _EXTRA used to pass your keywords 
     15; 
    916; @history 
    10 ; June 2005: Sebastien Masson => obsolete routine  
     17; June 2005: Sebastien Masson => obsolete routine 
     18; 
     19; @version $Id$ 
    1120; 
    1221;- 
  • trunk/SRC/Obsolete/ficdate.pro

    r114 r118  
    2121; @history Jerome Vialard  (jv\@lodyc.jussieu.fr) 
    2222;1/7/98 
     23; 
     24; @version $Id$ 
     25; 
    2326;- 
    2427;------------------------------------------------------------ 
     
    4649   rien = juldate(vdate) 
    4750; constitution de la date yymmdd 
    48    case 1 of  
     51   case 1 of 
    4952     year lt 10:                s_year ='0'+string(format='(i1)',year) 
    5053     year ge 10 and year lt 100:s_year=     string(format='(i2)',year) 
    51      year ge 100:BEGIN  
     54     year ge 100:BEGIN 
    5255        year = year-1900 
    5356        if year LT 10 then s_year ='0'+string(format='(i1)',year) $ 
  • trunk/SRC/Obsolete/fictype.pro

    r114 r118  
    55; @file_comments  gives fictype ('DA' 'MO' 'AN' 'SE' 'EX') corresponding 
    66;           to the given date in vairmer format ( yymmdd ) 
     7; 
     8; @param vdate {in}{required} 
     9; @param dim unused 
     10; 
    711; @uses common.pro 
    812; 
     
    1216; @history Jerome Vialard  (jv\@lodyc.jussieu.fr) 
    1317; 2/7/98 
     18; 
     19; @version $Id$ 
     20; 
    1421;- 
    1522;------------------------------------------------------------ 
  • trunk/SRC/Obsolete/imprime.pro

    r114 r118  
    33;------------------------------------------------------------ 
    44;+ 
    5 ; @file_comments  use printps instead... 
     5; @file_comments use printps instead... 
    66; @obsolete 
     7; 
     8; @param filename {in}{required} 
     9; 
     10; @keyword TRANS 
     11; 
     12; @keyword NB 
     13; 
     14; @uses printps 
    715; 
    816; @history 
    917; June 2005: Sebastien Masson, english version 
     18; 
     19; @version $Id$ 
     20; 
    1021;- 
    1122;------------------------------------------------------------ 
     
    1930; 
    2031  thisOS = strupcase(strmid(!version.os_family, 0, 3)) 
    21   CASE thisOS OF  
     32  CASE thisOS OF 
    2233    'MAC':return 
    2334    'WIN':return 
     
    3041    0:printps 
    3142    1:printps, filename 
    32     ELSE: BEGIN  
     43    ELSE: BEGIN 
    3344      ras = report('imprime accept only one element: psfilename') 
    3445      return 
    35     END  
    36   ENDCASE        
     46    END 
     47  ENDCASE 
    3748;------------------------------------------------------------ 
    3849  return 
  • trunk/SRC/Obsolete/jourdsmois.pro

    r114 r118  
    66; @obsolete 
    77; 
     8; @param mois {in}{optional} 
     9; @param annee {in}{optional} 
     10; 
    811; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 
    912; June 2005: Sebastien Masson, english version 
     13; 
     14; @version $Id$ 
     15; 
    1016;- 
    1117;------------------------------------------------------------ 
  • trunk/SRC/Obsolete/juldate.pro

    r114 r118  
    66; @obsolete 
    77; 
    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 
    1010;           sets month, day and year to the corresp values 
    1111; 
     
    1313; 
    1414; @param vvdate {in} date de la forme yymmdd ou yyyymmdd 
     15; 
     16; @keyword _EXTRA used to pass your keywords 
    1517; 
    1618; @uses common.pro vraidate 
     
    2022;     jour julien de l'annee -1. 
    2123;   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. 
    2325;   Attention, les variables globales year, month day sont attribuees 
    2426; 
    25 ; @examples  
     27; @examples 
    2628; IDL> date = juldate(930124) 
    2729; 
    2830; @history Jerome Vialard  (jv\@lodyc.jussieu.fr) 
    2931; 2/7/98 
     32; 
     33; @version $Id$ 
     34; 
    3035;- 
    3136;------------------------------------------------------------ 
     
    4752   mm = month & dd = day & yy = year 
    4853;------------------------------------------------------------ 
    49    ndate = n_elements(vdate)  
     54   ndate = n_elements(vdate) 
    5055   if total(mm EQ 0) EQ ndate then mm[*] = 6 
    5156   if total(dd EQ 0) EQ ndate then dd[*] = 15 
  • trunk/SRC/Obsolete/lec.pro

    r114 r118  
    2020; 
    2121; @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 
    2323; minuscule commencant par vo ou so). Dans cette methode on saute directement 
    2424; d'en-tete en en-tete jusqu'a trouver le bon fichier. 
     
    3131; 
    3232; 
    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 
    3434;             l'anomalie ('EX','AN','SE','MO','') 
    3535; 
     
    3737; ds ce cas en input on met seulement 'vo' ou 'so' la fonction retourne le 
    3838; nombre de fichiers lus. 
     39; 
     40; @keyword BOITE 
    3941; 
    4042; @keyword /EXPANOM {in} si on calcule l'anom par rapport a une exper 
     
    6062; celle du sous domaine definit par domdef (nx,ny,nz). avec /TOUT le 
    6163; 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: 
    6466;        http://www.ipsl.jussieu.fr/~smlod/sousdomaine.html 
    6567; Retourne -1 en cas d'erreur. 
    6668; 
    67 ; @uses common.pro isnumber.pro fivardate.pro   
     69; @uses common.pro isnumber.pro fivardate.pro 
    6870; 
    6971; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 
     
    7678;                       Sebastien Masson (decoupe pour les sous domaines...) 
    7779;                       2/99 
     80; 
     81; @version $Id$ 
     82; 
    7883;;- 
    7984;-------------------------------------------------------------- 
     
    141146; definition de la taille du fichier a aller chercher: jpidta,jpjdta,jpkdta... 
    142147;--------------------------------------------------------------------- 
    143    if n_elements(jpidta) EQ 0 THEN BEGIN  
     148   if n_elements(jpidta) EQ 0 THEN BEGIN 
    144149      if n_elements(ixmindta) EQ 0 OR n_elements(ixmaxdta) EQ 0 then $ 
    145150       jpidta = jpiglo else jpidta = ixmaxdta-ixmindta+1 
    146151   endif 
    147    if n_elements(jpjdta) EQ 0 THEN BEGIN  
     152   if n_elements(jpjdta) EQ 0 THEN BEGIN 
    148153      if n_elements(iymindta) EQ 0 OR n_elements(iymaxdta) EQ 0 then $ 
    149154       jpjdta = jpjglo else jpjdta = iymaxdta-iymindta+1 
    150155   endif 
    151    if n_elements(jpkdta) EQ 0 THEN BEGIN  
     156   if n_elements(jpkdta) EQ 0 THEN BEGIN 
    152157      if n_elements(izmindta) EQ 0 OR n_elements(izmaxdta) EQ 0 then $ 
    153158       jpkdta = jpkglo else jpkdta = izmaxdta-izmindta+1 
     
    267272; on shift z si key_shift est defininit 
    268273;--------------------------------------------------------------------- 
    269    if n_elements(key_shift) NE 0 THEN BEGIN  
     274   if n_elements(key_shift) NE 0 THEN BEGIN 
    270275      if dim EQ 'SO' then z = shift(z,key_shift, 0) $ 
    271276      ELSE z = shift(z,key_shift, 0, 0) 
     
    275280;  du zoom: nx,ny nz 
    276281;--------------------------------------------------------------------- 
    277    if NOT keyword_set(tout) then BEGIN  
     282   if NOT keyword_set(tout) then BEGIN 
    278283;------------------------------------------------------------- 
    279284; changement de domaine 
    280285;------------------------------------------------------------- 
    281       if keyword_set(boite) then BEGIN  
     286      if keyword_set(boite) then BEGIN 
    282287         Case 1 Of 
    283288            N_Elements(Boite) Eq 1:bte=[lon1, lon2, lat1, lat2, 0.,boite[0]] 
     
    290295         oldboite = [lon1, lon2, lat1, lat2, prof1, prof2] 
    291296         domdef, bte,GRILLE=vargrid 
    292       ENDIF  
     297      ENDIF 
    293298;------------------------------------------------------------- 
    294299      grille,mask,glam,gphi,gdep,nx,ny,nz,premierx,premiery,premierz,dernierx,derniery,dernierz 
     
    296301      if dim EQ 'SO' then z = z[premierx:dernierx, premiery:derniery]  $ 
    297302      ELSE z = z[premierx:dernierx, premiery:derniery, premierz:dernierz] 
    298    ENDIF ELSE BEGIN  
     303   ENDIF ELSE BEGIN 
    299304      case vargrid OF           ; on recupere le mask en entier ds le cas ou /TOUT 
    300305         'U':mask = umask()     ; n''est pas active et on le choisit en fontion 
    301306         'T':mask = tmask       ;  de la valeur de vargrid 
    302          'W':mask = tmask        
     307         'W':mask = tmask  
    303308         'V':mask = vmask() 
    304309         'F':mask = fmask() 
    305310      ENDCASE 
    306    ENDELSE  
     311   ENDELSE 
    307312;--------------------------------------------------------------------- 
    308313; calcul d'une anomalie si le keyword anom est active 
     
    311316      case anom of 
    312317         'EX' : adate = 0 
    313          'AN' : adate = floor(date/10000)*10000  
     318         'AN' : adate = floor(date/10000)*10000 
    314319         'SE' : adate = floor(date - floor(date/10000)*10000)/100 * 100 
    315320         'MO' : adate = floor(date/100)*100 
     
    327332;--------------------------------------------------------------------- 
    328333   IF n_elements(valmask) EQ 0 THEN valmask = 1e20 
    329    if dim EQ 'SO' then BEGIN  
     334   if dim EQ 'SO' then BEGIN 
    330335      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 
    333338      terre = where(mask[*,*,0] EQ 0) 
    334339      if terre[0] NE -1 then z[where(mask EQ 0)] = valmask 
    335    ENDELSE  
     340   ENDELSE 
    336341;--------------------------------------------------------------------- 
    337342   free_lun,numlec 
  • trunk/SRC/Obsolete/lect.pro

    r114 r118  
    33;------------------------------------------------------------ 
    44;+ 
    5 ;    
     5; 
    66; @file_comments 
    77;       lit les fichiers Vairmer de date1 a date2 et en sort un tableau 
     
    1919; IDL> resultat=lec('nom_Vairmer',date1,date2,'nomexp','direc',BOITE=boite) 
    2020; 
    21 ; @param nomchamp {in}{required} nom de champ Vairmer (chaine de 8 caracteres commencant 
    22 ;                   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 
     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 
    2525;             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 
    3837; (par defaut la meme que nomexp) 
    3938; 
     
    4241; @history Jerome Vialard (jv\@lodyc.jussieu.fr) 
    4342; 2/7/98 
     43; 
     44; @version $Id$ 
     45; 
    4446;- 
    4547;-------------------------------------------------------------- 
    4648;-------------------------------------------------------------- 
    4749;-------------------------------------------------------------- 
    48 ;                  1      2      3    4      5        
     50;                  1      2      3    4      5 
    4951function lect, nomchamp,date1,date2,nomexp,direc,BOITE=boite, $ 
    5052               ANOM=anom,EXPANOM=expanom 
     
    6870   if fictype(date1) ne fictype(date2) then $ 
    6971    return, report('Les deux dates doivent correspondre au meme type de fic vairmer') 
    70    fictyp = fictype(date1)  
     72   fictyp = fictype(date1) 
    7173;------------------------------------------------------------- 
    7274; creation du nom du fichier 
     
    7779   if (keyword_set(anom)) then ficname=ficname +'.'+anom 
    7880   if (keyword_set(expanom)) then ficname=ficname +'-'+expanom 
    79    case n_elements(boite) of  
     81   case n_elements(boite) of 
    8082      4 : box = strcompress(string(format='(i4,"_",i4,"_",i4,"_",i4)',boite) $ 
    8183                            ,/remove_all) 
     
    9092; 
    9193; 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) 
    9395; dimttab[0], dimttab[1], ... (valeur des dim ) 
    9496; time (axe des tps), ttab (tableau a lire) 
     
    101103      dimtableau = long(1) 
    102104      readu, numlec, jpt,dimtableau 
    103       case dimtableau of  
     105      case dimtableau of 
    104106         1 : begin 
    105107            n1 = long(1) 
     
    145147; changement de domaine 
    146148;------------------------------------------------------------- 
    147    if keyword_set(boite) then BEGIN  
     149   if keyword_set(boite) then BEGIN 
    148150      Case 1 Of 
    149151         N_Elements(Boite) Eq 1:bte=[lon1, lon2, lat1, lat2, 0.,boite[0]] 
     
    156158      oldboite = [lon1, lon2, lat1, lat2, prof1, prof2] 
    157159      domdef, bte 
    158    ENDIF  
    159 ;------------------------------------------------------------- 
    160 ;   Boucle de lecture des fichiers  
     160   ENDIF 
     161;------------------------------------------------------------- 
     162;   Boucle de lecture des fichiers 
    161163;------------------------------------------------------------- 
    162164   case fictyp of 
     
    190192; Moyenne du champs tab 
    191193;-------------------------------------------------------- 
    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 
    194196            if string(format='(a2)',nomchamp) EQ 'SO' then tab = reform(tab, nx, ny, /over) $ 
    195197            ELSE tab = reform(tab, nx, ny, nz, /over) 
    196          ENDIF  
     198         ENDIF 
    197199         tab = moyenne(tab,direc) 
    198200      endif 
    199201      if (jpt eq 1) then begin 
    200202         ttab = tab 
    201       endif else BEGIN  
     203      endif else BEGIN 
    202204         ttab = colle(ttab, tab, (size(tab))[0]+1 ) 
    203205      endelse 
     
    217219      endcase 
    218220;------------------------------------------------------------ 
    219 ;   Fin de boucle de lecture des fichiers  
     221;   Fin de boucle de lecture des fichiers 
    220222;------------------------------------------------------------ 
    221223      vdat=long(10000)*year+long(100)*month+day 
     
    229231   taille = size(ttab) 
    230232   writeu, numlec, long(jpt),long(taille[0]) 
    231    case taille[0] of  
     233   case taille[0] of 
    232234      1 : writeu, numlec,long(taille[1]) 
    233235      2 : writeu, numlec,long(taille[1]),long(taille[2]) 
     
    258260   close, /all 
    259261; 
    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 
    261263; 
    262264   return, ttab 
  • trunk/SRC/Obsolete/meshlec.pro

    r114 r118  
    44;+ 
    55; 
    6 ; @file_comments lecture du mask de des sorties d'OPA. les sources se trouvent ds les  
     6; @file_comments lecture du mask des sorties d'OPA. les sources se trouvent ds les 
    77; repertoires sur maia du type: 
    88;  /nom_exp/RESTARTS 
    99; 
    1010; @obsolete 
     11; 
    1112; @examples 
    1213; IDL> meshmask[,' nomfich'] 
    1314; 
    14 ; @param  nomfich {in}{required} string, c''est le nom du fichier a lire. Par defaut, c''est meshmask.  
    15 ; 
    16 ; 
    17 ;  @keyword GLAMBOUNDARY {in} un vecteur de 2 elements specifaint le min et le 
    18 ;    max qui doivent etre imposes en longitude (obligatoire si le 
    19 ;    tableau depasse 360 degres) 
     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) 
    2021; 
    2122; @keyword /pasblabla {in} pour suprimer les blablas 
     
    2425;        precision. ce Mot clef est maintenant active 
    2526;        automatiquement. 
     27; 
     28; @keyword GETDIMENSIONS 
    2629; 
    2730; @uses common.pro 
     
    3538; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 
    3639;     Marina Levy : lecture en double precision (cas calcul sur shine) 
     40; 
     41; @version $Id$ 
     42; 
    3743;- 
    3844;------------------------------------------------------------ 
     
    109115;       sizefile8:double = 1 
    110116;       sizefile4:double = 0 
    111 ;       ELSE:BEGIN  
     117;       ELSE:BEGIN 
    112118;          nothing = report('The OPA Mesh file as not the good size!') 
    113119;          free_lun,numlec 
    114120;          close, numlec 
    115 ;          return  
     121;          return 
    116122;       END 
    117123;    endcase 
     
    150156;------------------------------------------------------- 
    151157 
    152    IF NOT keyword_set(double) THEN BEGIN  
     158   IF NOT keyword_set(double) THEN BEGIN 
    153159      z3d    = fltarr(jpiglo, jpjglo, jpkglo) 
    154160      z2d    = fltarr(jpiglo, jpjglo) 
    155161      z1d    = fltarr(jpkglo) 
    156    ENDIF ELSE BEGIN  
     162   ENDIF ELSE BEGIN 
    157163      z3d    = dblarr(jpiglo, jpjglo, jpkglo) 
    158164      z2d    = dblarr(jpiglo, jpjglo) 
    159165      z1d    = dblarr(jpkglo) 
    160    ENDELSE  
     166   ENDELSE 
    161167 
    162168   if not keyword_set(pasblabla) then print, ' ' 
     
    175181 
    176182   if not keyword_set(pasblabla) then print, ' ' 
    177    readu, numlec, tab, z2d   
     183   readu, numlec, tab, z2d 
    178184   GPHIT=float(z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh]) 
    179185   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 
    181187   GPHIU=float(z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh]) 
    182188   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 
    184190   GPHIV=float(z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh]) 
    185191   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 
    187193   GPHIF=float(z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh]) 
    188194   if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur GPHIF[25,31]: ',z2d[25,31] 
     
    217223 
    218224   if not keyword_set(pasblabla) then print, ' ' 
    219    readu, numlec, tab, z3d      
     225   readu, numlec, tab, z3d 
    220226   TMASK=byte(z3d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh,izminmesh:izmaxmesh]) 
    221227   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 
    223229   UMASKred=byte(z3d[ixmaxmesh,iyminmesh:iymaxmesh,izminmesh:izmaxmesh]) 
    224230   umaskred = reform(umaskred) 
    225231   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 
    227233   VMASKred=byte(z3d[ixminmesh:ixmaxmesh,iymaxmesh,izminmesh:izmaxmesh]) 
    228234   vmaskred = reform(vmaskred) 
    229235   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 
    231237   fmaskredy=byte(z3d[ixmaxmesh,iyminmesh:iymaxmesh,izminmesh:izmaxmesh]) 
    232238   coast = where(fmaskredy NE 0 and fmaskredy NE 1) 
     
    240246; 
    241247   if not keyword_set(pasblabla) then print, ' ' 
    242    readu, numlec, tab, z2d    
     248   readu, numlec, tab, z2d 
    243249;FF=z2d[ixminmesh:ixmaxmesh,iyminmesh:iymaxmesh] 
    244250   if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur FF[25,5]: ',z2d[25,5] 
     
    246252   GDEPT=float(z1d[izminmesh:izmaxmesh]) 
    247253   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 
    249255   GDEPW=float(z1d[izminmesh:izmaxmesh]) 
    250256   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 
    252258   E3T=float(z1d[izminmesh:izmaxmesh]) 
    253259   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 
    255261   E3W=float(z1d[izminmesh:izmaxmesh]) 
    256262   if not keyword_set(pasblabla) then print, 'tableau: ',tab,' exemple de valeur E3W[3]: ',E3W[3] 
     
    260266;------------------------------------------------------- 
    261267;------------------------------------------------------- 
    262 ; bornes de glam qui ne doivent pas depasser 360 degres...  
     268; bornes de glam qui ne doivent pas depasser 360 degres.. 
    263269;------------------------------------------------------- 
    264270;    minglam = min(glamt, max = maxglam) 
     
    329335   if not keyword_set(pasblabla) then print,'lecture '+nomfich+' finie' 
    330336   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 
    332338 
    333339   return 
  • trunk/SRC/Obsolete/ncdf_meshlec.pro

    r114 r118  
    66; @obsolete 
    77; 
     8; @param filename 
     9; 
     10; @keyword _EXTRA used to pass your keywords 
     11; 
    812; @history 
    913; Aug. 2005, Sebastien Masson: switch to ncdf_meshread 
     14; 
     15; @version $Id$ 
     16; 
    1017;- 
    1118;------------------------------------------------------------ 
  • trunk/SRC/Obsolete/nlec.pro

    r114 r118  
    1515; 
    1616; 
    17 ;       @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 
     17; @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 
    2020;             effectuer les moyennes (si rien n'est donne on n'effectue pas de moyenne) 
    2121; 
    2222; @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 
    2424; 
    2525; @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 
    2834; 
    2935; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 
    3036;   14/8/98 
     37; 
     38; @version $Id$ 
     39; 
    3140;- 
    3241;------------------------------------------------------------ 
    3342;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 
    3645;REF.07.790101.grid.W.nc 
    3746;------------------------------------------------------------ 
     
    5160;------------------------------------------------------------ 
    5261   case n_params() of 
    53       1:BEGIN  
     62      1:BEGIN 
    5463         year=year+1900*(year ne 0 and year ne -1 and year lt 100) 
    5564         date = day+100*month+10000*year 
    5665      end 
    57       2:BEGIN  
     66      2:BEGIN 
    5867         if size(debut, /tname) EQ 'STRING' then begin 
    5968            prefix = strupcase(debut) 
     
    150159         datejul = 5*numsortie-1+julday(1, 1, year) 
    151160         if n_elements(date2) ne 0 then date2jul = 5*numsortie2-1+julday(1, 1, year2) $ 
    152          ELSE date2jul = datejul  
     161         ELSE date2jul = datejul 
    153162         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 
    155164            dejaserie = 1 
    156165            GOTO,  passerie 
    157          ENDIF  
     166         ENDIF 
    158167      ENDIF ELSE BEGIN 
    159168         passerie: 
     
    163172            dejapasserie = 1 
    164173            GOTO,  serie 
    165          ENDIF  
     174         ENDIF 
    166175      endelse 
    167176   ENDIF ELSE BEGIN 
     
    233242; moyenne eventuelle 
    234243;------------------------------------------------------------ 
    235    IF keyword_set(direc) THEN BEGIN  
     244   IF keyword_set(direc) THEN BEGIN 
    236245      IF jpt EQ 1 THEN res = moyenne(res,direc,BOITE=boite, nan = notan) $ 
    237246      ELSE res = grossemoyenne(res,direc,BOITE=boite, nan = notan) 
    238    ENDIF  
     247   ENDIF 
    239248;------------------------------------------------------------ 
    240249; mise en placer des parametres pour le trace... 
     
    248257   if keyword_set(structure) then res = {tab:res, grille:vargrid, unite:varunit, experience:varexp, nom:varname} 
    249258;------------------------------------------------------------ 
    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 
    251260   time = time-2l 
    252261   return,res 
     
    254263;------------------------------------------------------------ 
    255264end 
    256  
    257  
    258  
    259  
    260  
    261  
    262  
    263  
    264  
    265  
    266  
    267  
    268  
    269  
  • trunk/SRC/Obsolete/nlec5j.pro

    r114 r118  
    66; contiennent les sorties a 5j regroupees par type de grille par 
    77; paquets de 6 mois. 
    8 ; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS  
     8; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS 
    99; @obsolete 
    1010; 
     
    1616; @param nom {in}{required} nom du tableau contenant le champ que l'on veut (le 
    1717; 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 
    1921; fichier (compte a partir de 1) a partir de year 
     22; 
    2023; @param nbretps {in}{required} nombre de pas de temps a extraire 
    2124; 
    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) 
    2426; 
    2527; @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 
    2729; 
    2830; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 
    2931; 
    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; 
    3236; @uses common.pro 
    3337; 
     
    3539; 
    3640; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 
     41; 
     42; 
     43; @version $Id$ 
    3744; 
    3845;- 
     
    5158; gestion du nom du fichier 
    5259;------------------------------------------------------------ 
    53    if numsortie le 36 then BEGIN  
     60   if numsortie le 36 then BEGIN 
    5461      mmdd='0101' 
    5562      numsort=numsortie 
     
    8491          ras = report('LES FICHIERS: '+iodir+prefix+'.'+numfich+'.'+s_date+'.grid... n''existe pas.') 
    8592         return,  -1 
    86       ENDIF  
     93      ENDIF 
    8794      for i=0,n_elements(quelsfichiers)-1 do begin 
    8895         IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+quelsfichiers[i]+' > /dev/null' 
     
    94101               vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille 
    95102               goto, grilletrouvee 
    96             ENDIF  
     103            ENDIF 
    97104         endfor 
    98105         ncdf_close,cdfid 
     
    132139      premiery = 0 
    133140      premierz = 0 
    134    endif else BEGIN  
     141   endif else BEGIN 
    135142;------------------------------------------------------------ 
    136143; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 
    137144; + on recupere la dim du no9uveau domaine 
    138145;------------------------------------------------------------ 
    139       if keyword_set(boite) then BEGIN  
     146      if keyword_set(boite) then BEGIN 
    140147         Case  N_Elements(Boite) Of 
    141148            1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex 
     
    148155      ENDIF 
    149156      grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery, premierz, dernierx, derniery, dernierz 
    150    ENDELSE  
     157   ENDELSE 
    151158;------------------------------------------------------------ 
    152159; determination du nombre de pas de tps a extraire ds ce fichier 
     
    169176       res = [[[res]] $ 
    170177              ,[[nlec5j(nom,1+36*(mmdd eq '0101'),nbretps-nt,tout=tout,GRILLE=vargrid,BOITE=boite)]]] $ 
    171       ELSE BEGIN  
     178      ELSE BEGIN 
    172179       res=[res[*] $ 
    173180            ,(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 
    176183      ncdf_varget,cdfid,'time_counter',temps, offset = [numsort-1],count=[nt] 
    177184      time = [long(temps)+julday(1, 5, 1979), time] 
  • trunk/SRC/Obsolete/nlecan.pro

    r114 r118  
    44;+ 
    55; @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 
    77; 
    88; @obsolete 
    99; @categories lecture de NETCDF 
    1010; 
    11 ; @examples  
     11; @examples 
    1212; IDL> res=nlecan(nom, numsortie79, nbretps) 
    1313; 
     
    2121; 
    2222; @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 
    2424; 
    2525; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 
    2626; 
     27; @keyword _EXTRA used to pass your keywords 
     28; 
    2729; @returns tableau 2d (qd la serie ne fait que 1 pas de temps) ou 
    2830; 3d (valable ds 1 premier tps que pour les tableaux 2d) 
    29 ;  
     31; 
    3032; @uses common.pro 
    3133; 
     
    3335; 
    3436; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 
     37; 
     38; 
     39; @version $Id$ 
    3540; 
    3641;- 
     
    5560      vargrid=grille 
    5661      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' 
    5863      cdfid=ncdf_open(iodir+nomfich) 
    5964      varcontient=ncdf_varinq(cdfid,nom) 
     
    7277               vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille 
    7378               goto, grilletrouvee 
    74             ENDIF  
     79            ENDIF 
    7580         endfor 
    7681         ncdf_close,cdfid 
     
    121126      premiery = 0 
    122127      premierz = 0 
    123    endif else BEGIN  
     128   endif else BEGIN 
    124129;------------------------------------------------------------ 
    125130; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 
    126131; + on recupere la dim du no9uveau domaine 
    127132;------------------------------------------------------------ 
    128       if keyword_set(boite) then BEGIN  
     133      if keyword_set(boite) then BEGIN 
    129134         Case N_Elements(Boite) Of 
    130135            1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex 
     
    137142      ENDIF 
    138143      grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery, premierz 
    139    ENDELSE  
     144   ENDELSE 
    140145;------------------------------------------------------------ 
    141146   if varcontient.ndims eq 3 then $ 
  • trunk/SRC/Obsolete/nlecmois.pro

    r114 r118  
    1212; 
    1313; @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 
    1515; 
    1616; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 
     17; @keyword _EXTRA used to pass your keywords 
     18 
    1719; @uses common.pro 
    1820; 
    1921; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 
     22; 
     23; @version $Id$ 
    2024; 
    2125;- 
     
    5963               vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille 
    6064               goto, grilletrouvee 
    61             ENDIF  
     65            ENDIF 
    6266         endfor 
    6367         ncdf_close,cdfid 
     
    8589      premiery=0 
    8690      premierz=0 
    87    endif else BEGIN  
     91   endif else BEGIN 
    8892;------------------------------------------------------------ 
    8993; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 
    9094; + on recupere la dim du no9uveau domaine 
    9195;------------------------------------------------------------ 
    92       if keyword_set(boite) then BEGIN  
     96      if keyword_set(boite) then BEGIN 
    9397         Case N_Elements(Boite) Of 
    9498            1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex 
     
    101105      ENDIF 
    102106      grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery, premierz 
    103    ENDELSE  
     107   ENDELSE 
    104108;------------------------------------------------------------ 
    105109; determination du nombre de pas de tps a extraire ds ce fichier 
     
    121125       res = [[[res]] $ 
    122126              ,[[nlecmois(nom,numsortie79+nt,nbretps-nt,tout=tout,GRILLE=vargrid,BOITE=boite)]]] $ 
    123       else BEGIN  
     127      else BEGIN 
    124128         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 
    129133      ncdf_varget,cdfid,'time_counter',temps, offset = [numsortie-1],count=[nt] 
    130134      time = [long(temps)+julday(12, 31, 1978), time] 
  • trunk/SRC/Obsolete/nlecsaison.pro

    r114 r118  
    44;+ 
    55; @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 
    77; @obsolete 
    88; 
     
    2121; 
    2222; @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 
    2424; 
    2525; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 
    2626; 
     27; @keyword _EXTRA used to pass your keywords 
     28 
    2729; @returns tableau 2d (qd la serie ne fait que 1 pas de temps) ou 
    2830; 3d (valable ds 1 premier tps que pour les tableaux 2d) 
    29 ;  
     31; 
    3032; @uses common.pro 
    3133; 
     
    3335; 
    3436; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 
     37; 
     38; 
     39; @version $Id$ 
    3540; 
    3641;- 
     
    7277               vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille 
    7378               goto, grilletrouvee 
    74             ENDIF  
     79            ENDIF 
    7580         endfor 
    7681         ncdf_close,cdfid 
     
    118123      premiery = 0 
    119124      premierz = 0 
    120    endif else BEGIN  
     125   endif else BEGIN 
    121126;------------------------------------------------------------ 
    122127; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 
    123128; + on recupere la dim du no9uveau domaine 
    124129;------------------------------------------------------------ 
    125       if keyword_set(boite) then BEGIN  
     130      if keyword_set(boite) then BEGIN 
    126131         CASE N_Elements(Boite) Of 
    127132            1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex 
     
    134139      ENDIF 
    135140      grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery,premierz 
    136    ENDELSE  
     141   ENDELSE 
    137142;------------------------------------------------------------ 
    138143   if varcontient.ndims eq 3 then $ 
  • trunk/SRC/Obsolete/nlecserie.pro

    r114 r118  
    44;+ 
    55; @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 
    77; @obsolete 
    88; 
    99; @categories lecture de NETCDF 
    1010; 
    11 ; @examples  
     11; @examples 
    1212; IDL> res=nlecserie(nom, date1, date2) 
    1313; 
     
    2020; 
    2121; @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 
    2323; 
    2424; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef 
    2525; 
     26; @keyword _EXTRA used to pass your keywords 
     27; 
    2628; @returns tableau 2d (qd la serie ne fait que 1 pas de temps) ou 
    2729; 3d (valable ds 1 premier tps que pour les tableaux 2d) 
    28 ;  
     30; 
    2931; @uses common.pro 
    3032; 
     
    3234; 
    3335; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) 
     36; 
     37; @version $Id$ 
    3438; 
    3539;- 
     
    5660;      print, 'LES FICHIERS: '+iodir+prefix+'.'+nom+'.grid... n''existe pas.' 
    5761      return,  -1 
    58    ENDIF ELSE BEGIN  
     62   ENDIF ELSE BEGIN 
    5963      vargrid=strmid(quelsfichiers[0],strpos(quelsfichiers[0],'grid')+5,1) ;nom de grille 
    6064      IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+quelsfichiers[0]+' > /dev/null' 
    6165      cdfid=ncdf_open(quelsfichiers[0]) 
    62    ENDELSE  
     66   ENDELSE 
    6367;------------------------------------------------------------ 
    6468; lecture de certains attributs et de l'axe des temps 
     
    116120      premiery = 0 
    117121      premierz = 0 
    118    endif else BEGIN  
     122   endif else BEGIN 
    119123;------------------------------------------------------------ 
    120124; redefinition eventuelle du domaine ajuste a boite (a 6 elements) 
    121125; + on recupere la dim du no9uveau domaine 
    122126;------------------------------------------------------------ 
    123       if keyword_set(boite) then BEGIN  
     127      if keyword_set(boite) then BEGIN 
    124128         Case N_Elements(Boite) Of 
    125129            1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex 
  • trunk/SRC/Obsolete/vairdate.pro

    r114 r118  
    99; @file_comments you better use jul2date 
    1010; 
    11 ; @examples  
     11; @examples 
    1212; IDL> vdate = vairdate(1755087) 
    1313; 
    1414; @param jdate {in}{required} date en jours juliens 
    15 ; 
    1615; 
    1716; @keyword MENSUEL a activer si on veut pour  que les dates dont le jour est 
     
    2322;    0. par ex: 
    2423;    19990601   ->   19990000 
     24; 
     25; @keyword _EXTRA used to pass your keywords 
    2526; 
    2627; @returns vdate date vairmer 
     
    3637;                       d''utiliser vairdate avec des scalaire ou des 
    3738;                       tableaux 
     39; 
     40; @version $Id$ 
     41; 
    3842;- 
    3943;------------------------------------------------------------ 
     
    4953   index = where(year eq -1) 
    5054   if index[0] NE -1 then year[index] = 0 
    51    if keyword_set(mensuel) THEN BEGIN  
     55   if keyword_set(mensuel) THEN BEGIN 
    5256      index = where(day EQ 15) 
    5357      if index[0] NE -1 then day[index] = 0 
    5458   endif 
    55    if keyword_set(annuel) THEN BEGIN  
     59   if keyword_set(annuel) THEN BEGIN 
    5660      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 
    5862         day[index] = 0 
    5963         month[index] = 0 
  • trunk/SRC/Obsolete/vraidate.pro

    r114 r118  
    66; @obsolete 
    77; 
    8 ; 
    98; @param date {in}{required} une date du type yyyymmdd 
    109; 
    11 ;  @keyword /GRADS if  1 le year le 49 then year = 2000+year 
     10; @keyword /GRADS if  1 le year le 49 then year = 2000+year 
    1211;               if 50 le year le 99 then year = 1900+year 
     12; 
     13; @keyword _EXTRA used to pass your keywords 
    1314; 
    1415; @returns une date vairmer du type yyyymmdd 
    1516; si year est nulle ou egale a -1 ne change rien 
    1617; 
    17 ; @examples  
     18; @examples 
    1819; IDL> res=vraidate(date) 
    19 ; IDL> vraidate(980703)  
     20; IDL> vraidate(980703) 
    2021; donne 19980703 qui est un long 
    2122; 
     
    2324;   3/7/98 
    2425; remove automatic change from year 1 to 1901... Aug 2004 
     26; 
     27; @version $Id$ 
     28; 
    2529;- 
    2630;------------------------------------------------------------ 
Note: See TracChangeset for help on using the changeset viewer.