Changeset 118 for trunk/SRC/Obsolete/lec.pro
- Timestamp:
- 06/27/06 17:47:06 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Obsolete/lec.pro
r114 r118 20 20 ; 21 21 ; @param nomchamp {in}{required} 2 choix possibles: 22 ; 1) nom de champ Vairmer (chaine de 8 caracteres en majuscule ou 22 ; 1) nom de champ Vairmer (chaine de 8 caracteres en majuscule ou 23 23 ; minuscule commencant par vo ou so). Dans cette methode on saute directement 24 24 ; d'en-tete en en-tete jusqu'a trouver le bon fichier. … … 31 31 ; 32 32 ; 33 ; @keyword /ANOM {in} type du fichier vairmer par rapport auquel on doit calculer 33 ; @keyword /ANOM {in} type du fichier vairmer par rapport auquel on doit calculer 34 34 ; l'anomalie ('EX','AN','SE','MO','') 35 35 ; … … 37 37 ; ds ce cas en input on met seulement 'vo' ou 'so' la fonction retourne le 38 38 ; nombre de fichiers lus. 39 ; 40 ; @keyword BOITE 39 41 ; 40 42 ; @keyword /EXPANOM {in} si on calcule l'anom par rapport a une exper … … 60 62 ; celle du sous domaine definit par domdef (nx,ny,nz). avec /TOUT le 61 63 ; champ a la taille du domaine qui a etait selectionne pour la 62 ; cession en cours (jpi,jpj,jpk). 63 ; pour les sous domaines cf: 64 ; cession en cours (jpi,jpj,jpk). 65 ; pour les sous domaines cf: 64 66 ; http://www.ipsl.jussieu.fr/~smlod/sousdomaine.html 65 67 ; Retourne -1 en cas d'erreur. 66 68 ; 67 ; @uses common.pro isnumber.pro fivardate.pro 69 ; @uses common.pro isnumber.pro fivardate.pro 68 70 ; 69 71 ; @history Sebastien Masson (smasson\@lodyc.jussieu.fr) … … 76 78 ; Sebastien Masson (decoupe pour les sous domaines...) 77 79 ; 2/99 80 ; 81 ; @version $Id$ 82 ; 78 83 ;;- 79 84 ;-------------------------------------------------------------- … … 141 146 ; definition de la taille du fichier a aller chercher: jpidta,jpjdta,jpkdta... 142 147 ;--------------------------------------------------------------------- 143 if n_elements(jpidta) EQ 0 THEN BEGIN 148 if n_elements(jpidta) EQ 0 THEN BEGIN 144 149 if n_elements(ixmindta) EQ 0 OR n_elements(ixmaxdta) EQ 0 then $ 145 150 jpidta = jpiglo else jpidta = ixmaxdta-ixmindta+1 146 151 endif 147 if n_elements(jpjdta) EQ 0 THEN BEGIN 152 if n_elements(jpjdta) EQ 0 THEN BEGIN 148 153 if n_elements(iymindta) EQ 0 OR n_elements(iymaxdta) EQ 0 then $ 149 154 jpjdta = jpjglo else jpjdta = iymaxdta-iymindta+1 150 155 endif 151 if n_elements(jpkdta) EQ 0 THEN BEGIN 156 if n_elements(jpkdta) EQ 0 THEN BEGIN 152 157 if n_elements(izmindta) EQ 0 OR n_elements(izmaxdta) EQ 0 then $ 153 158 jpkdta = jpkglo else jpkdta = izmaxdta-izmindta+1 … … 267 272 ; on shift z si key_shift est defininit 268 273 ;--------------------------------------------------------------------- 269 if n_elements(key_shift) NE 0 THEN BEGIN 274 if n_elements(key_shift) NE 0 THEN BEGIN 270 275 if dim EQ 'SO' then z = shift(z,key_shift, 0) $ 271 276 ELSE z = shift(z,key_shift, 0, 0) … … 275 280 ; du zoom: nx,ny nz 276 281 ;--------------------------------------------------------------------- 277 if NOT keyword_set(tout) then BEGIN 282 if NOT keyword_set(tout) then BEGIN 278 283 ;------------------------------------------------------------- 279 284 ; changement de domaine 280 285 ;------------------------------------------------------------- 281 if keyword_set(boite) then BEGIN 286 if keyword_set(boite) then BEGIN 282 287 Case 1 Of 283 288 N_Elements(Boite) Eq 1:bte=[lon1, lon2, lat1, lat2, 0.,boite[0]] … … 290 295 oldboite = [lon1, lon2, lat1, lat2, prof1, prof2] 291 296 domdef, bte,GRILLE=vargrid 292 ENDIF 297 ENDIF 293 298 ;------------------------------------------------------------- 294 299 grille,mask,glam,gphi,gdep,nx,ny,nz,premierx,premiery,premierz,dernierx,derniery,dernierz … … 296 301 if dim EQ 'SO' then z = z[premierx:dernierx, premiery:derniery] $ 297 302 ELSE z = z[premierx:dernierx, premiery:derniery, premierz:dernierz] 298 ENDIF ELSE BEGIN 303 ENDIF ELSE BEGIN 299 304 case vargrid OF ; on recupere le mask en entier ds le cas ou /TOUT 300 305 'U':mask = umask() ; n''est pas active et on le choisit en fontion 301 306 'T':mask = tmask ; de la valeur de vargrid 302 'W':mask = tmask 307 'W':mask = tmask 303 308 'V':mask = vmask() 304 309 'F':mask = fmask() 305 310 ENDCASE 306 ENDELSE 311 ENDELSE 307 312 ;--------------------------------------------------------------------- 308 313 ; calcul d'une anomalie si le keyword anom est active … … 311 316 case anom of 312 317 'EX' : adate = 0 313 'AN' : adate = floor(date/10000)*10000 318 'AN' : adate = floor(date/10000)*10000 314 319 'SE' : adate = floor(date - floor(date/10000)*10000)/100 * 100 315 320 'MO' : adate = floor(date/100)*100 … … 327 332 ;--------------------------------------------------------------------- 328 333 IF n_elements(valmask) EQ 0 THEN valmask = 1e20 329 if dim EQ 'SO' then BEGIN 334 if dim EQ 'SO' then BEGIN 330 335 terre = where(mask[*,*,0] EQ 0) 331 if terre[0] NE -1 then z[terre] = valmask 332 ENDIF ELSE BEGIN 336 if terre[0] NE -1 then z[terre] = valmask 337 ENDIF ELSE BEGIN 333 338 terre = where(mask[*,*,0] EQ 0) 334 339 if terre[0] NE -1 then z[where(mask EQ 0)] = valmask 335 ENDELSE 340 ENDELSE 336 341 ;--------------------------------------------------------------------- 337 342 free_lun,numlec
Note: See TracChangeset
for help on using the changeset viewer.