source: trunk/SRC/Obsolete/nlecsaison.pro @ 114

Last change on this file since 114 was 114, checked in by smasson, 18 years ago

new compilation options (compile_opt idl2, strictarrsubs) in each routine
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 6.1 KB
Line 
1;------------------------------------------------------------
2;------------------------------------------------------------
3;------------------------------------------------------------
4;+
5; @file_comments lit les moyennes annuelles
6; sur maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS
7; @obsolete
8;
9; @categories lecture de NETCDF
10;
11; @examples
12; IDL> res=nlecsaison(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;REF.saisonnier.climato.grid.T.nc
39;------------------------------------------------------------
40;------------------------------------------------------------
41function nlecsaison, nom, numsortie79, nbretps,BOITE = boite, GRILLE=grille,TOUT=tout, _EXTRA = ex
42;
43  compile_opt idl2, strictarrsubs, obsolete
44;
45@common
46;------------------------------------------------------------
47;------------------------------------------------------------
48   res=-1
49   anneedepart=1979
50;------------------------------------------------------------
51; gestion du nom de la grille en testant les differentes possibilites
52; et ouverture du fichier
53;------------------------------------------------------------
54   if keyword_set(grille) then begin
55      vargrid=grille
56      nomfich=prefix+'.saisonnier.climato.grid.'+vargrid+'.nc'
57      IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+iodir+nomfich+' > /dev/null'
58      cdfid=ncdf_open(iodir+nomfich)
59      varcontient=ncdf_varinq(cdfid,nom)
60   endif else begin
61; liste des fichiers pouvant convenir
62      quelsfichiers = findfile(iodir+prefix+'.saisonnier.climato.grid.*')
63      IF quelsfichiers[0] EQ '' THEN $ ; liste vide
64         return, report('LES FICHIERS: '+iodir+prefix+'.saisonnier.climato.grid... n''existe pas.')
65      for i=0,n_elements(quelsfichiers)-1 do begin
66         IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+quelsfichiers[i]+' > /dev/null'
67         cdfid=ncdf_open(quelsfichiers[i])
68         contient=ncdf_inquire(cdfid)
69         for varid=0,contient.nvars-1 do BEGIN   ; ds les fichiers existants on
70            varcontient=ncdf_varinq(cdfid,varid) ; cherche le nom des variables
71            if varcontient.name eq nom then BEGIN
72               vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille
73               goto, grilletrouvee
74            ENDIF
75         endfor
76         ncdf_close,cdfid
77      endfor
78      return, report('La variable '+nom+' n''existe pas ds les fichiers'+iodir+prefix+'.saisonnier.climato.grid...')
79   endelse
80grilletrouvee:
81;------------------------------------------------------------
82; lecture de certains attributs et de l'axe des temps
83;------------------------------------------------------------
84   ncdf_attget,cdfid,nom,'title',value
85   varname=string(value)
86   ncdf_attget,cdfid,nom,'units',value
87   varunit=string(value)
88   if rstrpos(varname,'(') eq -1 then varname=varname+' ('+varunit+')'
89   varexp=prefix
90;------------------------------------------------------------
91; lecture de l'axe des temps
92;------------------------------------------------------------
93   ncdf_varget,cdfid,'time_counter',time, offset = [numsortie79-1],count=[nbretps]
94   time = long(time)+julday(12, 31, 1978)
95   jpt=nbretps
96   IF jpt EQ 1 THEN BEGIN
97      IF numsortie79 EQ 13 THEN vardate='CLIMATOLOGIE '+prefix $
98      ELSE begin
99         vardate='climato mensuelle '+strtrim(numsortie79,1)
100      endelse
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 $
140    ,offset=[premierx,premiery,numsortie79-1],count=[nx,ny,nbretps] else $
141    ncdf_varget,cdfid,nom,res $
142    ,offset=[premierx,premiery,premierz,numsortie79-1],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.