source: trunk/SRC/Obsolete/nlecmois.pro @ 118

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

add $ in Calendar, Grid, Interpolation, Obsolete and Postscript *.pro files, add svn:keywords Id to all these files, some improvements in header
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 6.9 KB
Line 
1;------------------------------------------------------------
2;------------------------------------------------------------
3;------------------------------------------------------------
4;+
5;
6; @file_comments lit les fichiers Net-Cdf  de moyenne mensuel  de
7; l'experience TOTEM/ECMWF qui sont sur
8; maia: /u/rech/eee/reee217/TOTEM/REF/OUTPUTS
9; @obsolete
10;
11; @keyword BOITE boite sur laquelle integrer (par defaut tt le domaine)
12;
13; @keyword GRILLE impose la grille a laquelle est rapporte le champ. rq permet
14; d'aller plus vite ds la lecture
15;
16; @keyword /TOUT oblige a lire le tableau entier en non pas celui reduit a domdef
17; @keyword _EXTRA used to pass your keywords
18
19; @uses common.pro
20;
21; @history Sebastien Masson (smasson\@lodyc.jussieu.fr)
22;
23; @version $Id$
24;
25;-
26;------------------------------------------------------------
27;------------------------------------------------------------
28;REF.moyenne.mensuelle.79-81.grid.T.nc
29;------------------------------------------------------------
30function nlecmois,nom,numsortie79,nbretps ,BOITE=boite,GRILLE=grille,TOUT=tout, _EXTRA = ex
31;
32  compile_opt idl2, strictarrsubs, obsolete
33;
34@common
35;------------------------------------------------------------
36   res =-1
37;------------------------------------------------------------
38; annee du nom du 1er fichier
39;------------------------------------------------------------
40   annee = floor(floor((numsortie79-1)/12.)/3.)*3+79
41;------------------------------------------------------------
42; gestion du nom de la grille en testant les differentes possibilites
43; et ouverture du fichier
44;------------------------------------------------------------
45   if keyword_set(grille) then begin
46      vargrid=grille
47      nomfich=prefix+'.moyenne.mensuelle.'+strtrim(annee, 1)+'-'+strtrim(annee+2, 1)+'.grid.'+vargrid+'.nc'
48      IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+iodir+nomfich+' > /dev/null'
49      cdfid=ncdf_open(iodir+nomfich)
50      varcontient=ncdf_varinq(cdfid,nom)
51   endif else begin
52; liste des fichiers pouvant convenir
53      quelsfichiers = findfile(iodir+prefix+'.moyenne.mensuelle.'+strtrim(annee, 1)+'-'+strtrim(annee+2, 1)+'.grid.*')
54      IF quelsfichiers[0] EQ '' THEN $ ; liste vide
55         return, report('LES FICHIERS: '+iodir+prefix+'.moyenne.mensuelle.'+strtrim(annee, 1)+'-'+strtrim(annee+2, 1)+'.grid... n''existe pas.')
56      for i=0,n_elements(quelsfichiers)-1 do begin
57         IF !version.OS_FAMILY EQ 'unix' THEN spawn, '\file '+quelsfichiers[i]+' > /dev/null'
58         cdfid=ncdf_open(quelsfichiers[i])
59         contient=ncdf_inquire(cdfid)
60         for varid=0,contient.nvars-1 do BEGIN ; ds les fichiers existants on
61            varcontient=ncdf_varinq(cdfid,varid) ; cherche le nom des variables
62            if varcontient.name eq nom then BEGIN
63               vargrid=strmid(quelsfichiers[i],strpos(quelsfichiers[i],'grid')+5,1) ;nom de grille
64               goto, grilletrouvee
65            ENDIF
66         endfor
67         ncdf_close,cdfid
68      endfor
69      return, report('La variable '+nom+' n''existe pas ds les fichiers'+iodir+prefix+'.moyenne.mensuelle.'+strtrim(annee, 1)+'-'+strtrim(annee+2, 1)+'grid...')
70   endelse
71grilletrouvee:
72;------------------------------------------------------------
73; lecture de certains attributs et de l'axe des temps
74;------------------------------------------------------------
75   ncdf_attget,cdfid,nom,'title',value
76   varname=string(value)
77   ncdf_attget,cdfid,nom,'units',value
78   varunit=string(value)
79   if rstrpos(varname,'(') eq -1 then varname=varname+' ('+varunit+')'
80   varexp=prefix
81;------------------------------------------------------------
82; determination du domaine geographique
83;------------------------------------------------------------
84   if keyword_set(tout) then begin
85      nx=jpi
86      ny=jpj
87      nz=jpk
88      premierx=0
89      premiery=0
90      premierz=0
91   endif else BEGIN
92;------------------------------------------------------------
93; redefinition eventuelle du domaine ajuste a boite (a 6 elements)
94; + on recupere la dim du no9uveau domaine
95;------------------------------------------------------------
96      if keyword_set(boite) then BEGIN
97         Case N_Elements(Boite) Of
98            1:Domdef, lon1, lon2, lat1, lat2, 0.,boite[0],GRILLE=vargrid, _EXTRA = ex
99            2:Domdef, lon1, lon2, lat1, lat2, boite[0],boite[1],GRILLE=vargrid, _EXTRA = ex
100            4:Domdef, [Boite, prof1, prof2],GRILLE=vargrid, _EXTRA = ex
101            5:Domdef, [Boite[0:3], 0, Boite[4]],GRILLE=vargrid, _EXTRA = ex
102            6:Domdef, Boite,GRILLE=vargrid, _EXTRA = ex
103            Else:return,  report('Mauvaise Definition de Boite')
104         endcase
105      ENDIF
106      grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery, premierz
107   ENDELSE
108;------------------------------------------------------------
109; determination du nombre de pas de tps a extraire ds ce fichier
110;------------------------------------------------------------
111   numsortie = numsortie79-12.*(annee-79)
112   if nbretps+numsortie-1 gt 36 then nt=36-numsortie+1 $
113   else nt=nbretps
114   numsortie = numsortie79-12.*(annee-79)
115;------------------------------------------------------------
116    if varcontient.ndims eq 3 then $
117    ncdf_varget,cdfid,nom,res,offset=[premierx,premiery,numsortie-1],count=[nx,ny,nt] else $
118    ncdf_varget,cdfid,nom,res,offset=[premierx,premiery,premierz,numsortie-1],count=[nx,ny,nz,nt]
119;------------------------------------------------------------
120; rappel en boucle de nlec si il faut ouvrir de nouveaux fichiers
121; pour constituer la serie temporelle
122;------------------------------------------------------------
123   if nbretps gt 36-numsortie+1 then begin
124      if varcontient.ndims eq 3 then $
125       res = [[[res]] $
126              ,[[nlecmois(nom,numsortie79+nt,nbretps-nt,tout=tout,GRILLE=vargrid,BOITE=boite)]]] $
127      else BEGIN
128         res=[res[*] $
129              ,(nlecmois(nom,numsortie79+nt,nbretps-nt,tout=tout,GRILLE=vargrid,BOITE=boite))[*]]
130
131         res=reform(res,nx,ny,nz,nbretps, /over)
132      ENDELSE
133      ncdf_varget,cdfid,'time_counter',temps, offset = [numsortie-1],count=[nt]
134      time = [long(temps)+julday(12, 31, 1978), time]
135      jpt=nt+jpt
136   endif else BEGIN
137      ncdf_varget,cdfid,'time_counter',temps, offset = [numsortie-1],count=[nt]
138      time = long(temps)+julday(12, 31, 1978)
139      jpt=nt
140   endelse
141;------------------------------------------------------------
142   ncdf_close,cdfid
143   IF n_elements(time) EQ 1 THEN BEGIN
144      caldat, time, m, d, y
145      if m lt 10 then m='0'+string(format='(i1)',m) $
146      else m=    string(format='(i2)',m)
147      if n_elements(langage) EQ 0 then langage = 'non definit'
148      if langage EQ 'gb' then vardate = strtrim(y, 1)+' '+string(format='(C(CMoA))',31*(m-1)) $
149      ELSE vardate = string(format='(C(CMoA))',31*(m-1))+' '+strtrim(y, 1)
150   endif
151;------------------------------------------------------------
152   return,res
153;------------------------------------------------------------
154;------------------------------------------------------------
155end
Note: See TracBrowser for help on using the repository browser.