- Timestamp:
- 08/26/11 17:37:13 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/ws_tropflux_1d_to_1m_ncdf.pro
r94 r95 50 50 ; IDL> ws_tropflux_1d_to_1m_ncdf 51 51 ; 52 ; KNOWN ISSUES 53 ; ============ 54 ; 55 ; test of existence of fullfilename_in not very efficient because 56 ; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 57 ; 52 58 ; TODO 53 59 ; ==== 54 60 ; 55 ; remove hard coded creation date 61 ; coding rules 62 ; 63 ; get rid of hard coded yyyymmdd and nb month 56 64 ; 57 65 ; include in the whole process … … 65 73 ; 66 74 ; $URL$ 75 ; 76 ; - fplod 20110826T151225Z aedon.locean-ipsl.upmc.fr (Darwin) 77 ; 78 ; * remove hard coded creation date 67 79 ; 68 80 ; - fplod 20110817T084001Z aedon.locean-ipsl.upmc.fr (Darwin) … … 118 130 ; 119 131 ; check if this file exists 132 msg='iii : looking for ' + filename 133 ras = report(msg) 120 134 fullfilename = isafile(iodirout + filename, NEW=0, /MUST_EXIST) 121 135 IF fullfilename[0] EQ '' THEN BEGIN … … 136 150 ; 137 151 initncdf, fullfilename 138 st=19890101 152 st=19890101 139 153 en=20101231 140 154 ws=read_ncdf("ws", st, en, file=fullfilename,/nostr) … … 146 160 ;; Monthly data 147 161 ; 148 jpt=22*12. 162 jpt=22*12. 149 163 time=timegen(jpt,start=julday(1,15,1989),units='Months') 150 164 caldat, tt,m,d,y … … 165 179 cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 166 180 time=time-julday(1,1,1950,00,00,00) 167 xlon=reform(glamt(firstxt:lastxt,0) )181 xlon=reform(glamt(firstxt:lastxt,0)) 168 182 ylat=reform(gphit(0,firstyt:lastyt)) 169 183 … … 175 189 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 176 190 ws_attr={units:'m/s',missing_value:1.e20,long_name:'Mean wind speed at 10m height',short_name:'ws',axis:'TYX'} 177 globattr={Source:'TropFlux - Air-Sea Fluxes for the Global Tropics', Metodology:'Praveen Kumar et al., Clim. Dyn 2011', Producer_Agency:'Joint research colloboration between IPSL, Paris and NIO, India', Time_range:cda0+' - '+cda1, Creation_date:'20110815', Website:'http://www.locean-ipsl.upmc.fr/tropflux/'} 191 192 creation_date=string(jul2date(SYSTIME( /JULIAN , /UTC )), format='(I8)') 193 194 globattr={Source:'TropFlux - Air-Sea Fluxes for the Global Tropics' $ 195 , Methodology:'Praveen Kumar et al., Clim. Dyn 2011' $ 196 , Producer_Agency:'Joint research colloboration between IPSL, Paris and NIO, India' $ 197 , Time_range:cda0+' - '+cda1 $ 198 , Creation_date:creation_date $ 199 , Website:'http://www.locean-ipsl.upmc.fr/tropflux/'} 178 200 179 201 … … 189 211 190 212 end 191
Note: See TracChangeset
for help on using the changeset viewer.