Changeset 25 for trunk/ToBeReviewed/CALCULS/level2mask.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/level2mask.pro (copied) (copied from trunk/CALCULS/level2mask.pro) (4 diffs)

trunk/ToBeReviewed/CALCULS/level2mask.pro

@@ -6 +6 @@
 ;
 ; PURPOSE: permet de passer d''un tableau 2d de niveau seuil au
-; tableau 3d de mask avec des 1 ds les niveaux au dessus (et sur) du
-; niveau seuil et des 0 en dessous.
+; tableau 3d de mask avec des 1 ds les niveaux au dessus du
+; niveau seuil et des 0 en dessous (et sur).
 ;
 ; CATEGORY: SANS BOUCLE
@@ -32 +32 @@
 ; MODIFICATION HISTORY: Sebastien Masson (smasson@lodyc.jussieu.fr)
 ;                       17/6/1999
+; Setp 2004: boundary level have 0 values and not 1 (as it was
+; explained before in the header). see:
+; print, array_equal(niveau, total(level2mask(niveau),3))
+;
 ;-
 ;------------------------------------------------------------
@@ -45 +49 @@
 ;------------------------------------------------------------
    niveaux = litchamp(tab)
-   grille,maskterre,glam,gphi,gdep,nx,ny,nz,premierx,premiery,premierz,dernierx,derniery,dernierz
+   grille,maskterre, -1, -1, -1,nx,ny,nz,firstx,firsty,firstz,lastx,lasty,lastz
 ;---------------------------------------------------------------
 ; verification de la coherence entre la taille du tableau et le domaine definit par domdef
@@ -53 +57 @@
    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:niveaux=niveaux[premierx:dernierx, premiery:derniery]
+      taille[1] eq jpi and taille[2] eq jpj:niveaux=niveaux[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.')