Changeset 25 for trunk/ToBeReviewed/CALCULS/hdyn.pro
- Timestamp:
- 05/02/06 14:59:12 (18 years ago)
- File:
-
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
trunk/ToBeReviewed/CALCULS/hdyn.pro
r23 r25 101 101 case 1 of 102 102 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, *] 105 105 end 106 106 tailles[1] eq nx and tailles[2] eq ny:BEGIN … … 116 116 e33d = replicate(1, nx*ny)#e3t 117 117 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) 119 119 if terre[0] NE -1 then vol[terre] = !values.f_nan 120 120 case level of … … 128 128 case 1 of 129 129 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, *, *] 132 132 end 133 133 tailles[1] eq nx and tailles[2] eq ny AND tailles[4] EQ jpt:BEGIN … … 144 144 e33d = e33d[*]#replicate(1, jpt) 145 145 e33d = reform(e33d, nx, ny, jpk, jpt, /over) 146 mask = tmask[ premierxt:dernierxt, premieryt:dernieryt, *]146 mask = tmask[firstxt:lastxt, firstyt:lastyt, *] 147 147 mask = mask[*]#replicate(1, jpt) 148 148 terre = where(mask EQ 0)
Note: See TracChangeset
for help on using the changeset viewer.