![(please configure the [header_logo] section in trac.ini)](/tropflux/chrome/site/your_project_logo.png){kind=logo}
Changeset 92
- Timestamp:
- 08/23/11 14:29:57 (13 years ago)
- Location:
- trunk/src
- Files:
-
- 12 edited
-
d2m_to_q2m_erai.pro (modified) (8 diffs)
-
interp_erai_dewt_1989_2009.pro (modified) (8 diffs)
-
interp_erai_lwr_1989_2009.pro (modified) (8 diffs)
-
interp_erai_msl_1989_2009.pro (modified) (5 diffs)
-
interp_erai_sst_1989_2009.pro (modified) (7 diffs)
-
interp_erai_t2m_1989_2009.pro (modified) (5 diffs)
-
interp_erai_ws_1989_2009.pro (modified) (9 diffs)
-
lwr_correction_ncdf.pro (modified) (2 diffs)
-
q2m_correction_ncdf.pro (modified) (2 diffs)
-
sst_correction_ncdf.pro (modified) (4 diffs)
-
t2m_correction_ncdf.pro (modified) (3 diffs)
-
ws_correction_ncdf.pro (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/d2m_to_q2m_erai.pro
r91 r92 9 9 ; Conversion of d2m to q2m on OAFLUX grid 10 10 ; 11 ; :file:`${PROJECT_OD}/erai_d2m_19890101_20091231_oafluxgrid.nc` have been 12 ; produced by :ref:`interp_erai_dewt_1989_2009.pro`. 13 ; 14 ; :file:`${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc` have been 15 ; produced by :ref:`interp_erai_t2m_1989_2009.pro`. 16 ; 17 ; :file:`${PROJECT_OD}/erai_msl_19890101_20091231_oafluxgrid.nc` have been 11 ; :file:`${PROJECT_OD}/erai_d2m_19890101_20091231_oafluxgrid.nc` 12 ; containing 13 ; ++ 14 ; has been produced by 15 ; :ref:`interp_erai_dewt_1989_2009.pro`. 16 ; 17 ; :file:`${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc` 18 ; containing 19 ; ++ 20 ; has been produced by 21 ; :ref:`interp_erai_t2m_1989_2009.pro`. 22 ; 23 ; :file:`${PROJECT_OD}/erai_msl_19890101_20091231_oafluxgrid.nc` 24 ; containing 25 ; ++ 26 ; has been 18 27 ; produced by :ref:`interp_erai_msl_1989_2009.pro`. 19 28 ; 20 ; q2m on OAFLUX grid is written in 29 ; q2m on OAFLUX grid 30 ; is written in 21 31 ; :file:`${PROJECT_OD}/erai_q2m_19890101_20091231_oafluxgrid.nc` 22 32 ; if this file not already exists. 23 33 ; 24 ; This file will be used by :ref:`q2m_correction_ncdf.pro`. 34 ; This file 35 ; :file:`${PROJECT_OD}/erai_q2m_19890101_20091231_oafluxgrid.nc` 36 ; will be used by :ref:`q2m_correction_ncdf.pro`. 37 ; 38 ; .. only:: man 39 ; 40 ; Figure is visible on PDF and HTML documents only. 41 ; 42 ; .. only:: html or latex 25 43 ; 26 44 ; .. graphviz:: … … 45 63 ; ======== 46 64 ; 65 ; :ref:`mooring_corrections` 66 ; 47 67 ; [Gill:AP:1982]_ 48 68 ; … … 50 70 ; 51 71 ; :ref:`project_profile.sh` 52 ;53 ; :ref:`mooring_corrections`54 72 ; 55 73 ; :ref:`interp_erai_dewt_1989_2009.pro` … … 94 112 ; 95 113 ; $Id$ 96 ; 114 ; 97 115 ; $URL$ 98 116 ; … … 174 192 ; 175 193 ; check if this file exists 194 msg='iii : looking for ' + filename_d2m 195 ras = report(msg) 176 196 fullfilename_d2m = isafile(iodirout + filename_d2m, NEW=0, /MUST_EXIST) 177 197 IF fullfilename_d2m[0] EQ '' THEN BEGIN … … 185 205 ; 186 206 ; check if this file exists 207 msg='iii : looking for ' + filename_t2m 208 ras = report(msg) 187 209 fullfilename_t2m = isafile(iodirout + filename_t2m, NEW=0, /MUST_EXIST) 188 210 IF fullfilename_t2m[0] EQ '' THEN BEGIN … … 196 218 ; 197 219 ; check if this file exists 220 msg='iii : looking for ' + filename_msl 221 ras = report(msg) 198 222 fullfilename_msl = isafile(iodirout + filename_msl, NEW=0, /MUST_EXIST) 199 223 IF fullfilename_msl[0] EQ '' THEN BEGIN … … 239 263 P=read_ncdf("msl", d1, d2, file=fullfilename_msl,/nostr) 240 264 help, P 241 ; ++ for debug ftp error 265 ; ++ for debug ftp error 242 266 ; ++print,P(0:1,0,0) 243 267 ; ++stop -
trunk/src/interp_erai_dewt_1989_2009.pro
r91 r92 12 12 ; Interpolation of d2 from ERA-I grid to OAFLUX grid 13 13 ; 14 ; :file:`${PROJECT_ID}/20c3m_erai_d2_TROP_1989_2009.nc` containing d2 from ERA-I have been produced 15 ; by :ref:`compute_erai_daily_region_2d.sh`. 16 ; 17 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` containing OAFLUX grid have been produced by :ref:`oaflux_mask_30N30S.pro`. 18 ; 19 ; Interpolated d2 is written in 20 ; :file:`${PROJECT_OD}/erai_d2m_19890101_20091231_oafluxgrid.nc` if this file not already exists. 21 ; 22 ; This output file :file:`${PROJECT_OD}/erai_d2m_19890101_20091231_oafluxgrid.nc` must be processed after by :ref:`d2m_to_q2m_erai.pro`. 14 ; :file:`${PROJECT_ID}/20c3m_erai_d2_TROP_1989_2009.nc` 15 ; containing 16 ; d2 from ERA-I 17 ; have been produced by 18 ; :ref:`compute_erai_daily_region_2d.sh`. 19 ; 20 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` 21 ; containing OAFLUX grid 22 ; have been produced by 23 ; :ref:`oaflux_mask_30N30S.pro`. 24 ; 25 ; Interpolated d2 26 ; is written in 27 ; :file:`${PROJECT_OD}/erai_d2m_19890101_20091231_oafluxgrid.nc` 28 ; if this file not already exists. 29 ; 30 ; This output file 31 ; :file:`${PROJECT_OD}/erai_d2m_19890101_20091231_oafluxgrid.nc` 32 ; must be processed after by 33 ; :ref:`d2m_to_q2m_erai.pro`. 34 ; 35 ; .. only:: man 36 ; 37 ; Figure is visible on PDF and HTML documents only. 38 ; 39 ; .. only:: html or latex 23 40 ; 24 41 ; .. graphviz:: … … 42 59 ; SEE ALSO 43 60 ; ======== 61 ; 62 ; :ref:`interpolate_data` 44 63 ; 45 64 ; :ref:`project_profile.sh` … … 69 88 ; make it work : pb on loholt1 idl8:: 70 89 ; 90 ; reach end but bad values (all NaN) :: 91 ; $ ncks -v d2m -d time,0,1 -d latitude,0,1 -d longitude,0,1 \ 92 ; $PROJECT_ID/erai_d2m_19890101_200 1231_oafluxgrid.nc 93 ; time[0]=14245.5 latitude[0]=-29.5 longitude[0]=30.5 d2m[0]=9.96921e+36 degK 94 ; time[0]=14245.5 latitude[0]=-29.5 longitude[1]=31.5 d2m[1]=9.96921e+36 degK 95 ; time[0]=14245.5 latitude[1]=-28.5 longitude[0]=30.5 d2m[350]=9.96921e+36 degK 96 ; time[0]=14245.5 latitude[1]=-28.5 longitude[1]=31.5 d2m[351]=9.96921e+36 degK 97 ; time[1]=14246.5 latitude[0]=-29.5 longitude[0]=30.5 d2m[21000]=9.96921e+36 degK 98 ; time[1]=14246.5 latitude[0]=-29.5 longitude[1]=31.5 d2m[21001]=9.96921e+36 degK 99 ; time[1]=14246.5 latitude[1]=-28.5 longitude[0]=30.5 d2m[21350]=9.96921e+36 degK 100 ; time[1]=14246.5 latitude[1]=-28.5 longitude[1]=31.5 d2m[21351]=9.96921e+36 degK 101 102 ; latitude[0]=-29.5 degrees_north 103 ; latitude[1]=-28.5 degrees_north 104 ; 105 ; longitude[0]=30.5 degrees_east 106 ; longitude[1]=31.5 degrees_east 107 108 ; time[0]=14245.5 days since 1950-01-01 00:00:00 109 ; time[1]=14246.5 days since 1950-01-01 00:00:00 110 ; 71 111 ; d2m[longitude,latitude,time]=d2mout:d2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr @ globattr 72 112 ; % NCDF_DIMDEF: 0 is not a valid cdfid. … … 82 122 ; /home/pinsard/tropflux_ws/src/interp_erai_dewt_1989_2009.pro 83 123 ; 124 ; no such trouble on idl64 on cratos 125 ; 84 126 ; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 85 127 ; … … 100 142 ; 101 143 ; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_dewt_1989_2009.pro $ 144 ; 145 ; - pinsard 2011-08-23T09:09:16Z loholt1.ipsl.polytechnique.fr (Linux) 146 ; 147 ; * retry on loholt1 with idl8 and idl6 102 148 ; 103 149 ; - fplod 20110808T094156Z cratos (Linux) … … 162 208 ; 163 209 ; check if this file exists 210 msg='iii : looking for ' + filename_msk 211 ras = report(msg) 164 212 fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 165 213 IF fullfilename_msk[0] EQ '' THEN BEGIN … … 173 221 ; 174 222 ; check if this file exists 223 msg='iii : looking for ' + filename 224 ras = report(msg) 175 225 fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST) 176 226 IF fullfilename[0] EQ '' THEN BEGIN … … 236 286 d2mout=fltarr(jpi,jpj,jptin) 237 287 for jt=0,jptin-1 do begin 238 print, 'Interpolation jt=',jt,' / ',jptin-1288 ; ++ print, 'Interpolation jt=',jt,' / ',jptin-1 239 289 tab=reform(d2min(*,*,jt)) 240 290 d2mout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $ -
trunk/src/interp_erai_lwr_1989_2009.pro
r91 r92 29 29 ; :file:`${PROJECT_OD}/erai_lwr_19890102_20091231_oafluxgrid.nc` 30 30 ; must be processed after by :ref:`add_198900101.sh`. 31 ; 32 ; .. only:: man 33 ; 34 ; Figure is visible on PDF and HTML documents only. 35 ; 36 ; .. only:: html or latex 31 37 ; 32 38 ; .. graphviz:: … … 50 56 ; SEE ALSO 51 57 ; ======== 58 ; 59 ; :ref:`interpolate_data` 52 60 ; 53 61 ; :ref:`project_profile.sh` … … 77 85 ; make it work : pb on loholt1 idl8:: 78 86 ; 79 ; OUTMASK_IND UNDEFINED = <Undefined>80 ; SET_OUTMSKVAL UNDEFINED = <Undefined>81 ;82 87 ; lwr[longitude,latitude,time]=lwrout:lwr_attr; longitude[]=xlon:lon_attr; latitude[]=ylat:lat_attr; tt[*time]=tt:time_attr @ globattr 83 88 ; % NCDF_DIMDEF: 0 is not a valid cdfid. … … 87 92 ; make it work : pb on loholt1 idl6:: 88 93 ; 89 ; OUTMASK_IND UNDEFINED = <Undefined>90 ; SET_OUTMSKVAL UNDEFINED = <Undefined>91 94 ; lwr[longitude,latitude,time]=lwrout:lwr_attr; longitude[]=xlon:lon_attr; latitude[]=ylat:lat_attr; tt[*time]=tt:time_attr @ globattr 92 95 ; % NCDF_VARPUT: Operation Failed, bad file (4) or variable (0) id ? … … 118 121 ; 119 122 ; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_lwr_1989_2009.pro $ 123 ; 124 ; - pinsard 2011-08-23T09:09:16Z loholt1.ipsl.polytechnique.fr (Linux) 125 ; 126 ; * retry on loholt1 with idl8 and idl6 120 127 ; 121 128 ; - fplod 20110819T084310Z aedon.locean-ipsl.upmc.fr (Darwin) … … 195 202 ; 196 203 ; check if this file exists 204 msg='iii : looking for ' + filename_msk 205 ras = report(msg) 197 206 fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 198 207 IF fullfilename_msk[0] EQ '' THEN BEGIN … … 208 217 ; 209 218 ; check if this file exists 219 msg='iii : looking for ' + filename 220 ras = report(msg) 210 221 fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST) 211 222 IF fullfilename[0] EQ '' THEN BEGIN … … 273 284 274 285 for jt=0,jptin-1 do begin 275 print, 'Interpolation jt=',jt,' / ',jptin-1286 ; ++print, 'Interpolation jt=',jt,' / ',jptin-1 276 287 277 288 tab=reform(lwrin(*,*,jt)) -
trunk/src/interp_erai_msl_1989_2009.pro
r91 r92 12 12 ; Interpolation of msl from ERA-I grid to OAFLUX grid 13 13 ; 14 ; :file:`${PROJECT_ID}/20c3m_erai_msl_TROP_1989_2009.nc` containing msl from ERA-I have been produced 15 ; by :ref:`compute_erai_daily_region_2d.sh`. 16 ; 17 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` containing OAFLUX grid have been produced by :ref:`oaflux_mask_30N30S.pro`. 18 ; 19 ; Interpolated msl is written in 20 ; :file:`${PROJECT_OD}/erai_msl_19890101_20091231_oafluxgrid.nc` if this file not already exists. 21 ; 22 ; This output file :file:`${PROJECT_OD}/erai_msl_19890101_20091231_oafluxgrid.nc` must be processed after by :ref:`d2m_to_q2m_erai.pro`. 23 ; 14 ; :file:`${PROJECT_ID}/20c3m_erai_msl_TROP_1989_2009.nc` 15 ; containing 16 ; msl from ERA-I 17 ; have been produced by 18 ; :ref:`compute_erai_daily_region_2d.sh`. 19 ; 20 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` 21 ; containing 22 ; OAFLUX grid 23 ; have been produced by 24 ; :ref:`oaflux_mask_30N30S.pro`. 25 ; 26 ; Interpolated msl 27 ; is written in 28 ; :file:`${PROJECT_OD}/erai_msl_19890101_20091231_oafluxgrid.nc` 29 ; if this file not already exists. 30 ; 31 ; This output file 32 ; :file:`${PROJECT_OD}/erai_msl_19890101_20091231_oafluxgrid.nc` 33 ; must be processed after by 34 ; :ref:`d2m_to_q2m_erai.pro`. 35 ; 36 ; .. only:: man 37 ; 38 ; Figure is visible on PDF and HTML documents only. 39 ; 40 ; .. only:: html or latex 24 41 ; 25 42 ; .. graphviz:: … … 44 61 ; ======== 45 62 ; 63 ; :ref:`interpolate_data` 64 ; 46 65 ; :ref:`project_profile.sh` 47 ;48 ; :ref:`interpolate_data`49 66 ; 50 67 ; :func:`report <saxo:report>` … … 70 87 ; ==== 71 88 ; 89 ; make it work loholt1 idl8:: 90 ; 91 ; msl[longitude,latitude,time]=mslout:msl_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr @ globattr 92 ; % NCDF_VARPUT: Operation Failed, bad file (6) or variable (0) id ? 93 ; (NC_ERROR=-31) 94 ; % Stop encountered: INTERP_ERAI_MSL_1989_2009 1 95 ; /home/pinsard/tropflux_ws/src/interp_erai_msl_1989_2009.pro 96 ; 97 ; use real output of :ref:`compute_erai_daily_region_2d.sh`. 98 ; 72 99 ; msl_attr={units:'milibars' .. millibar at least, hPa will be probably 73 100 ; more SI … … 89 116 ; 90 117 ; $URL$ 118 ; 119 ; - pinsard 2011-08-23T09:09:16Z loholt1.ipsl.polytechnique.fr (Linux) 120 ; 121 ; * retry on loholt1 with idl8 and idl6 91 122 ; 92 123 ; - fplod 20110811T113646Z aedon.locean-ipsl.upmc.fr (Darwin) … … 235 266 mslout=fltarr(jpi,jpj,jptin) 236 267 for jt=0,jptin-1 do begin 237 print, 'Interpolation jt=',jt,' / ',jptin-1268 ; ++ print, 'Interpolation jt=',jt,' / ',jptin-1 238 269 tab=reform(mslin(*,*,jt)) 239 270 mslout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $ -
trunk/src/interp_erai_sst_1989_2009.pro
r91 r92 12 12 ; Interpolation of sst from ERA-I grid to OAFLUX grid 13 13 ; 14 ; :file:`${PROJECT_ID}/20c3m_erai_sstk_TROP_1989_2009.nc` containing sst from ERA-I have been produced 15 ; by :ref:`compute_erai_daily_region_2d.sh`. 16 ; 17 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` containing OAFLUX grid have been produced by :ref:`oaflux_mask_30N30S.pro`. 18 ; 19 ; Interpolated sst is written in 20 ; :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc` if this file not already exists. 21 ; 22 ; This output file :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc` must be processed after by :ref:`sst_correction_ncdf.pro`. 23 ; 14 ; :file:`${PROJECT_ID}/20c3m_erai_sstk_TROP_1989_2009.nc` 15 ; containing sst from ERA-I 16 ; have been produced by 17 ; :ref:`compute_erai_daily_region_2d.sh`. 18 ; 19 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` 20 ; containing 21 ; OAFLUX grid 22 ; have been produced by 23 ; :ref:`oaflux_mask_30N30S.pro`. 24 ; 25 ; Interpolated sst 26 ; is written in 27 ; :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc` 28 ; if this file not already exists. 29 ; 30 ; This output file 31 ; :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc` 32 ; must be processed after by 33 ; :ref:`sst_correction_ncdf.pro`. 34 ; 35 ; .. only:: man 36 ; 37 ; Figure is visible on PDF and HTML documents only. 38 ; 39 ; .. only:: html or latex 24 40 ; 25 41 ; .. graphviz:: … … 43 59 ; SEE ALSO 44 60 ; ======== 61 ; 62 ; :ref:`interpolate_data` 45 63 ; 46 64 ; :ref:`project_profile.sh` … … 72 90 ; make it work : pb on loholt1 idl8:: 73 91 ; 74 ; sst[longitude,latitude,time]=sstkout:sstk_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr @ globattr 75 ; % NCDF_DIMDEF: 0 is not a valid cdfid. 76 ; % Stop encountered: INTERP_ERAI_SST_1989_2009 1 77 ; /home/pinsard/tropflux_ws/src/interp_erai_sst_1989_2009.pro 92 ; sst[longitude,latitude,time]=sstkout:sstk_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr @ globattr 93 ; % NCDF_VARPUT: Operation Failed, bad file (7) or variable (0) id ? 94 ; (NC_ERROR=-31) 95 ; % Stop encountered: INTERP_ERAI_SST_1989_2009 1 96 ; /home/pinsard/tropflux_ws/src/interp_erai_sst_1989_2009.pro 78 97 ; 79 98 ; but ok on loholt1 idl6 … … 97 116 ; 98 117 ; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_sst_1989_2009.pro $ 118 ; 119 ; - pinsard 2011-08-23T09:09:16Z loholt1.ipsl.polytechnique.fr (Linux) 120 ; 121 ; * retry on loholt1 with idl8 and idl6 99 122 ; 100 123 ; - fplod 20110808T094156Z cratos (Linux) … … 159 182 ; 160 183 ; check if this file exists 184 msg='iii : looking for ' + filename_msk 185 ras = report(msg) 161 186 fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 162 187 IF fullfilename_msk[0] EQ '' THEN BEGIN … … 170 195 ; 171 196 ; check if this file exists 197 msg='iii : looking for ' + filename 198 ras = report(msg) 172 199 fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST) 173 200 IF fullfilename[0] EQ '' THEN BEGIN … … 234 261 sstkout=fltarr(jpi,jpj,jptin) 235 262 for jt=0,jptin-1 do begin 236 print, 'Interpolation jt=',jt,' / ',jptin-1263 ; ++print, 'Interpolation jt=',jt,' / ',jptin-1 237 264 tab=reform(sstkin(*,*,jt)) 238 265 sstkout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $ -
trunk/src/interp_erai_t2m_1989_2009.pro
r91 r92 12 12 ; Interpolation of t2 from ERA-I grid to OAFLUX grid 13 13 ; 14 ; :file:`${PROJECT_ID}/20c3m_erai_t2_TROP_1989_2009.nc` containing t2 from ERA-I have been produced 15 ; by :ref:`compute_erai_daily_region_2d.sh`. 16 ; 17 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` containing OAFLUX grid have been produced by :ref:`oaflux_mask_30N30S.pro`. 18 ; 19 ; Interpolated t2 is written in 20 ; :file:`${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc` if this file not already exists. 21 ; 22 ; This output file :file:`${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc` must be processed after by :ref:`t2m_correction_ncdf.pro`. 14 ; :file:`${PROJECT_ID}/20c3m_erai_t2_TROP_1989_2009.nc` 15 ; containing 16 ; t2 from ERA-I 17 ; have been produced by 18 ; :ref:`compute_erai_daily_region_2d.sh`. 19 ; 20 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` 21 ; containing OAFLUX grid 22 ; have been produced by 23 ; :ref:`oaflux_mask_30N30S.pro`. 24 ; 25 ; Interpolated t2 26 ; is written in 27 ; :file:`${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc` 28 ; if this file not already exists. 29 ; 30 ; This output file 31 ; :file:`${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc` 32 ; must be processed after by 33 ; :ref:`t2m_correction_ncdf.pro`. 34 ; 35 ; .. only:: man 36 ; 37 ; Figure is visible on PDF and HTML documents only. 38 ; 39 ; .. only:: html or latex 23 40 ; 24 41 ; .. graphviz:: … … 77 94 ; make it work : pb on loholt1 idl8:: 78 95 ; 79 ; t2m[longitude,latitude,*time]=t2mout:t2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[]=timein:time_attr @globattr 80 ; % NCDF_DIMDEF: 0 is not a valid cdfid. 81 ; % Stop encountered: INTERP_ERAI_T2M_1989_2009 1 82 ; /home/pinsard/tropflux_ws/src/interp_erai_t2m_1989_2009.pro 83 ; 84 ; loholt1 idl6 : end without error but looking with ncview all brown ! 96 ; t2m[longitude,latitude,*time]=t2mout:t2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[]=timein:time_attr @globattr 97 ; % NCDF_VARPUT: Operation Failed, bad file (8) or variable (0) id ? 98 ; (NC_ERROR=-31) 99 ; % Stop encountered: INTERP_ERAI_T2M_1989_2009 1 100 ; /home/pinsard/tropflux_ws/src/interp_erai_t2m_1989_2009.pro 101 ; 102 ; loholt1 idl6 : end without error but looking with ncview all brown ! 85 103 ; 86 104 ; strange view (lat and lon shift with ncview) : check grid init … … 119 137 ; 120 138 ; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_t2m_1989_2009.pro $ 139 ; 140 ; - pinsard 2011-08-23T09:09:16Z loholt1.ipsl.polytechnique.fr (Linux) 141 ; 142 ; * retry on loholt1 with idl8 and idl6 121 143 ; 122 144 ; - pinsard 2011-08-08T16:11:48Z loholt1.ipsl.polytechnique.fr (Linux) … … 215 237 ; 216 238 ; check if this file exists 239 msg='iii : looking for ' + filename_msk 240 ras = report(msg) 217 241 fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 218 242 IF fullfilename_msk[0] EQ '' THEN BEGIN … … 226 250 ; 227 251 ; check if this file exists 252 msg='iii : looking for ' + filename_t2 253 ras = report(msg) 228 254 fullfilename_t2 = isafile(iodirin + filename_t2, NEW=0, /MUST_EXIST) 229 255 IF fullfilename_t2[0] EQ '' THEN BEGIN -
trunk/src/interp_erai_ws_1989_2009.pro
r91 r92 12 12 ; Interpolation of u10 and v10 from ERA-I grid to OAFLUX grid 13 13 ; 14 ; :file:`${PROJECT_ID}/20c3m_erai_u10_TROP_1989_2009.nc` containing u10 from ERA-I have been produced 15 ; by :ref:`compute_erai_daily_region_2d.sh`. 16 ; 17 ; :file:`${PROJECT_ID}/20c3m_erai_v10_TROP_1989_2009.nc` containing u10 from ERA-I have been produced 18 ; 19 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` containing OAFLUX grid have been produced by :ref:`oaflux_mask_30N30S.pro`. 20 ; 21 ; Interpolated u10 and v10 is written in 22 ; :file:`${PROJECT_OD}/erai_ws_19890101_20091231_oafluxgrid.nc` if this file not already exists. 23 ; 24 ; This output file :file:`${PROJECT_OD}/erai_ws_19890101_20091231_oafluxgrid.nc` must be processed after by :ref:`ws_correction_ncdf.pro`. 14 ; :file:`${PROJECT_ID}/20c3m_erai_u10_TROP_1989_2009.nc` 15 ; containing 16 ; u10 from ERA-I 17 ; have been produced by 18 ; :ref:`compute_erai_daily_region_2d.sh`. 19 ; 20 ; :file:`${PROJECT_ID}/20c3m_erai_v10_TROP_1989_2009.nc` 21 ; containing 22 ; v10 from ERA-I 23 ; have been produced by 24 ; :ref:`compute_erai_daily_region_2d.sh`. 25 ; 26 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` 27 ; containing OAFLUX grid 28 ; have been produced by 29 ; :ref:`oaflux_mask_30N30S.pro`. 30 ; 31 ; Interpolated u10 and v10 are written in 32 ; :file:`${PROJECT_OD}/erai_ws_19890101_20091231_oafluxgrid.nc` 33 ; if this file not already exists. 34 ; 35 ; This output file 36 ; :file:`${PROJECT_OD}/erai_ws_19890101_20091231_oafluxgrid.nc` 37 ; must be processed after by 38 ; :ref:`ws_correction_ncdf.pro`. 39 ; 40 ; .. only:: man 41 ; 42 ; Figure is visible on PDF and HTML documents only. 43 ; 44 ; .. only:: html or latex 25 45 ; 26 46 ; .. graphviz:: … … 45 65 ; SEE ALSO 46 66 ; ======== 67 ; 68 ; :ref:`interpolate_data` 47 69 ; 48 70 ; :ref:`project_profile.sh` … … 72 94 ; make it work : pb on loholt1 idl8:: 73 95 ; 74 ; Writing fields: 75 ; u10[longitude,latitude,time]=u10out:u10_attr; v10[longitude,latitude,time]=v10out:v10_attr; longitude[]=xlon:lon_attr; latitude[]=ylat:lat_attr; tt[*time]=tt:time_attr @ globattr 76 ; % NCDF_DIMDEF: 0 is not a valid cdfid. 77 ; % Stop encountered: INTERP_ERAI_WS_1989_2009 1 78 ; /home/pinsard/tropflux_ws/src/interp_erai_ws_1989_2009.pro 96 ; Writing fields: 97 ; u10[longitude,latitude,time]=u10out:u10_attr; v10[longitude,latitude,time]=v10out:v10_attr; longitude[]=xlon:lon_attr; latitude[]=ylat:lat_attr; tt[*time]=tt:time_attr @ globattr 98 ; % NCDF_VARPUT: Operation Failed, bad file (9) or variable (0) id ? 99 ; (NC_ERROR=-31) 100 ; % Stop encountered: INTERP_ERAI_WS_1989_2009 1 101 ; /home/pinsard/tropflux_ws/src/interp_erai_ws_1989_2009.pro 79 102 ; 80 103 ; make it work : pb on loholt1 idl6:: 81 104 ; 82 ; Writing fields: 105 ; Writing fields: 83 106 ; u10[longitude,latitude,time]=u10out:u10_attr; v10[longitude,latitude,time]=v10out:v10_attr; longitude[]=xlon:lon_attr; latitude[]=ylat:lat_attr; tt[*time]=tt:time_attr @ globattr 84 107 ; % NCDF_VARPUT: Operation Failed, bad file (6) or variable (0) id ? … … 104 127 ; 105 128 ; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_ws_1989_2009.pro $ 129 ; 130 ; - pinsard 2011-08-23T09:09:16Z loholt1.ipsl.polytechnique.fr (Linux) 131 ; 132 ; * retry on loholt1 with idl8 and idl6 106 133 ; 107 134 ; - fplod 20110808T094156Z cratos (Linux) … … 166 193 ; 167 194 ; check if this file exists 195 msg='iii : looking for ' + filename_msk 196 ras = report(msg) 168 197 fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 169 198 IF fullfilename_msk[0] EQ '' THEN BEGIN … … 179 208 ; 180 209 ; check if this file exists 210 msg='iii : looking for ' + filename 211 ras = report(msg) 181 212 fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST) 182 213 IF fullfilename[0] EQ '' THEN BEGIN … … 226 257 ; 227 258 ; check if this file exists 259 msg='iii : looking for ' + filename 260 ras = report(msg) 228 261 fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST) 229 262 IF fullfilename[0] EQ '' THEN BEGIN … … 259 292 260 293 for jt=0,jptin-1 do begin 261 print, 'Interpolation jt=',jt,' / ',jptin-1294 ; ++print, 'Interpolation jt=',jt,' / ',jptin-1 262 295 tab=reform(u10in(*,*,jt)) 263 296 u10out(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $ … … 293 326 u10_attr={units:'m/s',missing_value:valmask,long_name:'10 metre u wind component',short_name:'u10',axis:'TYX'} 294 327 v10_attr={units:'m/s',missing_value:valmask,long_name:'10 metre v wind component',short_name:'v10',axis:'TYX'} 295 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'}296 328 globattr={source:'Tropical ocean winds obtained from ERA Interim',timerange:cda0+' - '+cda1} 297 329 -
trunk/src/lwr_correction_ncdf.pro
r90 r92 9 9 ; Correction of lwr on OAFLUX grid 10 10 ; 11 ; :file:`${PROJECT_OD}/erai_lwr_19890101_20091231_oafluxgrid.nc` 11 ; :file:`${PROJECT_OD}/erai_lwr_19890101_20091231_oafluxgrid.nc` 12 12 ; containing 13 13 ; ++ 14 ; have been produced by 14 ; have been produced by 15 15 ; :ref:`add_19890101.sh`. 16 16 ; 17 ; Corrected lwr on OAFLUX grid 17 ; Corrected lwr on OAFLUX grid 18 18 ; is written in 19 19 ; :file:`${PROJECT_OD}/TropFlux_lwr_19890101_20091231.nc` 20 20 ; if this file not already exists. 21 21 ; 22 ; This file will be used by :ref:`TropFlux_19890101_20091231.pro`. 22 ; This output file 23 ; :file:`${PROJECT_OD}/TropFlux_lwr_19890101_20091231.nc` 24 ; will be used by :ref:`TropFlux_19890101_20091231.pro`. 25 ; 26 ; .. only:: man 27 ; 28 ; Figure is visible on PDF and HTML documents only. 29 ; 30 ; .. only:: html or latex 23 31 ; 24 32 ; .. graphviz:: … … 156 164 ; 157 165 ; check if this file exists 166 msg='iii : looking for ' + filename 167 ras = report(msg) 158 168 fullfilename = isafile(iodirout + filename, NEW=0, /MUST_EXIST) 159 169 IF fullfilename[0] EQ '' THEN BEGIN -
trunk/src/q2m_correction_ncdf.pro
r90 r92 12 12 ; Correction of q2m on OAFLUX grid 13 13 ; 14 ; :file:`${PROJECT_OD}/erai_q2m_19890101_20091231_oafluxgrid.nc` have been 15 ; produced by :ref:`d2m_to_q2m_erai.pro`. 16 ; 17 ; Corrected q2m on OAFLUX grid is written in 14 ; :file:`${PROJECT_OD}/erai_q2m_19890101_20091231_oafluxgrid.nc` 15 ; containing 16 ; ++ 17 ; has been produced by 18 ; :ref:`d2m_to_q2m_erai.pro`. 19 ; 20 ; Corrected q2m on OAFLUX grid 21 ; is written in 18 22 ; :file:`${PROJECT_OD}/TropFlux_q2m_19890101_20091231.nc` 19 23 ; if this file not already exists. 20 24 ; 21 25 ; This file will be used by :ref:`TropFlux_19890101_20091231.pro`. 26 ; 27 ; .. only:: man 28 ; 29 ; Figure is visible on PDF and HTML documents only. 30 ; 31 ; .. only:: html or latex 22 32 ; 23 33 ; .. graphviz:: … … 65 75 ; work on cratos idl7 even if NaNf values in erai_q2m_19890101_20091231_oafluxgrid.nc written by 66 76 ; d2m_to_q2m.pro. 67 ; 77 ; 68 78 ; No way ... NaNf also in this output !! 69 79 ; -
trunk/src/sst_correction_ncdf.pro
r91 r92 9 9 ; Correction of sst on OAFLUX grid 10 10 ; 11 ; :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc` have been 12 ; produced by :ref:`interp_erai_sst_1989_2009.pro`. 13 ; 14 ; Corrected sst on OAFLUX grid is written in 11 ; :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc` 12 ; containing 13 ; ++ 14 ; has been produced by 15 ; :ref:`interp_erai_sst_1989_2009.pro`. 16 ; 17 ; Corrected sst on OAFLUX grid 18 ; is written in 15 19 ; :file:`${PROJECT_OD}/TropFlux_sst_19890101_20091231.nc` 16 20 ; if this file not already exists. 17 21 ; 18 ; This file will be used by :ref:`cronin_gustiness_ncdf.pro` and 22 ; This output file 23 ; :file:`${PROJECT_OD}/TropFlux_sst_19890101_20091231.nc` 24 ; will be used by :ref:`cronin_gustiness_ncdf.pro` 25 ; and 19 26 ; :ref:`TropFlux_19890101_20091231.pro`. 27 ; 28 ; .. only:: man 29 ; 30 ; Figure is visible on PDF and HTML documents only. 31 ; 32 ; .. only:: html or latex 20 33 ; 21 34 ; .. graphviz:: … … 68 81 ; 9.9692100e+36 9.9692100e+36 9.9692100e+36 9.9692100e+36 2454616.5 69 82 ; 2454617.5 70 ; 83 ; 71 84 ; while :: 72 85 ; 73 86 ; $ ncks -H -C -v time -d time,7080,7090 ${PROJECT_ID}/20c3m_erai_sstk_TROP_1989_2009.nc 74 ; time[7080]=450432 75 ; time[7081]=450456 76 ; time[7082]=450480 77 ; time[7083]=450504 78 ; time[7084]=450528 79 ; time[7085]=450552 80 ; time[7086]=450576 81 ; time[7087]=450600 82 ; time[7088]=450624 83 ; time[7089]=450648 84 ; time[7090]=450672 85 ; 86 ; and 87 ; ncks -H -C -v time -d time,7080,7090 ${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc 88 ; time[7080]=21325 89 ; time[7081]=21326 90 ; time[7082]=21327 91 ; time[7083]=9.96920996839e+36 92 ; time[7084]=9.96920996839e+36 93 ; time[7085]=9.96920996839e+36 94 ; time[7086]=9.96920996839e+36 95 ; time[7087]=9.96920996839e+36 96 ; time[7088]=9.96920996839e+36 97 ; time[7089]=21334 98 ; time[7090]=21335 99 ; 87 ; time[7080]=450432 88 ; time[7081]=450456 89 ; time[7082]=450480 90 ; time[7083]=450504 91 ; time[7084]=450528 92 ; time[7085]=450552 93 ; time[7086]=450576 94 ; time[7087]=450600 95 ; time[7088]=450624 96 ; time[7089]=450648 97 ; time[7090]=450672 98 ; 99 ; and 100 ; ncks -H -C -v time -d time,7080,7090 ${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc 101 ; time[7080]=21325 102 ; time[7081]=21326 103 ; time[7082]=21327 104 ; time[7083]=9.96920996839e+36 105 ; time[7084]=9.96920996839e+36 106 ; time[7085]=9.96920996839e+36 107 ; time[7086]=9.96920996839e+36 108 ; time[7087]=9.96920996839e+36 109 ; time[7088]=9.96920996839e+36 110 ; time[7089]=21334 111 ; time[7090]=21335 112 ; 100 113 ; and :: 101 114 ; 102 115 ; fplod@cratos$ ncks -H -C -v time -d time,7080,7090 ${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc 103 ; time[7080]=21325 104 ; time[7081]=21326 105 ; time[7082]=21327 106 ; time[7083]=21328 107 ; time[7084]=21329 108 ; time[7085]=21330 109 ; time[7086]=21331 110 ; time[7087]=21332 111 ; time[7088]=21333 112 ; time[7089]=21334 113 ; time[7090]=21335 116 ; time[7080]=21325 117 ; time[7081]=21326 118 ; time[7082]=21327 119 ; time[7083]=21328 120 ; time[7084]=21329 121 ; time[7085]=21330 122 ; time[7086]=21331 123 ; time[7087]=21332 124 ; time[7088]=21333 125 ; time[7089]=21334 126 ; time[7090]=21335 114 127 115 128 ; coding rules … … 195 208 ; 196 209 ; check if this file exists 210 msg='iii : looking for ' + filename 211 ras = report(msg) 197 212 fullfilename = isafile(iodirout + filename, NEW=0, /MUST_EXIST) 198 213 IF fullfilename[0] EQ '' THEN BEGIN … … 221 236 tt=time & jpt=n_elements(time) 222 237 caldat, time,mon,day,yea 223 ; 238 ; 224 239 ; debug to understand some time value 9.9692100e+36 on idl 6 lohloht1 225 240 print, time[7080:7090] -
trunk/src/t2m_correction_ncdf.pro
r90 r92 11 11 ; 12 12 ; :file:`${PROJECT_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc` 13 ; containing 14 ; air temperature at 2 m height from ERA-I interpolated on OAFLUX grid 13 15 ; has been produced by 14 16 ; :ref:`interp_erai_t2m_1989_2009.pro`. 15 ; It contains air temperature at 2 m height from ERA-I interpolated on OAFLUX grid.16 17 ; 17 18 ; Corrected air temperature at 2 m height … … 20 21 ; if this file not already exists. 21 22 ; 22 ; This file will be used by :ref:`TropFlux_19890101_20091231.pro`. 23 ; This output file 24 ; :file:`${PROJECT_OD}/TropFlux_t2m_19890101_20091231.nc` 25 ; will be used by :ref:`TropFlux_19890101_20091231.pro`. 26 ; 27 ; .. only:: man 28 ; 29 ; Figure is visible on PDF and HTML documents only. 30 ; 31 ; .. only:: html or latex 23 32 ; 24 33 ; .. graphviz:: … … 195 204 ; 196 205 ; check if this file exists 206 msg='iii : looking for ' + filename_t2_uncor 207 ras = report(msg) 197 208 fullfilename_t2_uncor = isafile(iodirout + filename_t2_uncor, NEW=0, /MUST_EXIST) 198 209 IF fullfilename_t2_uncor[0] EQ '' THEN BEGIN -
trunk/src/ws_correction_ncdf.pro
r91 r92 9 9 ; Correction of ws on OAFLUX grid 10 10 ; 11 ; :file:`${PROJECT_OD}/erai_ws_19890101_20091231_oafluxgrid.nc` have been 12 ; produced by :ref:`interp_erai_ws_1989_2009.pro`. 11 ; :file:`${PROJECT_OD}/erai_ws_19890101_20091231_oafluxgrid.nc` 12 ; containing 13 ; ++ 14 ; have been produced by 15 ; :ref:`interp_erai_ws_1989_2009.pro`. 13 16 ; 14 17 ; Corrected ws on OAFLUX grid is written in … … 16 19 ; if this file not already exists. 17 20 ; 18 ; This file will be used by :ref:`TropFlux_19890101_20091231.pro`. 21 ; This output file 22 ; :file:`${PROJECT_OD}/TropFlux_ws_19890101_20091231.nc` 23 ; will be used by :ref:`TropFlux_19890101_20091231.pro`. 24 ; 25 ; .. only:: man 26 ; 27 ; Figure is visible on PDF and HTML documents only. 28 ; 29 ; .. only:: html or latex 19 30 ; 20 31 ; .. graphviz:: … … 141 152 ; 142 153 ; check if this file exists 154 msg='iii : looking for ' + filename 155 ras = report(msg) 143 156 fullfilename = isafile(iodirout + filename, NEW=0, /MUST_EXIST) 144 157 IF fullfilename[0] EQ '' THEN BEGIN
Note: See TracChangeset
for help on using the changeset viewer.
