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/v2t.pro

    r163 r231  
    1 ;------------------------------------------------------------ 
    2 ;------------------------------------------------------------ 
    3 ;------------------------------------------------------------ 
    41;+ 
    52; 
    63; @file_comments 
    7 ; Allows to pass a field referring to the grid V on the grid T,  
    8 ; thanks to the average:  
    9 ; res = 0.5*(res + shift(res, 1, 0))     
     4; Allows to pass a field referring to the grid V on the grid T, 
     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; 
     
    3734; $Id$ 
    3835; 
    39 ; @todo seb: nom du input? 
     36; @todo 
     37; seb: nom du input? 
    4038; 
    4139;- 
    42 ;------------------------------------------------------------ 
    43 ;------------------------------------------------------------ 
    44 ;------------------------------------------------------------ 
     40; 
    4541FUNCTION v2t, temp 
    46 ;--------------------------------------------------------- 
    4742; 
    4843  compile_opt idl2, strictarrsubs 
     
    7065; case on tha array's size and application 
    7166   taille = size(temp) 
    72    CASE taille[0] OF  
     67   CASE taille[0] OF 
    7368      1: res = -1 
    74       2: BEGIN  
     69      2: BEGIN 
    7570         case 1 of 
    7671            taille[1] eq nxv and taille[2] eq nyv: 
     
    8176         endcase 
    8277         mask = (vmask())[firstxv:lastxv, firstyv:lastyv, 0] 
    83          terre = where(mask EQ 0)  
     78         terre = where(mask EQ 0) 
    8479         IF terre[0] NE -1 THEN res[terre] = !values.f_nan 
    8580         res = 0.5*(res + shift(res, 0, +1)) 
    8681         res[*, 0] = !values.f_nan 
    8782         mask = tmask[firstxv:lastxv, firstyv:lastyv, 0] 
    88          terre = where(mask EQ 0)  
     83         terre = where(mask EQ 0) 
    8984         IF terre[0] NE -1 THEN res[terre] = valmask 
    9085      END 
    91       3: BEGIN  
     86      3: BEGIN 
    9287         case 1 of 
    9388            taille[1] eq nxv and taille[2] eq nyv AND taille[3] EQ nzt: 
     
    107102            mask = reform(mask, nxv, nyv, jpt, /over) 
    108103         ENDIF ELSE mask = (vmask())[firstxv:lastxv, firstyv:lastyv, firstzt:lastzt] 
    109          terre = where(temporary(mask) EQ 0)  
     104         terre = where(temporary(mask) EQ 0) 
    110105         IF terre[0] NE -1 THEN res[temporary(terre)] = !values.f_nan 
    111106         res = 0.5*(res + shift(res, 0, +1, 0)) 
     
    116111            mask = reform(mask, nxv, nyv, jpt, /over) 
    117112         ENDIF ELSE mask = tmask[firstxv:lastxv, firstyv:lastyv, firstzt:lastzt] 
    118          terre = where(temporary(mask) EQ 0)  
     113         terre = where(temporary(mask) EQ 0) 
    119114         IF terre[0] NE -1 THEN res[temporary(terre)] = valmask 
    120115      END 
    121       4: BEGIN  
     116      4: BEGIN 
    122117         case 1 of 
    123118            taille[1] eq nxv and taille[2] eq nyv AND taille[3] EQ nzt AND taille[4] EQ jpt: 
     
    132127         mask = temporary(mask[*])#replicate(1, jpt) 
    133128         mask = reform(mask, nxv, nyv, nzt, jpt, /over) 
    134          terre = where(temporary(mask) EQ 0)  
     129         terre = where(temporary(mask) EQ 0) 
    135130         IF terre[0] NE -1 THEN res[temporary(terre)] = !values.f_nan 
    136131         res = 0.5*(res + shift(res, 0, +1, 0, 0)) 
     
    139134         mask = temporary(mask[*])#replicate(1, jpt) 
    140135         mask = reform(mask, nxv, nyv, nzt, jpt, /over) 
    141          terre = where(temporary(mask) EQ 0)  
     136         terre = where(temporary(mask) EQ 0) 
    142137         IF terre[0] NE -1 THEN res[temporary(terre)] = valmask 
    143138      END 
     
    146141   IF NOT keyword_set(key_forgetold) THEN BEGIN 
    147142   @updateold 
    148    ENDIF  
     143   ENDIF 
    149144 
    150145    return, res 
    151146END 
    152  
Note: See TracChangeset for help on using the changeset viewer.