| 1 | ;+ |
|---|
| 2 | ; |
|---|
| 3 | ; make Z transport in Sv at T point |
|---|
| 4 | ; |
|---|
| 5 | ; @version |
|---|
| 6 | ; $Id$ |
|---|
| 7 | ; |
|---|
| 8 | ;- |
|---|
| 9 | FUNCTION make_ztran, file_name, ncdf_db $ |
|---|
| 10 | , BOXZOOM=boxzoom $ |
|---|
| 11 | , TIME_1=time_1 $ |
|---|
| 12 | , ALL_DATA=all_data $ |
|---|
| 13 | , TIME_2=time_2 $ |
|---|
| 14 | , ZMTYP=zmtyp |
|---|
| 15 | ; |
|---|
| 16 | compile_opt idl2, strictarrsubs |
|---|
| 17 | ; |
|---|
| 18 | @common |
|---|
| 19 | @com_eg |
|---|
| 20 | ; |
|---|
| 21 | ; |
|---|
| 22 | ; Read W |
|---|
| 23 | ; |
|---|
| 24 | file_nam = strmid(file_name, 0, strlen(file_name)-4) |
|---|
| 25 | ;old_boite = [lon1, lon2, lat1, lat2, prof1, prof2] |
|---|
| 26 | ;domdef |
|---|
| 27 | vargrid = 'W' |
|---|
| 28 | w = nc_read(file_nam+'W.nc','vovecrtz', ncdf_db, BOXZOOM = boxzoom, TIME_1 = time_1, TIME_2 = time_2, ALL_DATA = all_data) |
|---|
| 29 | |
|---|
| 30 | ; Transport along Z : trz = w*e1t*e2t (Sv) |
|---|
| 31 | |
|---|
| 32 | |
|---|
| 33 | e1t3d = reform(e1t, jpi*jpj) |
|---|
| 34 | e1t3d = reform(e1t3d#replicate(1, jpk), jpi, jpj, jpk) |
|---|
| 35 | |
|---|
| 36 | e2t3d = reform(e2t, jpi*jpj) |
|---|
| 37 | e2t3d = reform(e2t3d#replicate(1, jpk), jpi, jpj, jpk) |
|---|
| 38 | |
|---|
| 39 | trz = w.data*e1t3d*e2t3d |
|---|
| 40 | |
|---|
| 41 | trz = trz*1.e-6 |
|---|
| 42 | |
|---|
| 43 | ;domdef, old_boite |
|---|
| 44 | |
|---|
| 45 | field = {name: '', data: trz, legend: '', units: '', origin: '', dim: 0, direc: ''} |
|---|
| 46 | |
|---|
| 47 | field.origin = w.origin |
|---|
| 48 | field.dim = w.dim |
|---|
| 49 | field.direc = w.direc |
|---|
| 50 | |
|---|
| 51 | return, field |
|---|
| 52 | END |
|---|
