Changeset 100


Ignore:
Timestamp:
08/30/11 18:01:00 (13 years ago)
Author:
pinsard
Message:

start to homogenize (to be cont.)

Location:
trunk/src
Files:
24 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/TropFlux_19890101_20091231.pro

    r99 r100  
    6565; if this file not already exists. 
    6666; 
    67 ; This output file :file:`${PROJECT_OD}/TropFlux_19890101_20091231_coarev3.nc` 
    68 ; will be used by :ref:`TropFlux_NRT_ncdf.pro`. 
     67; This output file 
     68; :file:`${PROJECT_OD}/TropFlux_19890101_20091231_coarev3.nc` 
     69; will be used by 
     70; :ref:`TropFlux_NRT_ncdf.pro`. 
    6971; 
    7072;     .. graphviz:: 
    … …  
    9193;           {mask file_sst file_ws file_wg file_swr file_lwr file_t2m file_q2m} -> {tropflux_19890101_20091231} -> {file_out} 
    9294; 
    93 ;          } 
     95;        } 
    9496; 
    9597; SEE ALSO 
    … …  
    560562lon_attr={units:'degrees_east',long_name:'Longitude'} 
    561563lat_attr={units:'degrees_north',long_name:'Latitude'} 
     564time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     565globattr={source:'Fluxes for the Global Tropical Ocean - TropFlux',timerange:cda0+' - '+cda1} 
    562566swr_attr={units:'W/m2',missing_value:valmask,long_name:'Surface net shortwave radiation',short_name:'swr',axis:'TYX'} 
    563567lwr_attr={units:'W/m2',missing_value:valmask,long_name:'Surface net longwave radiation',short_name:'lwr',axis:'TYX'} 
    … …  
    567571wg_attr={units:'m/s',missing_value:valmask,long_name:'COARE convective gustiness',short_name:'wg',axis:'TYX'} 
    568572tau_attr={units:'N/m2',missing_value:valmask,long_name:'Wind stress magnitude',short_name:'tau',axis:'TYX'} 
    569 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
    570 ; 
    571 globattr={source:'Fluxes for the Global Tropical Ocean - TropFlux',timerange:cda0+' - '+cda1} 
    572573; 
    573574help, swr,lwr,lat,sen,time,xlon,ylat 
    574575; 
    575 ncfields = $ 
    576           + tau[longitude,latitude,time]=tau:tau_attr; ' $ 
    577           + wg[longitude,latitude,time]=wg:wg_attr; ' $ 
    578           + swr[longitude,latitude,time]=swr:swr_attr; ' $ 
    579 ;          +'lwr_coare[longitude,latitude,time]=lwr:lwr_attr; ' $ 
    580           +'lwr[longitude,latitude,time]=lwnet_clrk:lwr_clrk_attr; ' $ 
    581           +'lhf[longitude,latitude,time]=lat:lhf_attr; ' $ 
    582           +'shf[longitude,latitude,time]=sen:shf_attr; ' $ 
    583 ;          +'Ch[longitude,latitude,time]=Ch:Ch_attr; ' $ 
    584 ;          +'Ce[longitude,latitude,time]=Ce:Ce_attr; ' $ 
    585                       + 'longitude[]=xlon:lon_attr; ' $ 
    586                       + 'latitude[]=ylat:lat_attr; ' $ 
    587                       + 'time[*time]=time:time_attr ' $ 
    588                       + ' @ globattr' 
     576ncfields = 'tau[longitude,latitude,time]=tau:tau_attr; ' $ 
     577         + 'wg[longitude,latitude,time]=wg:wg_attr; ' $ 
     578         + 'swr[longitude,latitude,time]=swr:swr_attr; ' $ 
     579;         +'lwr_coare[longitude,latitude,time]=lwr:lwr_attr; ' $ 
     580         +'lwr[longitude,latitude,time]=lwnet_clrk:lwr_clrk_attr; ' $ 
     581         +'lhf[longitude,latitude,time]=lat:lhf_attr; ' $ 
     582         +'shf[longitude,latitude,time]=sen:shf_attr; ' $ 
     583;         +'Ch[longitude,latitude,time]=Ch:Ch_attr; ' $ 
     584;         +'Ce[longitude,latitude,time]=Ce:Ce_attr; ' $ 
     585         + 'longitude[]=xlon:lon_attr; ' $ 
     586         + 'latitude[]=ylat:lat_attr; ' $ 
     587         + 'time[*time]=time:time_attr ' $ 
     588         + ' @ globattr' 
    589589; 
    590590@ncdf_quickwrite 

  • trunk/src/TropFlux_NRT_ncdf.pro

    r99 r100  
    202202help, swr, lwr, lhf, shf 
    203203; 
    204 ;time=julday(1,1,1989)+lindgen(7670)+0.5 
    205 ;jpt=n_elements(time) 
    206 time=timegen(7670, units='days', start=julday(1,1,1989,0)) 
     204time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    207205jpt=n_elements(time) 
    208206; 
    … …  
    219217lat_attr={units:'degrees_north',long_name:'Latitude'} 
    220218time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     219globattr={Source:'TropFlux - Air-Sea Fluxes for the Global Tropics', Resolution:'1 x 1 degree and daily', IDL_program_name:'TropFlux_ncdf.pro',Timerange:cda0+' - '+cda1, More_Info:'ISCCP based SWR until 2007 and OLR based SWR after that with a 90 day linear transition between the two'} 
    221220sw_attr={units:'w/m^2',missing_value:1.e20,long_name:'Net Shortwave Radiation',short_name:'swr',axis:'TYX'} 
    222221lw_attr={units:'w/m^2',missing_value:1.e20,long_name:'Net Longwave Radiation',short_name:'lwr',axis:'TYX'} 
    223222sh_attr={units:'w/m^2',missing_value:1.e20,long_name:'Sensible Heat Flux',short_name:'shf',axis:'TYX'} 
    224223lh_attr={units:'w/m^2',missing_value:1.e20,long_name:'Latent Heat Flux',short_name:'lhf',axis:'TYX'} 
    225 ; 
    226 globattr={Source:'TropFlux - Air-Sea Fluxes for the Global Tropics', Resolution:'1 x 1 degree and daily', IDL_program_name:'TropFlux_ncdf.pro',Timerange:cda0+' - '+cda1, More_Info:'ISCCP based SWR until 2007 and OLR based SWR after that with a 90 day linear transition between the two'} 
    227 ; 
    228224; 
    229225ncfields = 'swr[longitude,latitude,time]=swr:sw_attr; ' $ 

  • trunk/src/TropFlux_swr_BLND_19890101_20091231.pro

    r97 r100  
    2828; if this file not already exists. 
    2929; 
    30 ; This file will be used by :ref:`TropFlux_NRT_ncdf.pro`. 
     30; This file will be used by 
     31; :ref:`TropFlux_NRT_ncdf.pro`. 
    3132; 
    3233;     .. graphviz:: 
    … …  
    274275swr_merged=swr_dt*0. 
    275276; 
    276 time=timegen(7670, units='days', start=julday(1,1,1989,0)) 
     277time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    277278jpt=n_elements(time) 
    278279; 
    … …  
    296297lat_attr={units:'degrees_north',long_name:'Latitude'} 
    297298time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     299globattr={source:'ISCCP based SWR until 2007, after 2007 timeseries is completed with OLR based SWR.  A linear transition is applied from ISCCP based SWR to OLR based SWR'} 
    298300sw_attr={units:'w/m^2',missing_value:1.e20,long_name:'Net Shortwave Radiation',short_name:'swr',axis:'TYX'} 
    299 globattr={source:'ISCCP based SWR until 2007, after 2007 timeseries is completed with OLR based SWR.  A linear transition is applied from ISCCP based SWR to OLR based SWR'} 
    300301; 
    301302ncfields = 'swr[longitude,latitude,time]=swr_merged:sw_attr; ' $ 

  • trunk/src/TropFlux_swr_DT_19890101_20071231.pro

    r98 r100  
    2020; if this file not already exists. 
    2121; 
    22 ; This file will be used by :ref:`TropFlux_swr_BLND_19890101_20091231.pro`. 
     22; This file will be used by 
     23; :ref:`TropFlux_swr_BLND_19890101_20091231.pro`. 
    2324; 
    2425;     .. graphviz:: 
    … …  
    200201; 
    201202swr=sw_m+swa-sw_bias       ; new swr, corrected for bias and variability 
    202 ;writing field 
    203 ;time=julday(1,1,1989)+lindgen(7670) 
     203; 
     204; writing field 
    204205cda0=string(jul2date(time(0)),format='(i8.8)') 
    205206cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    … …  
    214215lat_attr={units:'degrees_north',long_name:'Latitude'} 
    215216time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     217globattr={source:'Obtained from ISCCP.  Corrections has been done on mean bias and variability',timerange:cda0+' - '+cda1} 
    216218sw_attr={units:'w/m^2',missing_value:1.e20,long_name:'Net Shortwave Radiation',short_name:'swr',axis:'TYX'} 
    217 globattr={source:'Obtained from ISCCP.  Corrections has been done on mean bias and variability',timerange:cda0+' - '+cda1} 
    218219; 
    219220ncfields = 'swr[longitude,latitude,time]=swr:sw_attr; ' $ 

  • trunk/src/TropFlux_swr_NRT_19890101_20091231.pro

    r98 r100  
    2424; if this file not already exists. 
    2525; 
    26 ; This file will be used by :ref:`TropFlux_swr_BLND_19890101_20091231.pro`. 
     26; This file will be used by 
     27; :ref:`TropFlux_swr_BLND_19890101_20091231.pro`. 
    2728; 
    2829;     .. graphviz:: 
    … …  
    9697; 
    9798;   * replace tt by time 
     99;   * fill empty time_origin:wq 
    98100; 
    99101; - fplod 20110812T075406Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    246248olra=olr-olrm                       ; olr anomaly 
    247249rec=swm+(olra/0.760364)             ; reconstructed sw 
    248 ;writing field 
    249 time=julday(1,1,1989)+lindgen(7670) 
     250; 
     251; writing field 
     252time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    250253cda0=string(jul2date(time(0)),format='(i8.8)') 
    251254cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    … …  
    259262lon_attr={units:'degrees_east',long_name:'Longitude'} 
    260263lat_attr={units:'degrees_north',long_name:'Latitude'} 
    261 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:''} 
     264time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     265globattr={source:'Reconstructed from corrected OLR anomalies and TropFlux mean SWR climatology',timerange:cda0+' - '+cda1} 
    262266sw_attr={units:'w/m^2',missing_value:1.e20,long_name:'Net Shortwave Radiation',short_name:'swr',axis:'TYX'} 
    263 globattr={source:'Reconstructed from corrected OLR anomalies and TropFlux mean SWR climatology',timerange:cda0+' - '+cda1} 
    264267; 
    265268; 

  • trunk/src/TropFlux_wind_stress_19890101_20101231.pro

    r98 r100  
     1;+ 
     2; 
     3; TODO 
     4; ==== 
     5; 
     6; integration to the whole process 
     7; 
     8; EVOLUTIONS 
     9; ========== 
     10; 
     11; $Id$ 
     12; 
     13; $URL$ 
     14; 
     15; - fplod 20110830T153220Z cratos (Linux) 
     16; 
     17;   * remove multiple line statements 
     18; 
     19; - pk 20110830 
     20; 
     21;   * creation 
     22; 
     23;- 
    124pro TropFlux_wind_stress_19890101_20101231 
    225@common 
    326;------------------------------------------------------------ 
    4 st=19880101 & en=20101231 
    5 st=19880101 & en=20101231 
     27st=19880101 
     28en=20101231 
     29st=19880101 
     30en=20101231 
    631; 
    732file="/Volumes/PAYASAM/TropFlux/ERA_C_1989_2010/TropFlux_19890101_20101231.nc" 
    … …  
    2045wg=read_ncdf('wg',st,en,file=file_wg,/nostr) 
    2146s=sqrt((ws*ws)+(wg*wg)) 
    22 s=reform(s) & help, s 
     47s=reform(s) 
     48help, s 
    2349; 
    2450file="/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_ws_19890101_20101231.nc" 
    … …  
    2753v=read_ncdf('v10',st,en,file=file,/nostr) 
    2854; 
    29 u=reform(u) & v=reform(v) 
     55u=reform(u) 
     56v=reform(v) 
    3057help, u,v 
    3158taux=tau*u/s 
    … …  
    3461help, taux, tauy 
    3562; 
    36 time=timegen(8035, start=julday(1,1,1989,0), units='days') & jpt=n_elements(time) 
     63time=timegen(8035, start=julday(1,1,1989,0), units='days') 
     64jpt=n_elements(time) 
    3765;------------------------------------------------------------ 
    3866cda0=string(jul2date(time(0)),format='(i8.8)') 
    3967cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    4068time=time-julday(1,1,1950,00,00,00) 
    41 xlon=reform(glamt(*,0) ) & ylat=reform(gphit(0,*)) 
     69xlon=reform(glamt(*,0) ) 
     70ylat=reform(gphit(0,*)) 
    4271; 
    4372ncfile='!/Volumes/PAYASAM/TropFlux/ERA_C_1989_2010/tau_tropflux_1d_19890101_20101231.nc' 
    4473lon_attr={units:'degrees_east',long_name:'Longitude'} 
    4574lat_attr={units:'degrees_north',long_name:'Latitude'} 
     75time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     76globattr={Source:'TropFlux - Momentum Fluxes for the Global Tropics', Methodology:'Praveen Kumar et al., Clim. Dyn 2011', Producer_Agency:'Joint research colloboration between IPSL, Paris and NIO, India', Time_range:cda0+' - '+cda1, Creation_date:'20110707', Website:'http://www.locean-ipsl.upmc.fr/tropflux/'} 
    4677tau_attr={units:'N/m2',missing_value:valmask,long_name:'wind stress magnitude',short_name:'tau',axis:'TYX'} 
    4778taux_attr={units:'N/m2',missing_value:valmask,long_name:'zonal wind stress',short_name:'taux',axis:'TYX'} 
    4879tauy_attr={units:'N/m2',missing_value:valmask,long_name:'meridional wind stress',short_name:'tauy',axis:'TYX'} 
    49 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
    5080; 
    51 globattr={Source:'TropFlux - Momentum Fluxes for the Global Tropics', Methodology:'Praveen Kumar et al., Clim. Dyn 2011', Producer_Agency:'Joint research colloboration between IPSL, Paris and NIO, India', Time_range:cda0+' - '+cda1, Creation_date:'20110707', Website:'http://www.locean-ipsl.upmc.fr/tropflux/'} 
    5281; 
    5382help, xlon,ylat, taux, tauy, tau 

  • trunk/src/cronin_gustiness_ncdf.pro

    r98 r100  
    1717; if this file not already exists. 
    1818; 
    19 ; This file will be used by :ref:`TropFlux_swr_BLND_19890101_20091231.pro` and 
     19; This file will be used by 
     20; :ref:`TropFlux_swr_BLND_19890101_20091231.pro` 
     21; and 
    2022; :ref:`TropFlux_19890101_20091231.pro`. 
    2123; 
    … …  
    181183; 
    182184; writing field 
    183 time=timegen(7670, units='days', start=julday(1,1,1989,0)) 
     185time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    184186jpt=n_elements(time) 
    185187; 
    … …  
    193195lon_attr={units:'degrees_east',long_name:'Longitude'} 
    194196lat_attr={units:'degrees_north',long_name:'Latitude'} 
    195 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
     197time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
    196198gust_attr={units:'m/s',missing_value:valmask,long_name:'Climatological gustiness',short_name:'wg',axis:'TYX'} 
    197199globattr={source:'Climatological wind gustiness obtained by fitting cronins gustiness values against sst ',timerange:cda0+' - '+cda1} 

  • trunk/src/d2m_to_q2m_erai.pro

    r98 r100  
    2525; ++ 
    2626; has been 
    27 ; produced by :ref:`interp_erai_msl_1989_2009.pro`. 
     27; produced by 
     28; :ref:`interp_erai_msl_1989_2009.pro`. 
    2829; 
    2930; q2m on OAFLUX grid 
    … …  
    3435; This file 
    3536; :file:`${PROJECT_OD}/erai_q2m_19890101_20091231_oafluxgrid.nc` 
    36 ; will be used by :ref:`q2m_correction_ncdf.pro`. 
     37; will be used by 
     38; :ref:`q2m_correction_ncdf.pro`. 
    3739; 
    3840; .. only:: man 
    … …  
    5860; 
    5961;           {file_d2m file_t2m file_msl} -> {d2m_to_q2m_erai} -> {file_q2m} 
     62; 
    6063;        } 
    6164; 
    … …  
    271274;writing field 
    272275ncfile='!' + fullfilename_out 
    273 time=timegen(7670, units='days', start=julday(1,1,1989)) 
     276time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    274277jpt=n_elements(time) 
    275278cda0=string(jul2date(time(0)),format='(i8.8)') 
    … …  
    284287lat_attr={units:'degrees_north',long_name:'Latitude'} 
    285288q2m_attr={units:'kg/kg',missing_value:valmask,long_name:'Surface specific humidity at 2m ',short_name:'q2m',axis:'TYX'} 
    286 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
     289time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
    287290globattr={source:'Surface specific humidity at 2m calculated from dew point temperature',timerange:cda0+' - '+cda1} 
    288291; 

  • trunk/src/interp_erai_dewt_1989_2009.pro

    r99 r100  
    241241d2min=read_ncdf("d2",19881201,20101231,file=fullfilename,/nostr) 
    242242; 
    243 timein=time 
    244243jptin=jpt 
    245244tab=d2min(*,*,0) 
    … …  
    256255help, d2min,lonin,latin,mskin,lonout,latout,mskout 
    257256; 
    258 si=size(d2min) 
    259257d2mout=fltarr(jpi,jpj,jptin) 
    260258for jt=0,jptin-1 do begin 
    … …  
    267265endfor 
    268266; 
    269 jptin=jpt 
    270 ; 
    271267initncdf, fullfilename_msk 
    272 cda0=string(jul2date(timein(0)),format='(i8.8)') 
    273 cda1=string(jul2date(timein(jpt-1)),format='(i8.8)') 
    274 ; 
    275 time=timegen(7670, units='days', start=julday(1,1,1989)) 
     268; 
     269time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    276270jpt=n_elements(time) 
    277271cda0=string(jul2date(time(0)),format='(i8.8)') 
    278272cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    279 timein=time-julday(1,1,1950,00,00) 
     273time=time-julday(1,1,1950,00,00) 
    280274; 
    281275lat=latout 
    … …  
    284278lon_attr={units:'degrees_east',long_name:'Longitude'} 
    285279lat_attr={units:'degrees_north',long_name:'Latitude'} 
    286 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
     280time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     281globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    287282d2m_attr={units:'degK',missing_value:1e20,long_name:'Dew Point Temperature at 2m',short_name:'d2m',axis:'TYX'} 
    288 globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    289283; 
    290284ncfields = 'd2m[longitude,latitude,time]=d2mout:d2m_attr; ' $ 
    291285                      + 'longitude[]=lon:lon_attr; ' $ 
    292286                      + 'latitude[]=lat:lat_attr; ' $ 
    293                       + 'time[*time]=timein:time_attr ' $ 
     287                      + 'time[*time]=time:time_attr ' $ 
    294288                      + ' @ globattr' 
    295289; 

  • trunk/src/interp_erai_lwr_1989_2009.pro

    r98 r100  
    1414; :file:`${PROJECT_ID}/20c3m_erai_str_TROP_1989_2009.nc` 
    1515; containing str from ERA-I 
    16 ; has been produced 
    17 ; by :ref:`compute_erai_daily_region_2d.sh`. 
     16; has been produced by 
     17; :ref:`compute_erai_daily_region_2d.sh`. 
    1818; 
    1919; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` 
    … …  
    2222; :ref:`oaflux_mask_30N30S.pro`. 
    2323; 
    24 ; Interpolated str is written in 
     24; Interpolated str 
     25; is written in 
    2526; :file:`${PROJECT_OD}/erai_lwr_19890102_20091231_oafluxgrid.nc` 
    2627; if this file not already exists. 
    … …  
    2829; This output file 
    2930; :file:`${PROJECT_OD}/erai_lwr_19890102_20091231_oafluxgrid.nc` 
    30 ; must be processed after by :ref:`add_198900101.sh`. 
     31; must be processed after by 
     32; :ref:`add_198900101.sh`. 
    3133; 
    3234; .. only:: man 
    … …  
    5254;           {file_in mask} -> {interp_erai_lwr_1989_2009} -> {file_out} 
    5355; 
    54 ;         } 
     56;        } 
    5557; 
    5658; SEE ALSO 
    … …  
    118120; 
    119121; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_lwr_1989_2009.pro $ 
     122; 
     123; - fplod 20110830T145010Z cratos (Linux) 
     124; 
     125;   * replace xlon by lon and xlat by lat; suppress extra initcdf call 
    120126; 
    121127; - fplod 20110830T135030Z cratos (Linux) 
    … …  
    272278print, 'lon grid ',min(lonin),max(lonin),lonin(1)-lonin(0) 
    273279; 
    274 timein=time 
    275280jptin=jpt 
    276281tab=lwrin(*,*,0) 
    … …  
    279284initncdf, fullfilename_msk 
    280285domdef 
    281 latout=reform(gphit(0,*)) lonout=reform(glamt(*,0)) 
     286latout=reform(gphit(0,*)) 
     287lonout=reform(glamt(*,0)) 
    282288print, 'lat grid ',min(latout),max(latout),latout(1)-latout(0) 
    283289print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0) 
    284290mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 
    285291; 
    286 si=size(lwrin) 
    287292lwrout=fltarr(jpi,jpj,jptin) 
    288293; 
    … …  
    300305time=timegen(7669, start=julday(1,2,1989,0), units='days') 
    301306jpt=n_elements(time) 
    302 jptin=jpt 
    303307cda0=string(jul2date(time(0)),format='(i8.8)') 
    304308cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    305309time=time-julday(1,1,1950,00,00,00) 
    306 xlon=reform(glamt(*,0) ) 
    307 ylat=reform(gphit(0,*)) 
    308 ; 
    309 initncdf, fullfilename_msk 
    310 valmask=1.e20 
    311 ; 
    312 ylat=latout 
    313 xlon=lonout 
     310; 
     311lat=latout 
     312lon=lonout 
    314313ncfile='!' + fullfilename_out 
    315314lon_attr={units:'degrees_east',long_name:'Longitude'} 
    316315lat_attr={units:'degrees_north',long_name:'Latitude'} 
    317 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
    318 ; 
     316time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     317globattr={source:'Tropical ocean fluxes obtained from ERA Interim',timerange:cda0+' - '+cda1} 
    319318lwr_attr={units:'W/m2',missing_value:valmask,long_name:'Surface net longwave radiation ',short_name:'lwr',axis:'TYX'} 
    320 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
    321 globattr={source:'Tropical ocean fluxes obtained from ERA Interim',timerange:cda0+' - '+cda1} 
    322319; 
    323320ncfields = 'lwr[longitude,latitude,time]=lwrout:lwr_attr; ' $ 
    324                       + 'longitude[]=xlon:lon_attr; ' $ 
    325                       + 'latitude[]=ylat:lat_attr; ' $ 
     321                      + 'longitude[]=lon:lon_attr; ' $ 
     322                      + 'latitude[]=lat:lat_attr; ' $ 
    326323                      + 'time[*time]=time:time_attr ' $ 
    327324                      + ' @ globattr' 

  • trunk/src/interp_erai_msl_1989_2009.pro

    r99 r100  
    5656;           {file_in mask} -> {interp_erai_msl_1989_2009} -> {file_out} 
    5757; 
    58 ;          } 
     58;        } 
    5959; 
    6060; SEE ALSO 
    … …  
    258258help, mslin,lonin,latin,mskin,lonout,latout,mskout 
    259259; 
    260 si=size(mslin) 
    261260mslout=fltarr(jpi,jpj,jptin) 
    262261for jt=0,jptin-1 do begin 
    … …  
    269268endfor 
    270269; 
    271 jptin=jpt 
    272 ; 
    273270initncdf, fullfilename_msk 
    274 ;time=julday(1,2,1989)+lindgen(7516) 
    275 time=timegen(7670, units='days', start=julday(1,1,1989)) 
    276 jpt=n_elements(time) 
     271time=timegen(jptin, start=julday(1,1,1989,0), units='days') 
    277272cda0=string(jul2date(time(0)),format='(i8.8)') 
    278273cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    … …  
    284279lon_attr={units:'degrees_east',long_name:'Longitude'} 
    285280lat_attr={units:'degrees_north',long_name:'Latitude'} 
    286 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
     281time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     282globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    287283msl_attr={units:'milibars',missing_value:1e20,long_name:'Mean Sea leval pressure',short_name:'msl',axis:'TYX'} 
    288 globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    289284; 
    290285ncfields = 'msl[longitude,latitude,time]=mslout:msl_attr; ' $ 

  • trunk/src/interp_erai_sst_1989_2009.pro

    r99 r100  
    5555;           {file_in mask} -> {interp_erai_sst_1989_2009} -> {file_out} 
    5656; 
    57 ;          } 
     57;        } 
    5858; 
    5959; SEE ALSO 
    … …  
    252252help, sstkin,lonin,latin,mskin,lonout,latout,mskout 
    253253; 
    254 si=size(sstkin) 
    255254sstkout=fltarr(jpi,jpj,jptin) 
    256255for jt=0,jptin-1 do begin 
    … …  
    263262endfor 
    264263; 
    265 jptin=jpt 
    266 ; 
    267264initncdf, fullfilename_msk 
    268 time=timegen(7670, units='days', start=julday(1,1,1989)) 
     265time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    269266jpt=n_elements(time) 
    270267cda0=string(jul2date(time(0)),format='(i8.8)') 
    … …  
    277274lon_attr={units:'degrees_east',long_name:'Longitude'} 
    278275lat_attr={units:'degrees_north',long_name:'Latitude'} 
    279 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
     276time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     277globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    280278sstk_attr={units:'degK',missing_value:1e20,long_name:'Sea Surface Temperature',short_name:'sst',axis:'TYX'} 
    281 globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    282279; 
    283280ncfields = 'sst[longitude,latitude,time]=sstkout:sstk_attr; ' $ 

  • trunk/src/interp_erai_t2m_1989_2009.pro

    r99 r100  
    5050;           fontname=Courier, 
    5151;           color=blue, 
    52 ;           URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/interp_erai_t2m_1989_2009.pro"", 
     52;           URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/interp_erai_t2m_1989_2009.pro", 
    5353;           label="${PROJECT}/src/interp_erai_t2m_1989_2009.pro"]; 
    5454; 
    5555;           {file_in mask} -> {interp_erai_t2m_1989_2009} -> {file_out} 
    5656; 
    57 ;          } 
     57;        } 
    5858; 
    5959; SEE ALSO 
    … …  
    373373ncdf_getatt, fullfilename_t2, 'time', long_name=long_name 
    374374time_attr={units:units, long_name:long_name, time_origin:' 1957-JAN-01 00:00:00'} 
     375globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    375376ncdf_getatt, fullfilename_t2, 't2', units=units 
    376377ncdf_getatt, fullfilename_t2, 't2', long_name=long_name 
    377378t2m_attr={units:units, long_name:long_name, missing_value:1e20, short_name:'t2m',axis:'TYX'} 
    378 globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    379379; 
    380380help, t2mout 

  • trunk/src/interp_erai_ws_1989_2009.pro

    r98 r100  
    6161;           {file_u10 file_v10 mask} -> {interp_erai_ws_1989_2009} -> {file_out} 
    6262; 
    63 ;         } 
     63;        } 
    6464; 
    6565; SEE ALSO 
    … …  
    306306help, u10in, v10in,lonin,latin,mskin,lonout,latout,mskout 
    307307; 
    308 si=size(u10in) 
    309308u10out=fltarr(jpi,jpj,jptin) 
    310309v10out=fltarr(jpi,jpj,jptin) 
    … …  
    326325endfor 
    327326; 
    328 time=timegen(7670, units='days', start=julday(1,1,1989)) 
     327time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    329328jpt=n_elements(time) 
    330 jptin=jpt 
    331329cda0=string(jul2date(time(0)),format='(i8.8)') 
    332330cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    … …  
    343341lon_attr={units:'degrees_east',long_name:'Longitude'} 
    344342lat_attr={units:'degrees_north',long_name:'Latitude'} 
    345 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
    346 ; 
     343time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     344globattr={source:'Tropical ocean winds obtained from ERA Interim',timerange:cda0+' - '+cda1} 
    347345u10_attr={units:'m/s',missing_value:valmask,long_name:'10 metre u wind component',short_name:'u10',axis:'TYX'} 
    348346v10_attr={units:'m/s',missing_value:valmask,long_name:'10 metre v wind component',short_name:'v10',axis:'TYX'} 
    349 globattr={source:'Tropical ocean winds obtained from ERA Interim',timerange:cda0+' - '+cda1} 
    350347; 
    351348help, u10out,v10out,time,xlon,ylat 

  • trunk/src/interp_olr_30n30s_1989_2009.pro

    r99 r100  
    4343;           {file_in mask} -> {interp_olr_30n30s_1989_2009} -> {file_out} 
    4444; 
    45 ;         } 
     45;        } 
    4646; 
    4747; SEE ALSO 
    … …  
    261261endfor 
    262262; 
    263 jptin=jpt 
    264263timein=timein-julday(1,1,1950,00,00) 
    265264; 
    … …  
    269268lon_attr={units:'degrees_east',long_name:'Longitude'} 
    270269lat_attr={units:'degrees_north',long_name:'Latitude'} 
    271 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'} 
     270time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     271globattr={source:'NOAA OLR, cf http://www.cdc.noaa.gov/cdc/data.interp_OLR.html, interpolated to 1x1 degree grid and with oaflux mask'} 
    272272olr_attr={units:'w/m^2',missing_value:1e20,long_name:'Outgoing Longwave Radiation',short_name:'olr',axis:'TYX'} 
    273 olra_attr={units:'w/m^2',missing_value:1e20,long_name:'Outgoing Longwave Radiation Anomaly',short_name:'olra',axis:'TYX'} 
    274 globattr={source:'NOAA OLR, cf http://www.cdc.noaa.gov/cdc/data.interp_OLR.html, interpolated to 1x1 degree grid and with oaflux mask'} 
    275273; 
    276274ncfields = 'olr[longitude,latitude,time]=olrout:olr_attr; ' $ 

  • trunk/src/lwr_correction_ncdf.pro

    r98 r100  
    2222; This output file 
    2323; :file:`${PROJECT_OD}/TropFlux_lwr_19890101_20091231.nc` 
    24 ; will be used by :ref:`TropFlux_19890101_20091231.pro`. 
     24; will be used by 
     25; :ref:`TropFlux_19890101_20091231.pro`. 
    2526; 
    2627; .. only:: man 
    … …  
    220221lat=reform(gphit(0,0:jpj-1)) 
    221222lon=reform(glamt(0:jpi-1,0)) 
    222 time=timegen(7670, units='days', start=julday(1,1,1989)) 
     223time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    223224jpt=n_elements(time) 
    224225; 
    … …  
    233234lat_attr={units:'degrees_north',long_name:'Latitude'} 
    234235time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     236globattr={source:'Basic data obtained from ERAI.  Mean correction for  bias and correction for variability are applied',timerange:cda0+' - '+cda1} 
    235237lwr_attr={units:'degK',missing_value:1.e20,long_name:'Air Temperature at 2m',short_name:'lwr',axis:'TYX'} 
    236 globattr={source:'Basic data obtained from ERAI.  Mean correction for  bias and correction for variability are applied',timerange:cda0+' - '+cda1} 
    237238; 
    238239ncfields = 'lwr[longitude,latitude,time]=lwr_new:lwr_attr; ' $ 

  • trunk/src/paper01/fig1/fig1_no_obs_map.pro

    r97 r100  
    1111; :file:`${PROJECT_OD}/mask_oaflux_30N30S.nc` 
    1212; containing OAFLUX grid 
    13 ; has been produced by :ref:`oaflux_mask_30N30S.pro`. 
     13; has been produced by 
     14; :ref:`oaflux_mask_30N30S.pro`. 
    1415; 
    1516; Figure of 

  • trunk/src/paper01/fig11/swr_isccp_tropflux_correlation_fig10.pro

    r97 r100  
    1212; containing 
    1313; sst corrected on OAFLUX grid 
    14 ; has been produced by :ref:`sst_correction_ncdf.pro`. 
     14; has been produced by 
     15; :ref:`sst_correction_ncdf.pro`. 
    1516; 
    1617; :file:`${PROJECT_OD}/correlation_isccp_olr.idl` 
    1718; containing 
    1819; ++ 
    19 ; has been produced by :ref:`swr_isccp_tropflux_new_v1.pro` 
     20; has been produced by 
     21; :ref:`swr_isccp_tropflux_new_v1.pro` 
    2022; 
    2123; .. graphviz:: 

  • trunk/src/paper01/fig11/swr_isccp_tropflux_new_v1.pro

    r97 r100  
    2121; This output file 
    2222; :file:`${PROJECT_OD}/correlation_isccp_olr.idl` 
    23 ; will be used by :ref:`swr_isccp_tropflux_correlation_fig10.pro`. 
     23; will be used by 
     24; :ref:`swr_isccp_tropflux_correlation_fig10.pro`. 
    2425; 
    2526; .. graphviz:: 

  • trunk/src/q2m_correction_ncdf.pro

    r98 r100  
    2323; if this file not already exists. 
    2424; 
    25 ; This file will be used by :ref:`TropFlux_19890101_20091231.pro`. 
     25; This file will be used by 
     26; :ref:`TropFlux_19890101_20091231.pro`. 
    2627; 
    2728; .. only:: man 
    … …  
    204205lat=reform(gphit(0,0:jpj-1)) 
    205206lon=reform(glamt(0:jpi-1,0)) 
    206 time=timegen(7670, units='days', start=julday(1,1,1989,0)) 
     207time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    207208jpt=n_elements(time) 
    208209; 
    … …  
    217218lat_attr={units:'degrees_north',long_name:'Latitude'} 
    218219time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     220globattr={source:'Basic data obtained from ERAI.  Mean correction for bias and correction for variability are applied',timerange:cda0+' - '+cda1} 
    219221q2m_attr={units:'g/kg',missing_value:1.e20,long_name:'Specific humidity at 2m height',short_name:'q2m',axis:'TYX'} 
    220 globattr={source:'Basic data obtained from ERAI.  Mean correction for bias and correction for variability are applied',timerange:cda0+' - '+cda1} 
    221222; 
    222223ncfields = 'q2m[longitude,latitude,time]=q2m_new:q2m_attr; ' $ 

  • trunk/src/sst_correction_ncdf.pro

    r99 r100  
    2222; This output file 
    2323; :file:`${PROJECT_OD}/TropFlux_sst_19890101_20091231.nc` 
    24 ; will be used by :ref:`cronin_gustiness_ncdf.pro` 
     24; will be used by 
     25; :ref:`cronin_gustiness_ncdf.pro` 
    2526; and 
    2627; :ref:`TropFlux_19890101_20091231.pro`. 
    … …  
    279280lat=reform(gphit(0,0:jpj-1)) 
    280281lon=reform(glamt(0:jpi-1,0)) 
    281 time=timegen(7670, units='days', start=julday(1,1,1989)) 
     282time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    282283jpt=n_elements(time) 
    283284; 
    … …  
    292293lat_attr={units:'degrees_north',long_name:'Latitude'} 
    293294time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     295globattr={source:'Basic data obtained from ERAI.  Bias and variability correction are applied',timerange:cda0+' - '+cda1} 
    294296sst_attr={units:'degK',missing_value:1.e20,long_name:'Sea Surface Temperature',short_name:'sst',axis:'TYX'} 
    295 globattr={source:'Basic data obtained from ERAI.  Bias and variability correction are applied',timerange:cda0+' - '+cda1} 
    296297; 
    297298ncfields = 'sst[longitude,latitude,time]=sst_new:sst_attr; ' $ 

  • trunk/src/t2m_correction_ncdf.pro

    r97 r100  
    2323; This output file 
    2424; :file:`${PROJECT_OD}/TropFlux_t2m_19890101_20091231.nc` 
    25 ; will be used by :ref:`TropFlux_19890101_20091231.pro`. 
     25; will be used by 
     26; :ref:`TropFlux_19890101_20091231.pro`. 
    2627; 
    2728; .. only:: man 
    … …  
    294295cda1=string(da2) 
    295296; 
    296 % conversion for hours since 19570101 to days from 19500101 
     297; conversion for hours since 19570101 to days from 19500101 
    297298time=julday(01, 01, 1957,timein)-julday(1,1,1950) 
    298299jpt=n_elements(time) 
    … …  
    306307ncdf_getatt, fullfilename_t2_uncor, 'latitude', long_name=long_name 
    307308lat_attr={units:units, long_name:long_name} 
    308 ; 
    309309time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
    310 ; 
     310globattr={source:'Basic data obtained from ERAI.  Mean correction for air temperautre bias and correction for variability are applied',timerange:cda0+' - '+cda1} 
    311311ncdf_getatt, fullfilename_t2_uncor, 't2m', units=units 
    312312ncdf_getatt, fullfilename_t2_uncor, 't2m', long_name=long_name 
    313313t2m_attr={units:units,missing_value:1.e20,long_name:long_name,short_name:'t2m',axis:'TYX'} 
    314 globattr={source:'Basic data obtained from ERAI.  Mean correction for air temperautre bias and correction for variability are applied',timerange:cda0+' - '+cda1} 
    315314; 
    316315ncfields = 't2m[longitude,latitude,time]=t2m_new:t2m_attr; ' $ 

  • trunk/src/ws_correction_ncdf.pro

    r98 r100  
    1515; :ref:`interp_erai_ws_1989_2009.pro`. 
    1616; 
    17 ; Corrected ws on OAFLUX grid is written in 
     17; Corrected ws on OAFLUX grid 
     18; is written in 
    1819; :file:`${PROJECT_OD}/TropFlux_ws_19890101_20091231.nc` 
    1920; if this file not already exists. 
    … …  
    2122; This output file 
    2223; :file:`${PROJECT_OD}/TropFlux_ws_19890101_20091231.nc` 
    23 ; will be used by :ref:`TropFlux_19890101_20091231.pro`. 
     24; will be used by 
     25; :ref:`TropFlux_19890101_20091231.pro`. 
    2426; 
    2527; .. only:: man 
    … …  
    4446;           {file_in} -> {ws_correction_ncdf} -> {file_out} 
    4547; 
    46 ;          } 
     48;        } 
    4749; 
    4850; SEE ALSO 
    … …  
    211213lat=reform(gphit(0,0:jpj-1)) 
    212214lon=reform(glamt(0:jpi-1,0)) 
    213 time=timegen(7670, units='days', start=julday(1,1,1989,0)) 
     215time=timegen(7670, start=julday(1,1,1989,0), units='days') 
    214216jpt=n_elements(time) 
    215217; 
    … …  
    224226lat_attr={units:'degrees_north',long_name:'Latitude'} 
    225227time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'} 
     228globattr={source:'Basic data obtained from ERAI. Bias and variability correction are applied',timerange:cda0+' - '+cda1} 
    226229w_attr={units:'m/s',missing_value:1.e20,long_name:'mean wind speed',short_name:'w',axis:'TYX'} 
    227 globattr={source:'Basic data obtained from ERAI. Bias and variability correction are applied',timerange:cda0+' - '+cda1} 
    228230; 
    229231ncfields = 'ws[longitude,latitude,time]=w_new:w_attr; ' $ 

  • trunk/src/ws_tropflux_1d_to_1m_ncdf.pro

    r99 r100  
    3636;           {file_1d} -> {ws_tropflux_1d_to_1m_ncdf} -> {file_1m} 
    3737; 
    38 ;         } 
     38;        } 
    3939; 
    4040; SEE ALSO 
    … …  
    159159; 
    160160jpt=22*12. 
    161 time=timegen(jpt,start=julday(1,15,1989),units='Months') 
    162 caldat, time,m,d,y 
    163 caldat, time, mon,day,yea 
     161time=timegen(jpt, start=julday(1,15,1989,0), units='Months') 
     162caldat, time, m, d, y 
     163caldat, time, mon, day, yea 
    164164; 
    165165wsm=fltarr(nxt,nyt,jpt) 
Note: See TracChangeset for help on using the changeset viewer.