source: trunk/SRC/Utilities/fitintobox.pro @ 327

Last change on this file since 327 was 327, checked in by pinsard, 17 years ago

modification of headers : mainly blanks around = sign for keywords in declaration of function and pro
  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 11.8 KB
RevLine 
[2]1;+
2;
[224]3; @file_comments
[133]4; Check that the input array has size and dimensions
[11]5; compatible with the domain that was defined with the previous call
[254]6; of <pro>domdef</pro>.
[2]7;
[224]8; @categories
[157]9; Compatibility
[224]10;
11; @param sz {in}{required}
12;
[133]13; @param nx {in}{required}
[2]14;
[133]15; @param ny {in}{required}
[2]16;
[133]17; @param nz {in}{required}
[2]18;
[133]19; @param jpi {in}{required}
[2]20;
[133]21; @param jpj {in}{required}
[2]22;
[133]23; @param jpk {in}{required}
[224]24;
[133]25; @param jpt {in}{required}
[2]26;
[224]27; @returns
28; -1
[11]29;
[224]30; @history
31; Sebastien Masson (smasson\@lodyc.jussieu.fr)
[133]32;          10 juin 2000.
33;          June 2005: S. Masson rewrite all.
34;
[224]35; @version
36; $Id$
[133]37;
38; @hidden
39;
[2]40;-
[11]41FUNCTION err_mess, sz, jpi, nx, jpj, ny, jpk, nz, jpt
[114]42;
43  compile_opt idl2, strictarrsubs
44;
[11]45  IF n_elements(sz EQ 1) THEN $
46    RETURN, report(['Error: ' $
47                  , 'the vector size (' + tostr(sz) + ') is incompatible' $
[238]48                  , 'with the domain dimensions ' $
[11]49                  , '[jpi/nx, jpj/ny, jpk/nz, jpt] = [' $
50                  + strtrim(jpi, 1) + '/' + strtrim(nx, 1) $
51                  + ', ' + strtrim(jpj, 1) + '/' + strtrim(ny, 1) $
52                  + ', ' + strtrim(jpk, 1) + '/' + strtrim(nz, 1) $
53                  + ', ' + strtrim(jpt, 1) + ']'], /simple)
54  RETURN, report(['Error: ' $
55                  , 'the array dimensions ' + tostr(sz) + ' are incompatible' $
[238]56                  , 'with the domain dimensions ' $
[11]57                  , '[jpi/nx, jpj/ny, jpk/nz, jpt] = [' $
58                  + strtrim(jpi, 1) + '/' + strtrim(nx, 1) $
59                  + ', ' + strtrim(jpj, 1) + '/' + strtrim(ny, 1) $
60                  + ', ' + strtrim(jpk, 1) + '/' + strtrim(nz, 1) $
61                  + ', ' + strtrim(jpt, 1) + ']'], /simple)
62END
[133]63;
64;+
65;
[224]66; @categories
[157]67; Compatibility
[133]68;
[163]69; @param field {in}{required}{type=array or struct}
[224]70; an array or a structure that can be read by the function
[231]71; <pro>litchamp</pro>
[224]72;
[231]73; @param nx {in}{optional}{default=define by <pro>grille</pro>}
[133]74;
[231]75; @param ny {in}{optional}{default=define by <pro>grille</pro>}
[133]76;
[231]77; @param nz {in}{optional}{default=define by <pro>grille</pro>}
[133]78;
[231]79; @param firstx {in}{optional}{default=define by <pro>grille</pro>}
[133]80;
[231]81; @param firsty {in}{optional}{default=define by <pro>grille</pro>}
[133]82;
[231]83; @param firstz {in}{optional}{default=define by <pro>grille</pro>}
[133]84;
[231]85; @param lastx {in}{optional}{default=define by <pro>grille</pro>}
[133]86;
[231]87; @param lasty {in}{optional}{default=define by <pro>grille</pro>}
[133]88;
[231]89; @param lastz {in}{optional}{default=define by <pro>grille</pro>}
[133]90;
[224]91; @keyword WDEPTH
[136]92; To specify that we are at W level
[133]93;
[224]94; @returns
[136]95; an array with dimensions matching the domain
96; or -1 if there is an error...
[133]97;
[238]98; @uses
[231]99; cm_4mesh
100; cm_4cal
[133]101;
[224]102; @examples
[136]103; IDL> help, fitintobox(findgen(jpi,jpj))
104; <Expression>    FLOAT     = Array[41, 3]
105; IDL> help, fitintobox(findgen(jpi,jpj,78))
[224]106; Error:
[136]107; the array dimensions [180,148,78] are incompatible
[238]108; with the domain dimensions
[136]109; [jpi/nx, jpj/ny, jpk/nz, jpt] = [180/41, 148/3, 31/31, 1]
110; <Expression>    INT       =       -1
[133]111;
[224]112; @history
113; Sebastien Masson (smasson\@lodyc.jussieu.fr)
[133]114;          10 juin 2000.
115;          June 2005: S. Masson rewrite all.
116;
[224]117; @version
118; $Id$
[238]119;
[133]120;-
[11]121FUNCTION fitintobox, field, nx, ny, nz, firstx, firsty $
[327]122                   , firstz, lastx, lasty, lastz, WDEPTH=wdepth
[114]123;
124  compile_opt idl2, strictarrsubs
125;
[11]126@cm_4mesh
127@cm_4cal
128  IF NOT keyword_set(key_forgetold) THEN BEGIN
129@updatenew
130  ENDIF
[2]131;---------------------
[11]132  arr = litchamp(field)
133  IF n_params() EQ 1 THEN grille, -1, -1, -1, -1, nx, ny, nz $
134    , firstx, firsty, firstz, lastx, lasty, lastz, WDEPTH = wdepth
135;--------------------------------------------------------------
136;--------------------------------------------------------------
137; case according the number of dimensions of the array
138;--------------------------------------------------------------
139;--------------------------------------------------------------
140  sz = size(arr)
141  case sz[0] of
142;--------------------------------------------------------------
143    0:BEGIN                     ; scalar
144;--------------------------------------------------------------
145      return, report('Error: scalar value = ' + strtrim(arr, 1), /simple)
146    END
147;--------------------------------------------------------------
148    1:BEGIN                     ; 1D arrays
149;--------------------------------------------------------------
150      CASE 1 OF
151; x arrays
[224]152        sz[1] EQ jpi                                                   :arr = (temporary(arr))[firstx:lastx                               ]
153        sz[1] EQ  nx                                                   :
[11]154; y arrays
[224]155        sz[1] EQ jpj                                                   :arr = (temporary(arr))[              firsty:lasty                 ]
156        sz[1] EQ  ny                                                   :
[11]157; z arrays
[224]158        sz[1] EQ jpk                                                   :arr = (temporary(arr))[                            firstz:lastz   ]
159        sz[1] EQ  nz                                                   :
[11]160; t arrays
161        sz[1] EQ jpt                                                   :
[224]162        ELSE:return, err_mess(sz[1], jpi, nx, jpj, ny, jpk, nz, jpt)
[11]163      ENDCASE
164    END
165;--------------------------------------------------------------
166    2:BEGIN                     ; 2D arrays
167;--------------------------------------------------------------
168      CASE 1 OF
169; xy arrays
[224]170        sz[1] EQ jpi AND sz[2] EQ jpj                                  :arr = (temporary(arr))[firstx:lastx, firsty:lasty                 ]
[11]171        sz[1] EQ jpi AND sz[2] EQ  ny                                  :arr = (temporary(arr))[firstx:lastx,            *                 ]
172        sz[1] EQ  nx AND sz[2] EQ jpj                                  :arr = (temporary(arr))[           *, firsty:lasty                 ]
173        sz[1] EQ  nx AND sz[2] EQ  ny                                  :arr = (temporary(arr))[           *,            *                 ]
174; x(y)z arrays
[224]175        sz[1] EQ jpi AND ny EQ 1      AND sz[2] EQ jpk                 :arr = (temporary(arr))[firstx:lastx, firstz:lastz   ]
[11]176        sz[1] EQ jpi AND ny EQ 1      AND sz[2] EQ  nz                 :arr = (temporary(arr))[firstx:lastx,            *   ]
177        sz[1] EQ  nx AND ny EQ 1      AND sz[2] EQ jpk                 :arr = (temporary(arr))[           *, firstz:lastz   ]
178        sz[1] EQ  nx AND ny EQ 1      AND sz[2] EQ  nz                 :
179; (x)yz arrays
[224]180        nx EQ 1      AND sz[1] EQ jpj AND sz[2] EQ jpk                 :arr = (temporary(arr))[              firsty:lasty, firstz:lastz   ]
[11]181        nx EQ 1      AND sz[1] EQ jpj AND sz[2] EQ  nz                 :arr = (temporary(arr))[              firsty:lasty,            *   ]
182        nx EQ 1      AND sz[1] EQ  ny AND sz[2] EQ jpk                 :arr = (temporary(arr))[                         *, firstz:lastz   ]
183        nx EQ 1      AND sz[1] EQ  ny AND sz[2] EQ  nz                 :
184; xt arrays
[224]185        sz[1] EQ jpi                                   AND sz[2] EQ jpt:arr = (temporary(arr))[firstx:lastx                            , *]
[11]186        sz[1] EQ  nx                                   AND sz[2] EQ jpt:
187; yt arrays
188                         sz[1] EQ jpj                  AND sz[2] EQ jpt:arr = (temporary(arr))[              firsty:lasty              , *]
189                         sz[1] EQ  ny                  AND sz[2] EQ jpt:
190; zt arrays
[224]191                                          sz[1] EQ jpk AND sz[2] EQ jpt:arr = (temporary(arr))[                            firstz:lastz, *]
[11]192                                          sz[1] EQ  nz AND sz[2] EQ jpt:
193        ELSE:return, err_mess(sz[1:2], jpi, nx, jpj, ny, jpk, nz, jpt)
194      ENDCASE
195    END
196;--------------------------------------------------------------
197    3:BEGIN                     ; 3D arrays
198;--------------------------------------------------------------
199      CASE 1 OF
200; xyz arrays
[224]201        sz[1] EQ jpi AND sz[2] EQ jpj AND sz[3] EQ jpk                 :arr = (temporary(arr))[firstx:lastx, firsty:lasty, firstz:lastz   ]
[11]202        sz[1] EQ jpi AND sz[2] EQ  ny AND sz[3] EQ jpk                 :arr = (temporary(arr))[firstx:lastx,            *, firstz:lastz   ]
203        sz[1] EQ  nx AND sz[2] EQ jpj AND sz[3] EQ jpk                 :arr = (temporary(arr))[           *, firsty:lasty, firstz:lastz   ]
204        sz[1] EQ  nx AND sz[2] EQ  ny AND sz[3] EQ jpk                 :arr = (temporary(arr))[           *,            *, firstz:lastz   ]
[224]205        sz[1] EQ jpi AND sz[2] EQ jpj AND sz[3] EQ  nz                 :arr = (temporary(arr))[firstx:lastx, firsty:lasty,            *   ]
[11]206        sz[1] EQ jpi AND sz[2] EQ  ny AND sz[3] EQ  nz                 :arr = (temporary(arr))[firstx:lastx,            *,            *   ]
207        sz[1] EQ  nx AND sz[2] EQ jpj AND sz[3] EQ  nz                 :arr = (temporary(arr))[           *, firsty:lasty,            *   ]
208        sz[1] EQ  nx AND sz[2] EQ  ny AND sz[3] EQ  nz                 :
209; xyt arrays
[224]210        sz[1] EQ jpi AND sz[2] EQ jpj                  AND sz[3] EQ jpt:arr = (temporary(arr))[firstx:lastx, firsty:lasty,               *]
[11]211        sz[1] EQ jpi AND sz[2] EQ  ny                  AND sz[3] EQ jpt:arr = (temporary(arr))[firstx:lastx,            *,               *]
212        sz[1] EQ  nx AND sz[2] EQ jpj                  AND sz[3] EQ jpt:arr = (temporary(arr))[           *, firsty:lasty,               *]
213        sz[1] EQ  nx AND sz[2] EQ  ny                  AND sz[3] EQ jpt:
214; (x)yzt arrays
[224]215        nx EQ 1      AND sz[1] EQ jpj AND sz[2] EQ jpk AND sz[3] EQ jpt:arr = (temporary(arr))[              firsty:lasty, firstz:lastz, *]
[11]216        nx EQ 1      AND sz[1] EQ jpj AND sz[2] EQ  nz AND sz[3] EQ jpt:arr = (temporary(arr))[              firsty:lasty,            *, *]
217        nx EQ 1      AND sz[1] EQ  ny AND sz[2] EQ jpk AND sz[3] EQ jpt:arr = (temporary(arr))[                         *, firstz:lastz, *]
218        nx EQ 1      AND sz[1] EQ  ny AND sz[2] EQ  nz AND sz[3] EQ jpt:
219; x(y)zt arrays
[224]220        sz[1] EQ jpi AND ny EQ 1      AND sz[2] EQ jpk AND sz[3] EQ jpt:arr = (temporary(arr))[firstx:lastx,               firstz:lastz, *]
[11]221        sz[1] EQ jpi AND ny EQ 1      AND sz[2] EQ  nz AND sz[3] EQ jpt:arr = (temporary(arr))[firstx:lastx,                          *, *]
222        sz[1] EQ  nx AND ny EQ 1      AND sz[2] EQ jpk AND sz[3] EQ jpt:arr = (temporary(arr))[           *,               firstz:lastz, *]
223        sz[1] EQ  nx AND ny EQ 1      AND sz[2] EQ  nz AND sz[3] EQ jpt:
224        ELSE:return, err_mess(sz[1:3], jpi, nx, jpj, ny, jpk, nz, jpt)
225      ENDCASE
226    END
227;--------------------------------------------------------------
228    4:BEGIN                     ; 4D arrays
229;--------------------------------------------------------------
230      CASE 1 OF
231; xyzt arrays
[224]232        sz[1] EQ jpi AND sz[2] EQ jpj AND sz[3] EQ jpk AND sz[4] EQ jpt:arr = (temporary(arr))[firstx:lastx, firsty:lasty, firstz:lastz, *]
[11]233        sz[1] EQ jpi AND sz[2] EQ  ny AND sz[3] EQ jpk AND sz[4] EQ jpt:arr = (temporary(arr))[firstx:lastx,            *, firstz:lastz, *]
234        sz[1] EQ  nx AND sz[2] EQ jpj AND sz[3] EQ jpk AND sz[4] EQ jpt:arr = (temporary(arr))[           *, firsty:lasty, firstz:lastz, *]
235        sz[1] EQ  nx AND sz[2] EQ  ny AND sz[3] EQ jpk AND sz[4] EQ jpt:arr = (temporary(arr))[           *,            *, firstz:lastz, *]
[224]236        sz[1] EQ jpi AND sz[2] EQ jpj AND sz[3] EQ  nz AND sz[4] EQ jpt:arr = (temporary(arr))[firstx:lastx, firsty:lasty,            *, *]
[11]237        sz[1] EQ jpi AND sz[2] EQ  ny AND sz[3] EQ  nz AND sz[4] EQ jpt:arr = (temporary(arr))[firstx:lastx,            *,            *, *]
238        sz[1] EQ  nx AND sz[2] EQ jpj AND sz[3] EQ  nz AND sz[4] EQ jpt:arr = (temporary(arr))[           *, firsty:lasty,            *, *]
239        sz[1] EQ  nx AND sz[2] EQ  ny AND sz[3] EQ  nz AND sz[4] EQ jpt:
240        ELSE:return, err_mess(sz[1:4], jpi, nx, jpj, ny, jpk, nz, jpt)
241      ENDCASE
242    END
243    ELSE:return, report('Error: fitintobox is managing arrays with a maximum of 4 dimensions', /simple)
244  ENDCASE
[2]245
[11]246  return, arr
[2]247end
Note: See TracBrowser for help on using the repository browser.