Changeset 73


Ignore:
Timestamp:
08/09/11 12:56:34 (13 years ago)
Author:
pinsard
Message:

try (and failed) correction on climserv

Location:
trunk/src
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/interp_erai_dewt_1989_2009.pro

    r71 r73  
    7070; ==== 
    7171; 
    72 ; make it work : pb on loholt1:: 
     72; make it work : pb on loholt1 idl8:: 
    7373; 
    7474;   d2m[longitude,latitude,time]=d2mout:d2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr  @ globattr 
    … …  
    7676;   % Stop encountered: INTERP_ERAI_DEWT_1989_2009    1 
    7777;      /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 
    7886; 
    7987; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 

  • trunk/src/interp_erai_lwr_1989_2009.pro

    r71 r73  
    7070; ==== 
    7171; 
    72 ; make it work : pb on loholt1:: 
     72; make it work : pb on loholt1 idl8:: 
    7373; 
    7474;  OUTMASK_IND     UNDEFINED = <Undefined> 
    … …  
    8080;     /home/pinsard/tropflux_ws/src/interp_erai_lwr_1989_2009.pro 
    8181; 
     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; 
    8292; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 
    8393; 
    … …  
    101111; 
    102112; $URL$ 
     113; 
     114; - pinsard 2011-08-08T16:03:48Z loholt1.ipsl.polytechnique.fr (Linux) 
     115; 
     116;   * remove return statement 
    103117; 
    104118; - fplod 20110808T094156Z cratos (Linux) 
    … …  
    288302@ncdf_quickwrite 
    289303 
    290 return 
    291304end 

  • trunk/src/interp_erai_sst_1989_2009.pro

    r72 r73  
    7373; on cratos strange red data looking with ncview 
    7474; 
    75 ; make it work : pb on loholt1:: 
     75; make it work : pb on loholt1 idl8:: 
    7676; 
    7777;  sst[longitude,latitude,time]=sstkout:sstk_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr  @ globattr 
    … …  
    7979;  % Stop encountered: INTERP_ERAI_SST_1989_2009    1 
    8080;     /home/pinsard/tropflux_ws/src/interp_erai_sst_1989_2009.pro 
     81; 
     82; but ok on loholt1 idl6 
    8183; 
    8284; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 

  • trunk/src/interp_erai_t2m_1989_2009.pro

    r71 r73  
    7878; ==== 
    7979; 
    80 ; make it work : pb on loholt1:: 
     80; make it work : pb on loholt1 idl8:: 
    8181; 
    8282;  t2m[longitude,latitude,*time]=t2mout:t2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[]=timein:time_attr  @globattr 
    … …  
    8585;     /home/pinsard/tropflux_ws/src/interp_erai_t2m_1989_2009.pro 
    8686; 
     87; loholt1 idl6 : end without error but looking with ncview all brown !  
     88; 
    8789; strange view (lat and lon shift with ncview) : check grid init 
    8890; 
    89 ; remove useless netcdf_read 
     91; remove useless netcdf_read or choose between read_ncdf and ncdf_lec 
    9092; 
    9193; check time values 
    … …  
    108110; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 
    109111; 
    110 ; remove ncdf_quickwrite exercices 
    111 ; 
    112112; KNOWN ISSUES 
    113113; ============ 
    … …  
    122122; 
    123123; $URL$ 
     124; 
     125; - pinsard 2011-08-08T16:11:48Z loholt1.ipsl.polytechnique.fr (Linux) 
     126; 
     127;   * remove ncdf_quickwrite exercices 
    124128; 
    125129; - fplod 20110808T094156Z cratos (Linux) 
    … …  
    349353@ncdf_quickwrite 
    350354 
    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 one 
    359 nlat=73 
    360 nlon=96 
    361 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_quickwrite 
    370  
    371 timepb=findgen(2) 
    372 t2m=findgen(350,60,2) 
    373 help, lat 
    374 help, lon 
    375 help, t2m 
    376 help, time 
    377 help, timein 
    378 help, timepb 
    379 ncfile='!' + fullfilename_out 
    380 ; 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 ; ok 
    384 ;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_quickwrite 
    392  
    393355end 

  • trunk/src/interp_erai_ws_1989_2009.pro

    r71 r73  
    7373; ==== 
    7474; 
    75 ; make it work : pb on loholt1:: 
     75; make it work : pb on loholt1 idl8:: 
    7676; 
    7777;  Writing fields:  
    … …  
    8080;  % Stop encountered: INTERP_ERAI_WS_1989_2009    1 
    8181;     /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 
    8290; 
    8391; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 

  • trunk/src/sst_correction_ncdf.pro

    r72 r73  
    6565; ==== 
    6666; 
     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 
    67118; coding rules 
    68119; 
    … …  
    169220tt=time & jpt=n_elements(time) 
    170221caldat, time,mon,day,yea 
     222;  
     223; debug to understand some time value 9.9692100e+36 on idl 6 lohloht1 
     224print, time[7080:7090] 
     225; 
    171226sst_m=sst*0. 
    172227 
Note: See TracChangeset for help on using the changeset viewer.