Changeset 48 for trunk/divfred.pro


Ignore:
Timestamp:
03/16/14 20:38:39 (10 years ago)
Author:
pinsard
Message:

fix thanks to coding rules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/divfred.pro

    r41 r48  
    2525; - Requires SAXO tools 
    2626; - les matrices u et v peuvent de 2 a 4 dimensions. 
    27 ; attention pour distinguer les différentes configurations de u et v 
     27; attention pour distinguer les différentes configurations de u et v 
    2828; (xy, xyz, xyt, xyzt), on regarde la variable du common 
    2929;        -time qui contient le calendrier en jour julien d''IDL auquel 
    3030;        se rapportent u et v ansi que la variable 
    31 ;        -jpt qui est le nombre de pas de temps à considérer dans time. 
    32 ; les tableaux u et v sont découpés sur le même domaine 
    33 ; géographique. A cause du décalage des grilles T, U, V et F il est 
    34 ; possible que ces 2 tableaux n''aient pas la même taille et se 
    35 ; reportent à des indices différents. Si tel est le cas les tableaux 
    36 ; sont redécoupés sur les indices qu'ils ont en commun et le dommaine 
    37 ; est redéfini pour qu'il colle à ces indices communs. 
    38 ; pour éviter ces redécoupes utiliser le mot clé /memeindice dans 
     31;        -jpt qui est le nombre de pas de temps à considérer dans time. 
     32; les tableaux u et v sont découpés sur le même domaine 
     33; géographique. A cause du décalage des grilles T, U, V et F il est 
     34; possible que ces 2 tableaux n''aient pas la même taille et se 
     35; reportent à des indices différents. Si tel est le cas les tableaux 
     36; sont redécoupés sur les indices qu'ils ont en commun et le dommaine 
     37; est redéfini pour qu'il colle à ces indices communs. 
     38; pour éviter ces redécoupes utiliser le mot clé /memeindice dans 
    3939; <pro>domdef</pro> 
    4040; 
    41 ; les points sur le bord du dessin sont mis à !values.f_nan 
     41; les points sur le bord du dessin sont mis à !values.f_nan 
    4242; 
    4343; @todo 
    44 ; ++ pas fini de comprendre, tester (compatibilité saxo), adapter, commenter 
     44; ++ pas fini de comprendre, tester (compatibilité saxo), adapter, commenter 
    4545; nettoyer 
    4646; ++ a comparer et merger avec SAXO_DIR/ToBeReviewed/CALCULS/div.pro 
     
    111111               v = v[indicex[0]:indicex[0]+nx-1,indicey[0]:indicey[0]+ny-1, *] 
    112112            END 
    113             ELSE:BEGIN 
     113            ELSE :BEGIN 
    114114               zdiv = -1 
    115115               GOTO, sortie 
     
    142142; mise a !values.f_nan de la bordure 
    143143;------------------------------------------------------------ 
    144          if  NOT keyword_set(key_periodic) OR nx NE jpi then begin 
     144         if NOT keyword_set(key_periodic) OR nx NE jpi then begin 
    145145            zdiv(0, *, *) = !values.f_nan 
    146146            zdiv(nx-1, *, *) = !values.f_nan 
     
    161161         varunits = '1e6*s-1' 
    162162         domdef, (glamt[indice2d])[0, 0], (glamu[indice2d])[nx-1, 0],(gphit[indice2d])[0, 0], (gphiv[indice2d])[0, ny-1], vert1, vert2, grille = ['t'] 
    163          if keyword_set(direc) then  zdiv = moyenne(zdiv,direc,/nan, boite = boite) 
     163         if keyword_set(direc) then zdiv = moyenne(zdiv,direc,/nan, boite = boite) 
    164164      END 
    165165;---------------------------------------------------------------------------- 
     
    189189               v = v[indicex[0]:indicex[0]+nx-1,indicey[0]:indicey[0]+ny-1, *] 
    190190            END 
    191             ELSE:return,  -1 
     191            ELSE : return,  -1 
    192192         endcase 
    193193;------------------------------------------------------------ 
     
    216216; mise a !values.f_nan de la bordure 
    217217;------------------------------------------------------------ 
    218          if  NOT keyword_set(key_periodic) OR nx NE jpi then begin 
     218         if NOT keyword_set(key_periodic) OR nx NE jpi then begin 
    219219            zdiv(0, *, *) = !values.f_nan 
    220220            zdiv(nx-1, *, *) = !values.f_nan 
     
    232232         varunits = '1e6*s-1' 
    233233         domdef, (glamt[indice2d])[0, 0], (glamu[indice2d])[nx-1, 0],(gphit[indice2d])[0, 0], (gphiv[indice2d])[0, ny-1], vert1, vert2, grille = ['t'] 
    234          if keyword_set(direc) then  zdiv = grossemoyenne(zdiv,direc,/nan, boite = boite) 
     234         if keyword_set(direc) then zdiv = grossemoyenne(zdiv,direc,/nan, boite = boite) 
    235235      END 
    236236;---------------------------------------------------------------------------- 
     
    247247;---------------------------------------------------------------------------- 
    248248;---------------------------------------------------------------------------- 
    249       ELSE:BEGIN                ;xy 
     249      ELSE : BEGIN                ;xy 
    250250         indice3d = lindgen(jpi, jpj, jpk) 
    251251         indice3d = indice3d[indicex[0]:indicex[0]+nx-1,indicey[0]:indicey[0]+ny-1, firstzt] 
     
    273273               v = v[indice2d] 
    274274            END 
    275             ELSE:return,  -1 
     275            ELSE : return,  -1 
    276276         endcase 
    277277;------------------------------------------------------------ 
     
    289289; mise a !values.f_nan de la bordure 
    290290;------------------------------------------------------------ 
    291          if  NOT keyword_set(key_periodic) OR nx NE jpi then begin 
     291         if NOT keyword_set(key_periodic) OR nx NE jpi then begin 
    292292            zdiv(0, *) = !values.f_nan 
    293293            zdiv(nx-1, *) = !values.f_nan 
     
    308308         varunits = '1e6*s-1' 
    309309         domdef, (glamt[indice2d])[0, 0], (glamu[indice2d])[nx-1, 0],(gphit[indice2d])[0, 0], (gphiv[indice2d])[0, ny-1], vert1, vert2, grille = ['t'] 
    310          if keyword_set(direc) then  zdiv = moyenne(zdiv,direc,/nan, boite = boite) 
     310         if keyword_set(direc) then zdiv = moyenne(zdiv,direc,/nan, boite = boite) 
    311311 
    312312      END 
Note: See TracChangeset for help on using the changeset viewer.