Changeset 594

Timestamp:

07/23/12 15:06:54 (12 years ago)

Author:

lelod

Message:

adaptation de cresamsu a amsub

Location:

trunk/src

Files:

• cresamsu.pro (modified) (3 diffs)
• idl_amsu_netcdf.pro (modified) (7 diffs)

trunk/src/cresamsu.pro

@@ -234 +234 @@
    hdeb=0
    hfin=24                      ; on prend toutes les orbites de la journee (regule par rtemp)  
-   if rtemp eq 0.5 then pastemp=0.5
    if rtemp gt 0.5 then pastemp=1
-   if rtemp lt 0.5 then pastemp=0.25
+   if rtemp le 0.5 then pastemp=0.5
+   if rtemp le 0.175 then pastemp=0.25
 endif
 pasthor=pastemp*24.
@@ -510 +510 @@
 
 if nbfile GT 0 THEN BEGIN
-    good=where (kont1 ge 1,nb)
+    good=where (kont1 ge 3,nb)
     msg = report(['iii : ' + routine + ' : points valides = ' + string(nb)])
     if (nb ne 0) then tb1fin[good]=stb1[good]/z1[good]
@@ -519 +519 @@
     msg = report(['iii : ' + routine + ' : extremes du compteur = ' + string(min(kont1)) + ' ' + string(max(kont1))])
 
-    result=idl_amsu_netcdf(numch, temp,rtemp,rhdeg, tb1fin, yylat, xxlon, jjours)
+    result=idl_amsu_netcdf(numch, temp,rtemp,rhdeg, tb1fin,z1, yylat, xxlon, jjours)
     result_type=size(result,/type)
     IF (result_type EQ 3) THEN BEGIN

trunk/src/idl_amsu_netcdf.pro

@@ -163 +163 @@
 ;
 ;-
-function idl_amsu_netcdf, numch, temp, rtemp, rhdeg, datain, lats, lons, timesin
+function idl_amsu_netcdf, numch, temp, rtemp, rhdeg, datain, weight, lats, lons, timesin
 ;
 compile_opt idl2, strictarrsubs
@@ -220 +220 @@
 ; dans les fichiers NetCDF
 moyenne_tb = fltarr(nblon, nblat, nbtime)
+poids = fltarr(nblon, nblat, nbtime)
 xlon = fltarr(nblon)
 ylat = fltarr(nblat)
@@ -225 +226 @@
 ;
 moyenne_tb[*,*,*] = datain[*,*,*]
+poids[*,*,*] = weight[*,*,*]
 xlon[*] = lons
 ylat[*] = lats
@@ -233 +235 @@
 if nnbad ne 0 then begin
    moyenne_tb[ind]=-9999.
+   poids[ind]=-9999.
 endif
 msg = report(['iii : ' + routine + ' : nb valeurs manquantes = ' + string(nnbad)])
@@ -268 +271 @@
 ;
 ; declaration des variables
-varid = lonarr(4)
+varid = lonarr(5)
 varid[0] = ncdf_vardef(id, 'xlon', [dimidx], /float)
 varid[1] = ncdf_vardef(id, 'ylat', [dimidy], /float)
 varid[2] = ncdf_vardef(id, 'time', [dimidt], /float)
 varid[3] = ncdf_vardef(id, 'moyenne_tb',[dimidx,dimidy,dimidt], /float)
+varid[4] = ncdf_vardef(id, 'poids',[dimidx,dimidy,dimidt], /float)
 ;
 ; attributs de chaque variable
@@ -291 +295 @@
 ncdf_attput, id, varid[3], 'unit', 'degree_kelvin'
 ;
+ncdf_attput, id, varid[4], 'short_name', 'poids'
+ncdf_attput, id, varid[4], 'missing_value', -9999.00
+ncdf_attput, id, varid[4], 'unit', 'valeur entre 0 et 1'
+
 ; remplissage de l'entete
 ncdf_control, id, /endef
@@ -299 +307 @@
 ; remplissage des variables
 ncdf_varput, id, varid[3], moyenne_tb
+ncdf_varput, id, varid[4], poids
 ;
 ; fermeture du fichier