Ignore:
Timestamp:
03/19/07 18:15:51 (17 years ago)
Author:
pinsard
Message:

improvements/corrections of some *.pro headers

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SRC/ToBeReviewed/GRILLE/f2v.pro

    r163 r231  
    1 ;------------------------------------------------------------ 
    2 ;------------------------------------------------------------ 
    3 ;------------------------------------------------------------ 
    41;+ 
    52; 
    63; @file_comments 
    7 ; Allows to pass a field referring to the grid F on the grid V,  
    8 ; thanks to the average:  
    9 ; res = 0.5*(res + shift(res, 1, 0))     
     4; Allows to pass a field referring to the grid F on the grid V, 
     5; thanks to the average: 
     6; res = 0.5*(res + shift(res, 1, 0)) 
    107; 
    11 ; @categories  
     8; @categories 
    129; Grid 
    1310; 
     
    2118; common.pro 
    2219; 
    23 ; @restrictions  
    24 ; Force parameters of the zoom on the grid V to be the same as  
     20; @restrictions 
     21; Force parameters of the zoom on the grid V to be the same as 
    2522; those on the grid T. 
    2623; 
    27 ; @restrictions  
    28 ; Points which can not be calculated are put at the value NaN  
     24; @restrictions 
     25; Points which can not be calculated are put at the value NaN 
    2926; consecrated by IDL: !values.f_nan 
    3027; 
     
    3835; 
    3936;- 
    40 ;------------------------------------------------------------ 
    41 ;------------------------------------------------------------ 
    42 ;------------------------------------------------------------ 
     37; 
    4338FUNCTION f2v, temp 
    44 ;--------------------------------------------------------- 
    4539; 
    4640  compile_opt idl2, strictarrsubs 
     
    6862; case on tha array's size and application 
    6963   taille = size(temp) 
    70    CASE taille[0] OF  
     64   CASE taille[0] OF 
    7165      1: res = -1 
    72       2: BEGIN  
     66      2: BEGIN 
    7367         case 1 of 
    7468            taille[1] eq nxf and taille[2] eq nyf: 
     
    7973         endcase 
    8074         mask = (fmask())[firstxf:lastxf, firstyf:lastyf, 0] 
    81          terre = where(mask EQ 0)  
     75         terre = where(mask EQ 0) 
    8276         IF terre[0] NE -1 THEN res[terre] = !values.f_nan 
    8377         res = 0.5*(res + shift(res, 1, 0)) 
    8478         if NOT (keyword_set(key_periodic) AND nxf EQ jpi) then res[0, *] = !values.f_nan 
    8579         mask = (vmask())[firstxf:lastxf, firstyf:lastyf, 0] 
    86          terre = where(mask EQ 0)  
     80         terre = where(mask EQ 0) 
    8781         IF terre[0] NE -1 THEN res[terre] = valmask 
    8882      END 
    89       3: BEGIN  
     83      3: BEGIN 
    9084         case 1 of 
    9185            taille[1] eq nxf and taille[2] eq nyf AND taille[3] EQ nzt: 
     
    10599            mask = reform(mask, nxf, nyf, jpt, /over) 
    106100         ENDIF ELSE mask = (fmask())[firstxf:lastxf, firstyf:lastyf, firstzt:lastzt] 
    107          terre = where(temporary(mask) EQ 0)  
     101         terre = where(temporary(mask) EQ 0) 
    108102         IF terre[0] NE -1 THEN res[temporary(terre)] = !values.f_nan 
    109103         res = 0.5*(res + shift(res, 1, 0, 0)) 
     
    114108            mask = reform(mask, nxf, nyf, jpt, /over) 
    115109         ENDIF ELSE mask = (vmask())[firstxf:lastxf, firstyf:lastyf, firstzt:lastzt] 
    116          terre = where(temporary(mask) EQ 0)  
     110         terre = where(temporary(mask) EQ 0) 
    117111         IF terre[0] NE -1 THEN res[temporary(terre)] = valmask 
    118112      END 
    119       4: BEGIN  
     113      4: BEGIN 
    120114         case 1 of 
    121115            taille[1] eq nxf and taille[2] eq nyf AND taille[3] EQ nzt AND taille[4] EQ jpt: 
     
    130124         mask = temporary(mask[*])#replicate(1, jpt) 
    131125         mask = reform(mask, nxf, nyf, nzt, jpt, /over) 
    132          terre = where(temporary(mask) EQ 0)  
     126         terre = where(temporary(mask) EQ 0) 
    133127         IF terre[0] NE -1 THEN res[temporary(terre)] = !values.f_nan 
    134128         res = 0.5*(res + shift(res, 1, 0, 0, 0)) 
     
    137131         mask = temporary(mask[*])#replicate(1, jpt) 
    138132         mask = reform(mask, nxf, nyf, nzt, jpt, /over) 
    139          terre = where(temporary(mask) EQ 0)  
     133         terre = where(temporary(mask) EQ 0) 
    140134         IF terre[0] NE -1 THEN res[temporary(terre)] = valmask 
    141135      END 
     
    144138   IF NOT keyword_set(key_forgetold) THEN BEGIN 
    145139   @updateold 
    146    ENDIF  
    147     
     140   ENDIF 
     141 
    148142   return, res 
    149143END 
    150  
Note: See TracChangeset for help on using the changeset viewer.