Changeset 48 for trunk/extrap.pro


Ignore:
Timestamp:
03/16/14 20:38:39 (10 years ago)
Author:
pinsard
Message:

fix thanks to coding rules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/extrap.pro

    r41 r48  
    22; NAME: extrap.pro (based on remplit.pro) 
    33; 
    4 ; PURPOSE: Extrapolates ocean data values on land points  
     4; PURPOSE: Extrapolates ocean data values on land points 
    55; 
    66; CATEGORY: Subroutine 
     
    88; CALLING SEQUENCE: extrap, zdata, zmask, it, val 
    99; 
    10 ; INPUTS:  
     10; INPUTS: 
    1111;          zdata : field to extrapolate 
    1212;          zmask : field's mask 
    1313;          it    : iteration 
    14 ;          val   : field value on masked points   
     14;          val   : field value on masked points 
    1515; 
    1616; KEYWORD PARAMETERS: None 
    1717; 
    18 ; OUTPUTS:  
    19 ;          zdata : extrapolated field  
     18; OUTPUTS: 
     19;          zdata : extrapolated field 
    2020;          zmask : mask after extrapolation 
    2121; 
    2222; COMMON BLOCKS: 
    23 ;       common_interp.pro  
     23;          common_interp.pro 
    2424; 
    2525; 
     
    3838; IF keyword_set(key_performance) THEN print, systime(1)-temp 
    3939; 
    40 ; 2. Extrapolating when needed  
     40; 2. Extrapolating when needed 
    4141; ============================ 
    4242; 
     
    7171; IF keyword_set(key_performance) THEN print, systime(1)-temp 
    7272; IF keyword_set(key_performance) THEN temp = systime(1) 
    73 ; IF n_elements(ocean2) NE jpiatm*(jpjatm+4) THEN BEGIN  
     73; IF n_elements(ocean2) NE jpiatm*(jpjatm+4) THEN BEGIN 
    7474; land = where(zmask EQ 0.) 
    7575; zremplit(land) = val 
    76 ; ENDIF  
     76; ENDIF 
    7777; IF keyword_set(key_performance) THEN print, systime(1)-temp 
    7878; 
     
    9595   terre = where(mmmask EQ 0) 
    9696   if terre[0] EQ -1 then GOTO, fini 
    97 ; les points du bord du cadre doivent maintenan etre dans la terre 
     97; les points du bord du cadre doivent maintenant etre dans la terre 
    9898   mmmask[0, *] = 0 
    9999   mmmask[jpiatm+1, *] = 0 
     
    119119    +1./sqrt(2)*(z[(jpiatm+2)+1+cote]+z[(jpiatm+2)-1+cote]+z[-(jpiatm+2)+1+cote]+z[-(jpiatm+2)-1+cote]) 
    120120   poids = voisins[cote] 
    121        
     121 
    122122   z[cote] = zcote/poids 
    123123;--------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.