- Timestamp:
- 08/09/11 12:56:34 (13 years ago)
- Location:
- trunk/src
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/interp_erai_dewt_1989_2009.pro
r71 r73 70 70 ; ==== 71 71 ; 72 ; make it work : pb on loholt1 ::72 ; make it work : pb on loholt1 idl8:: 73 73 ; 74 74 ; d2m[longitude,latitude,time]=d2mout:d2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr @ globattr … … 76 76 ; % Stop encountered: INTERP_ERAI_DEWT_1989_2009 1 77 77 ; /home/pinsard/tropflux_ws/src/interp_erai_dewt_1989_2009.pro 78 ; 79 ; make it work : pb on loholt1 idl6:: 80 ; 81 ; d2m[longitude,latitude,time]=d2mout:d2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr @ globattr 82 ; % NCDF_VARPUT: Operation Failed, bad file (4) or variable (0) id ? 83 ; (NC_ERROR=-31) 84 ; % Stop encountered: INTERP_ERAI_DEWT_1989_2009 1 85 ; /home/pinsard/tropflux_ws/src/interp_erai_dewt_1989_2009.pro 78 86 ; 79 87 ; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call -
trunk/src/interp_erai_lwr_1989_2009.pro
r71 r73 70 70 ; ==== 71 71 ; 72 ; make it work : pb on loholt1 ::72 ; make it work : pb on loholt1 idl8:: 73 73 ; 74 74 ; OUTMASK_IND UNDEFINED = <Undefined> … … 80 80 ; /home/pinsard/tropflux_ws/src/interp_erai_lwr_1989_2009.pro 81 81 ; 82 ; make it work : pb on loholt1 idl6:: 83 ; 84 ; OUTMASK_IND UNDEFINED = <Undefined> 85 ; SET_OUTMSKVAL UNDEFINED = <Undefined> 86 ; lwr[longitude,latitude,time]=lwrout:lwr_attr; longitude[]=xlon:lon_attr; latitude[]=ylat:lat_attr; tt[*time]=tt:time_attr @ globattr 87 ; % NCDF_VARPUT: Operation Failed, bad file (4) or variable (0) id ? 88 ; (NC_ERROR=-31) 89 ; % Stop encountered: INTERP_ERAI_LWR_1989_2009 1 90 ; /home/pinsard/tropflux_ws/src/interp_erai_lwr_1989_2009.pro 91 ; 82 92 ; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 83 93 ; … … 101 111 ; 102 112 ; $URL$ 113 ; 114 ; - pinsard 2011-08-08T16:03:48Z loholt1.ipsl.polytechnique.fr (Linux) 115 ; 116 ; * remove return statement 103 117 ; 104 118 ; - fplod 20110808T094156Z cratos (Linux) … … 288 302 @ncdf_quickwrite 289 303 290 return291 304 end -
trunk/src/interp_erai_sst_1989_2009.pro
r72 r73 73 73 ; on cratos strange red data looking with ncview 74 74 ; 75 ; make it work : pb on loholt1 ::75 ; make it work : pb on loholt1 idl8:: 76 76 ; 77 77 ; sst[longitude,latitude,time]=sstkout:sstk_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr @ globattr … … 79 79 ; % Stop encountered: INTERP_ERAI_SST_1989_2009 1 80 80 ; /home/pinsard/tropflux_ws/src/interp_erai_sst_1989_2009.pro 81 ; 82 ; but ok on loholt1 idl6 81 83 ; 82 84 ; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call -
trunk/src/interp_erai_t2m_1989_2009.pro
r71 r73 78 78 ; ==== 79 79 ; 80 ; make it work : pb on loholt1 ::80 ; make it work : pb on loholt1 idl8:: 81 81 ; 82 82 ; t2m[longitude,latitude,*time]=t2mout:t2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[]=timein:time_attr @globattr … … 85 85 ; /home/pinsard/tropflux_ws/src/interp_erai_t2m_1989_2009.pro 86 86 ; 87 ; loholt1 idl6 : end without error but looking with ncview all brown ! 88 ; 87 89 ; strange view (lat and lon shift with ncview) : check grid init 88 90 ; 89 ; remove useless netcdf_read 91 ; remove useless netcdf_read or choose between read_ncdf and ncdf_lec 90 92 ; 91 93 ; check time values … … 108 110 ; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 109 111 ; 110 ; remove ncdf_quickwrite exercices111 ;112 112 ; KNOWN ISSUES 113 113 ; ============ … … 122 122 ; 123 123 ; $URL$ 124 ; 125 ; - pinsard 2011-08-08T16:11:48Z loholt1.ipsl.polytechnique.fr (Linux) 126 ; 127 ; * remove ncdf_quickwrite exercices 124 128 ; 125 129 ; - fplod 20110808T094156Z cratos (Linux) … … 349 353 @ncdf_quickwrite 350 354 351 ; ok ncfields = 'longitude[]=lon:lon_attr;latitude[]=lat:lat_attr'352 ; ok ncfields = 'longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[]=timein:time_attr'353 ; ok ncfields = 'longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[*]=timein:time_attr'354 ; ok ncfields = 'longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[*]=timein:time_attr@globattr'355 ; pas ok ncfields = 't2m[longitude,latitude,*time]=t2mout:t2m_attr;longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[*]=timein:time_attr@globattr'356 ; pas ok ncfields = 'longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[*]=timein:time_attr;t2m[longitude,latitude,*time]=t2mout:t2m_attr@globattr'357 358 ; the simpliest one359 nlat=73360 nlon=96361 latitudetp=findgen(nlat)362 longitudetp=findgen(nlon)363 timetp=findgen(2)364 pressure3d=findgen(nlon,nlat,2)365 patts = {units:'pascal',missing_value:1e10}366 globatts={source:'my program',version:2}367 NCFIELDS = 'pressure3d[longitudetp,latitudetp,*timetp]=pressure3d:patts;longitudetp[];latitudetp[];timetp[]@globatts'368 ncfile='my4.nc'369 @ncdf_quickwrite370 371 timepb=findgen(2)372 t2m=findgen(350,60,2)373 help, lat374 help, lon375 help, t2m376 help, time377 help, timein378 help, timepb379 ncfile='!' + fullfilename_out380 ; pas ok ncfields = 'longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[]=timein:time_attr;t2m[longitude,latitude,*time]'381 ; pas ok ncfields = 't2m[longitude,latitude,*time];longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[]=timein:time_attr'382 ; ok avec timepb=findgen(2) ncfields = 't2m[longitude,latitude,*time];longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[]=timepb:time_attr'383 ; ok384 ;timepb=findgen(2)385 ;timepb[0]=280512.386 ;timepb[1]=280536.387 ;ncfields = 't2m[longitude,latitude,*time];longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[]=timepb:time_attr'388 389 ;++timepb=findgen(7670)390 ;++ncfields = 't2m[longitude,latitude,*time];longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[]=timepb:time_attr'391 ;++@ncdf_quickwrite392 393 355 end -
trunk/src/interp_erai_ws_1989_2009.pro
r71 r73 73 73 ; ==== 74 74 ; 75 ; make it work : pb on loholt1 ::75 ; make it work : pb on loholt1 idl8:: 76 76 ; 77 77 ; Writing fields: … … 80 80 ; % Stop encountered: INTERP_ERAI_WS_1989_2009 1 81 81 ; /home/pinsard/tropflux_ws/src/interp_erai_ws_1989_2009.pro 82 ; 83 ; make it work : pb on loholt1 idl6:: 84 ; 85 ; Writing fields: 86 ; 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 87 ; % NCDF_VARPUT: Operation Failed, bad file (6) or variable (0) id ? 88 ; (NC_ERROR=-31) 89 ; % Stop encountered: INTERP_ERAI_WS_1989_2009 1 82 90 ; 83 91 ; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call -
trunk/src/sst_correction_ncdf.pro
r72 r73 65 65 ; ==== 66 66 ; 67 ; very strange value of time on lohlot1 idl6 but not always and not always at the same indices !!:: 68 ; 69 ; print, time[7080:7090] 70 ; 2454607.5 2454608.5 2454609.5 9.9692100e+36 9.9692100e+36 71 ; 9.9692100e+36 9.9692100e+36 9.9692100e+36 9.9692100e+36 2454616.5 72 ; 2454617.5 73 ; 74 ; while :: 75 ; 76 ; $ ncks -H -C -v time -d time,7080,7090 ${PROJECT_ID}/20c3m_erai_sstk_TROP_1989_2009.nc 77 ; time[7080]=450432 78 ; time[7081]=450456 79 ; time[7082]=450480 80 ; time[7083]=450504 81 ; time[7084]=450528 82 ; time[7085]=450552 83 ; time[7086]=450576 84 ; time[7087]=450600 85 ; time[7088]=450624 86 ; time[7089]=450648 87 ; time[7090]=450672 88 ; 89 ; and 90 ; ncks -H -C -v time -d time,7080,7090 ${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc 91 ; time[7080]=21325 92 ; time[7081]=21326 93 ; time[7082]=21327 94 ; time[7083]=9.96920996839e+36 95 ; time[7084]=9.96920996839e+36 96 ; time[7085]=9.96920996839e+36 97 ; time[7086]=9.96920996839e+36 98 ; time[7087]=9.96920996839e+36 99 ; time[7088]=9.96920996839e+36 100 ; time[7089]=21334 101 ; time[7090]=21335 102 ; 103 ; and :: 104 ; 105 ; fplod@cratos$ ncks -H -C -v time -d time,7080,7090 ${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc 106 ; time[7080]=21325 107 ; time[7081]=21326 108 ; time[7082]=21327 109 ; time[7083]=21328 110 ; time[7084]=21329 111 ; time[7085]=21330 112 ; time[7086]=21331 113 ; time[7087]=21332 114 ; time[7088]=21333 115 ; time[7089]=21334 116 ; time[7090]=21335 117 67 118 ; coding rules 68 119 ; … … 169 220 tt=time & jpt=n_elements(time) 170 221 caldat, time,mon,day,yea 222 ; 223 ; debug to understand some time value 9.9692100e+36 on idl 6 lohloht1 224 print, time[7080:7090] 225 ; 171 226 sst_m=sst*0. 172 227
Note: See TracChangeset
for help on using the changeset viewer.