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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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 
Note: See TracChangeset for help on using the changeset viewer.