Ignore:
Timestamp:
03/16/07 10:22:26 (17 years ago)
Author:
pinsard
Message:

corrections of some misspellings in some *.pro

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SRC/Computation/curl.pro

    r168 r226  
    99; @categories 
    1010; Calculation 
    11 ;  
     11; 
    1212; @param UU 
    1313; Matrix representing the zonal coordinates (U point) of a field of vectors 
    1414; A 2D (xy), 3D (xyz or yt) or a structure readable by litchamp and containing 
    1515; a 2D (xy), 3D (xyz or yt) array (4d case is not coded yet). 
    16 ; note that the dimension of the arry must suit the domain dimension. 
    17 ; 
    18 ; @param VV  
     16; note that the dimension of the array must suit the domain dimension. 
     17; 
     18; @param VV 
    1919; Matrix representing the meridional coordinates (V point) of a field of vectors 
    2020; A 2D (xy), 3D (xyz or yt) or a structure readable by litchamp and containing 
    2121; a 2D (xy), 3D (xyz or yt) array (4d case is not coded yet). 
    22 ; note that the dimension of the arry must suit the domain dimension. 
     22; note that the dimension of the array must suit the domain dimension. 
    2323; 
    2424; @keyword DIREC {type=scalar string} 
     
    3030; 
    3131; @uses 
    32 ; cm_4cal, cm_4data, cm_4mmesh  
     32; cm_4cal, cm_4data, cm_4mmesh 
    3333; 
    3434; @restrictions 
    3535; 
    36 ; - Works only for Arakawa C-grid.  
     36; - Works only for Arakawa C-grid. 
    3737; - UU must be on U grid, VV must be on V grid 
    3838; - 4d case is not coded yet 
    3939; - the common variable jpt is used to differ xyz (jpt=1) and xyt (jpt\=1) cases. 
    40 ; - U and V arrays are cut in the same geographic domain. Because of the shift between  
    41 ;   T, U, V and F grids, it is possible that these two arrays do not have the same  
    42 ;   size and refer to different indexes. In this case, arrays are re-cut on  
     40; - U and V arrays are cut in the same geographic domain. Because of the shift between 
     41;   T, U, V and F grids, it is possible that these two arrays do not have the same 
     42;   size and refer to different indexes. In this case, arrays are re-cut on 
    4343;   common indexes. To avoid these re-cuts, use the keyword /memeindice in domdef.pro 
    4444; - When computing the divergence, we update, vargrid, varname, varunits and the 
    4545;   grid position parameters (firstxf, lastxf, nxf, firstyf, lastyf, nyf). 
    46 ; - points that cannot be computed (domain bondaries, coastline) are set to NaN 
     46; - points that cannot be computed (domain boundaries, coastline) are set to NaN 
    4747; 
    4848; @examples 
     
    5050; IDL> plt, curl(dist(jpi,jpj), dist(jpi,jpj)) 
    5151; 
    52 ; @history  
     52; @history 
    5353; Guillaume Roullet (grlod\@ipsl.jussieu.fr) 
    5454; Sebastien Masson (smasson\@lodyc.jussieu.fr) 
    5555; adaptation to work with a reduce domain 
    5656; 21/5/1999: missing values at !values.f_nan 
    57 ;  
     57; 
    5858; @version 
    5959; $Id$ 
    6060; 
    61 ; @todo  
     61; @todo 
    6262; code the 4d case 
    6363;- 
     
    7171@cm_4cal                        ; for jpt 
    7272@cm_4data                       ; for varname, vargrid, vardate, varunit, valmask 
    73 @cm_4mesh  
     73@cm_4mesh 
    7474; 
    7575  tempsun = systime(1)          ; To key_performance 
     
    8484  szu = size(u) 
    8585  szv = size(v) 
    86    
     86 
    8787  if szu[0] NE szv[0] then return, report('U and V input data must have the same number of dimensions') 
    8888 
     
    9696  indiceyv = (lindgen(jpj))[firstyv:firstyv+nyv-1] 
    9797  indicey = inter(indiceyu, indiceyv) 
    98   nx = n_elements(indicex)  
     98  nx = n_elements(indicex) 
    9999  ny = n_elements(indicey) 
    100100  indice2d = lindgen(jpi, jpj) 
     
    105105  varunits = 's-1' 
    106106  if n_elements(valmask) EQ 0 THEN valmask = 1e20 
    107   firstxt = indicex[0] & lastxt = indicex[0]+nx-1 & nxt = nx  
     107  firstxt = indicex[0] & lastxt = indicex[0]+nx-1 & nxt = nx 
    108108  firstyt = indicey[0] & lastyt = indicey[0]+ny-1 & nyt = ny 
    109109;---------------------------------------------------------------------------- 
     
    113113;xyz 
    114114;---------------------------------------------------------------------------- 
    115     szu[0] EQ 3 AND jpt EQ 1:BEGIN  
     115    szu[0] EQ 3 AND jpt EQ 1:BEGIN 
    116116;------------------------------------------------------------ 
    117117; extraction of U and V on the appropriated domain 
     
    119119      case 1 of 
    120120        szu[1] EQ nxu AND szu[2] EQ nyu AND $ 
    121            szv[1] EQ nxv AND szv[2] EQ nyv:BEGIN  
     121           szv[1] EQ nxv AND szv[2] EQ nyv:BEGIN 
    122122          case 1 of 
    123123            nxu NE nx:if indicex[0] EQ firstxu then u = u[0:nx-1, *, *] ELSE u = u[1: nx, *, *] 
     
    129129        END 
    130130        szu[1] EQ jpi AND szu[2] EQ jpj AND $ 
    131            szv[1] EQ jpi AND szv[2] EQ jpj:BEGIN  
     131           szv[1] EQ jpi AND szv[2] EQ jpj:BEGIN 
    132132          u = u[indicex[0]:indicex[0]+nx-1, indicey[0]:indicey[0]+ny-1, *] 
    133133          v = v[indicex[0]:indicex[0]+nx-1, indicey[0]:indicey[0]+ny-1, *] 
     
    142142      landu = where(coefu EQ 0) 
    143143      if landu[0] NE -1 then coefu[temporary(landu)] = !values.f_nan 
    144        
     144 
    145145      coefv = ((e2v[indice2d])[*]#replicate(1., nzt)) $ 
    146146              *(vmask())[indicex[0]:indicex[0]+nx-1, indicey[0]:indicey[0]+ny-1, firstzt:lastzt] 
     
    175175;---------------------------------------------------------------------------- 
    176176;---------------------------------------------------------------------------- 
    177     szu[0] EQ 3 AND jpt GT 1:BEGIN  
     177    szu[0] EQ 3 AND jpt GT 1:BEGIN 
    178178;------------------------------------------------------------ 
    179179; extraction of U and V on the appropriated domain 
     
    181181      case 1 of 
    182182        szu[1] EQ nxu AND szu[2] EQ nyu AND $ 
    183            szv[1] EQ nxv AND szv[2] EQ nyv:BEGIN  
     183           szv[1] EQ nxv AND szv[2] EQ nyv:BEGIN 
    184184          if nxu NE nx then $ 
    185              if indicex[0] EQ firstxu then u = u[0:nx-1, *, *] ELSE u = u[1: nx, *, *]  
     185             if indicex[0] EQ firstxu then u = u[0:nx-1, *, *] ELSE u = u[1: nx, *, *] 
    186186          IF nxv NE nx THEN $ 
    187187             if indicex[0] EQ firstxv then v = v[0:nx-1, *, *] ELSE v = v[1: nx, *, *] 
    188188          IF nyu NE ny THEN $ 
    189              if indicey[0] EQ firstyu then u = u[*, 0:ny-1, *] ELSE u = u[*, 1: ny, *]  
     189             if indicey[0] EQ firstyu then u = u[*, 0:ny-1, *] ELSE u = u[*, 1: ny, *] 
    190190          IF nyv NE ny THEN $ 
    191191             if indicey[0] EQ firstyv then v = v[*, 0:ny-1, *] ELSE v = v[*, 1: ny, *] 
    192192        END 
    193193        szu[1] EQ jpi AND szu[2] EQ jpj AND $ 
    194            szv[1] EQ jpi AND szv[2] EQ jpj:BEGIN  
     194           szv[1] EQ jpi AND szv[2] EQ jpj:BEGIN 
    195195          u = u[indicex[0]:indicex[0]+nx-1, indicey[0]:indicey[0]+ny-1, *] 
    196196          v = v[indicex[0]:indicex[0]+nx-1, indicey[0]:indicey[0]+ny-1, *] 
    197197        END 
    198         ELSE:BEGIN  
     198        ELSE:BEGIN 
    199199          print, 'problemes d''adequation entre la taille du domaine et la taille des matrices necessaires a tracer des vecteurs' 
    200200          return, -1 
     
    240240;---------------------------------------------------------------------------- 
    241241;---------------------------------------------------------------------------- 
    242     szu[0] EQ 4:BEGIN  
     242    szu[0] EQ 4:BEGIN 
    243243      return, report('Case not coded contact saxo team or make a do loop!') 
    244244    END 
     
    248248;---------------------------------------------------------------------------- 
    249249;---------------------------------------------------------------------------- 
    250     szu[0] EQ 2:BEGIN  
     250    szu[0] EQ 2:BEGIN 
    251251;------------------------------------------------------------ 
    252252;------------------------------------------------------------ 
    253253      case 1 of 
    254254        szu[1] EQ nxu AND szu[2] EQ nyu AND $ 
    255            szv[1] EQ nxv AND szv[2] EQ nyv:BEGIN  
     255           szv[1] EQ nxv AND szv[2] EQ nyv:BEGIN 
    256256          if nxu NE nx then $ 
    257              if indicex[0] EQ firstxu then u = u[0:nx-1, *] ELSE u = u[1: nx, *]  
     257             if indicex[0] EQ firstxu then u = u[0:nx-1, *] ELSE u = u[1: nx, *] 
    258258          IF nxv NE nx THEN $ 
    259259             if indicex[0] EQ firstxv then v = v[0:nx-1, *] ELSE v = v[1: nx, *] 
    260260          IF nyu NE ny THEN $ 
    261              if indicey[0] EQ firstyu then u = u[*, 0:ny-1] ELSE u = u[*, 1: ny]  
     261             if indicey[0] EQ firstyu then u = u[*, 0:ny-1] ELSE u = u[*, 1: ny] 
    262262          IF nyv NE ny THEN $ 
    263263             if indicey[0] EQ firstyv then v = v[*, 0:ny-1] ELSE v = v[*, 1: ny] 
    264264        END 
    265265        szu[1] EQ jpi AND szu[2] EQ jpj AND $ 
    266            szv[1] EQ jpi AND szv[2] EQ jpj:BEGIN  
     266           szv[1] EQ jpi AND szv[2] EQ jpj:BEGIN 
    267267          u = u[indice2d] 
    268268          v = v[indice2d] 
     
    305305  ENDCASE 
    306306;------------------------------------------------------------ 
    307   if keyword_set(key_performance) THEN print, 'temps curl', systime(1)-tempsun  
     307  if keyword_set(key_performance) THEN print, 'temps curl', systime(1)-tempsun 
    308308 
    309309  return, psi 
Note: See TracChangeset for help on using the changeset viewer.