Changeset 97 for trunk/src/paper01/fig2/read_lw.pro
- Timestamp:
- 08/30/11 14:38:29 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/paper01/fig2/read_lw.pro
r94 r97 76 76 ; 77 77 @cm_project 78 79 ; ;DEFINE THE OUTPUT TIME AXIS78 ; 79 ; DEFINE THE OUTPUT TIME AXIS 80 80 ; 81 81 jda1=date2jul(date1) … … 83 83 jpt=(jda2-jda1+1l) 84 84 time=jda1+dindgen(jpt) 85 86 85 ; 87 ;; FIRST READ ALL DATA FROM SITE88 ;; - if file exists, extract correct time axis89 ;; - if not, fill variable with missing values90 86 ; 91 92 ;; LWR 87 ; FIRST READ ALL DATA FROM SITE 88 ; - if file exists, extract correct time axis 89 ; - if not, fill variable with missing values 90 ; 91 ; 92 ; LWR 93 93 fi=project_id_env+'lwnet'+csite+'_dy.cdf' 94 94 f=file_test(fi) … … 105 105 print, fi, ' not found' 106 106 endelse 107 108 109 107 ; 110 ;; Replace missing values by "NaN" 108 ; 109 ; 110 ; Replace missing values by "NaN" 111 111 ; 112 112 ;tsvars=['at','bf','bp','dyn','emp','evap','heat','iso','lw','lwnet','qlat','qsen','sw', $ … … 114 114 ;depvars=['u','v','d','t','s'] 115 115 vars=['lw'] 116 116 ; 117 117 nn=n_elements(vars) 118 118 for n=0,nn-1 do begin 119 119 var=vars(n) 120 com='ind=where('+var+' ge 1.e20) 120 com='ind=where('+var+' ge 1.e20) 121 121 if (ind(0) ne -1) then '+var+'(ind)=!values.f_nan' 122 122 r=execute(com) 123 123 endfor 124 124 ; 125 125 ;++nsmooth=nsmooth 126 126 ;++lw=smooth(lw,nsmooth,/nan) 127 127 ; 128 128 end
Note: See TracChangeset
for help on using the changeset viewer.