source: trunk/SRC/Obsolete/nlecan.pro @ 97

Last change on this file since 97 was 97, checked in by pinsard, 18 years ago

start to modify headers of Obsolete *.pro files for better idldoc output
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 6.2 KB
Line 
1;------------------------------------------------------------
2;------------------------------------------------------------
3;------------------------------------------------------------
4;+
5; @file_comments lit les moyennes annuelles
6; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS
7;
8; @obsolete
9; @categories lecture de NETCDF
10;
11; @examples
12; IDL> res=nlecan(nom, numsortie79, nbretps)
13;
14; @param nom {in}{required} nom du tableau contenant le champ que l'on veut (le
15; trouver avec ncdflec,fichier,/var
16; @param numsortie79  {in}{required} le numero du pas de temps que l'on veut sortir du
17; fichier (compte a partir de 1) a partir de 79
18; @param nbretps  {in}{required} nombre de pas de temps a extraire
19;
20; @keyword BOITE  boite sur laquelle integrer (par defaut tt le domaine)
21;
22; @keyword GRILLE impose la grille a laquelle est rapporte le champ. rq permet
23; d'aller plus vite ds la lecture
24;
25; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef
26;
27; @returns tableau 2d (qd la serie ne fait que 1 pas de temps) ou
28; 3d (valable ds 1 premier tps que pour les tableaux 2d)
29;
30; @uses common.pro
31;
32; @restrictions appele par nlec
33;
34; @history Sebastien Masson (smasson\@lodyc.jussieu.fr)
35;
36;-
37;------------------------------------------------------------
38;------------------------------------------------------------
39;REF.moyenne.annuelle.grid.T.nc
40;------------------------------------------------------------
41function nlecan, nom, numsortie79, nbretps,BOITE = boite,GRILLE=grille,TOUT=tout, _EXTRA = ex
42@common
43;------------------------------------------------------------
44;------------------------------------------------------------
45   res=-1
46   anneedepart=1979
47;------------------------------------------------------------
48; gestion du nom de la grille en testant les differentes possibilites
49; et ouverture du fichier
50;------------------------------------------------------------
51   if keyword_set(grille) then begin
52      vargrid=grille
53      nomfich=prefix+'.moyenne.annuelle.grid.'+vargrid+'.nc'
54      IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+iodir+nomfich+' > /dev/null'     
55      cdfid=ncdf_open(iodir+nomfich)
56      varcontient=ncdf_varinq(cdfid,nom)
57   endif else begin
58; liste des fichiers pouvant convenir
59      quelsfichiers = findfile(iodir+prefix+'.moyenne.annuelle.grid*')
60      IF quelsfichiers[0] EQ '' THEN $ ; liste vide
61         return,  report('LES FICHIERS: '+iodir+prefix+'.moyenne.annuelle.grid... n''existe pas.')
62      for i=0,n_elements(quelsfichiers)-1 do begin
63         IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+quelsfichiers[i]+' > /dev/null'
64         cdfid=ncdf_open(quelsfichiers[i])
65         contient=ncdf_inquire(cdfid)
66         for varid=0,contient.nvars-1 do BEGIN   ; ds les fichiers existants on
67            varcontient=ncdf_varinq(cdfid,varid) ; cherche le nom des variables
68            if varcontient.name eq nom then BEGIN
69               vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille
70               goto, grilletrouvee
71            ENDIF
72         endfor
73         ncdf_close,cdfid
74      endfor
75      return, report('La variable '+nom+' n''existe pas ds les fichiers'+iodir+prefix+'.moyenne.annuelle.grid...')
76   endelse
77grilletrouvee:
78;------------------------------------------------------------
79; lecture de certains attributs et de l'axe des temps
80;------------------------------------------------------------
81   ncdf_attget,cdfid,nom,'title',value
82   varname=string(value)
83   ncdf_attget,cdfid,nom,'units',value
84   varunit=string(value)
85   if rstrpos(varname,'(') eq -1 then varname=varname+' ('+varunit+')'
86   varexp=prefix
87;------------------------------------------------------------
88; lecture de l'axe des temps
89;------------------------------------------------------------
90   ncdf_varget,cdfid,'time_counter',time, offset = [numsortie79-1],count=[nbretps]
91   time = long(time)+julday(12, 31, 1978)
92   jpt=nbretps
93   IF jpt EQ 1 THEN BEGIN
94      caldat, time[0], month, day, year
95      case 1 of
96         year lt 10:                 s_year ='0'+string(format='(i1)',year)
97         year lt 100 and year ge 10 :s_year =    string(format='(i2)',year)
98         year ge 100:                s_year =    string(format='(i2)',year-1900)
99      endcase
100      vardate='annee: '+s_year
101   endif
102;------------------------------------------------------------
103;------------------------------------------------------------
104; extraction du tableau qui nous interesse
105;------------------------------------------------------------
106;------------------------------------------------------------
107;------------------------------------------------------------
108; determination du domaine geographique
109;------------------------------------------------------------
110   if keyword_set(tout) then begin
111      nx=jpi
112      ny=jpj
113      nz=jpk
114      glam=[1]
115      gphi=[1]
116      gdep=[1]
117      premierx = 0
118      premiery = 0
119      premierz = 0
120   endif else BEGIN
121;------------------------------------------------------------
122; redefinition eventuelle du domaine ajuste a boite (a 6 elements)
123; + on recupere la dim du no9uveau domaine
124;------------------------------------------------------------
125      if keyword_set(boite) then BEGIN
126         Case N_Elements(Boite) Of
127            1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex
128            2:Domdef, lon1, lon2, lat1, lat2, boite[0],boite[1],GRILLE=vargrid, _EXTRA = ex
129            4:Domdef, [Boite,prof1, prof2],GRILLE=vargrid, _EXTRA = ex
130            5:Domdef, [Boite[0:3], 0, Boite[4]],GRILLE=vargrid, _EXTRA = ex
131            6:Domdef, Boite,GRILLE=vargrid, _EXTRA = ex
132            Else:return, report('Mauvaise Definition de Boite')
133         endcase
134      ENDIF
135      grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery, premierz
136   ENDELSE
137;------------------------------------------------------------
138   if varcontient.ndims eq 3 then $
139    ncdf_varget,cdfid,nom,res,offset=[premierx,premiery,numsortie79-1] $
140    ,count=[nx,ny,nbretps] else $
141    ncdf_varget,cdfid,nom,res,offset=[premierx,premiery, premierz,numsortie79-1] $
142    ,count=[nx,ny,nz,nbretps]
143;------------------------------------------------------------
144   ncdf_close,cdfid
145   return,res
146;------------------------------------------------------------
147;------------------------------------------------------------
148end
Note: See TracBrowser for help on using the repository browser.