- Timestamp:
- 08/09/11 16:12:48 (13 years ago)
- Location:
- trunk/src
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/TropFlux_19890101_20091231.pro
r50 r74 14 14 ; gustiness values. 15 15 ; 16 ; 16 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` containing OAFLUX grid have been produced by :ref:`oaflux_mask_30N30S.pro`. 17 ; 18 ; :file:`${PROJECT_OD}/TropFlux_sst_19890101_20091231.nc` containing sst corrected on OAFLUX grid have been produced by :ref:`sst_correction_ncdf.pro`. 19 ; 20 ; :file:`${PROJECT_OD}/TropFlux_ws_19890101_20091231.nc` containing ws corrected on OAFLUX grid have been produced by :ref:`ws_correction_ncdf.pro`. 21 ; 22 ; :file:`${PROJECT_OD}/TropFlux_gustiness_19890101_20091231.nc` containing ++ 23 ; have been produced by :ref:` cronin_gustiness_ncdf.pro`. 24 ; 25 ; :file:`${PROJECT_OD}/TropFlux_swr_19890101_20091231_BLND.nc` containing ws corrected on OAFLUX grid 26 ; have been produced by :ref:`TropFlux_swr_BLND_19890101_20091231.pro`. 27 ; 28 ; :file:`${PROJECT_OD}/TropFlux_lwr_19890101_20091231.nc` containing lwr corrected on OAFLUX grid have been produced by :ref:`lwr_correction_ncdf.pro`. 29 ; 30 ; :file:`${PROJECT_OD}/TropFlux_t2m_19890101_20091231.nc` containing t2m corrected on OAFLUX grid have been produced by :ref:`t2m_correction_ncdf.pro`. 31 ; 32 ; :file:`${PROJECT_OD}/TropFlux_q2m_19890101_20091231.nc` containing q2m corrected on OAFLUX grid have been produced by :ref:`d2m_to_q2m_erai.pro`. 33 ; 34 ; net heat flux components are written 35 ; in :file:`${PROJECT_OD}/TropFlux_19890101_20091231.nc` 36 ; if this file not already exists. 37 ; 17 38 ; .. graphviz:: 18 39 ; 19 ; digraph TropFlux_19890101_20091231 {40 ; digraph tropflux_19890101_20091231 { 20 41 ; graph [ 21 42 ; rankdir="LR", 22 43 ; ] 44 ; 23 45 ; mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/mask_oaflux_30N30S.nc"]; 24 ; file_sst [shape=ellipse,fontname=Courier,label=" /Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_sst_19890101_20091231_v20.nc"];25 ; file_w ind [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc"];26 ; file_wg [shape=ellipse,fontname=Courier,label=" /Users/pkb/data/TropFlux/TropFlux_gustiness_19890101_20091231_v50.nc"];27 ; file_sw [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20091231_BLND_v50.nc"];28 ; file_lw [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_lwr_19890101_20091231_v2.nc"];29 ; file_ air [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_t2m_19890101_20091231_v50.nc"];30 ; file_q [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc"];31 ; 32 ; ncfile [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_19890101_20091231_v51.nc"];33 ; 34 ; TropFlux_19890101_20091231 [shape=box,46 ; file_sst [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_sst_19890101_20091231.nc"]; 47 ; file_ws [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_ws_19890101_20091231.nc"]; 48 ; file_wg [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_gustiness_19890101_20091231.nc"]; 49 ; file_swr [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_swr_19890101_20091231_BLND.nc"]; 50 ; file_lwr [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_lwr_19890101_20091231.nc"]; 51 ; file_t2m [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_t2m_19890101_20091231.nc"]; 52 ; file_q2m [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_q2m_19890101_20091231.nc"]; 53 ; 54 ; file_out[shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_19890101_20091231.nc"]; 55 ; 56 ; tropflux_19890101_20091231 [shape=box, 35 57 ; fontname=Courier, 36 58 ; color=blue, … … 38 60 ; label="${PROJECT}/src/TropFlux_19890101_20091231.pro"]; 39 61 ; 40 ; {mask file_sst file_w ind file_wg file_sw file_lw file_air file_q} -> {TropFlux_19890101_20091231} -> {ncfile}62 ; {mask file_sst file_ws file_wg file_swr file_lwr file_t2m file_q2m} -> {tropflux_19890101_20091231} -> {file_out} 41 63 ; 42 64 ; } … … 64 86 ; :: 65 87 ; 88 ; IDL> .compile TropFlux_19890101_20091231 66 89 ; IDL> tropflux_19890101_20091231 67 90 ; … … 69 92 ; ==== 70 93 ; 71 ; hard coded directory - usage of ${PROJECT_ID} 94 ; I (fp) do not know how 95 ; ${PROJECT_OD}/TropFlux_swr_19890101_20091231_BLND.nc is produced. 96 ; 97 ; avoid mix lower/uppercase in pro name to avoid compile 72 98 ; 73 99 ; coding rules … … 82 108 ; ========== 83 109 ; 110 ; - fplod 20110809T110911Z aedon.locean-ipsl.upmc.fr (Darwin) 111 ; 112 ; * complete descritption 113 ; * remove v* from filenames (in and out) 114 ; * usage of ${PROJECT_OD} 115 ; * remove return statement 116 ; * add test on IO files 117 ; 84 118 ; - fplod 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin) 85 119 ; … … 101 135 ; 102 136 pro TropFlux_19890101_20091231 103 @common 137 ; 138 @cm_4cal 139 @cm_4data 140 @cm_4mesh 141 @cm_4data 104 142 @cm_project 105 143 ; … … 132 170 ; 133 171 ; check if this file exists 172 msg='iii : looking for ' + filename_msk 173 ras = report(msg) 134 174 fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 135 175 IF fullfilename_msk[0] EQ '' THEN BEGIN … … 138 178 STOP 139 179 ENDIF 140 180 ; 181 ; test if ${PROJECT_OD} defined 182 CASE project_od_env OF 183 '' : BEGIN 184 msg = 'eee : ${PROJECT_OD} is not defined' 185 ras = report(msg) 186 STOP 187 END 188 ELSE: BEGIN 189 msg = 'iii : ${PROJECT_OD} is ' + project_od_env 190 ras = report(msg) 191 END 192 ENDCASE 193 ; 194 ; check if output data will be possible 195 iodirout = isadirectory(project_od_env) 196 ; 197 ; existence and protection for reading 198 IF (FILE_TEST(iodirout, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 199 msg = 'eee : the directory' + iodirout + ' is not accessible.' 200 ras = report(msg) 201 STOP 202 ENDIF 203 ; 204 ; existence and protection for writing 205 IF (FILE_TEST(iodirout, /DIRECTORY, /WRITE) EQ 0) THEN BEGIN 206 msg = 'eee : the directory' + iodirout + ' was not found.' 207 ras = report(msg) 208 STOP 209 ENDIF 210 ; 211 ; build sst data filename 212 filename_sst='TropFlux_sst_19890101_20091231.nc' 213 ; 214 ; check if this file exists 215 msg='iii : looking for ' + filename_sst 216 ras = report(msg) 217 fullfilename_sst = isafile(iodirout + filename_sst, NEW=0, /MUST_EXIST) 218 IF fullfilename_sst[0] EQ '' THEN BEGIN 219 msg = 'eee : the file ' + fullfilename_sst + ' was not found.' 220 ras = report(msg) 221 STOP 222 ENDIF 223 ; 224 ; build ws data filename 225 filename_ws='TropFlux_ws_19890101_20091231.nc' 226 ; 227 ; check if this file exists 228 msg='iii : looking for ' + filename_ws 229 ras = report(msg) 230 fullfilename_ws = isafile(iodirout + filename_ws, NEW=0, /MUST_EXIST) 231 IF fullfilename_ws[0] EQ '' THEN BEGIN 232 msg = 'eee : the file ' + fullfilename_ws + ' was not found.' 233 ras = report(msg) 234 STOP 235 ENDIF 236 ; 237 ; build swr data filename 238 filename_swr='TropFlux_swr_19890101_20091231_BLND.nc' 239 ; 240 ; check if this file exists 241 msg='iii : looking for ' + filename_swr 242 ras = report(msg) 243 fullfilename_swr = isafile(iodirout + filename_swr, NEW=0, /MUST_EXIST) 244 IF fullfilename_swr[0] EQ '' THEN BEGIN 245 msg = 'eee : the file ' + fullfilename_swr + ' was not found.' 246 ras = report(msg) 247 STOP 248 ENDIF 249 ; 250 ; build lwr data filename 251 filename_lwr='TropFlux_lwr_19890101_20091231.nc' 252 ; 253 ; check if this file exists 254 msg='iii : looking for ' + filename_lwr 255 ras = report(msg) 256 fullfilename_lwr = isafile(iodirout + filename_lwr, NEW=0, /MUST_EXIST) 257 IF fullfilename_lwr[0] EQ '' THEN BEGIN 258 msg = 'eee : the file ' + fullfilename_lwr + ' was not found.' 259 ras = report(msg) 260 STOP 261 ENDIF 262 ; 263 ; build t2m data filename 264 filename_t2m='TropFlux_t2m_19890101_20091231.nc' 265 ; 266 ; check if this file exists 267 msg='iii : looking for ' + filename_t2m 268 ras = report(msg) 269 fullfilename_t2m = isafile(iodirout + filename_t2m, NEW=0, /MUST_EXIST) 270 IF fullfilename_t2m[0] EQ '' THEN BEGIN 271 msg = 'eee : the file ' + fullfilename_t2m + ' was not found.' 272 ras = report(msg) 273 STOP 274 ENDIF 275 ; 276 ; build q2m data filename 277 filename_q2m='TropFlux_q2m_19890101_20091231.nc' 278 ; 279 ; check if this file exists 280 msg='iii : looking for ' + filename_q2m 281 ras = report(msg) 282 fullfilename_q2m = isafile(iodirout + filename_q2m, NEW=0, /MUST_EXIST) 283 IF fullfilename_q2m[0] EQ '' THEN BEGIN 284 msg = 'eee : the file ' + fullfilename_q2m + ' was not found.' 285 ras = report(msg) 286 STOP 287 ENDIF 288 ; 289 ; build wg data filename 290 filename_wg='TropFlux_gustiness_19890101_20091231.nc' 291 ; 292 ; check if this file exists 293 msg='iii : looking for ' + filename_wg 294 ras = report(msg) 295 fullfilename_wg = isafile(iodirout + filename_wg, NEW=0, /MUST_EXIST) 296 IF fullfilename_wg[0] EQ '' THEN BEGIN 297 msg = 'eee : the file ' + fullfilename_wg + ' was not found.' 298 ras = report(msg) 299 STOP 300 ENDIF 301 ; 141 302 da1=19880101 & da2=20101231 142 303 ; … … 144 305 msk=ncdf_lec(fullfilename_msk,var='msk') 145 306 146 dir='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/' 147 dir1='/Users/pkb/data/TropFlux/' 148 149 file_sst=dir+'TropFlux_sst_19890101_20091231_v20.nc' 150 file_wind=dir+'TropFlux_ws_19890101_20091231_v20.nc' 151 file_sw=dir1+'TropFlux_swr_19890101_20091231_BLND_v50.nc' 152 file_lw=dir+'TropFlux_lwr_19890101_20091231_v2.nc' 153 file_air=dir1+'TropFlux_t2m_19890101_20091231_v50.nc' 154 file_q=dir+'TropFlux_q2m_19890101_20091231_v20.nc' 155 file_wg=dir1+'TropFlux_gustiness_19890101_20091231_v50.nc' 156 157 initncdf, file_sst 158 159 ws=read_ncdf('ws',da1,da2,file=file_wind,/nostr) 160 wg=read_ncdf('wg',da1-1,da2,file=file_wg,/nostr) 307 initncdf, fullfilename_sst 308 309 ws=read_ncdf('ws',da1,da2,file=fullfilename_ws,/nostr) 310 wg=read_ncdf('wg',da1-1,da2,file=fullfilename_wg,/nostr) 161 311 tt=time & jpt=n_elements(time) 162 sst=read_ncdf('sst',da1,da2,file=f ile_sst,/nostr)312 sst=read_ncdf('sst',da1,da2,file=fullfilename_sst,/nostr) 163 313 sst=reform(sst-273.15) 164 swd=read_ncdf('swr',da1,da2,file=f ile_sw,/nostr)165 lw=read_ncdf('lwr',da1,da2,file=f ile_lw,/nostr)314 swd=read_ncdf('swr',da1,da2,file=fullfilename_swr,/nostr) 315 lw=read_ncdf('lwr',da1,da2,file=fullfilename_lwr,/nostr) 166 316 swd=swd/0.94 ;; converting from net swr to downward swr 167 317 168 t2m=read_ncdf('t2m',da1,da2,file=f ile_air,/nostr)-273.15 ; in C169 q2m=read_ncdf('q2m',da1,da2,file=f ile_q,/nostr) ; in g/kg318 t2m=read_ncdf('t2m',da1,da2,file=fullfilename_t2m,/nostr)-273.15 ; in C 319 q2m=read_ncdf('q2m',da1,da2,file=fullfilename_q2m,/nostr) ; in g/kg 170 320 171 321 w=sqrt(ws*ws+wg*wg) ;; wind corrected for gustiness … … 242 392 xlon=reform(glamt(*,0) ) & ylat=reform(gphit(0,*)) 243 393 244 ncfile='! /Users/pkb/data/TropFlux/TropFlux_19890101_20091231_v51.nc'394 ncfile='!${PROJECT_OD}/TropFlux_19890101_20091231.nc' 245 395 lon_attr={units:'degrees_east',long_name:'Longitude'} 246 396 lat_attr={units:'degrees_north',long_name:'Latitude'} … … 272 422 @ncdf_quickwrite 273 423 274 return275 424 end -
trunk/src/TropFlux_NRT_ncdf.pro
r50 r74 1 1 ;+ 2 2 ; 3 ; .. _TropFlux_NRT_ncdf _v50.pro:3 ; .. _TropFlux_NRT_ncdf.pro: 4 4 ; 5 ; ===================== ====6 ; TropFlux_NRT_ncdf _v50.pro7 ; ===================== ====5 ; ===================== 6 ; TropFlux_NRT_ncdf.pro 7 ; ===================== 8 8 ; 9 9 ; .. graphviz:: 10 10 ; 11 ; digraph TropFlux_NRT_ncdf_v50{11 ; digraph tropflux_nrt_ncdf { 12 12 ; graph [ 13 13 ; rankdir="LR", … … 17 17 ; file3 [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_19890101_20091231_v51.nc"]; 18 18 ; ncfile [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_19890101_20091231.nc"]; 19 ; TropFlux_NRT_ncdf_v50 [shape=box, 19 ; 20 ; tropflux_nrt_ncdf [shape=box, 20 21 ; fontname=Courier, 21 22 ; color=blue, 22 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/TropFlux_NRT_ncdf_v50.pro", 23 ; label="${PROJECT}/src/TropFlux_NRT_ncdf_v50.pro" ]; 24 ; {file1 file2 file3} -> {TropFlux_NRT_ncdf_v50} -> {ncfile} 23 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/TropFlux_NRT_ncdf.pro", 24 ; label="${PROJECT}/src/TropFlux_NRT_ncdf.pro" ]; 25 ; 26 ; {file1 file2 file3} -> {tropflux_nrt_ncdf} -> {ncfile} 25 27 ; } 26 28 ; … … 39 41 ; :: 40 42 ; 41 ; IDL> tropflux_nrt_ncdf _v5043 ; IDL> tropflux_nrt_ncdf 42 44 ; 43 45 ; TODO … … 51 53 ; ========== 52 54 ; 55 ; - fplod 20110809T133815Z cratos (Linux) 56 ; 57 ; * remove v50 from pro name and in the graph 58 ; 53 59 ; - fplod 20101214T112526Z aedon.locean-ipsl.upmc.fr (Darwin) 54 60 ; 55 61 ; * add graph 62 ; 56 63 ; - fplod 20101214T093615Z aedon.locean-ipsl.upmc.fr (Darwin) 57 64 ; … … 63 70 ; 64 71 ;- 65 pro TropFlux_NRT_ncdf _v5072 pro TropFlux_NRT_ncdf 66 73 @common 67 74 -
trunk/src/TropFlux_swr_BLND_19890101_20091231.pro
r50 r74 7 7 ; ======================================= 8 8 ; 9 ; :file:`${PROJECT_OD}/TropFlux_swr_19890101_20071231_DT.nc` 10 ; containing ++ 11 ; have been produced by :ref:`TropFlux_swr_DT_19890101_20071231.pro`. 12 ; 13 ; :file:`${PROJECT_OD}/TropFlux_swr_19890101_20091231_NRT.nc` 14 ; containing ++ 15 ; have been produced by :ref:`TropFlux_swr_NRT_19890101_20091231.pro`. 16 ; 17 ; :file:`${PROJECT_OD}/TropFlux_gustiness_19890101_20091231.nc` containing ++ 18 ; have been produced by :ref:` cronin_gustiness_ncdf.pro`. 19 ; 20 ; ++ are written 21 ; in :file:`${PROJECT_OD}/TropFlux_swr_19890101_20091231_BLND.nc` 22 ; if this file not already exists. 23 ; 24 ; This file will be used by :ref:`TropFlux_19890101_20091231.pro`. 25 ; 9 26 ; .. graphviz:: 10 27 ; 11 ; digraph TropFlux_swr_BLND_19890101_20091231 {28 ; digraph tropflux_swr_blnd_19890101_20091231 { 12 29 ; graph [ 13 30 ; rankdir="LR", 14 31 ; ] 15 ; file_dt [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20071231_DT_v50.nc"];16 ; file_ nrt [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20091231_NRT_v50.nc"];17 ; file_ gustiness [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_gustiness_19890101_20091231_v50.nc"];18 ; 19 ; ncfile [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20091231_BLND.nc"];20 ; 21 ; TropFlux_swr_BLND_19890101_20091231 [shape=box,32 ; 33 ; file_dt [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_swr_19890101_20071231_DT.nc"]; 34 ; file_nrt [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_swr_19890101_20091231_NRT.nc"]; 35 ; file_gustiness [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_gustiness_19890101_20091231.nc"]; 36 ; file_out [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_swr_19890101_20091231_blnd.nc"]; 37 ; 38 ; tropflux_swr_blnd_19890101_20091231 [shape=box, 22 39 ; fontname=Courier, 23 40 ; color=blue, 24 41 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/TropFlux_swr_BLND_19890101_20091231.pro", 25 42 ; label="${PROJECT}/src/TropFlux_swr_BLND_19890101_20091231.pro"]; 26 ; {file_dt file_nrt file_gustiness} -> {TropFlux_swr_BLND_19890101_20091231} ->{ ncfile} 43 ; 44 ; {file_dt file_nrt file_gustiness} -> {tropflux_swr_blnd_19890101_20091231} -> {file_out} 27 45 ; 28 46 ; } … … 41 59 ; :: 42 60 ; 61 ; IDL> .compile TropFlux_swr_BLND_19890101_20091231 43 62 ; IDL> tropflux_swr_blnd_19890101_20091231 44 63 ; … … 47 66 ; ==== 48 67 ; 49 ; hard coded directory - usage of ${PROJECT_ID}68 ; avoid mix lower/uppercase in pro name to avoid compile 50 69 ; 51 70 ; coding rules 71 ; 72 ; complete description 73 ; 74 ; why different date1 and date2 when readind dt and nrt 75 ; 76 ; why two read_ncdf for nrt 77 ; 78 ; KNOWN ISSUES 79 ; ============ 80 ; 81 ; test of existence of fullfilename_in not very efficient because 82 ; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 52 83 ; 53 84 ; EVOLUTIONS 54 85 ; ========== 55 86 ; 87 ; - fplod 20110809T115747Z aedon.locean-ipsl.upmc.fr (Darwin) 88 ; 89 ; * usage of ${PROJECT_OD} 90 ; * remove v* from filenames (in and out) 91 ; * add test on IO files 92 ; 56 93 ; - fplod 20101215T085754Z aedon.locean-ipsl.upmc.fr (Darwin) 57 94 ; … … 68 105 ;- 69 106 pro TropFlux_swr_BLND_19890101_20091231 70 @common 71 dir="/Users/pkb/data/TropFlux/" 72 73 file=dir+"TropFlux_swr_19890101_20071231_DT_v50.nc" 74 initncdf, file 75 dt=read_ncdf("swr", 19890101, 20071231, file=file,/nostr) 107 ; 108 @cm_4cal 109 @cm_4data 110 @cm_4mesh 111 @cm_4data 112 @cm_project 113 ; 114 ; 115 ; test if ${PROJECT_OD} defined 116 CASE project_od_env OF 117 '' : BEGIN 118 msg = 'eee : ${PROJECT_OD} is not defined' 119 ras = report(msg) 120 STOP 121 END 122 ELSE: BEGIN 123 msg = 'iii : ${PROJECT_OD} is ' + project_od_env 124 ras = report(msg) 125 END 126 ENDCASE 127 ; 128 ; check if output data will be possible 129 iodirout = isadirectory(project_od_env) 130 ; 131 ; existence and protection for reading 132 IF (FILE_TEST(iodirout, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 133 msg = 'eee : the directory' + iodirout + ' is not accessible.' 134 ras = report(msg) 135 STOP 136 ENDIF 137 ; 138 ; existence and protection for writing 139 IF (FILE_TEST(iodirout, /DIRECTORY, /WRITE) EQ 0) THEN BEGIN 140 msg = 'eee : the directory' + iodirout + ' was not found.' 141 ras = report(msg) 142 STOP 143 ENDIF 144 ; 145 ; build dt data filename 146 filename_dt="TropFlux_swr_19890101_20071231_DT.nc" 147 ; 148 ; check if this file exists 149 msg='iii : looking for ' + filename_dt 150 ras = report(msg) 151 fullfilename_dt = isafile(iodirout + filename_dt, NEW=0, /MUST_EXIST) 152 IF fullfilename_dt[0] EQ '' THEN BEGIN 153 msg = 'eee : the file ' + fullfilename_dt + ' was not found.' 154 ras = report(msg) 155 STOP 156 ENDIF 157 158 ; build nrt data filename 159 filename_nrt="TropFlux_swr_19890101_20091231_NRT.nc" 160 ; 161 ; check if this file exists 162 msg='iii : looking for ' + filename_nrt 163 ras = report(msg) 164 fullfilename_nrt = isafile(iodirout + filename_nrt, NEW=0, /MUST_EXIST) 165 IF fullfilename_nrt[0] EQ '' THEN BEGIN 166 msg = 'eee : the file ' + fullfilename_nrt + ' was not found.' 167 ras = report(msg) 168 STOP 169 ENDIF 170 ; 171 ; build wg data filename 172 filename_wg="TropFlux_gustiness_19890101_20091231.nc" 173 ; 174 ; check if this file exists 175 msg='iii : looking for ' + filename_wg 176 ras = report(msg) 177 fullfilename_wg = isafile(iodirout + filename_wg, NEW=0, /MUST_EXIST) 178 IF fullfilename_wg[0] EQ '' THEN BEGIN 179 msg = 'eee : the file ' + fullfilename_wg + ' was not found.' 180 ras = report(msg) 181 STOP 182 ENDIF 183 ; 184 ; build output filename 185 filename_out = 'TropFlux_swr_19890101_20091231_BLND.nc' 186 fullfilename_out = iodirout + filename_out 187 ; in order to avoid unexpected overwritten 188 IF (FILE_TEST(fullfilename_out) EQ 1) THEN BEGIN 189 msg = 'eee : the file ' + fullfilename_out + ' already exists.' 190 ras = report(msg) 191 STOP 192 ENDIF 193 ; 194 initncdf, fullfilename_dt 195 dt=read_ncdf("swr", 19890101, 20071231, file=fullfilename_dt,/nostr) 76 196 help, dt 77 197 78 file=dir+"TropFlux_swr_19890101_20091231_NRT_v50.nc" 79 initncdf, file 80 nrt=read_ncdf("sw", 20080100, 20100112, file=file,/nostr) 198 initncdf, fullfilename_nrt 199 nrt=read_ncdf("sw", 20080100, 20100112, file=fullfilename_nrt,/nostr) 81 200 help, nrt 82 201 83 swr_nrt=read_ncdf("sw", 19880101, 20100112, file=f ile,/nostr)202 swr_nrt=read_ncdf("sw", 19880101, 20100112, file=fullfilename_nrt,/nostr) 84 203 85 204 swr_dt=[[[dt]],[[nrt]]] … … 94 213 for jt=0,jpt-1 do swr_merged(*,*,jt)=swr_dt(*,*,jt)*a(jt)+(1-a(jt))*swr_nrt(*,*,jt) 95 214 96 file=dir+"TropFlux_gustiness_19890101_20091231_v50.nc" 97 initncdf, file 215 initncdf, file_wg 98 216 time=timegen(7670, units='days', start=julday(1,1,1989,0)) & jpt=n_elements(time) 99 217 lat=reform(gphit(0,0:jpj-1)) … … 103 221 104 222 time=time-julday(1,1,1950,00,00) 105 ncfile='! /Users/pkb/data/TropFlux/TropFlux_swr_19890101_20091231_BLND.nc223 ncfile='!' + fullfilename_out 106 224 lon_attr={units:'degrees_east',long_name:'Longitude'} 107 225 lat_attr={units:'degrees_north',long_name:'Latitude'} -
trunk/src/TropFlux_swr_DT_19890101_20071231.pro
r50 r74 10 10 ; .. graphviz:: 11 11 ; 12 ; digraph TropFlux_swr_DT_19890101_20071231 {12 ; digraph tropflux_swr_dt_19890101_20071231 { 13 13 ; graph [ 14 14 ; rankdir="LR", … … 19 19 ; ncfile [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20071231_DT_v50.nc"]; 20 20 ; 21 ; TropFlux_swr_DT_19890101_20071231 [shape=box,21 ; tropflux_swr_dt_19890101_20071231 [shape=box, 22 22 ; fontname=Courier, 23 23 ; color=blue, … … 25 25 ; label="${PROJECT}/src/TropFlux_swr_DT_19890101_20071231.pro"]; 26 26 ; 27 ; {file_swr}-> { TropFlux_swr_DT_19890101_20071231} -> {ncfile}27 ; {file_swr}-> {tropflux_swr_dt_19890101_20071231} -> {ncfile} 28 28 ; 29 29 ; } -
trunk/src/TropFlux_swr_NRT_19890101_20091231.pro
r50 r74 10 10 ; .. graphviz:: 11 11 ; 12 ; digraph TropFlux_swr_NRT_19890101_200912311 {12 ; digraph tropflux_swr_nrt_19890101_200912311 { 13 13 ; graph [ 14 14 ; rankdir="LR", … … 17 17 ; file_olr [shape=ellipse,fontname=Courier,label="/Users/pkb/data/heat_budget/NOAA_OLR/olr_oafluxgrid_30n30s_19890101_20091231.nc"]; 18 18 ; ncfile [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20091231_NRT_v50.nc"]; 19 ; TropFlux_swr_NRT_19890101_200912311 [shape=box, 19 ; 20 ; tropflux_swr_nrt_19890101_200912311 [shape=box, 20 21 ; fontname=Courier, 21 22 ; color=blue, 22 23 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/TropFlux_swr_NRT_19890101_200912311.pro", 23 24 ; label="${PROJECT}/src/TropFlux_swr_NRT_19890101_200912311.pro"]; 24 ; {file_sw file_olr}-> {TropFlux_swr_NRT_19890101_200912311} -> {ncfile} 25 ; 26 ; {file_sw file_olr}-> {tropflux_swr_nrt_19890101_200912311} -> {ncfile} 25 27 ; } 26 28 ; -
trunk/src/cor30a.pro
r50 r74 67 67 ; ======== 68 68 ; 69 ; :ref:`TropFlux_19890101_20091231.pro`69 ; used by :ref:`TropFlux_19890101_20091231.pro` 70 70 ; 71 71 ; TODO 72 72 ; ==== 73 ;74 ; hard coded directory - usage of ${PROJECT_ID}75 73 ; 76 74 ; coding rules
Note: See TracChangeset
for help on using the changeset viewer.