| [85] | 1 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 2 | OPENR, 22, 'input_coord' & READF, 22, lonu & READF, 22, latu & READF, 22, lsu & READF, 22, lctu & CLOSE, 22 |
|---|
| 3 | OPENR, 23, 'input_more' & READF, 23, hgtu, tsurfu & CLOSE, 23 |
|---|
| 4 | domain='d01' & filesWRF = FindFile('wrfout_'+domain+'_????-??-??_??:??:??') & nf=n_elements(filesWRF) |
|---|
| 5 | id=ncdf_open(filesWRF(0)) |
|---|
| 6 | NCDF_DIMINQ, id, NCDF_DIMID(id, 'west_east' ), toto, nx & NCDF_DIMINQ, id, NCDF_DIMID(id, 'south_north' ), toto, ny |
|---|
| 7 | NCDF_DIMINQ, id, NCDF_DIMID(id, 'bottom_top' ), toto, nz & NCDF_DIMINQ, id, NCDF_DIMID(id, 'Time' ), toto, nt |
|---|
| 8 | NCDF_CLOSE, id |
|---|
| 9 | id=ncdf_open(filesWRF(nf-1)) ;; for interrupted runs |
|---|
| 10 | NCDF_DIMINQ, id, NCDF_DIMID(id, 'Time' ), toto, ntlast |
|---|
| 11 | NCDF_CLOSE, id |
|---|
| 12 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 13 | yeye = 0 & nttot = (nf-1)*nt + ntlast & localtime = lctu + history_interval_s*findgen(nttot)/3700. |
|---|
| 14 | for loop = 0, nf-1 do begin |
|---|
| 15 | timetime = SYSTIME(1) |
|---|
| 16 | if (loop ne nf-1) then nloop2=nt else nloop2=ntlast |
|---|
| 17 | for loop2 = 0, nloop2-1 do begin |
|---|
| 18 | ;wprime = getget(filesWRF(loop), 'W', count=[0,0,0,1], offset=[0,0,0,loop2]) |
|---|
| 19 | |
|---|
| 20 | ;if (localtime(yeye) eq 11.) then begin |
|---|
| 21 | if (localtime(yeye) eq 14.) then begin |
|---|
| 22 | wprime = getget(filesWRF(loop), 'W', count=[0,0,0,1], offset=[0,0,0,loop2]) |
|---|
| 23 | print, 'save !! ', localtime(yeye) |
|---|
| 24 | save, wprime, h, filename='savew' |
|---|
| 25 | stop |
|---|
| 26 | uprime = getget(filesWRF(loop), 'U', count=[0,0,0,1], offset=[0,0,0,loop2]) |
|---|
| 27 | vprime = getget(filesWRF(loop), 'V', count=[0,0,0,1], offset=[0,0,0,loop2]) |
|---|
| 28 | help, uprime |
|---|
| 29 | help, vprime |
|---|
| 30 | veltot = sqrt(uprime^2+vprime^2) |
|---|
| 31 | print, 'save !! ', localtime(yeye) |
|---|
| 32 | save, veltot, h, filename='saveu' |
|---|
| 33 | stop |
|---|
| 34 | endif |
|---|
| 35 | |
|---|
| 36 | yeye = TEMPORARY(yeye) + 1 |
|---|
| 37 | endfor |
|---|
| 38 | endfor |
|---|
| 39 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 40 | |
|---|
