Changeset 226 for trunk/SRC/Grid/romsdepth.pro
- Timestamp:
- 03/16/07 10:22:26 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Grid/romsdepth.pro
r192 r226 16 16 ; common variable (cm_4mesh) romszinfos must be correctly defined 17 17 ; 18 ; @history 18 ; @history 19 19 ; Sept 2006 Sebastien Masson (smasson\@lodyc.jussieu.fr) 20 20 ; … … 41 41 grille, -1, -1, -1, -1, nx, ny, nz, firstx, firsty, firstz, lastx, lasty, lastz 42 42 hroms = hroms[firstx:lastx, firsty:lasty] 43 nt = n_elements(zeta)/nx/ny 43 nt = n_elements(zeta)/nx/ny 44 44 ; 45 45 cff1 = 1./sinh(theta_s) … … 47 47 ; 48 48 IF type EQ 'W' THEN BEGIN 49 sc = (findgen(jpk)-jpk)/jpk 50 ; sc = (dindgen(jpk+1)-jpk)/jpk 49 sc = (findgen(jpk)-jpk)/jpk 50 ; sc = (dindgen(jpk+1)-jpk)/jpk 51 51 ; jpk = jpk+1 52 52 ENDIF ELSE BEGIN 53 sc = (findgen(jpk)-jpk-0.5+1)/jpk 53 sc = (findgen(jpk)-jpk-0.5+1)/jpk 54 54 ENDELSE 55 55 ; … … 60 60 hinv = 1./hroms 61 61 hinv = hinv[*]#replicate(1., jpk) 62 ; put a z dimens tion to zeta62 ; put a z dimension to zeta 63 63 zeta = transpose(temporary(zeta)) 64 64 zeta = reform((temporary(zeta))[*]#replicate(1., jpk), nt, ny, nx, jpk, /overwrite)
Note: See TracChangeset
for help on using the changeset viewer.