00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117
00118
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140
00141
00142
00143
00144
00145
00146
00147
00148
00149
00150
00151
00152
00153
00154
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207
00208
00209
00210
00211
00212
00213
00214
00215
00216 #ifdef HAVE_CONFIG_H
00217 # include <config.h>
00218 #endif
00219
00220
00227
00228
00229
00230
00231
00232 #include <uves_mbias_impl.h>
00233
00234 #include <uves_utils.h>
00235 #include <uves_corrbadpix.h>
00236 #include <uves_parameters.h>
00237 #include <uves.h>
00238 #include <uves_dfs.h>
00239 #include <uves_pfits.h>
00240 #include <uves_qclog.h>
00241 #include <uves_recipe.h>
00242 #include <uves_utils_wrappers.h>
00243 #include <uves_error.h>
00244 #include <irplib_mkmaster.h>
00245 #include <uves_msg.h>
00246
00247
00248 #include <cpl.h>
00249 #include <float.h>
00250 #include <ctype.h>
00251 #include <string.h>
00252
00253
00254
00255
00256 static void uves_mbias_qclog(const cpl_imagelist* raw_imgs,
00257 uves_propertylist **raw_headers,
00258 enum uves_chip chip,
00259 const cpl_image* mbia,
00260
00261
00262 cpl_table* qclog
00263 );
00264
00265 static void
00266 uves_mbias_qc_ron_raw(const cpl_image* rbia,
00267 enum uves_chip chip,
00268 const int x_cent_s,
00269 const int x_cent_e,
00270 const int y_cent_s,
00271 const int y_cent_e,
00272 cpl_table* qclog);
00273
00274 static int
00275 uves_mbias_define_parameters(cpl_parameterlist *parameters);
00276
00277
00278
00279
00280 #define cpl_plugin_get_info uves_mbias_get_info
00281 UVES_RECIPE_DEFINE(
00282 UVES_MBIAS_ID, UVES_MBIAS_DOM,
00283
00284
00285 uves_mbias_define_parameters,
00286 "Jonas M. Larsen", "cpl@eso.org",
00287 "Creates the master bias frame",
00288 "This recipe creates a master bias frame by computing the median of all input\n"
00289 "bias frames. All input frames must have same tag and size and must be either\n"
00290 "BIAS_BLUE or BIAS_RED.\n"
00291 "On blue input the recipe computes one master bias frame; on red input a \n"
00292 "master bias frame for each chip is produced. The average, standard deviation\n"
00293 "and median of the master bias image(s) are written to the FITS header(s)");
00294
00295
00296
00297
00298
00301
00308
00309 int uves_mbias_define_parameters_body(cpl_parameterlist *parameters,
00310 const char *recipe_id)
00311 {
00312
00313
00314
00315
00316 if (uves_define_global_parameters(parameters) != CPL_ERROR_NONE)
00317 {
00318 return -1;
00319 }
00320
00321
00322
00323
00324
00325 if (uves_corr_traps_define_parameters(parameters,recipe_id)
00326 != CPL_ERROR_NONE)
00327 {
00328 return -1;
00329 }
00330
00331
00332
00333
00334
00335 if (uves_master_stack_define_parameters(parameters,recipe_id)
00336 != CPL_ERROR_NONE)
00337 {
00338 return -1;
00339 }
00340
00341 return (cpl_error_get_code() != CPL_ERROR_NONE);
00342 }
00343
00344
00345
00346
00352
00353 static int
00354 uves_mbias_define_parameters(cpl_parameterlist *parameters)
00355 {
00356 return uves_mbias_define_parameters_body(parameters, make_str(UVES_MBIAS_ID));
00357 }
00358
00359
00360
00361
00362
00363
00380
00381 static cpl_image *
00382 uves_mbias_process_chip(const cpl_imagelist *raw_images,
00383 uves_propertylist **raw_headers,
00384 uves_propertylist *mbias_header,
00385 int binx, int biny,
00386 enum uves_chip chip,
00387 bool CLEAN_TRAPS,
00388 const char* STACK_METHOD,
00389 const double STACK_KLOW,
00390 const double STACK_KHIGH,
00391 const int STACK_NITER)
00392 {
00393 cpl_image *master_bias = NULL;
00394 double exposure_time = 0;
00395 int badpixels_cleaned;
00396 int i;
00397 int nraw=0;
00398 bool red_ccd_is_new=false;
00399 cpl_vector* bias_levels=NULL;
00400 double bias_mean=0;
00401
00402 uves_msg("Calculating master bias...");
00403
00404 check_nomsg(red_ccd_is_new=uves_ccd_is_new(raw_headers[0]));
00405
00406 if(strcmp(STACK_METHOD,"MEDIAN")==0) {
00407 uves_msg("method median");
00408 master_bias=irplib_mkmaster_median(raw_images,5.,5,1.e-5);
00409 } else {
00410 uves_msg("method mean");
00411 master_bias=irplib_mkmaster_mean(raw_images,5.,5,1.e-5,STACK_KLOW,STACK_KHIGH,STACK_NITER);
00412 }
00413
00414
00415 exposure_time = 0;
00416 nraw=cpl_imagelist_get_size(raw_images);
00417 for (i = 0; i < nraw; i++)
00418 {
00419 check( exposure_time += uves_pfits_get_exptime(raw_headers[i]),
00420 "Error reading exposure time");
00421 }
00422 exposure_time /= nraw;
00423
00424 check( uves_pfits_set_exptime(mbias_header, exposure_time),
00425 "Error setting master bias exposure time");
00426 if(CLEAN_TRAPS) {
00427 check( badpixels_cleaned =
00428 uves_correct_badpix_all(master_bias, mbias_header,
00429 chip, binx, biny, false,red_ccd_is_new),
00430 "Error replacing bad pixels");
00431
00432 uves_msg("%d bad pixels replaced", badpixels_cleaned);
00433 }
00434
00435 cleanup:
00436 if (cpl_error_get_code() != CPL_ERROR_NONE)
00437 {
00438 uves_free_image(&master_bias);
00439 }
00440
00441 return master_bias;
00442 }
00443
00444
00451
00452 static void
00453 UVES_CONCAT2X(UVES_MBIAS_ID,exe)(cpl_frameset *frames,
00454 const cpl_parameterlist *parameters,
00455 const char *starttime)
00456 {
00457 uves_mbias_exe_body(frames, parameters, starttime, make_str(UVES_MBIAS_ID));
00458 return;
00459 }
00460
00461
00472
00473 void
00474 uves_mbias_exe_body(cpl_frameset *frames,
00475 const cpl_parameterlist *parameters,
00476 const char *starttime,
00477 const char *recipe_id)
00478 {
00479
00480 cpl_imagelist *raw_images[2] = {NULL, NULL};
00481 uves_propertylist **raw_headers[2] = {NULL, NULL};
00482
00483 cpl_table* qclog[2] = {NULL, NULL};
00484
00485
00486 uves_propertylist *product_header[2] = {NULL, NULL};
00487 cpl_image *master_bias = NULL;
00488 cpl_stats *mbias_stats = NULL;
00489
00490
00491 char *product_filename = NULL;
00492 bool blue;
00493 enum uves_chip chip;
00494 int binx, biny;
00495 const char* PROCESS_CHIP=NULL;
00496 bool CLEAN_TRAPS;
00497 int raw_index = 0;
00498 int i=0;
00499
00500 const char* STACK_METHOD=NULL;
00501 double STACK_KLOW=0;
00502 double STACK_KHIGH=0;
00503 int STACK_NITER=0;
00504
00505
00506
00507 check( uves_load_raw_imagelist(frames,
00508 false,
00509 UVES_BIAS(true), UVES_BIAS(false),
00510 CPL_TYPE_DOUBLE,
00511 raw_images, raw_headers, product_header,
00512 &blue), "Error loading raw frames");
00513
00514
00515 check( binx = uves_pfits_get_binx(raw_headers[0][0]),
00516 "Could not get raw frame x-binning");
00517 check( biny = uves_pfits_get_biny(raw_headers[0][0]),
00518 "Could not get raw frame y-binning");
00519 check( uves_get_parameter(parameters, NULL, "uves", "process_chip", CPL_TYPE_STRING, &PROCESS_CHIP),
00520 "Could not read parameter");
00521 uves_string_toupper((char*)PROCESS_CHIP);
00522
00523 check( uves_get_parameter(parameters, NULL, recipe_id, "clean_traps", CPL_TYPE_BOOL, &CLEAN_TRAPS),
00524 "Could not read parameter");
00525
00526 check( uves_get_parameter(parameters, NULL, recipe_id, "stack_method", CPL_TYPE_STRING, &STACK_METHOD),
00527 "Could not read parameter");
00528 uves_string_toupper((char*)STACK_METHOD);
00529
00530 check( uves_get_parameter(parameters, NULL, recipe_id, "klow", CPL_TYPE_DOUBLE, &STACK_KLOW),
00531 "Could not read parameter");
00532 check( uves_get_parameter(parameters, NULL, recipe_id, "khigh", CPL_TYPE_DOUBLE, &STACK_KHIGH),
00533 "Could not read parameter");
00534 check( uves_get_parameter(parameters, NULL, recipe_id, "niter", CPL_TYPE_INT, &STACK_NITER),
00535 "Could not read parameter");
00536
00537
00538 for (chip = uves_chip_get_first(blue);
00539 chip != UVES_CHIP_INVALID;
00540 chip = uves_chip_get_next(chip))
00541 {
00542 if(strcmp(PROCESS_CHIP,"REDU") == 0) {
00543 chip = uves_chip_get_next(chip);
00544 }
00545 raw_index = uves_chip_get_index(chip);
00546
00547 uves_msg("Processing %s chip",
00548 uves_chip_tostring_upper(chip));
00549
00550 uves_msg_debug("Binning = %dx%d", binx, biny);
00551
00552
00553 uves_free_image(&master_bias);
00554 check( master_bias = uves_mbias_process_chip(raw_images[raw_index],
00555 raw_headers[raw_index],
00556 product_header[raw_index],
00557 binx, biny,
00558 chip,CLEAN_TRAPS,
00559 STACK_METHOD,
00560 STACK_KLOW,
00561 STACK_KHIGH,
00562 STACK_NITER),
00563 "Error processing chip");
00564
00565
00566 cpl_free(product_filename);
00567 check( product_filename = uves_masterbias_filename(chip),
00568 "Error getting filename");
00569
00570
00571 uves_msg("Calculating QC parameters");
00572 uves_qclog_delete(&qclog[0]);
00573 qclog[0] = uves_qclog_init(raw_headers[raw_index][0], chip);
00574 check(uves_mbias_qclog(raw_images[raw_index],
00575 raw_headers[raw_index],
00576 chip,
00577 master_bias,
00578
00579 qclog[0]),"error computing qclog");
00580
00581
00582 uves_msg("Saving product...");
00583
00584 check( uves_frameset_insert(frames,
00585 master_bias,
00586 CPL_FRAME_GROUP_PRODUCT,
00587 CPL_FRAME_TYPE_IMAGE,
00588 CPL_FRAME_LEVEL_INTERMEDIATE,
00589 product_filename,
00590 UVES_MASTER_BIAS(chip),
00591 raw_headers[raw_index][0],
00592 product_header[raw_index],
00593 NULL,
00594 parameters,
00595 recipe_id,
00596 PACKAGE "/" PACKAGE_VERSION,qclog,
00597 starttime, true,
00598 UVES_ALL_STATS),
00599 "Could not add master bias %s to frameset", product_filename);
00600 uves_qclog_delete(&qclog[0]);
00601 uves_msg("Master bias '%s' added to frameset", product_filename);
00602
00603 if(strcmp(PROCESS_CHIP,"REDL") == 0) {
00604 chip = uves_chip_get_next(chip);
00605 }
00606
00607
00608 }
00609
00610 cleanup:
00611
00612 if (raw_images[0] != NULL)
00613 {
00614
00615 for (i = 0; i < cpl_imagelist_get_size(raw_images[0]); i++)
00616 {
00617 if (raw_headers[0] != NULL) uves_free_propertylist(&raw_headers[0][i]);
00618 if (raw_headers[1] != NULL) uves_free_propertylist(&raw_headers[1][i]);
00619 }
00620 cpl_free(raw_headers[0]); raw_headers[0] = NULL;
00621 cpl_free(raw_headers[1]); raw_headers[1] = NULL;
00622 }
00623 uves_free_imagelist(&raw_images[0]);
00624 uves_free_imagelist(&raw_images[1]);
00625
00626
00627 uves_qclog_delete(&qclog[0]);
00628 uves_free_image(&master_bias);
00629 uves_free_propertylist(&product_header[0]);
00630 uves_free_propertylist(&product_header[1]);
00631 cpl_free(product_filename);
00632 uves_free_stats(&mbias_stats);
00633
00634 return;
00635 }
00636
00637
00638
00639 static int
00640 count_good(const cpl_image *image)
00641 {
00642 return
00643 cpl_image_get_size_x(image) * cpl_image_get_size_y(image) -
00644 cpl_image_count_rejected(image);
00645 }
00646
00653
00654 static void
00655 reject_lo_hi(cpl_image *image, double min, double max)
00656 {
00657 cpl_mask *mask_lo = NULL;
00658 cpl_mask *mask_hi = NULL;
00659
00660 mask_lo = cpl_mask_threshold_image_create(image, -DBL_MAX, min);
00661 mask_hi = cpl_mask_threshold_image_create(image, max, DBL_MAX);
00662 assure_mem( mask_lo );
00663 assure_mem( mask_hi );
00664
00665 cpl_mask_or(mask_lo, mask_hi);
00666
00667 cpl_image_reject_from_mask(image, mask_lo);
00668
00669 cleanup:
00670 uves_free_mask(&mask_lo);
00671 uves_free_mask(&mask_hi);
00672 return;
00673 }
00674
00685 static void uves_mbias_qclog(const cpl_imagelist* raw_imgs,
00686 uves_propertylist **raw_headers,
00687 enum uves_chip chip,
00688 const cpl_image* mbia,
00689
00690
00691 cpl_table* qclog
00692 )
00693 {
00694 int nx_pix= 0;
00695 int ny_pix= 0;
00696
00697 int sample_x= 100;
00698 int sample_y= 100;
00699 int x_cent_s= 0;
00700 int x_cent_e= 0;
00701 int y_cent_s= 0;
00702 int y_cent_e= 0;
00703
00704
00705
00706
00707
00708 double upp_threshold= 0.0;
00709 double low_threshold= 0.0;
00710 double extra=0.1;
00711 double qc_ron_master= 0.0;
00712
00713 double master_median=0.0;
00714 int pn= 0;
00715
00716
00717 double min=0.0;
00718 double max=0.0;
00719 double struct_col=0.0;
00720 double struct_row=0.0;
00721
00722 double time_s=+9999999.0;
00723 double time_e=-9999999.0;
00724 int nraw=0;
00725 double qc_duty_cycle=0.;
00726 double exposure_time=0;
00727 int i=0;
00728 char key_name[80];
00729
00730 const cpl_image* rbia=NULL;
00731 cpl_image* tima=NULL;
00732 cpl_image* avg_col=NULL;
00733 cpl_image* avg_row=NULL;
00734
00735 uves_qclog_add_string(qclog,
00736 "QC TEST1 ID",
00737 "Test-on-Master-Bias",
00738 "Name of QC test",
00739 "%s");
00740
00741 uves_msg("Computing duty cycle...");
00742
00743
00744 exposure_time = 0;
00745 nraw = cpl_imagelist_get_size(raw_imgs);
00746 check_nomsg(uves_qclog_add_int(qclog,
00747 "PRO DATANCOM",
00748 nraw,
00749 "Number of frames combined",
00750 "%d"));
00751
00752
00753 for (i = 0; i < nraw; i++)
00754 {
00755 check( exposure_time = uves_pfits_get_mjdobs(raw_headers[i]),
00756 "Error reading exposure time");
00757 if(exposure_time >= time_e) time_e = exposure_time;
00758 if(exposure_time <= time_s) time_s = exposure_time;
00759 }
00760 if(nraw > 1) {
00761 qc_duty_cycle = (time_e-time_s)/ (nraw-1);
00762 }
00763 else
00764 {
00765 qc_duty_cycle = 0;
00766 }
00767
00768 check_nomsg(uves_qclog_add_double(qclog,
00769 "QC DUTYCYCL",
00770 qc_duty_cycle,
00771 "Time to store a frame",
00772 "%.5e"));
00773
00774
00775
00776
00777
00778
00779
00780
00781
00782
00783
00784 nx_pix = cpl_image_get_size_x(mbia);
00785 ny_pix = cpl_image_get_size_y(mbia);
00786
00787
00788 x_cent_s = (nx_pix - sample_x)/2;
00789 x_cent_e = (nx_pix + sample_x)/2;
00790 y_cent_s = (ny_pix - sample_y)/2;
00791 y_cent_e = (ny_pix + sample_y)/2;
00792
00793
00794 check_nomsg(upp_threshold =
00795 cpl_image_get_median_window(mbia,
00796 x_cent_s,
00797 y_cent_s,
00798 x_cent_e,
00799 y_cent_e)*(1 + extra));
00800 check_nomsg(low_threshold =
00801 cpl_image_get_median_window(mbia,
00802 x_cent_s,
00803 y_cent_s,
00804 x_cent_e,
00805 y_cent_e)*(1 - extra));
00806
00807
00808
00809
00810 check_nomsg(pn = PORT_ID(chip));
00811 uves_msg_debug("Port number = %d", pn);
00812
00813 rbia = cpl_imagelist_get_const(raw_imgs,0);
00814 check_nomsg(uves_mbias_qc_ron_raw(rbia, chip,
00815 x_cent_s,x_cent_e,y_cent_s,y_cent_e,qclog));
00816
00817
00818
00819
00820
00821 check_nomsg(tima=cpl_image_duplicate(mbia));
00822
00823 check_nomsg( reject_lo_hi(tima, low_threshold, upp_threshold) );
00824 if (count_good(tima) >= 2)
00825 {
00826 check_nomsg(master_median = cpl_image_get_median(tima));
00827 check_nomsg(qc_ron_master = cpl_image_get_stdev(tima));
00828 }
00829 else
00830 {
00831 master_median = -1;
00832 qc_ron_master = -1;
00833 uves_msg_warning("Only %d good pixels in image. Setting QC parameters to -1",
00834 count_good(tima));
00835 }
00836 uves_free_image(&tima);
00837
00838 check_nomsg(uves_qclog_add_double(qclog,
00839 "PRO DATAMED",
00840 master_median,
00841 "Median of pixel values",
00842 "%7.3f"));
00843
00844 sprintf(key_name, "QC OUT%d RON MASTER", pn);
00845 check_nomsg(uves_qclog_add_double(qclog,
00846 key_name,
00847 qc_ron_master,
00848 "Read noise frame in ADU",
00849 "%8.4f"));
00850
00851
00852
00853
00854
00855
00856
00857
00858
00859
00860
00861
00862
00863
00864 check_nomsg(tima=cpl_image_duplicate(mbia));
00865 if (chip != UVES_CHIP_BLUE) {
00866
00867
00868
00869 check_nomsg(cpl_image_threshold(tima,
00870 -DBL_MAX,300,
00871 -DBL_MAX,300));
00872 }
00873
00874
00875 check_nomsg(avg_col = cpl_image_collapse_create(tima,1));
00876 check_nomsg(cpl_image_divide_scalar(avg_col,cpl_image_get_size_x(tima)));
00877
00878
00879 min = cpl_image_get_mean(avg_col) - 2;
00880 max = cpl_image_get_mean(avg_col) + 2;
00881
00882
00883
00884
00885 check_nomsg( reject_lo_hi(avg_col, min, max) );
00886 if (count_good(avg_col) >= 2)
00887 {
00888 check_nomsg(struct_col = cpl_image_get_stdev(avg_col));
00889 }
00890 else
00891 {
00892 struct_col = -1;
00893 uves_msg_warning("Only %d good pixels in image. Setting QC parameter to -1",
00894 count_good(avg_col));
00895 }
00896
00897 sprintf(key_name,"%s%d%s","QC OUT",pn," STRUCTY");
00898 check_nomsg(uves_qclog_add_double(qclog,
00899 key_name,
00900 struct_col,
00901 "structure in Y (bias slope)",
00902 "%8.4f"));
00903
00904
00905
00906 check_nomsg(avg_row = cpl_image_collapse_create(tima,0));
00907 check_nomsg(cpl_image_divide_scalar(avg_row,cpl_image_get_size_y(tima)));
00908
00909
00910 min = cpl_image_get_mean(avg_row) - 2;
00911 max = cpl_image_get_mean(avg_row) + 2;
00912
00913
00914
00915
00916 check_nomsg( reject_lo_hi(avg_row, min, max) );
00917 if (count_good(avg_row) >= 2)
00918 {
00919 check_nomsg(struct_row = cpl_image_get_stdev(avg_row));
00920 }
00921 else
00922 {
00923 struct_row = -1;
00924 uves_msg_warning("Only %d good pixels in image. Setting QC parameter to -1",
00925 count_good(avg_row));
00926 }
00927
00928
00929 sprintf(key_name,"%s%d%s","QC OUT",pn," STRUCTX");
00930 check_nomsg(uves_qclog_add_double(qclog,
00931 key_name,
00932 struct_row,
00933 "structure in X (bias slope)",
00934 "%8.4f"));
00935
00936
00937
00938
00939 cleanup:
00940 uves_free_image(&avg_col);
00941 uves_free_image(&avg_row);
00942 uves_free_image(&tima);
00943
00944 return;
00945
00946 }
00960 static void
00961 uves_mbias_qc_ron_raw(const cpl_image* rbia,
00962 enum uves_chip chip,
00963 const int x_cent_s,
00964 const int x_cent_e,
00965 const int y_cent_s,
00966 const int y_cent_e,
00967 cpl_table* qclog)
00968 {
00969
00970 double qc_ron_raw=0.0;
00971 double upp_threshold=0.0;
00972 double low_threshold=0.0;
00973 double extra=0.1;
00974 char key_name[80];
00975 int pn=0;
00976 cpl_image* tima=NULL;
00977
00978
00979
00980
00981
00982 check_nomsg(upp_threshold =
00983 cpl_image_get_median_window(rbia,
00984 x_cent_s,
00985 y_cent_s,
00986 x_cent_e,
00987 y_cent_e)*(1 + extra));
00988
00989 check_nomsg(low_threshold =
00990 cpl_image_get_median_window(rbia,
00991 x_cent_s,
00992 y_cent_s,
00993 x_cent_e,
00994 y_cent_e)*(1 - extra));
00995
00996
00997
00998
00999
01000
01001 check_nomsg(tima=cpl_image_duplicate(rbia));
01002
01003 check_nomsg( reject_lo_hi(tima, low_threshold, upp_threshold) );
01004 if (count_good(tima) >= 2)
01005 {
01006 check_nomsg(qc_ron_raw = cpl_image_get_stdev(tima));
01007 }
01008 else
01009 {
01010 qc_ron_raw = -1;
01011 uves_msg_warning("Only %d good pixels in image. Setting QC parameter to -1",
01012 count_good(tima));
01013 }
01014
01015
01016
01017
01018
01019 check_nomsg(pn = PORT_ID(chip));
01020
01021
01022 sprintf(key_name,"%s%d%s","QC OUT",pn," RON RAW");
01023 check_nomsg(uves_qclog_add_double(qclog,
01024 key_name,
01025 qc_ron_raw,
01026 "Read noise frame in ADU",
01027 "%8.4f"));
01028
01029 cleanup:
01030 uves_free_image(&tima);
01031 return;
01032 }
01033