Changeset 477 for trunk

Timestamp:

01/04/12 11:42:48 (12 years ago)

Author:

lelod

Message:

ca marche???

Location:

trunk/src

Files:

• amsu2ncdf.pro (modified) (1 diff)
• cresamsu.pro (modified) (7 diffs)
• extract_amsua.pro (modified) (2 diffs)
• idl_amsu_netcdf.pro (modified) (10 diffs)

trunk/src/amsu2ncdf.pro

@@ -249 +249 @@
 ;domaine -25 0 W, 0 - 20 N
 ; old grid filemsgref=project_id_env + '/MSG/2006/08/01/200608010000_msg-tb108_map_15min.nc'
-filemsgref=project_id_env + '/MSG/2006/06/01/msg_IR108_20060601.nc'
+filemsgref=project_id_env + '/EPSAT/extracted-surf-rr_epsat-sg_multi-sat_010d_30min_20060801_v3.1-02c.nc'
 print, 'Lecture de la grille reguliere ',filemsgref
 reggrid=read_regular_grid(filemsgref)

trunk/src/cresamsu.pro

@@ -56 +56 @@
 ; aligner lal/lon de l'exemple sur les autres modules amsu
 ;
-; lelod 2011/12/15
-; modifier la creation du ncdf pour ajouter la variable "mask"
-;
 ; lelod 2011/10/24
 ; actuellement codes en dur
@@ -147 +144 @@
 
 ; rayons de cressmann
-rhor=60. ; en gros deux fois la resol de AMSUB, et proche de celle d'AMSUA
+rhor=120. ; en gros deux fois la resol de AMSUA
 rtemp=0.5 ;en jours rayon de Cressmann (0.25, 0.5, 1, 2,... - nb entiers!)
 r2=rhor*rhor
@@ -244 +241 @@
 ; boucle sur les points du fichier
 
-         for i=0,nn do begin
+         for i=0,nn-1 do begin
 
 ; recup de tb1,hour,lon,lat
@@ -256 +253 @@
                if ((jour ge hrdeb-rtemphor) and (jour le hrfin+rtemphor) and (hour[i] le hfin) and (hour[i] ge hdeb) ) then begin
                                 ; print,'on cherche les donnees dans la fenetre temporelle'
-                                ; utilisation seulement des orbites du
-                                ; matin - projet heat low
+                                ; utilisation seulement des orbites du matin - projet heat low
                                 ; grilles temporelle et spatiale utile
                   jdp=fix((jour-hrdeb-rtemphor)/pasthor)>0
@@ -291 +287 @@
                               rnorm=2.
                               poids=(rnorm-dist-disj)/(rnorm+dist+disj)
-                              stb1[k,l,nj]=stb1[k,l,nj]+tb1*poids
+                              stb1[k,l,nj]=stb1[k,l,nj]+tb1[i]*poids
                               z1[k,l,nj]=z1[k,l,nj]+poids
                               kont1[k,l,nj]=kont1[k,l,nj]+1
@@ -301 +297 @@
                endif
             endif
+; fin test sur donnees valides (inutil normalement)
          endfor
          
@@ -321 +318 @@
 print,'extremes du compteur',min(kont1),max(kont1)
 
-; recreation du masque en 2D. On pourrait le deplacer dans idl_amsu_netcdf?
-file=project_id_env+'/MASK/ETOPO1_Ice_g_gmt4.nc'
-domdef,min(xxlon), max(xxlon), min(yylat), max(yylat)
-bate = ncdf_lec(file, var = 'z') GT 0
-xxe = float(ncdf_lec(file, var = 'lon'))
-yye = float(ncdf_lec(file, var = 'lat'))
-jpie = n_elements(xxe)
-jpje = n_elements(yye)
-maskfin=fltarr(nblon,nblat)
-maskfin=interpolate(bate,indgen(nblon),indgen(nblat),/grid)
-ind=where(maskfin gt 0.5,nland)
-if nland ne 0 then begin
-   maskfin(ind)=1
-endif
-ind=where(maskfin le 0.5,nsea)
-if nsea ne 0 then begin
-   maskfin(ind)=0
-endif
 ;toto=reform(tb1fin[*,*,0])
                                 ;plt,xxlon,yylat,toto
-idl_amsu_netcdf, numch, tb1fin, yylat, xxlon, jjours, maskfin
-
-
-
+idl_amsu_netcdf, numch, tb1fin, yylat, xxlon, jjours
+
+
+print,'fin cresamsu'
 end

trunk/src/extract_amsua.pro

@@ -422 +422 @@
 ; selection des taches au sol situees dans la zone d'interet
          zone=where((longrid ge lon_min) and (longrid le lon_max) $
-                    and (latgrid ge lat_min) and (latgrid le lat_max) and (tbgrid gt tbmin), npt)
+                    and (latgrid ge lat_min) and (latgrid le lat_max) and (tbgrid gt tbmin) and (tbgrid lt tbmax), npt)
          if npt ne 0 then begin
             if (n_elements(chint) EQ 0) then begin
@@ -450 +450 @@
 
    nn=n_elements(chint)
- ;  for i=0L, npt-1L do begin
- ;     if (( chint[i] lt tbmin) or (chint[i] gt tbmax) ) then begin
- ;        chint[i]=!values.f_nan
- ;     endif
- ;  endfor
-
 
    print,'nb points dans la zone',nn

trunk/src/idl_amsu_netcdf.pro

@@ -27 +27 @@
 ;           lats [shape=tab, label="latitude 1d array"]
 ;           timesin [shape=tab, label="time 1d array"]
-;           maskin [shape=tab, label="landseamask 3d array lon, lat, time"]
 ;
 ;           amsu_t5 [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/AMSU/numch_yyyymmdd_yyyymmdd_geomin_geomax.nc"];
@@ -51 +50 @@
 ;   IDL> timesin = indgen(12) + julday(8,1,2006L)
 ;   IDL> datain = findgen(size(lons,/dim), size(lats,/dim), size(timesin,/dim))
-;   IDL> maskin = findgen(size(lons,/dim), size(lats,/dim))
-;   IDL> idl_amsu_netcdf, numch , datain, lats, lons, timesin, maskin
+;   IDL> idl_amsu_netcdf, numch , datain, lats, lons, timesin
 ;
 ; SEE ALSO
@@ -100 +98 @@
 ; - fplod 20111216T160934Z cratos (Linux)
 ;
-;   * add maskin in the example and graph
 ;   * add a print to see filename on log
 ;   * fix landseamask unit
@@ -130 +127 @@
 ;
 ;-
-pro idl_amsu_netcdf, numch, datain, lats, lons, timesin, maskin
+pro idl_amsu_netcdf, numch, datain, lats, lons, timesin
 ;
 compile_opt idl2, strictarrsubs
@@ -161 +158 @@
 ; definition des tableaux qui vont contenir les valeurs des fichiers ascii
 moyenne_tb = fltarr(nblon, nblat, nbtime)
-landseamask= intarr(nblon, nblat)
 xlon = fltarr(nblon)
 ylat = fltarr(nblat)
@@ -167 +163 @@
 ;
 moyenne_tb[*,*,*] = datain[*,*,*]
-landseamask[*,*] = maskin[*,*]
 xlon[*] = lons
 ylat[*] = lats
@@ -173 +168 @@
 ;
 ; masque des valeurs aberrantes
-ind = where((moyenne_tb le 0L) OR (moyenne_tb GT 350L), nnbad)
+ind = where(finite(moyenne_tb,/nan), nnbad)
+if nnbad ne 0 then begin
+   moyenne_tb[ind]=-9999.
+endif
 print,"nb valeurs manquantes", nnbad
 ;
@@ -192 +190 @@
 ncdf_attput, id, 'File_Name', cdfout, /GLOBAL
 ncdf_attput, id, 'IDL_Program_Name', 'idl_amsu_netcdf', /GLOBAL
-creation_date = string(jul2date(SYSTIME( /JULIAN , /UTC )), format='(I8)')
-ncdf_attput, id, 'Creation_date', creation_date, /GLOBAL
+;print,'creation date'
+;creation_date = string(jul2date(SYSTIME(/JULIAN, /UTC )), format='(I8)')
+;ncdf_attput, id, 'Creation_date', creation_date, /GLOBAL
+;print,creation_date
 ;
 ; declaration des variables
-varid = lonarr(5)
+varid = lonarr(4)
 varid[0] = ncdf_vardef(id, 'xlon', [dimidx], /double)
 varid[1] = ncdf_vardef(id, 'ylat', [dimidy], /double)
 varid[2] = ncdf_vardef(id, 'time', [dimidt], /double)
 varid[3] = ncdf_vardef(id, 'moyenne_tb',[dimidx,dimidy,dimidt], /float)
-varid[4] = ncdf_vardef(id, 'landseamask',[dimidx,dimidy,dimidt], /short)
 ;
 ; attributs de chaque variable
@@ -212 +211 @@
 ncdf_attput, id, varid[3], 'time_origin', '2000-01-01 00:00:00'
 ncdf_attput, id, varid[3], 'short_name', 'amsu_tb'
-ncdf_attput, id, varid[3], 'missing_value', '-999'
+ncdf_attput, id, varid[3], 'missing_value', '-9999.00'
 ncdf_attput, id, varid[3], 'unit', 'degree_kelvin'
-ncdf_attput, id, varid[4], 'short_name', 'mask'
-ncdf_attput, id, varid[4], 'unit', 'number (0 ocean, 1 land)'
 ;
 ; remplissage de l'entete
@@ -225 +222 @@
 ; remplissage des variables
 ncdf_varput, id, varid[3], moyenne_tb
-ncdf_varput, id, varid[4], landseamask
 ;
 ; fermeture du fichier