- Timestamp:
- 12/14/10 12:36:07 (13 years ago)
- Location:
- trunk/src
- Files:
-
- 14 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/TropFlux_19890101_20091231.pro
r6 r7 14 14 ; gustiness values. 15 15 ; 16 ; 17 ; .. only:: man 18 ; 19 ; Figure is visible on PDF and HTML documents only. 20 ; 21 ; .. only:: html 22 ; 23 ; .. graphviz:: 24 ; 25 ; digraph TropFlux_19890101_20091231 { 26 ; graph [ 27 ; rankdir="LR", 28 ; ] 29 ; mask [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc"]; 30 ; file_sst [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_sst_19890101_20091231_v20.nc"]; 31 ; file_wind [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc"]; 32 ; file_wg [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_gustiness_19890101_20091231_v50.nc"]; 33 ; file_sw [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20091231_BLND_v50.nc"]; 34 ; file_lw [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_lwr_19890101_20091231_v2.nc"]; 35 ; file_air [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_t2m_19890101_20091231_v50.nc"]; 36 ; file_q [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc"]; 37 ; 38 ; ncfile [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_19890101_20091231_v51.nc"]; 39 ; 40 ; TropFlux_19890101_20091231 [shape=box, 41 ; fontname=Courier, 42 ; color=blue, 43 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/TropFlux_19890101_20091231.pro", 44 ; label="${TROPFLUX}/src/TropFlux_19890101_20091231.pro"]; 45 ; 46 ; {mask file_sst file_wind file_wg file_sw file_lw file_air file_q} -> {TropFlux_19890101_20091231} -> {ncfile} 47 ; 48 ; } 49 ; 16 50 ; SEE ALSO 17 51 ; ======== 18 52 ; 53 ; :func:`cor30a` 54 ; 19 55 ; EXAMPLES 20 56 ; ======== … … 34 70 ; ========== 35 71 ; 72 ; - fplod 20101214T112131Z aedon.locean-ipsl.upmc.fr (Darwin) 73 ; 74 ; * add graph 75 ; 36 76 ; - fplod 20101214T093615Z aedon.locean-ipsl.upmc.fr (Darwin) 37 77 ; … … 46 86 pro TropFlux_19890101_20091231 47 87 @common 48 ;------------------------------------------------------------49 ;50 88 51 89 da1=19880101 & da2=20101231 … … 146 184 tt=time 147 185 time=timegen(7670, start=julday(1,1,1989,0), units='days') & jpt=n_elements(time) 148 ; ------------------------------------------------------------186 ; 149 187 cda0=string(jul2date(time(0)),format='(i8.8)') 150 188 cda1=string(jul2date(time(jpt-1)),format='(i8.8)') -
trunk/src/TropFlux_NRT_ncdf.pro
r6 r7 6 6 ; TropFlux_NRT_ncdf_v50.pro 7 7 ; ========================= 8 ; 9 ; 10 ; .. only:: man 11 ; 12 ; Figure is visible on PDF and HTML documents only. 13 ; 14 ; .. only:: html 15 ; 16 ; .. graphviz:: 17 ; 18 ; digraph TropFlux_NRT_ncdf_v50 { 19 ; graph [ 20 ; rankdir="LR", 21 ; ] 22 ; file1 [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20091231_BLND.nc"]; 23 ; file2 [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_lwr_19890101_20091231_v50.nc"]; 24 ; file3 [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_19890101_20091231_v51.nc"]; 25 ; ncfile [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_19890101_20091231.nc"]; 26 ; TropFlux_NRT_ncdf_v50 [shape=box, 27 ; fontname=Courier, 28 ; color=blue, 29 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/TropFlux_NRT_ncdf_v50.pro", 30 ; label="${TROPFLUX}/src/TropFlux_NRT_ncdf_v50.pro" ]; 31 ; {file1 file2 file3} -> {TropFlux_NRT_ncdf_v50} -> {ncfile} 32 ; } 8 33 ; 9 34 ; SEE ALSO … … 15 40 ; :: 16 41 ; 17 ; IDL> trop Flux_nrt_ncdf_v5042 ; IDL> tropflux_nrt_ncdf_v50 18 43 ; 19 44 ; TODO … … 27 52 ; ========== 28 53 ; 54 ; - fplod 20101214T112526Z aedon.locean-ipsl.upmc.fr (Darwin) 55 ; 56 ; * add graph 29 57 ; - fplod 20101214T093615Z aedon.locean-ipsl.upmc.fr (Darwin) 30 58 ; … … 38 66 pro TropFlux_NRT_ncdf_v50 39 67 @common 40 ;----------------------------------------------------------------- 68 41 69 dir="/Users/pkb/data/TropFlux/" 42 70 file=dir+"TropFlux_swr_19890101_20091231_BLND.nc" … … 64 92 65 93 time=time-julday(1,1,1950,00,00) 66 ; -----------------------------------------------------------------94 ; 67 95 lat=reform(gphit(0,0:jpj-1)) 68 96 lon=reform(glamt(0:jpi-1,0)) -
trunk/src/TropFlux_swr_BLND_19890101_20091231.pro
r6 r7 40 40 pro TropFlux_swr_BLND_19890101_20091231 41 41 @common 42 ;------------------------------------------------43 42 dir="/Users/pkb/data/TropFlux/" 44 43 … … 68 67 file=dir+"TropFlux_gustiness_19890101_20091231_v50.nc" 69 68 initncdf, file 70 ;------------------------------------------------71 69 time=timegen(7670, units='days', start=julday(1,1,1989,0)) & jpt=n_elements(time) 72 70 lat=reform(gphit(0,0:jpj-1)) -
trunk/src/TropFlux_swr_DT_19890101_20071231.pro
r6 r7 39 39 pro TropFlux_swr_DT_19890101_20071231 40 40 @common 41 ;----------------------------------------------------------------42 41 file='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/swr_oafluxgrid_1985_2007.nc' 43 42 initncdf, file 44 43 45 44 sw=read_ncdf("swr", 19890101, 20071231, file=file,/nostr) 46 ;----------------------------------------------------------------47 45 help, sw 48 46 … … 50 48 help, sw_mean 51 49 52 ;----------------------------------------------------------------53 50 tt=time & jpt=n_elements(time) 54 51 caldat, time,mon,day,yea … … 67 64 68 65 swr=sw_m+swa-sw_bias ;; new swr, corrected for bias and variability 69 ;----------------------------------------------------------------70 66 ;writing field 71 67 ;time=julday(1,1,1989)+lindgen(7670) -
trunk/src/TropFlux_swr_NRT_19890101_20091231.pro
r6 r7 39 39 pro TropFlux_swr_NRT_19890101_20091231 40 40 @common 41 ;----------------------------------------------------------------42 41 file='/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20071231_DT_v50.nc' 43 42 initncdf, file … … 45 44 46 45 swi=sw 47 ;----------------------------------------------------------------48 46 file_olr='/Users/pkb/data/heat_budget/NOAA_OLR/olr_oafluxgrid_30n30s_19890101_20091231.nc' 49 47 … … 79 77 swi=swi[*,*,365:365+364] & olri=olri[*,*,365:365+364] 80 78 help, swi 81 ;----------------------------------------------------------------82 79 83 80 initncdf, file_olr … … 96 93 olra=olr-olrm ; olr anomaly 97 94 rec=swm+(olra/0.760364) ; reconstructed sw 98 ;----------------------------------------------------------------99 95 ;writing field 100 96 time=julday(1,1,1989)+lindgen(7670) -
trunk/src/cronin_gustiness_ncdf.pro
r6 r7 38 38 pro cronin_gustiness_ncdf 39 39 @common 40 ;---------------------------------------------------------------------------41 40 da1=19880101 & da2=20101231 42 41 date1=19890101 & date2=20081231 … … 60 59 61 60 ;; writing field 62 ;---------------------------------------------------------------------------63 61 time=timegen(7670, units='days', start=julday(1,1,1989,0)) & jpt=n_elements(time) 64 62 -
trunk/src/d2m_to_q2m_erai.pro
r6 r7 9 9 ; SEE ALSO 10 10 ; ======== 11 ; 12 ; :func:`rh_to_spechum(rh,t,P)` 11 13 ; 12 14 ; EXAMPLES … … 59 61 rh=(e/es)*100 60 62 61 ;----------------------------------------------------------------------62 63 ;; conversion of relative humidity to sp. humidity 63 64 ;; Ref.1 - Gill, Appendix 4, … … 72 73 q2m=rh_to_spechum(rh,t,P) 73 74 74 ;----------------------------------------------------------------------75 75 ;;writing field 76 76 ncfile='!/Volumes/PRAVEEN/TropFlux/input_uncor/q2m_erai_19890101_20091231_oafluxgrid.nc' -
trunk/src/dew_temp_airt_rel_humidity_jv.pro
r6 r7 44 44 ;- 45 45 function dew_temp_airt_rel_humidity_jv, t,td 46 ;----------------------------------------------------------------------47 46 48 47 if (not keyword_set(t)) then t=28 ; air temperature in degC … … 50 49 es0=6.11 ; reference saturation vapor pressure (hPa) 51 50 52 ;----------------------------------------------------------------------53 51 ;; calculating fields. First saturation vapor pressure (Es) and 54 52 ;; actual vapor pressure are calculated. The ratio (E/Es)*100 -
trunk/src/interp_erai_t2m_1989_2009.pro
r6 r7 76 76 initncdf, gridout 77 77 help, t2mout 78 ;------------------------------------------------------------79 78 tt=timein 80 79 time=julday(1,1,1989)+lindgen(7670)+0.5 & jpt=n_elements(time) -
trunk/src/lwr_correction_ncdf.pro
r6 r7 39 39 pro lwr_correction_ncdf 40 40 @common 41 ;----------------------------------------------------------------------42 41 da1=19880101 & da2=20091231 43 42 … … 49 48 lwr_mean=grossemoyenne(lwr,'t',/nan) 50 49 help,lwr_mean 51 ;----------------------------------------------------------------------52 50 53 51 tt=time & jpt=n_elements(time) … … 77 75 help, lwr_new 78 76 79 ;----------------------------------------------------------------------80 77 ;writing field 81 78 lat=reform(gphit(0,0:jpj-1)) -
trunk/src/q2m_correction_ncdf.pro
r6 r7 38 38 pro q2m_correction_ncdf 39 39 @common 40 ;----------------------------------------------------------------------41 40 da1=19880101 & da2=20101231 42 41 … … 48 47 q2m_mean=grossemoyenne(q2m,'t',/nan) 49 48 help, q2m_mean 50 ;----------------------------------------------------------------------51 49 52 50 tt=time & jpt=n_elements(time) … … 68 66 help, q2m_ano 69 67 70 ;; applying the correction for var yability based on the scatter68 ;; applying the correction for variability based on the scatter 71 69 ;q2m_ano=q2m_ano*(1/0.919333) ;; (2000-2008) 72 70 q2m_ano=q2m_ano*(1/0.924674) ;; (2000-2009) … … 75 73 help, q2m_new 76 74 77 ;----------------------------------------------------------------------78 75 ;writing field 79 76 lat=reform(gphit(0,0:jpj-1)) -
trunk/src/sst_correction_ncdf.pro
r6 r7 39 39 pro sst_correction_ncdf 40 40 @common 41 ;----------------------------------------------------------------------42 41 da1=19880101 & da2=20091231 43 42 … … 49 48 sst_mean=grossemoyenne(sst,'t',/nan) 50 49 help, sst_mean 51 ;----------------------------------------------------------------------52 50 53 51 tt=time & jpt=n_elements(time) … … 77 75 help, sst_new 78 76 79 ;----------------------------------------------------------------------80 77 ;writing field 81 78 lat=reform(gphit(0,0:jpj-1)) -
trunk/src/t2m_correction_ncdf.pro
r6 r7 38 38 pro t2m_correction_ncdf 39 39 @common 40 ;----------------------------------------------------------------------41 40 da1=19890101 & da2=20091231 42 41 … … 48 47 t2m_mean=grossemoyenne(t2m,'t',/nan) 49 48 help, w_mean 50 ;----------------------------------------------------------------------51 49 52 50 tt=time & jpt=n_elements(time) … … 82 80 help, t2m_new 83 81 84 ;----------------------------------------------------------------------85 82 ;writing field 86 83 lat=reform(gphit(0,0:jpj-1)) -
trunk/src/ws_correction_ncdf.pro
r6 r7 39 39 pro ws_correction_ncdf 40 40 @common 41 ;----------------------------------------------------------------------42 41 da1=19880101 & da2=20101231 43 42 … … 52 51 w_mean=grossemoyenne(w,'t',/nan) 53 52 help, w_mean 54 ;----------------------------------------------------------------------55 53 56 54 tt=time & jpt=n_elements(time) … … 80 78 help, w_new 81 79 82 ;----------------------------------------------------------------------83 80 ;writing field 84 81 lat=reform(gphit(0,0:jpj-1))
Note: See TracChangeset
for help on using the changeset viewer.