Changeset 494 for trunk

Timestamp:

02/23/12 12:32:08 (12 years ago)

Author:

lelod

Message:

correction bug dans extract

File:

• trunk/src/extract_amsua.pro (modified) (3 diffs)

trunk/src/extract_amsua.pro

@@ -92 +92 @@
 ; EVOLUTIONS
 ; ==========
+; lelod 20120223
+; modif lecture du fichier bathy (ETOPO) pour limiter le temps
+; necessaire. Introduction de coordonnees limite dans initncdf, puis
+; domdef
+; Attention, si on change la zone, il faut changer les coordonnees!!!
 ;
 ; lelod 20111215
@@ -249 +254 @@
 ; lecture fichier land - sea (S. Masson)
 file=project_id_env+'/MASK/ETOPO1_Ice_g_gmt4.nc'
-initncdf,file, zaxisname='toto',xaxisname='lon',yaxisname='lat',XMINMESH=8900,XMAXMESH=14000, YMINMESH=4500,YMAXMESH=8000 
+;attention: zone prise en compte limitee par le choix des nos en x et
+;y dans initncdf. La region traitee ne doit pas exceder lon_min=-25 &
+;lon_max=40 & lat_min=-5 & lat_max=40
+
+initncdf,file, zaxisname='toto',xaxisname='lon',yaxisname='lat',XMINMESH=8000,XMAXMESH=15000, YMINMESH=4500,YMAXMESH=8500 
 domdef,lon_min-15, lon_max+15, lat_min-15, lat_max+15
 bathy=read_ncdf('z',file=file)
 bate=bathy.arr gt 0
-xxe=reform(glamt[*,0])
-yye=reform(gphit[0,*])
+xxe1=reform(glamt[*,0])
+zonx=where(xxe1 ge lon_min-15 and xxe1 le lon_max+15)
+xxe=xxe1[zonx]
+yye1=reform(gphit[0,*])
+zony=where(yye1 ge lat_min-15 and yye1 le lat_max+15)
+yye=yye1[zony]
 jpie = n_elements(xxe)
 jpje = n_elements(yye)
+;print,jpie,jpje,size(bate)
 ;PRINT, 'lecture bathy terminee',SYSTIME() 
 ; appel au ssprgm qui interpole les fonctions de correction (inutile
@@ -272 +286 @@
 a = STRARR(nb_file)
    ; filea = STRARR(nb_filea)
-print, 'ouverture pour lecture de ', list_file
+print, 'iii ouverture pour lecture de ', list_file
+print, 'iii traitement du jour ',yyyy,mm,dd
+
 openr, lun_a, list_file, /get_lun
 filea=''