Changeset 559


Ignore:
Timestamp:
05/22/12 10:01:43 (12 years ago)
Author:
llmlod
Message:

fic orbites separees

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/create_amsu_netcdf.pro

    r513 r559  
    8686;- 
    8787pro create_amsu_netcdf $ 
    88     , file_ncdf, yyyy, mm, dd, hh, mn, ss, lons, lats, datatab, masktab 
     88    , file_ncdf, yyyy, mm, dd, hh, mn, ss, lons, lats, datatab, masktab, val_err, numch   
    8989; 
    9090compile_opt idl2, strictarrsubs 
     
    151151; Variable : data : 
    152152; ----------------- 
    153 data_id = NCDF_VARDEF(id, 'tbch5b', [dlong_id,dlatd_id,dtime_id], /float) 
     153data_id = NCDF_VARDEF(id, 'tb'+numch, [dlong_id,dlatd_id,dtime_id], /float) 
    154154NCDF_ATTPUT, id, data_id, 'units', "Kelvin" 
    155 NCDF_ATTPUT, id, data_id, 'long_name', "AMSU-B Channel 5 Brightness Temperature" 
     155NCDF_ATTPUT, id, data_id, 'long_name', "AMSU-"+strupcase(strmid(numch,0,1))+" Channel "+strmid(numch,1,1)+" Brightness Temperature" 
    156156NCDF_ATTPUT, id, data_id, 'valid_range', [100,500] 
    157157if (debug eq 1) then begin 
    158158   print, 'dans create_amsu_netcdf' 
    159159   help, datatab 
    160    datatab_novm=where(datatab ne -999.) 
     160   datatab_novm=where(datatab ne val_err) 
    161161   help, datatab_novm 
    162162endif 
    163 NCDF_ATTPUT, id, data_id, 'actual_range', [min(datatab[where(datatab ne -999.)]),max(datatab[where(datatab ne -999.)])] 
    164 NCDF_ATTPUT, id, data_id, 'missing_value', -999. 
    165 NCDF_ATTPUT, id, data_id, '_FillValue', -999. 
     163NCDF_ATTPUT, id, data_id, 'actual_range', [min(datatab[where(datatab ne val_err)]),max(datatab[where(datatab ne val_err)])] 
     164NCDF_ATTPUT, id, data_id, 'missing_value', val_err 
     165NCDF_ATTPUT, id, data_id, '_FillValue', val_err 
    166166; 
    167167; Variable : Mask : 
     
    171171NCDF_ATTPUT, id, mask_id, 'long_name', "land-sea mask, 1 all land, 0 all sea" 
    172172NCDF_ATTPUT, id, mask_id, 'valid_range', [0,1] 
    173 NCDF_ATTPUT, id, mask_id, 'actual_range', [min(masktab[where(datatab ne -999.)]),max(masktab[where(datatab ne -999.)])] 
    174 NCDF_ATTPUT, id, mask_id, 'missing_value', -999. 
    175 NCDF_ATTPUT, id, mask_id, '_FillValue', -999. 
     173NCDF_ATTPUT, id, mask_id, 'missing_value', val_err 
     174NCDF_ATTPUT, id, mask_id, '_FillValue', val_err 
    176175 
    177176NCDF_CONTROL, id, /ENDEF 
Note: See TracChangeset for help on using the changeset viewer.