Changeset 25 for trunk/ToBeReviewed/CALCULS/depth2level.pro

Timestamp:

05/02/06 14:59:12 (18 years ago)

Author:

pinsard

Message:

upgrade of CALCULS according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/ : files

File:

• trunk/ToBeReviewed/CALCULS/depth2level.pro (copied) (copied from trunk/CALCULS/depth2level.pro) (4 diffs)

trunk/ToBeReviewed/CALCULS/depth2level.pro

@@ -60 +60 @@
 ;------------------------------------------------------------
    in = litchamp(tab)
-   grille,mask,glam,gphi,gdep,nx,ny,nz,premierx,premiery,premierz,dernierx,derniery,dernierz
-   glam = 1
-   gphi = 1
+   grille,mask,-1,-1,gdep,nx,ny,nz,firstx,firsty,firstz,lastx,lasty,lastz
 ;---------------------------------------------------------------
 ; verification de la coherence entre la taille du tableau et le domaine definit par domdef
@@ -70 +68 @@
    if taille[0] NE 2 then return, report('le champ en entree doit contenir un tableau 2d')
    case 1 of
-      taille[1] eq jpi and taille[2] eq jpj:in=in[premierx:dernierx, premiery:derniery]
+      taille[1] eq jpi and taille[2] eq jpj:in=in[firstx:lastx, firsty:lasty]
       taille[1] eq  nx and taille[2] eq  ny:
       else:return, report('Probleme d''adequation entre les tailles du domaine et celle du champ.')
@@ -85 +83 @@
 ;------------------------------------------------------------
 ; on passe en tableaux qui ont la taille des tableaux 3d
-   prof=replicate(1,nx*ny)#gdep[premierz:dernierz] 
+   prof=replicate(1,nx*ny)#gdep[firstz:lastz] 
    in = in[*]#replicate(1, nz)
 ;
@@ -94 +92 @@
    if keyword_set(upper) then begin
       levels = levels-1
-      notvalid = where(levels EQ 0)
+      notvalid = where(levels EQ -1)
    ENDIF ELSE notvalid = where(levels EQ nz)
    IF notvalid[0] NE -1 THEN levels[notvalid] = !values.f_nan