Ignore:
Timestamp:
05/02/06 14:59:12 (18 years ago)
Author:
pinsard
Message:

upgrade of CALCULS according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/ : files

File:
1 copied

Legend:

Unmodified
Added
Removed
  • trunk/ToBeReviewed/CALCULS/hdyn.pro

    r23 r25  
    101101         case 1 of 
    102102            tailles[1] eq jpi and tailles[2] eq jpj: BEGIN 
    103                sn = tabsn[premierxt:dernierxt, premieryt:dernieryt, *] 
    104                tn = tabtn[premierxt:dernierxt, premieryt:dernieryt, *] 
     103               sn = tabsn[firstxt:lastxt, firstyt:lastyt, *] 
     104               tn = tabtn[firstxt:lastxt, firstyt:lastyt, *] 
    105105            end 
    106106            tailles[1] eq  nx and tailles[2] eq  ny:BEGIN  
     
    116116         e33d = replicate(1, nx*ny)#e3t 
    117117         e33d = reform(e33d, nx, ny, jpk, /over) 
    118          terre = where(tmask[premierxt:dernierxt, premieryt:dernieryt, *] EQ 0) 
     118         terre = where(tmask[firstxt:lastxt, firstyt:lastyt, *] EQ 0) 
    119119         if terre[0] NE -1 then vol[terre] = !values.f_nan 
    120120         case level of 
     
    128128         case 1 of 
    129129            tailles[1] eq jpi and tailles[2] eq jpj AND tailles[4] EQ jpt: BEGIN 
    130                sn = tabsn[premierxt:dernierxt, premieryt:dernieryt, *, *] 
    131                tn = tabtn[premierxt:dernierxt, premieryt:dernieryt, *, *] 
     130               sn = tabsn[firstxt:lastxt, firstyt:lastyt, *, *] 
     131               tn = tabtn[firstxt:lastxt, firstyt:lastyt, *, *] 
    132132            end 
    133133            tailles[1] eq  nx and tailles[2] eq  ny AND tailles[4] EQ jpt:BEGIN  
     
    144144         e33d = e33d[*]#replicate(1, jpt) 
    145145         e33d = reform(e33d, nx, ny, jpk, jpt, /over) 
    146          mask = tmask[premierxt:dernierxt, premieryt:dernieryt, *] 
     146         mask = tmask[firstxt:lastxt, firstyt:lastyt, *] 
    147147         mask = mask[*]#replicate(1, jpt) 
    148148         terre = where(mask EQ 0) 
Note: See TracChangeset for help on using the changeset viewer.