- Timestamp:
- 01/24/11 17:50:19 (13 years ago)
- Location:
- trunk/src
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/sst_correction_ncdf.pro
r20 r31 1 1 ;+ 2 2 ; 3 ; .. sst_correction_ncdf.pro:3 ; .. _sst_correction_ncdf.pro: 4 4 ; 5 5 ; ======================= -
trunk/src/test/read.jnl
r30 r31 19 19 ! ======== 20 20 ! 21 ! ferret> CANCEL DATA/ALL 22 ! ferret> GO read 21 ! :: 22 ! 23 ! ferret> GO read 24 ! 25 ! or:: 26 ! 27 ! $ ferret -nojnl -script read.jnl 23 28 ! 24 29 ! SEE ALSO … … 63 68 !- 64 69 ! 70 CANCEL DATA/ALL 71 CANCEL SYMBOL/ALL 72 ! 65 73 ! pour creer des metafile convertible en ps 66 74 set mode metafile … … 68 76 go bold 69 77 70 LET yyyymmdd_b = "19890101" 71 LET yyyymmdd_e = "20091231" 78 !++ http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2010/msg00690.html 79 !let mm = 1 80 !say `mm,zw=2` 81 ! ++ http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2009/msg00771.html 82 !repeat/range=1:2/name=m (use "http://dodsserver/sst_`m`.nc"; show dat) 83 !repeat/range=1:12/name=mm ( show sym mm; say `mm,zw=2`) 84 !repeat/range=1:12/name=mm ( show sym mm; LET CMM B`mm,zw=2`A; LIST CMM) 72 85 73 ! ++ pb floatstr si yyyymmdd_b est un entier LET fullfilename = "/usr/zeus/incas/fplod/tropflux_d"+"TropFlux_" + floatstr(yyyymmdd_b,"(I8.8)") + "_" + floatstr(yyyymmdd_e,"(I8.8)") + ".nc" 74 LET filename = "TropFlux_" + yyyymmdd_b + "_" + yyyymmdd_e + ".nc" 75 LET fullfilename = "/usr/zeus/incas/fplod/tropflux_d/" + filename 76 LIST fullfilename 86 ! pour les log 87 ! http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2009/msg00036.html 88 show symbol session* 89 let date = "($session_date)" 90 let datetime = "($session_date) ($session_time)" 91 92 ! DEFINE VARIABLE yyyymmdd_b = "19890101" 93 ! DEFINE VARIABLE yyyymmdd_e = "20091231" 94 !++ DEFINE VARIABLE cb = {spawn:printf '%8.8d' `b` "} 95 96 DEFINE SYMBOL yyyymmdd_b = 19890101 97 DEFINE SYMBOL yyyymmdd_e = 20091231 98 99 ! ++ pb floatstr si yyyymmdd_b est un entier DEFINE VARIABLE fullfilename = "/usr/zeus/incas/fplod/tropflux_d"+"TropFlux_" + floatstr(yyyymmdd_b,"(I8.8)") + "_" + floatstr(yyyymmdd_e,"(I8.8)") + ".nc" 100 !DEFINE VARIABLE filename = "TropFlux_" + yyyymmdd_b + "_" + yyyymmdd_e + ".nc" 101 DEFINE SYMBOL filename = TropFlux_($yyyymmdd_b)_($yyyymmdd_e).nc 102 ! ++ $TROPFLUX_D defini ? 103 !DEFINE SYMBOL tropflux_d = /usr/zeus/incas/fplod/tropflux_d/ 104 !DEFINE SYMBOL fullfilename = ($tropflux_d)($filename) 105 DEFINE VARIABLE tropflux_d = {spawn:"echo $TROPFLUX_ID"} 106 DEFINE SYMBOL fullfilename = "$tropflux_d($filename)" 107 CANCEL SYMBOL filename 77 108 ! open 109 ! test if the file exists 110 !++DEFINE VARIABLE a = {"dummy", SPAWN:"ls `fullfilename`"} 111 !++DEFINE VARIABLE nfiles = `a,RETURN=IEND` 112 !++CANCEL a 113 !++LIST nfiles 114 !++IF `nfiles EQ 2` THEN EXIT 115 !++CANCEL nfiles 116 78 117 ! ++ use fullfilename 79 use "/usr/zeus/incas/fplod/tropflux_d/TropFlux_19890101_20091231.nc" 118 !use "/usr/zeus/incas/fplod/tropflux_d/TropFlux_19890101_20091231.nc" 119 use "`fullfilename`" 80 120 ! 81 121 ! plot the first timestep of swr … … 86 126 ! 87 127 ! create a Postscript file 88 !++sp Fprint -l cps -p portrait -o TropFlux_19890101_20091231.ferret.ps metafile.plt128 sp Fprint -l cps -p portrait -o TropFlux_19890101_20091231.ferret.ps metafile.plt 89 129 ! 90 ! 130 ! create a GIF file 91 131 frame/file=TropFlux_19890101_20091231.ferret.gif 92 132 ! -
trunk/src/test/read.m
r29 r31 35 35 % TODO 36 36 % ==== 37 % 38 % 39 % ++ see http://www.esrl.noaa.gov/psd/data/gridded/tools.html 37 40 % 38 41 % make it usable for real end-user -
trunk/src/test/read.py
r29 r31 7 7 from Numeric import * 8 8 from Scientific.IO.NetCDF import * 9 from pylab import * 10 9 11 import time 10 12 … … 46 48 ========== 47 49 50 - fplod 20110120T164029Z aedon.locean-ipsl.upmc.fr (Darwin) 51 52 * thanks to http://www.unidata.ucar.edu/software/netcdf/examples/programs/sfc_pres_temp_rd.py 53 54 - fplod 20110107T134657Z aedon.locean-ipsl.upmc.fr (Darwin) 55 56 * add a main level for sphinx referencement (not yet ok) 57 48 58 - fplod 20110106T133620Z aedon.locean-ipsl.upmc.fr (Darwin) 49 59 … … 53 63 """ 54 64 55 yyyymmdd_b=19890101; 56 yyyymmdd_e=20091231; 57 # 58 filename = 'TropFlux_' + str(yyyymmdd_b) + '_' + str(yyyymmdd_e) + '.nc' 59 fullfilename = '/usr/zeus/incas/fplod/tropflux_d/' + filename 60 # 61 file = NetCDFFile(fullfilename, 'r') 62 print file.dimensions 63 print file.variables 64 # 65 # close NetCDF file 66 file.close() 65 if __name__ == '__main__': 67 66 68 67 yyyymmdd_b=19890101; 68 yyyymmdd_e=20091231; 69 # 70 filename = 'TropFlux_' + str(yyyymmdd_b) + '_' + str(yyyymmdd_e) + '.nc' 71 fullfilename = '/usr/zeus/incas/fplod/tropflux_d/' + filename 72 # 73 file = NetCDFFile(fullfilename, 'r') 74 # controle print 75 print file.dimensions 76 print file.variables 77 # 78 # get latitudes 79 ydata=file.variables['latitude'] 80 # controle print 81 print ydata.units 82 # get longitudes 83 xdata=file.variables['longitude'] 84 # get swr variable 85 zdata=file.variables['swr'] 86 # controle print 87 print zdata.units 88 # plot 89 xlabel('xlabel+todo') 90 ylabel('ylabel+todo') 91 # close NetCDF file 92 file.close() -
trunk/src/tropflux_profile.sh
r2 r31 44 44 # ======== 45 45 # 46 # for IDL/GDL environment :47 #48 # :ref:`tropflux_startup.pro`49 46 # 50 47 # EXAMPLES
Note: See TracChangeset
for help on using the changeset viewer.