Changeset 98 for trunk/src


Ignore:
Timestamp:
08/30/11 16:10:19 (13 years ago)
Author:
pinsard
Message:

replace tt by time

Location:
trunk/src
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/TropFlux_19890101_20091231.pro

    r97 r98  
    127127; TODO 
    128128; ==== 
     129; 
     130; describe usage of tau 
    129131; 
    130132; make it work :: 
    … …  
    192194; 
    193195; $URL$ 
     196; 
     197; - fplod 20110830T135832Z cratos (Linux) 
     198; 
     199;   * replace tt by time 
    194200; 
    195201; - fplod 20110830T084129Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    545551  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,8)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & wg(*,*,jt)=tab*msk+valmask*(1-msk) 
    546552endfor 
    547 tt=time 
    548553time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    549554jpt=n_elements(time) 
    … …  
    551556cda0=string(jul2date(time(0)),format='(i8.8)') 
    552557cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    553 tt=time-julday(1,1,1950,00,00,00) 
     558time=time-julday(1,1,1950,00,00,00) 
    554559xlon=reform(glamt(*,0) ) 
    555560ylat=reform(gphit(0,*)) 
    … …  
    569574globattr={source:'Fluxes for the Global Tropical Ocean - TropFlux',timerange:cda0+' - '+cda1} 
    570575; 
    571 help, swr,lwr,lat,sen,tt,xlon,ylat 
     576help, swr,lwr,lat,sen,time,xlon,ylat 
    572577; 
    573578ncfields = $ 
    … …  
    583588                      + 'longitude[]=xlon:lon_attr; ' $ 
    584589                      + 'latitude[]=ylat:lat_attr; ' $ 
    585                       + 'tt[*time]=tt:time_attr ' $ 
     590                      + 'time[*time]=time:time_attr ' $ 
    586591                      + ' @ globattr' 
    587592; 

  • trunk/src/TropFlux_swr_DT_19890101_20071231.pro

    r97 r98  
    9090; 
    9191; $URL$ 
     92; 
     93; - fplod 20110830T135743Z cratos (Linux) 
     94; 
     95;   * replace tt by time 
    9296; 
    9397; - fplod 20110812T073138Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    179183help, sw_mean 
    180184; 
    181 tt=time 
    182185jpt=n_elements(time) 
    183186caldat, time,mon,day,yea 
    … …  
    199202;writing field 
    200203;time=julday(1,1,1989)+lindgen(7670) 
    201 cda0=string(jul2date(tt(0)),format='(i8.8)') 
    202 cda1=string(jul2date(tt(jpt-1)),format='(i8.8)') 
     204cda0=string(jul2date(time(0)),format='(i8.8)') 
     205cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    203206; 
    204207lat=reform(gphit(0,0:jpj-1)) 

  • trunk/src/TropFlux_swr_NRT_19890101_20091231.pro

    r97 r98  
    9292; 
    9393; $URL$ 
     94; 
     95; - fplod 20110830T135938Z cratos (Linux) 
     96; 
     97;   * replace tt by time 
    9498; 
    9599; - fplod 20110812T075406Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    228232initncdf, fullfilename_olr 
    229233olr=read_ncdf('olr', 19890100, 20091231, box=box, file=fullfilename_olr,/nostr) 
    230 tt=time 
    231234jpt=n_elements(time) 
    232235caldat, time,mon,day,yea 
    … …  
    245248;writing field 
    246249time=julday(1,1,1989)+lindgen(7670) 
    247 cda0=string(jul2date(tt(0)),format='(i8.8)') 
    248 cda1=string(jul2date(tt(jpt-1)),format='(i8.8)') 
     250cda0=string(jul2date(time(0)),format='(i8.8)') 
     251cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    249252; 
    250253lat=reform(gphit(0,0:jpj-1)) 

  • trunk/src/TropFlux_wind_stress_19890101_20101231.pro

    r97 r98  
    3838cda0=string(jul2date(time(0)),format='(i8.8)') 
    3939cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    40 tt=time-julday(1,1,1950,00,00,00) 
     40time=time-julday(1,1,1950,00,00,00) 
    4141xlon=reform(glamt(*,0) ) & ylat=reform(gphit(0,*)) 
    4242; 
    … …  
    5858                      + 'longitude[]=xlon:lon_attr; ' $ 
    5959                      + 'latitude[]=ylat:lat_attr; ' $ 
    60                       + 'tt[*time]=tt:time_attr ' $ 
     60                      + 'time[*time]=time:time_attr ' $ 
    6161                      + ' @ globattr' 
    6262; 

  • trunk/src/cronin_gustiness_ncdf.pro

    r97 r98  
    7676; 
    7777; $URL$ 
     78; 
     79; - fplod 20110830T140029Z cratos (Linux) 
     80; 
     81;   * replace tt by time 
    7882; 
    7983; - fplod 20110808T144208Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    182186cda0=string(jul2date(time(0)),format='(i8.8)') 
    183187cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    184 tt=time-julday(1,1,1950,00,00,00) 
     188time=time-julday(1,1,1950,00,00,00) 
    185189xlon=reform(glamt(*,0) ) 
    186190ylat=reform(gphit(0,*)) 
    … …  
    196200                      + 'longitude[]=xlon:lon_attr; ' $ 
    197201                      + 'latitude[]=ylat:lat_attr; ' $ 
    198                       + 'tt[*time]=tt:time_attr ' $ 
     202                      + 'time[*time]=time:time_attr ' $ 
    199203                      + ' @ globattr' 
    200204; 

  • trunk/src/d2m_to_q2m_erai.pro

    r97 r98  
    245245initncdf, fullfilename_t2m 
    246246t=read_ncdf("t2m", d1,d2, file=fullfilename_t2m,/nostr) 
    247 tt=time 
    248247jpt=n_elements(time) 
    249248help, t, td 
    … …  
    277276cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    278277; 
    279 tt=time-julday(1,1,1950,00,00,00) 
     278time=time-julday(1,1,1950,00,00,00) 
    280279xlon=reform(glamt(*,0) ) 
    281280ylat=reform(gphit(0,*)) 
    … …  
    291290                      + 'longitude[]=xlon:lon_attr; ' $ 
    292291                      + 'latitude[]=ylat:lat_attr; ' $ 
    293                       + 'tt[*time]=tt:time_attr ' $ 
     292                      + 'time[*time]=time:time_attr ' $ 
    294293                      + ' @ globattr' 
    295294; 

  • trunk/src/interp_erai_lwr_1989_2009.pro

    r97 r98  
    118118; 
    119119; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_lwr_1989_2009.pro $ 
     120; 
     121; - fplod 20110830T135030Z cratos (Linux) 
     122; 
     123;   * replace tt by time 
    120124; 
    121125; - pinsard 2011-08-23T09:09:16Z loholt1.ipsl.polytechnique.fr (Linux) 
    … …  
    296300time=timegen(7669, start=julday(1,2,1989,0), units='days') 
    297301jpt=n_elements(time) 
    298 tt=time 
    299302jptin=jpt 
    300 cda0=string(jul2date(tt(0)),format='(i8.8)') 
    301 cda1=string(jul2date(tt(jpt-1)),format='(i8.8)') 
    302 tt=tt-julday(1,1,1950,00,00,00) 
     303cda0=string(jul2date(time(0)),format='(i8.8)') 
     304cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
     305time=time-julday(1,1,1950,00,00,00) 
    303306xlon=reform(glamt(*,0) ) 
    304307ylat=reform(gphit(0,*)) 
    … …  
    321324                      + 'longitude[]=xlon:lon_attr; ' $ 
    322325                      + 'latitude[]=ylat:lat_attr; ' $ 
    323                       + 'tt[*time]=tt:time_attr ' $ 
     326                      + 'time[*time]=time:time_attr ' $ 
    324327                      + ' @ globattr' 
    325328; 

  • trunk/src/interp_erai_sst_1989_2009.pro

    r97 r98  
    241241ind=where((tab eq 1.e9) or finite(tab,/nan)) 
    242242mskin(ind)=0. 
    243 ;mskin=glamt*0.+1. 
    244243; 
    245244initncdf, fullfilename_msk 
    … …  
    270269time=timegen(7670, units='days', start=julday(1,1,1989)) 
    271270jpt=n_elements(time) 
    272 tt=time 
    273271cda0=string(jul2date(time(0)),format='(i8.8)') 
    274272cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 

  • trunk/src/interp_erai_ws_1989_2009.pro

    r97 r98  
    138138; 
    139139; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_ws_1989_2009.pro $ 
     140; 
     141; - fplod 20110830T135212Z cratos (Linux) 
     142; 
     143;   * replace tt by time 
    140144; 
    141145; - pinsard 2011-08-23T09:09:16Z loholt1.ipsl.polytechnique.fr (Linux) 
    … …  
    324328time=timegen(7670, units='days', start=julday(1,1,1989)) 
    325329jpt=n_elements(time) 
    326 tt=time 
    327330jptin=jpt 
    328 cda0=string(jul2date(tt(0)),format='(i8.8)') 
    329 cda1=string(jul2date(tt(jpt-1)),format='(i8.8)') 
    330 tt=tt-julday(1,1,1950,00,00,00) 
     331cda0=string(jul2date(time(0)),format='(i8.8)') 
     332cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
     333time=time-julday(1,1,1950,00,00,00) 
    331334xlon=reform(glamt(*,0) ) 
    332335ylat=reform(gphit(0,*)) 
    … …  
    346349globattr={source:'Tropical ocean winds obtained from ERA Interim',timerange:cda0+' - '+cda1} 
    347350; 
    348 help, u10out,v10out,tt,xlon,ylat 
     351help, u10out,v10out,time,xlon,ylat 
    349352; 
    350353ncfields = 'u10[longitude,latitude,time]=u10out:u10_attr; ' $ 
    … …  
    352355                      + 'longitude[]=xlon:lon_attr; ' $ 
    353356                      + 'latitude[]=ylat:lat_attr; ' $ 
    354                       + 'tt[*time]=tt:time_attr ' $ 
     357                      + 'time[*time]=time:time_attr ' $ 
    355358                      + ' @ globattr' 
    356359; 

  • trunk/src/interp_olr_30n30s_1989_2009.pro

    r97 r98  
    233233olrin=read_ncdf("olr",19881231, 20091231, file=fullfilename,/nostr) 
    234234; 
    235 mskin=glamt*0.+1. 
    236235; 
    237236timein=time 
    238237jptin=jpt 
     238mskin=glamt*0.+1. 
    239239; 
    240240initncdf, fullfilename_msk 

  • trunk/src/lwr_correction_ncdf.pro

    r97 r98  
    191191help,lwr_mean 
    192192; 
    193 tt=time 
    194193jpt=n_elements(time) 
    195194caldat, time,mon,day,yea 

  • trunk/src/paper01/fig13/net_flux_validation_scatter_2000_2007.pro

    r97 r98  
    135135jpt=(jda2-jda1+1l) 
    136136time=jda1+dindgen(jpt) 
    137 tt=time 
    138137; 
    139138ocean='net_flux' 

  • trunk/src/paper01/fig3/read_variables_v2.pro

    r97 r98  
    6565; ==== 
    6666; 
    67 ; externalize function 
    68 ; 
    6967; make it work (missing data) 
    7068; 
    … …  
    7573; 
    7674; $URL$ 
     75; 
     76; - fplod 20110830T140518Z cratos (Linux) 
     77; 
     78;   * remove function time_lec (already externalized) 
    7779; 
    7880; - fplod 20110412T110534Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    249251; 
    250252end 
    251 ; 
    252 function time_lec, fi 
    253 tt=ncdf_lec(fi,var='time') 
    254 fid=ncdf_open(fi) 
    255 vid=ncdf_varid(fid,'time') 
    256 ncdf_attget, fid,vid,'units',orig 
    257 orig=string(orig) 
    258 yy=long(strmid(orig,11,4)) 
    259 mm=long(strmid(orig,16,2)) 
    260 dd=long(strmid(orig,19,2)) 
    261 tt=julday(mm,dd,yy,12,00)+tt 
    262 ; 
    263 return, tt 
    264 ; 
    265 end 

  • trunk/src/q2m_correction_ncdf.pro

    r97 r98  
    176176help, q2m_mean 
    177177; 
    178 tt=time 
    179178jpt=n_elements(time) 
    180179caldat, time,mon,day,yea 

  • trunk/src/sst_correction_ncdf.pro

    r97 r98  
    247247help, sst_mean 
    248248; 
    249 tt=time 
    250249jpt=n_elements(time) 
    251250caldat, time,mon,day,yea 

  • trunk/src/ws_correction_ncdf.pro

    r97 r98  
    182182help, w_mean 
    183183; 
    184 tt=time 
    185184jpt=n_elements(time) 
    186185caldat, time,mon,day,yea 

  • trunk/src/ws_tropflux_1d_to_1m_ncdf.pro

    r97 r98  
    156156help, ws 
    157157; 
    158 tt=time 
    159 ; 
    160158; Monthly data 
    161159; 
    162160jpt=22*12. 
    163161time=timegen(jpt,start=julday(1,15,1989),units='Months') 
    164 caldat, tt,m,d,y 
     162caldat, time,m,d,y 
    165163caldat, time, mon,day,yea 
    166164; 
Note: See TracChangeset for help on using the changeset viewer.