source: trunk/Obsolete/nlecmois.pro @ 44

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

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