Ignore:
Timestamp:
Jan 8, 2017, 4:44:10 PM (8 years ago)
Author:
oboucher
Message:

Changes to StratAer? module, including setting up a new output file
hence dimension of output metadata is increased to 10 items.

Location:
LMDZ5/trunk/libf/phylmd/Dust
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/Dust/spla_output_dat.h

    r2630 r2752  
    33
    44!Dust emission module
    5   type(ctrl_out),save :: o_m1dflux      = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    6   'm1dflux','m1dflux','', (/ ('', i=1, 9) /))
    7   type(ctrl_out),save :: o_m2dflux      = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    8   'm2dflux','m2dflux','', (/ ('', i=1, 9) /))
    9   type(ctrl_out),save :: o_m3dflux      = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    10   'm3dflux','m3dflux','', (/ ('', i=1, 9) /))
     5  type(ctrl_out),save :: o_m1dflux      = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     6  'm1dflux','m1dflux','', (/ ('', i=1, 10) /))
     7  type(ctrl_out),save :: o_m2dflux      = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     8  'm2dflux','m2dflux','', (/ ('', i=1, 10) /))
     9  type(ctrl_out),save :: o_m3dflux      = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     10  'm3dflux','m3dflux','', (/ ('', i=1, 10) /))
    1111
    1212! traceur_spl
    13   type(ctrl_out),save :: o_taue550    = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    14   'taue550','Tau ext 550','', (/ ('', i=1, 9) /))
    15   type(ctrl_out),save :: o_taue670     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    16   'taue670','Tau ext 670','', (/ ('', i=1, 9) /))
    17   type(ctrl_out),save :: o_taue865     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    18   'taue865','Tau ext 865','', (/ ('', i=1, 9) /))
    19   type(ctrl_out),save :: o_taue550_tr2     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    20   'taue550_tr2','Tau ext 550tr2','', (/ ('', i=1, 9) /))
    21   type(ctrl_out),save :: o_taue670_tr2     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    22   'taue670_tr2','Tau ext 670tr2','', (/ ('', i=1, 9) /))
    23   type(ctrl_out),save :: o_taue865_tr2     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    24   'taue865_tr2','Tau ext 865tr2','', (/ ('', i=1, 9) /))
    25   type(ctrl_out),save :: o_taue550_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    26   'taue550_ss','Tau ext 550ss','', (/ ('', i=1, 9) /))
    27   type(ctrl_out),save :: o_taue670_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    28   'taue670_ss','Tau ext 670ss','', (/ ('', i=1, 9) /))
    29   type(ctrl_out),save :: o_taue865_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    30   'taue865_ss','Tau ext 865ss','', (/ ('', i=1, 9) /))
    31   type(ctrl_out),save :: o_taue550_dust     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    32   'taue550_dust','Tau ext 550dust','', (/ ('', i=1, 9) /))
    33   type(ctrl_out),save :: o_taue670_dust     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    34   'taue670_dust','Tau ext 670dust','', (/ ('', i=1, 9) /))
    35   type(ctrl_out),save :: o_taue865_dust     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    36   'taue865_dust','Tau ext 865dust','', (/ ('', i=1, 9) /))
    37   type(ctrl_out),save :: o_taue550_dustsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    38   'taue550_dustsco','Tau ext 550dustsco','', (/ ('', i=1, 9) /))
    39   type(ctrl_out),save :: o_taue670_dustsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    40   'taue670_dustsco','Tau ext 670dustsco','', (/ ('', i=1, 9) /))
    41   type(ctrl_out),save :: o_taue865_dustsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    42   'taue865_dustsco','Tau ext 865dustsco','', (/ ('', i=1, 9) /))
    43 
    44 
    45   type(ctrl_out),save :: o_taue550_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     13  type(ctrl_out),save :: o_taue550    = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     14  'taue550','Tau ext 550','', (/ ('', i=1, 10) /))
     15  type(ctrl_out),save :: o_taue670     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     16  'taue670','Tau ext 670','', (/ ('', i=1, 10) /))
     17  type(ctrl_out),save :: o_taue865     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     18  'taue865','Tau ext 865','', (/ ('', i=1, 10) /))
     19  type(ctrl_out),save :: o_taue550_tr2     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     20  'taue550_tr2','Tau ext 550tr2','', (/ ('', i=1, 10) /))
     21  type(ctrl_out),save :: o_taue670_tr2     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     22  'taue670_tr2','Tau ext 670tr2','', (/ ('', i=1, 10) /))
     23  type(ctrl_out),save :: o_taue865_tr2     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     24  'taue865_tr2','Tau ext 865tr2','', (/ ('', i=1, 10) /))
     25  type(ctrl_out),save :: o_taue550_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     26  'taue550_ss','Tau ext 550ss','', (/ ('', i=1, 10) /))
     27  type(ctrl_out),save :: o_taue670_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     28  'taue670_ss','Tau ext 670ss','', (/ ('', i=1, 10) /))
     29  type(ctrl_out),save :: o_taue865_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     30  'taue865_ss','Tau ext 865ss','', (/ ('', i=1, 10) /))
     31  type(ctrl_out),save :: o_taue550_dust     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     32  'taue550_dust','Tau ext 550dust','', (/ ('', i=1, 10) /))
     33  type(ctrl_out),save :: o_taue670_dust     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     34  'taue670_dust','Tau ext 670dust','', (/ ('', i=1, 10) /))
     35  type(ctrl_out),save :: o_taue865_dust     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     36  'taue865_dust','Tau ext 865dust','', (/ ('', i=1, 10) /))
     37  type(ctrl_out),save :: o_taue550_dustsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     38  'taue550_dustsco','Tau ext 550dustsco','', (/ ('', i=1, 10) /))
     39  type(ctrl_out),save :: o_taue670_dustsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     40  'taue670_dustsco','Tau ext 670dustsco','', (/ ('', i=1, 10) /))
     41  type(ctrl_out),save :: o_taue865_dustsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     42  'taue865_dustsco','Tau ext 865dustsco','', (/ ('', i=1, 10) /))
     43
     44  type(ctrl_out),save :: o_taue550_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    4645  'taue550_terra','Tau ext 550 terra','', &
    47       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    48          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    49   type(ctrl_out),save :: o_taue550_fine_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     46      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     47         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     48  type(ctrl_out),save :: o_taue550_fine_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    5049  'taue550_fine_terra','Tau ext fine 550 terra','', &
    51       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    52          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    53   type(ctrl_out),save :: o_taue550_coss_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     50      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     51         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     52  type(ctrl_out),save :: o_taue550_coss_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    5453  'taue550_coss_terra','Tau ext coss 550 terra','', &
    55       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    56          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    57   type(ctrl_out),save :: o_taue550_codu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     54      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     55         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     56  type(ctrl_out),save :: o_taue550_codu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    5857  'taue550_codu_terra','Tau ext codu 550 terra','', &
    59       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    60          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    61   type(ctrl_out),save :: o_taue550_scdu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     58      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     59         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     60  type(ctrl_out),save :: o_taue550_scdu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    6261  'taue550_scdu_terra','Tau ext scdu 550 terra','', &
    63       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    64          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    65 
    66   type(ctrl_out),save :: o_taue670_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     62      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     63         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     64
     65  type(ctrl_out),save :: o_taue670_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    6766  'taue670_terra','Tau ext 670 terra','', &
    68       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    69          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    70   type(ctrl_out),save :: o_taue670_fine_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     67      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     68         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     69  type(ctrl_out),save :: o_taue670_fine_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    7170  'taue670_fine_terra','Tau ext fine 670 terra','', &
    72       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    73          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    74   type(ctrl_out),save :: o_taue670_coss_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     71      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     72         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     73  type(ctrl_out),save :: o_taue670_coss_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    7574  'taue670_coss_terra','Tau ext coss 670 terra','', &
    76       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    77          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    78   type(ctrl_out),save :: o_taue670_codu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     75      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     76         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     77  type(ctrl_out),save :: o_taue670_codu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    7978  'taue670_codu_terra','Tau ext codu 670 terra','', &
    80       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    81          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    82   type(ctrl_out),save :: o_taue670_scdu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     79      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     80         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     81  type(ctrl_out),save :: o_taue670_scdu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    8382  'taue670_scdu_terra','Tau ext scdu 670 terra','', &
    84       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    85          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    86 
    87   type(ctrl_out),save :: o_taue865_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     83      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     84         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     85
     86  type(ctrl_out),save :: o_taue865_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    8887  'taue865_terra','Tau ext 865 terra','', &
    89       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    90          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    91   type(ctrl_out),save :: o_taue865_fine_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     88      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     89         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     90  type(ctrl_out),save :: o_taue865_fine_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    9291  'taue865_fine_terra','Tau ext fine 865 terra','', &
    93       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    94          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    95   type(ctrl_out),save :: o_taue865_coss_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     92      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     93         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     94  type(ctrl_out),save :: o_taue865_coss_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    9695  'taue865_coss_terra','Tau ext coss 865 terra','', &
    97       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    98          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    99   type(ctrl_out),save :: o_taue865_codu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     96      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     97         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     98  type(ctrl_out),save :: o_taue865_codu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    10099  'taue865_codu_terra','Tau ext codu 865 terra','', &
    101       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    102          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    103   type(ctrl_out),save :: o_taue865_scdu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     100      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     101         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     102  type(ctrl_out),save :: o_taue865_scdu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    104103  'taue865_scdu_terra','Tau ext scdu 865 terra','', &
    105       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    106          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    107 
    108   type(ctrl_out),save :: o_taue550_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     104      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     105         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     106
     107  type(ctrl_out),save :: o_taue550_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    109108  'taue550_aqua','Tau ext 550 aqua','', &
    110       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    111          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    112   type(ctrl_out),save :: o_taue550_fine_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     109      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     110         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     111  type(ctrl_out),save :: o_taue550_fine_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    113112  'taue550_fine_aqua','Tau ext fine 550 aqua','', &
    114       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    115          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    116   type(ctrl_out),save :: o_taue550_coss_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     113      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     114         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     115  type(ctrl_out),save :: o_taue550_coss_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    117116  'taue550_coss_aqua','Tau ext coss 550 aqua','', &
    118       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    119          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    120   type(ctrl_out),save :: o_taue550_codu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     117      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     118         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     119  type(ctrl_out),save :: o_taue550_codu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    121120  'taue550_codu_aqua','Tau ext codu 550 aqua','', &
    122       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    123          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    124   type(ctrl_out),save :: o_taue550_scdu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     121      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     122         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     123  type(ctrl_out),save :: o_taue550_scdu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    125124  'taue550_scdu_aqua','Tau ext scdu 550 aqua','', &
    126       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    127          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    128 
    129   type(ctrl_out),save :: o_taue670_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     125      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     126         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     127
     128  type(ctrl_out),save :: o_taue670_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    130129  'taue670_aqua','Tau ext 670 aqua','', &
    131       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    132          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    133   type(ctrl_out),save :: o_taue670_fine_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     130      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     131         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     132  type(ctrl_out),save :: o_taue670_fine_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    134133  'taue670_fine_aqua','Tau ext fine 670 aqua','', &
    135       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    136          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    137   type(ctrl_out),save :: o_taue670_coss_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     134      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     135         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     136  type(ctrl_out),save :: o_taue670_coss_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    138137  'taue670_coss_aqua','Tau ext coss 670 aqua','', &
    139       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    140          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    141   type(ctrl_out),save :: o_taue670_codu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     138      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     139         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     140  type(ctrl_out),save :: o_taue670_codu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    142141  'taue670_codu_aqua','Tau ext codu 670 aqua','', &
    143       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    144          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    145   type(ctrl_out),save :: o_taue670_scdu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     142      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     143         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     144  type(ctrl_out),save :: o_taue670_scdu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    146145  'taue670_scdu_aqua','Tau ext scdu 670 aqua','', &
    147       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    148          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    149 
    150   type(ctrl_out),save :: o_taue865_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     146      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     147         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     148
     149  type(ctrl_out),save :: o_taue865_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    151150  'taue865_aqua','Tau ext 865 aqua','', &
    152       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    153          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    154   type(ctrl_out),save :: o_taue865_fine_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     151      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     152         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     153  type(ctrl_out),save :: o_taue865_fine_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    155154  'taue865_fine_aqua','Tau ext fine 865 aqua','', &
    156       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    157          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    158   type(ctrl_out),save :: o_taue865_coss_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     155      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     156         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     157  type(ctrl_out),save :: o_taue865_coss_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    159158  'taue865_coss_aqua','Tau ext coss 865 aqua','', &
    160       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    161          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    162   type(ctrl_out),save :: o_taue865_codu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     159      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     160         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     161  type(ctrl_out),save :: o_taue865_codu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    163162  'taue865_codu_aqua','Tau ext codu 865 aqua','', &
    164       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    165          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    166   type(ctrl_out),save :: o_taue865_scdu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     163      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     164         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     165  type(ctrl_out),save :: o_taue865_scdu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    167166  'taue865_scdu_aqua','Tau ext scdu 865 aqua','', &
    168       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    169          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    170 
    171 
    172 
    173 !  type(ctrl_out),save :: o_taue550_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     167      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
     168         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     169
     170!  type(ctrl_out),save :: o_taue550_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    174171!  'taue550_terra','Tau ext 550 terra','', &
    175 !      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     172!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
    176173!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    177174!
    178 !  type(ctrl_out),save :: o_taue670_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     175!  type(ctrl_out),save :: o_taue670_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    179176!  'taue670_aqua','Tau ext 670 aqua','', &
    180 !      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     177!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
    181178!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    182179!
    183 !  type(ctrl_out),save :: o_taue670_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     180!  type(ctrl_out),save :: o_taue670_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    184181!  'taue670_terra','Tau ext 670 terra','', &
    185 !      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     182!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
    186183!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    187184!
    188 !  type(ctrl_out),save :: o_taue865_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     185!  type(ctrl_out),save :: o_taue865_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    189186!  'taue865_aqua','Tau ext 865 aqua','', &
    190 !      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     187!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
    191188!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    192189!
    193 !  type(ctrl_out),save :: o_taue865_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     190!  type(ctrl_out),save :: o_taue865_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
    194191!  'taue865_terra','Tau ext 865 terra','', &
    195 !      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     192!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)',  &
    196193!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    197 
    198194       
    199   type(ctrl_out),save :: o_trm01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    200   'trm01','Burden PREC','', (/ ('', i=1, 9) /))
    201   type(ctrl_out),save :: o_trm02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    202   'trm02','Burden FINE','', (/ ('', i=1, 9) /))
    203   type(ctrl_out),save :: o_trm03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    204   'trm03','Burden COSS','', (/ ('', i=1, 9) /))
    205   type(ctrl_out),save :: o_trm04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    206   'trm04','Burden CODU','', (/ ('', i=1, 9) /))
    207   type(ctrl_out),save :: o_trm05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    208   'trm05','Burden SCDU','', (/ ('', i=1, 9) /))
    209 
    210   type(ctrl_out),save :: o_sconc01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    211   'sconc01','Surf. Conc. PREC','', (/ ('', i=1, 9) /))
    212   type(ctrl_out),save :: o_sconc02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    213   'sconc02','Surf. Conc. FINE','', (/ ('', i=1, 9) /))
    214   type(ctrl_out),save :: o_sconc03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    215   'sconc03','Surf. Conc. COSS','', (/ ('', i=1, 9) /))
    216   type(ctrl_out),save :: o_sconc04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    217   'sconc04','Surf. Conc. CODU','', (/ ('', i=1, 9) /))
    218   type(ctrl_out),save :: o_sconc05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    219   'sconc05','Surf. Conc. SCDU','', (/ ('', i=1, 9) /))
    220 
    221 
    222 
     195  type(ctrl_out),save :: o_trm01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     196  'trm01','Burden PREC','', (/ ('', i=1, 10) /))
     197  type(ctrl_out),save :: o_trm02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     198  'trm02','Burden FINE','', (/ ('', i=1, 10) /))
     199  type(ctrl_out),save :: o_trm03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     200  'trm03','Burden COSS','', (/ ('', i=1, 10) /))
     201  type(ctrl_out),save :: o_trm04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     202  'trm04','Burden CODU','', (/ ('', i=1, 10) /))
     203  type(ctrl_out),save :: o_trm05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     204  'trm05','Burden SCDU','', (/ ('', i=1, 10) /))
     205
     206  type(ctrl_out),save :: o_sconc01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     207  'sconc01','Surf. Conc. PREC','', (/ ('', i=1, 10) /))
     208  type(ctrl_out),save :: o_sconc02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     209  'sconc02','Surf. Conc. FINE','', (/ ('', i=1, 10) /))
     210  type(ctrl_out),save :: o_sconc03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     211  'sconc03','Surf. Conc. COSS','', (/ ('', i=1, 10) /))
     212  type(ctrl_out),save :: o_sconc04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     213  'sconc04','Surf. Conc. CODU','', (/ ('', i=1, 10) /))
     214  type(ctrl_out),save :: o_sconc05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     215  'sconc05','Surf. Conc. SCDU','', (/ ('', i=1, 10) /))
    223216
    224217!lessivage
    225218
    226   type(ctrl_out),save :: o_flux01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    227   'flux01','emission PREC','', (/ ('', i=1, 9) /))
    228   type(ctrl_out),save :: o_flux02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    229   'flux02','emission FINE','', (/ ('', i=1, 9) /))
    230   type(ctrl_out),save :: o_flux03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    231   'flux03','emission COSS','', (/ ('', i=1, 9) /))
    232   type(ctrl_out),save :: o_flux04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    233   'flux04','emission CODU','', (/ ('', i=1, 9) /))
    234   type(ctrl_out),save :: o_flux05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    235   'flux05','emission SCDU','', (/ ('', i=1, 9) /))
    236 
    237 
    238   type(ctrl_out),save :: o_ds01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    239   'ds01','Depot sec PREC','', (/ ('', i=1, 9) /))
    240   type(ctrl_out),save :: o_ds02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    241   'ds02','Depot sec FINE','', (/ ('', i=1, 9) /))
    242   type(ctrl_out),save :: o_ds03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    243   'ds03','Depot sec COSS','', (/ ('', i=1, 9) /))
    244   type(ctrl_out),save :: o_ds04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    245   'ds04','Depot sec CODU','', (/ ('', i=1, 9) /))
    246   type(ctrl_out),save :: o_ds05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    247   'ds05','Depot sec SCDU','', (/ ('', i=1, 9) /))
    248 
    249   type(ctrl_out),save :: o_dh01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    250   'dh01','Depot hum PREC','', (/ ('', i=1, 9) /))
    251   type(ctrl_out),save :: o_dh02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    252   'dh02','Depot hum FINE','', (/ ('', i=1, 9) /))
    253   type(ctrl_out),save :: o_dh03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    254   'dh03','Depot hum COSS','', (/ ('', i=1, 9) /))
    255   type(ctrl_out),save :: o_dh04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    256   'dh04','Depot hum CODU','', (/ ('', i=1, 9) /))
    257   type(ctrl_out),save :: o_dh05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    258   'dh05','Depot hum SCDU','', (/ ('', i=1, 9) /))
    259 
    260   type(ctrl_out),save :: o_dtrconv01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    261   'dtrconv01','Tiedke convective PREC','', (/ ('', i=1, 9) /))
    262   type(ctrl_out),save :: o_dtrconv02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    263   'dtrconv02','Tiedke convective FINE','', (/ ('', i=1, 9) /))
    264   type(ctrl_out),save :: o_dtrconv03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    265   'dtrconv03','Tiedke convective COSS','', (/ ('', i=1, 9) /))
    266   type(ctrl_out),save :: o_dtrconv04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    267   'dtrconv04','Tiedke convective CODU','', (/ ('', i=1, 9) /))
    268   type(ctrl_out),save :: o_dtrconv05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    269   'dtrconv05','Tiedke convective SCDU','', (/ ('', i=1, 9) /))
    270 
    271   type(ctrl_out),save :: o_dtherm01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    272   'dtherm01','Thermals dtracer PREC','', (/ ('', i=1, 9) /))
    273   type(ctrl_out),save :: o_dtherm02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    274   'dtherm02','Thermals dtracer FINE','', (/ ('', i=1, 9) /))
    275   type(ctrl_out),save :: o_dtherm03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    276   'dtherm03','Thermals dtracer COSS','', (/ ('', i=1, 9) /))
    277   type(ctrl_out),save :: o_dtherm04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    278   'dtherm04','Thermals dtracer CODU','', (/ ('', i=1, 9) /))
    279   type(ctrl_out),save :: o_dtherm05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    280   'dtherm05','Thermals dtracer SCDU','', (/ ('', i=1, 9) /))
    281 
    282   type(ctrl_out),save :: o_dhkecv01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    283   'dhkecv01','KE dep hum convective PREC','', (/ ('', i=1, 9) /))
    284   type(ctrl_out),save :: o_dhkecv02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    285   'dhkecv02','KE dep hum convective FINE','', (/ ('', i=1, 9) /))
    286   type(ctrl_out),save :: o_dhkecv03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    287   'dhkecv03','KE dep hum convective COSS','', (/ ('', i=1, 9) /))
    288   type(ctrl_out),save :: o_dhkecv04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    289   'dhkecv04','KE dep hum convective CODU','', (/ ('', i=1, 9) /))
    290   type(ctrl_out),save :: o_dhkecv05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    291   'dhkecv05','KE dep hum convective SCDU','', (/ ('', i=1, 9) /))
     219  type(ctrl_out),save :: o_flux01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     220  'flux01','emission PREC','', (/ ('', i=1, 10) /))
     221  type(ctrl_out),save :: o_flux02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     222  'flux02','emission FINE','', (/ ('', i=1, 10) /))
     223  type(ctrl_out),save :: o_flux03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     224  'flux03','emission COSS','', (/ ('', i=1, 10) /))
     225  type(ctrl_out),save :: o_flux04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     226  'flux04','emission CODU','', (/ ('', i=1, 10) /))
     227  type(ctrl_out),save :: o_flux05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     228  'flux05','emission SCDU','', (/ ('', i=1, 10) /))
     229
     230  type(ctrl_out),save :: o_ds01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     231  'ds01','Depot sec PREC','', (/ ('', i=1, 10) /))
     232  type(ctrl_out),save :: o_ds02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     233  'ds02','Depot sec FINE','', (/ ('', i=1, 10) /))
     234  type(ctrl_out),save :: o_ds03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     235  'ds03','Depot sec COSS','', (/ ('', i=1, 10) /))
     236  type(ctrl_out),save :: o_ds04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     237  'ds04','Depot sec CODU','', (/ ('', i=1, 10) /))
     238  type(ctrl_out),save :: o_ds05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     239  'ds05','Depot sec SCDU','', (/ ('', i=1, 10) /))
     240
     241  type(ctrl_out),save :: o_dh01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     242  'dh01','Depot hum PREC','', (/ ('', i=1, 10) /))
     243  type(ctrl_out),save :: o_dh02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     244  'dh02','Depot hum FINE','', (/ ('', i=1, 10) /))
     245  type(ctrl_out),save :: o_dh03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     246  'dh03','Depot hum COSS','', (/ ('', i=1, 10) /))
     247  type(ctrl_out),save :: o_dh04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     248  'dh04','Depot hum CODU','', (/ ('', i=1, 10) /))
     249  type(ctrl_out),save :: o_dh05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     250  'dh05','Depot hum SCDU','', (/ ('', i=1, 10) /))
     251
     252  type(ctrl_out),save :: o_dtrconv01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     253  'dtrconv01','Tiedke convective PREC','', (/ ('', i=1, 10) /))
     254  type(ctrl_out),save :: o_dtrconv02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     255  'dtrconv02','Tiedke convective FINE','', (/ ('', i=1, 10) /))
     256  type(ctrl_out),save :: o_dtrconv03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     257  'dtrconv03','Tiedke convective COSS','', (/ ('', i=1, 10) /))
     258  type(ctrl_out),save :: o_dtrconv04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     259  'dtrconv04','Tiedke convective CODU','', (/ ('', i=1, 10) /))
     260  type(ctrl_out),save :: o_dtrconv05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     261  'dtrconv05','Tiedke convective SCDU','', (/ ('', i=1, 10) /))
     262
     263  type(ctrl_out),save :: o_dtherm01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     264  'dtherm01','Thermals dtracer PREC','', (/ ('', i=1, 10) /))
     265  type(ctrl_out),save :: o_dtherm02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     266  'dtherm02','Thermals dtracer FINE','', (/ ('', i=1, 10) /))
     267  type(ctrl_out),save :: o_dtherm03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     268  'dtherm03','Thermals dtracer COSS','', (/ ('', i=1, 10) /))
     269  type(ctrl_out),save :: o_dtherm04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     270  'dtherm04','Thermals dtracer CODU','', (/ ('', i=1, 10) /))
     271  type(ctrl_out),save :: o_dtherm05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     272  'dtherm05','Thermals dtracer SCDU','', (/ ('', i=1, 10) /))
     273
     274  type(ctrl_out),save :: o_dhkecv01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     275  'dhkecv01','KE dep hum convective PREC','', (/ ('', i=1, 10) /))
     276  type(ctrl_out),save :: o_dhkecv02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     277  'dhkecv02','KE dep hum convective FINE','', (/ ('', i=1, 10) /))
     278  type(ctrl_out),save :: o_dhkecv03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     279  'dhkecv03','KE dep hum convective COSS','', (/ ('', i=1, 10) /))
     280  type(ctrl_out),save :: o_dhkecv04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     281  'dhkecv04','KE dep hum convective CODU','', (/ ('', i=1, 10) /))
     282  type(ctrl_out),save :: o_dhkecv05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     283  'dhkecv05','KE dep hum convective SCDU','', (/ ('', i=1, 10) /))
    292284
    293285!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    294   type(ctrl_out),save :: o_d_tr_ds01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    295   'd_tr_ds01','Tendance dep sec  PREC','', (/ ('', i=1, 9) /))
    296   type(ctrl_out),save :: o_d_tr_ds02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    297   'd_tr_ds02','Tendance dep sec FINE','', (/ ('', i=1, 9) /))
    298   type(ctrl_out),save :: o_d_tr_ds03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    299   'd_tr_ds03','Tendance dep sec COSS','', (/ ('', i=1, 9) /))
    300   type(ctrl_out),save :: o_d_tr_ds04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    301   'd_tr_ds04','Tendance depot sec CODU','', (/ ('', i=1, 9) /))
    302   type(ctrl_out),save :: o_d_tr_ds05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    303   'd_tr_ds05','Tendance dep sec SCDU','', (/ ('', i=1, 9) /))
     286  type(ctrl_out),save :: o_d_tr_ds01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     287  'd_tr_ds01','Tendance dep sec  PREC','', (/ ('', i=1, 10) /))
     288  type(ctrl_out),save :: o_d_tr_ds02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     289  'd_tr_ds02','Tendance dep sec FINE','', (/ ('', i=1, 10) /))
     290  type(ctrl_out),save :: o_d_tr_ds03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     291  'd_tr_ds03','Tendance dep sec COSS','', (/ ('', i=1, 10) /))
     292  type(ctrl_out),save :: o_d_tr_ds04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     293  'd_tr_ds04','Tendance depot sec CODU','', (/ ('', i=1, 10) /))
     294  type(ctrl_out),save :: o_d_tr_ds05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     295  'd_tr_ds05','Tendance dep sec SCDU','', (/ ('', i=1, 10) /))
    304296!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    305297
    306   type(ctrl_out),save :: o_dhkelsc01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    307   'dhkelsc01','KE dep hum large scale PREC','', (/ ('', i=1, 9) /))
    308   type(ctrl_out),save :: o_dhkelsc02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    309   'dhkelsc02','KE dep hum large scale FINE','', (/ ('', i=1, 9) /))
    310   type(ctrl_out),save :: o_dhkelsc03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    311   'dhkelsc03','KE dep hum large scale COSS','', (/ ('', i=1, 9) /))
    312   type(ctrl_out),save :: o_dhkelsc04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    313   'dhkelsc04','KE dep hum large scale CODU','', (/ ('', i=1, 9) /))
    314   type(ctrl_out),save :: o_dhkelsc05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    315   'dhkelsc05','KE dep hum large scale SCDU','', (/ ('', i=1, 9) /))
    316 
    317 
    318   type(ctrl_out),save :: o_d_tr_cv01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    319   'd_tr_cv01','cvltr d_tr_cv PREC','', (/ ('', i=1, 9) /))
    320   type(ctrl_out),save :: o_d_tr_cv02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    321   'd_tr_cv02','cvltr d_tr_cv FINE','', (/ ('', i=1, 9) /))
    322   type(ctrl_out),save :: o_d_tr_cv03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    323   'd_tr_cv03','cvltr d_tr_cv COSS','', (/ ('', i=1, 9) /))
    324   type(ctrl_out),save :: o_d_tr_cv04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    325   'd_tr_cv04','cvltr d_tr_cv CODU','', (/ ('', i=1, 9) /))
    326   type(ctrl_out),save :: o_d_tr_cv05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    327   'd_tr_cv05','cvltr d_tr_cv SCDU','', (/ ('', i=1, 9) /))
    328 
    329 
    330   type(ctrl_out),save :: o_d_tr_trsp01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    331   'd_tr_trsp01','cvltr d_tr_trsp PREC','', (/ ('', i=1, 9) /))
    332   type(ctrl_out),save :: o_d_tr_trsp02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    333   'd_tr_trsp02','cvltr d_tr_trsp FINE','', (/ ('', i=1, 9) /))
    334   type(ctrl_out),save :: o_d_tr_trsp03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    335   'd_tr_trsp03','cvltr d_tr_trsp COSS','', (/ ('', i=1, 9) /))
    336   type(ctrl_out),save :: o_d_tr_trsp04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    337   'd_tr_trsp04','cvltr d_tr_trsp CODU','', (/ ('', i=1, 9) /))
    338   type(ctrl_out),save :: o_d_tr_trsp05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    339   'd_tr_trsp05','cvltr d_tr_trsp SCDU','', (/ ('', i=1, 9) /))
    340 
    341 
    342   type(ctrl_out),save :: o_d_tr_sscav01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    343   'd_tr_sscav01','cvltr d_tr_sscav PREC','', (/ ('', i=1, 9) /))
    344   type(ctrl_out),save :: o_d_tr_sscav02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    345   'd_tr_sscav02','cvltr d_tr_sscav FINE','', (/ ('', i=1, 9) /))
    346   type(ctrl_out),save :: o_d_tr_sscav03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    347   'd_tr_sscav03','cvltr d_tr_sscav COSS','', (/ ('', i=1, 9) /))
    348   type(ctrl_out),save :: o_d_tr_sscav04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    349   'd_tr_sscav04','cvltr d_tr_sscav CODU','', (/ ('', i=1, 9) /))
    350   type(ctrl_out),save :: o_d_tr_sscav05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    351   'd_tr_sscav05','cvltr d_tr_sscav SCDU','', (/ ('', i=1, 9) /))
    352 
    353   type(ctrl_out),save :: o_d_tr_sat01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    354   'd_tr_sat01','cvltr d_tr_sat PREC','', (/ ('', i=1, 9) /))
    355   type(ctrl_out),save :: o_d_tr_sat02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    356   'd_tr_sat02','cvltr d_tr_sat FINE','', (/ ('', i=1, 9) /))
    357   type(ctrl_out),save :: o_d_tr_sat03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    358   'd_tr_sat03','cvltr d_tr_sat COSS','', (/ ('', i=1, 9) /))
    359   type(ctrl_out),save :: o_d_tr_sat04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    360   'd_tr_sat04','cvltr d_tr_sat CODU','', (/ ('', i=1, 9) /))
    361   type(ctrl_out),save :: o_d_tr_sat05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    362   'd_tr_sat05','cvltr d_tr_sat SCDU','', (/ ('', i=1, 9) /))
    363 
    364   type(ctrl_out),save :: o_d_tr_uscav01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    365   'd_tr_uscav01','cvltr d_tr_uscav PREC','', (/ ('', i=1, 9) /))
    366   type(ctrl_out),save :: o_d_tr_uscav02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    367   'd_tr_uscav02','cvltr d_tr_uscav FINE','', (/ ('', i=1, 9) /))
    368   type(ctrl_out),save :: o_d_tr_uscav03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    369   'd_tr_uscav03','cvltr d_tr_uscav COSS','', (/ ('', i=1, 9) /))
    370   type(ctrl_out),save :: o_d_tr_uscav04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    371   'd_tr_uscav04','cvltr d_tr_uscav CODU','', (/ ('', i=1, 9) /))
    372   type(ctrl_out),save :: o_d_tr_uscav05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    373   'd_tr_uscav05','cvltr d_tr_uscav SCDU','', (/ ('', i=1, 9) /))
     298  type(ctrl_out),save :: o_dhkelsc01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     299  'dhkelsc01','KE dep hum large scale PREC','', (/ ('', i=1, 10) /))
     300  type(ctrl_out),save :: o_dhkelsc02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     301  'dhkelsc02','KE dep hum large scale FINE','', (/ ('', i=1, 10) /))
     302  type(ctrl_out),save :: o_dhkelsc03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     303  'dhkelsc03','KE dep hum large scale COSS','', (/ ('', i=1, 10) /))
     304  type(ctrl_out),save :: o_dhkelsc04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     305  'dhkelsc04','KE dep hum large scale CODU','', (/ ('', i=1, 10) /))
     306  type(ctrl_out),save :: o_dhkelsc05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     307  'dhkelsc05','KE dep hum large scale SCDU','', (/ ('', i=1, 10) /))
     308
     309  type(ctrl_out),save :: o_d_tr_cv01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     310  'd_tr_cv01','cvltr d_tr_cv PREC','', (/ ('', i=1, 10) /))
     311  type(ctrl_out),save :: o_d_tr_cv02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     312  'd_tr_cv02','cvltr d_tr_cv FINE','', (/ ('', i=1, 10) /))
     313  type(ctrl_out),save :: o_d_tr_cv03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     314  'd_tr_cv03','cvltr d_tr_cv COSS','', (/ ('', i=1, 10) /))
     315  type(ctrl_out),save :: o_d_tr_cv04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     316  'd_tr_cv04','cvltr d_tr_cv CODU','', (/ ('', i=1, 10) /))
     317  type(ctrl_out),save :: o_d_tr_cv05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     318  'd_tr_cv05','cvltr d_tr_cv SCDU','', (/ ('', i=1, 10) /))
     319
     320  type(ctrl_out),save :: o_d_tr_trsp01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     321  'd_tr_trsp01','cvltr d_tr_trsp PREC','', (/ ('', i=1, 10) /))
     322  type(ctrl_out),save :: o_d_tr_trsp02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     323  'd_tr_trsp02','cvltr d_tr_trsp FINE','', (/ ('', i=1, 10) /))
     324  type(ctrl_out),save :: o_d_tr_trsp03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     325  'd_tr_trsp03','cvltr d_tr_trsp COSS','', (/ ('', i=1, 10) /))
     326  type(ctrl_out),save :: o_d_tr_trsp04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     327  'd_tr_trsp04','cvltr d_tr_trsp CODU','', (/ ('', i=1, 10) /))
     328  type(ctrl_out),save :: o_d_tr_trsp05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     329  'd_tr_trsp05','cvltr d_tr_trsp SCDU','', (/ ('', i=1, 10) /))
     330
     331  type(ctrl_out),save :: o_d_tr_sscav01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     332  'd_tr_sscav01','cvltr d_tr_sscav PREC','', (/ ('', i=1, 10) /))
     333  type(ctrl_out),save :: o_d_tr_sscav02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     334  'd_tr_sscav02','cvltr d_tr_sscav FINE','', (/ ('', i=1, 10) /))
     335  type(ctrl_out),save :: o_d_tr_sscav03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     336  'd_tr_sscav03','cvltr d_tr_sscav COSS','', (/ ('', i=1, 10) /))
     337  type(ctrl_out),save :: o_d_tr_sscav04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     338  'd_tr_sscav04','cvltr d_tr_sscav CODU','', (/ ('', i=1, 10) /))
     339  type(ctrl_out),save :: o_d_tr_sscav05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     340  'd_tr_sscav05','cvltr d_tr_sscav SCDU','', (/ ('', i=1, 10) /))
     341
     342  type(ctrl_out),save :: o_d_tr_sat01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     343  'd_tr_sat01','cvltr d_tr_sat PREC','', (/ ('', i=1, 10) /))
     344  type(ctrl_out),save :: o_d_tr_sat02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     345  'd_tr_sat02','cvltr d_tr_sat FINE','', (/ ('', i=1, 10) /))
     346  type(ctrl_out),save :: o_d_tr_sat03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     347  'd_tr_sat03','cvltr d_tr_sat COSS','', (/ ('', i=1, 10) /))
     348  type(ctrl_out),save :: o_d_tr_sat04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     349  'd_tr_sat04','cvltr d_tr_sat CODU','', (/ ('', i=1, 10) /))
     350  type(ctrl_out),save :: o_d_tr_sat05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     351  'd_tr_sat05','cvltr d_tr_sat SCDU','', (/ ('', i=1, 10) /))
     352
     353  type(ctrl_out),save :: o_d_tr_uscav01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     354  'd_tr_uscav01','cvltr d_tr_uscav PREC','', (/ ('', i=1, 10) /))
     355  type(ctrl_out),save :: o_d_tr_uscav02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     356  'd_tr_uscav02','cvltr d_tr_uscav FINE','', (/ ('', i=1, 10) /))
     357  type(ctrl_out),save :: o_d_tr_uscav03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     358  'd_tr_uscav03','cvltr d_tr_uscav COSS','', (/ ('', i=1, 10) /))
     359  type(ctrl_out),save :: o_d_tr_uscav04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     360  'd_tr_uscav04','cvltr d_tr_uscav CODU','', (/ ('', i=1, 10) /))
     361  type(ctrl_out),save :: o_d_tr_uscav05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     362  'd_tr_uscav05','cvltr d_tr_uscav SCDU','', (/ ('', i=1, 10) /))
    374363
    375364!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     365  type(ctrl_out),save :: o_d_tr_insc01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     366  'd_tr_insc01','large-scale d_tr_insc PREC','', (/ ('', i=1, 10) /))
     367  type(ctrl_out),save :: o_d_tr_insc02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     368  'd_tr_insc02','large-scale d_tr_insc FINE','', (/ ('', i=1, 10) /))
     369  type(ctrl_out),save :: o_d_tr_insc03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     370  'd_tr_insc03','large-scale d_tr_insc COSS','', (/ ('', i=1, 10) /))
     371  type(ctrl_out),save :: o_d_tr_insc04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     372  'd_tr_insc04','large-scale d_tr_insc CODU','', (/ ('', i=1, 10) /))
     373  type(ctrl_out),save :: o_d_tr_insc05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     374  'd_tr_insc05','large-scale d_tr_insc SCDU','', (/ ('', i=1, 10) /))
    376375                !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    377   type(ctrl_out),save :: o_d_tr_insc01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    378   'd_tr_insc01','large-scale d_tr_insc PREC','', (/ ('', i=1, 9) /))
    379   type(ctrl_out),save :: o_d_tr_insc02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    380   'd_tr_insc02','large-scale d_tr_insc FINE','', (/ ('', i=1, 9) /))
    381   type(ctrl_out),save :: o_d_tr_insc03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    382   'd_tr_insc03','large-scale d_tr_insc COSS','', (/ ('', i=1, 9) /))
    383   type(ctrl_out),save :: o_d_tr_insc04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    384   'd_tr_insc04','large-scale d_tr_insc CODU','', (/ ('', i=1, 9) /))
    385   type(ctrl_out),save :: o_d_tr_insc05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    386   'd_tr_insc05','large-scale d_tr_insc SCDU','', (/ ('', i=1, 9) /))
     376  type(ctrl_out),save :: o_d_tr_bcscav01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     377  'd_tr_bcscav01','large-scale d_tr_bcscav PREC','', (/ ('', i=1, 10) /))
     378  type(ctrl_out),save :: o_d_tr_bcscav02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     379  'd_tr_bcscav02','large-scale d_tr_bcscav FINE','', (/ ('', i=1, 10) /))
     380  type(ctrl_out),save :: o_d_tr_bcscav03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     381  'd_tr_bcscav03','large-scale d_tr_bcscav COSS','', (/ ('', i=1, 10) /))
     382  type(ctrl_out),save :: o_d_tr_bcscav04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     383  'd_tr_bcscav04','large-scale d_tr_bcscav CODU','', (/ ('', i=1, 10) /))
     384  type(ctrl_out),save :: o_d_tr_bcscav05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     385  'd_tr_bcscav05','large-scale d_tr_bcscav SCDU','', (/ ('', i=1, 10) /))
    387386                !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    388   type(ctrl_out),save :: o_d_tr_bcscav01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    389   'd_tr_bcscav01','large-scale d_tr_bcscav PREC','', (/ ('', i=1, 9) /))
    390   type(ctrl_out),save :: o_d_tr_bcscav02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    391   'd_tr_bcscav02','large-scale d_tr_bcscav FINE','', (/ ('', i=1, 9) /))
    392   type(ctrl_out),save :: o_d_tr_bcscav03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    393   'd_tr_bcscav03','large-scale d_tr_bcscav COSS','', (/ ('', i=1, 9) /))
    394   type(ctrl_out),save :: o_d_tr_bcscav04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    395   'd_tr_bcscav04','large-scale d_tr_bcscav CODU','', (/ ('', i=1, 9) /))
    396   type(ctrl_out),save :: o_d_tr_bcscav05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    397   'd_tr_bcscav05','large-scale d_tr_bcscav SCDU','', (/ ('', i=1, 9) /))
     387  type(ctrl_out),save :: o_d_tr_evapls01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     388  'd_tr_evapls01','large-scale d_tr_evapls PREC','', (/ ('', i=1, 10) /))
     389  type(ctrl_out),save :: o_d_tr_evapls02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     390  'd_tr_evapls02','large-scale d_tr_evapls FINE','', (/ ('', i=1, 10) /))
     391  type(ctrl_out),save :: o_d_tr_evapls03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     392  'd_tr_evapls03','large-scale d_tr_evapls COSS','', (/ ('', i=1, 10) /))
     393  type(ctrl_out),save :: o_d_tr_evapls04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     394  'd_tr_evapls04','large-scale d_tr_evapls CODU','', (/ ('', i=1, 10) /))
     395  type(ctrl_out),save :: o_d_tr_evapls05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     396  'd_tr_evapls05','large-scale d_tr_evapls SCDU','', (/ ('', i=1, 10) /))
    398397                !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    399   type(ctrl_out),save :: o_d_tr_evapls01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    400   'd_tr_evapls01','large-scale d_tr_evapls PREC','', (/ ('', i=1, 9) /))
    401   type(ctrl_out),save :: o_d_tr_evapls02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    402   'd_tr_evapls02','large-scale d_tr_evapls FINE','', (/ ('', i=1, 9) /))
    403   type(ctrl_out),save :: o_d_tr_evapls03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    404   'd_tr_evapls03','large-scale d_tr_evapls COSS','', (/ ('', i=1, 9) /))
    405   type(ctrl_out),save :: o_d_tr_evapls04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    406   'd_tr_evapls04','large-scale d_tr_evapls CODU','', (/ ('', i=1, 9) /))
    407   type(ctrl_out),save :: o_d_tr_evapls05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    408   'd_tr_evapls05','large-scale d_tr_evapls SCDU','', (/ ('', i=1, 9) /))
     398  type(ctrl_out),save :: o_d_tr_ls01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     399  'd_tr_ls01','large-scale d_tr_ls PREC','', (/ ('', i=1, 10) /))
     400  type(ctrl_out),save :: o_d_tr_ls02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     401  'd_tr_ls02','large-scale d_tr_ls FINE','', (/ ('', i=1, 10) /))
     402  type(ctrl_out),save :: o_d_tr_ls03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     403  'd_tr_ls03','large-scale d_tr_ls COSS','', (/ ('', i=1, 10) /))
     404  type(ctrl_out),save :: o_d_tr_ls04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     405  'd_tr_ls04','large-scale d_tr_ls CODU','', (/ ('', i=1, 10) /))
     406  type(ctrl_out),save :: o_d_tr_ls05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     407  'd_tr_ls05','large-scale d_tr_ls SCDU','', (/ ('', i=1, 10) /))
    409408                !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    410   type(ctrl_out),save :: o_d_tr_ls01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    411   'd_tr_ls01','large-scale d_tr_ls PREC','', (/ ('', i=1, 9) /))
    412   type(ctrl_out),save :: o_d_tr_ls02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    413   'd_tr_ls02','large-scale d_tr_ls FINE','', (/ ('', i=1, 9) /))
    414   type(ctrl_out),save :: o_d_tr_ls03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    415   'd_tr_ls03','large-scale d_tr_ls COSS','', (/ ('', i=1, 9) /))
    416   type(ctrl_out),save :: o_d_tr_ls04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    417   'd_tr_ls04','large-scale d_tr_ls CODU','', (/ ('', i=1, 9) /))
    418   type(ctrl_out),save :: o_d_tr_ls05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    419   'd_tr_ls05','large-scale d_tr_ls SCDU','', (/ ('', i=1, 9) /))
     409  type(ctrl_out),save :: o_d_tr_dyn01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     410  'd_tr_dyn01','cvltr d_tr_dyn PREC','', (/ ('', i=1, 10) /))
     411  type(ctrl_out),save :: o_d_tr_dyn02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     412  'd_tr_dyn02','cvltr d_tr_dyn FINE','', (/ ('', i=1, 10) /))
     413  type(ctrl_out),save :: o_d_tr_dyn03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     414  'd_tr_dyn03','cvltr d_tr_dyn COSS','', (/ ('', i=1, 10) /))
     415  type(ctrl_out),save :: o_d_tr_dyn04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     416  'd_tr_dyn04','cvltr d_tr_dyn CODU','', (/ ('', i=1, 10) /))
     417  type(ctrl_out),save :: o_d_tr_dyn05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     418  'd_tr_dyn05','cvltr d_tr_dyn SCDU','', (/ ('', i=1, 10) /))
    420419                !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    421 
    422 
     420  type(ctrl_out),save :: o_d_tr_cl01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     421  'd_tr_cl01','cvltr d_tr_cl PREC','', (/ ('', i=1, 10) /))
     422  type(ctrl_out),save :: o_d_tr_cl02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     423  'd_tr_cl02','cvltr d_tr_cl FINE','', (/ ('', i=1, 10) /))
     424  type(ctrl_out),save :: o_d_tr_cl03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     425  'd_tr_cl03','cvltr d_tr_cl COSS','', (/ ('', i=1, 10) /))
     426  type(ctrl_out),save :: o_d_tr_cl04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     427  'd_tr_cl04','cvltr d_tr_cl CODU','', (/ ('', i=1, 10) /))
     428  type(ctrl_out),save :: o_d_tr_cl05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     429  'd_tr_cl05','cvltr d_tr_cl SCDU','', (/ ('', i=1, 10) /))
    423430                !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    424   type(ctrl_out),save :: o_d_tr_dyn01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    425   'd_tr_dyn01','cvltr d_tr_dyn PREC','', (/ ('', i=1, 9) /))
    426   type(ctrl_out),save :: o_d_tr_dyn02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    427   'd_tr_dyn02','cvltr d_tr_dyn FINE','', (/ ('', i=1, 9) /))
    428   type(ctrl_out),save :: o_d_tr_dyn03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    429   'd_tr_dyn03','cvltr d_tr_dyn COSS','', (/ ('', i=1, 9) /))
    430   type(ctrl_out),save :: o_d_tr_dyn04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    431   'd_tr_dyn04','cvltr d_tr_dyn CODU','', (/ ('', i=1, 9) /))
    432   type(ctrl_out),save :: o_d_tr_dyn05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    433   'd_tr_dyn05','cvltr d_tr_dyn SCDU','', (/ ('', i=1, 9) /))
    434                 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    435                 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    436   type(ctrl_out),save :: o_d_tr_cl01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    437   'd_tr_cl01','cvltr d_tr_cl PREC','', (/ ('', i=1, 9) /))
    438   type(ctrl_out),save :: o_d_tr_cl02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    439   'd_tr_cl02','cvltr d_tr_cl FINE','', (/ ('', i=1, 9) /))
    440   type(ctrl_out),save :: o_d_tr_cl03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    441   'd_tr_cl03','cvltr d_tr_cl COSS','', (/ ('', i=1, 9) /))
    442   type(ctrl_out),save :: o_d_tr_cl04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    443   'd_tr_cl04','cvltr d_tr_cl CODU','', (/ ('', i=1, 9) /))
    444   type(ctrl_out),save :: o_d_tr_cl05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    445   'd_tr_cl05','cvltr d_tr_cl SCDU','', (/ ('', i=1, 9) /))
    446                 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    447   type(ctrl_out),save :: o_d_tr_th01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    448   'd_tr_th01','cvltr d_tr_th PREC','', (/ ('', i=1, 9) /))
    449   type(ctrl_out),save :: o_d_tr_th02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    450   'd_tr_th02','cvltr d_tr_th FINE','', (/ ('', i=1, 9) /))
    451   type(ctrl_out),save :: o_d_tr_th03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    452   'd_tr_th03','cvltr d_tr_th COSS','', (/ ('', i=1, 9) /))
    453   type(ctrl_out),save :: o_d_tr_th04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    454   'd_tr_th04','cvltr d_tr_th CODU','', (/ ('', i=1, 9) /))
    455   type(ctrl_out),save :: o_d_tr_th05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    456   'd_tr_th05','cvltr d_tr_th SCDU','', (/ ('', i=1, 9) /))
    457                 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    458 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     431  type(ctrl_out),save :: o_d_tr_th01     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     432  'd_tr_th01','cvltr d_tr_th PREC','', (/ ('', i=1, 10) /))
     433  type(ctrl_out),save :: o_d_tr_th02     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     434  'd_tr_th02','cvltr d_tr_th FINE','', (/ ('', i=1, 10) /))
     435  type(ctrl_out),save :: o_d_tr_th03     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     436  'd_tr_th03','cvltr d_tr_th COSS','', (/ ('', i=1, 10) /))
     437  type(ctrl_out),save :: o_d_tr_th04     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     438  'd_tr_th04','cvltr d_tr_th CODU','', (/ ('', i=1, 10) /))
     439  type(ctrl_out),save :: o_d_tr_th05     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     440  'd_tr_th05','cvltr d_tr_th SCDU','', (/ ('', i=1, 10) /))
    459441!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    460   type(ctrl_out),save :: o_sed_ss3D     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    461   'sed_ss3D','Tendance Sedmet. COSS','', (/ ('', i=1, 9) /))
    462   type(ctrl_out),save :: o_sed_dust3D     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    463   'sed_dust3D','Tendance Sedmet. CODU','', (/ ('', i=1, 9) /))
    464   type(ctrl_out),save :: o_sed_dustsco3D     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    465   'sed_dustsco3D','Tendance Sedmet. SCDU','', (/ ('', i=1, 9) /))
     442  type(ctrl_out),save :: o_sed_ss3D     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     443  'sed_ss3D','Tendance Sedmet. COSS','', (/ ('', i=1, 10) /))
     444  type(ctrl_out),save :: o_sed_dust3D     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     445  'sed_dust3D','Tendance Sedmet. CODU','', (/ ('', i=1, 10) /))
     446  type(ctrl_out),save :: o_sed_dustsco3D     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     447  'sed_dustsco3D','Tendance Sedmet. SCDU','', (/ ('', i=1, 10) /))
    466448!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    467 
    468   type(ctrl_out),save :: o_sed_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    469   'sed_ss','Sedmet. COSS','', (/ ('', i=1, 9) /))
    470   type(ctrl_out),save :: o_sed_dust     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    471   'sed_dust','Sedmet. CODU','', (/ ('', i=1, 9) /))
    472   type(ctrl_out),save :: o_sed_dustsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    473   'sed_dustsco','Sedmet. SCDU','', (/ ('', i=1, 9) /))
    474   type(ctrl_out),save :: o_g2p_gas     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    475   'g2p_gas','Gas2particle gas sink','', (/ ('', i=1, 9) /))
    476   type(ctrl_out),save :: o_g2p_aer     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    477   'g2p_aer','Gas2particle tr2 src','', (/ ('', i=1, 9) /))
    478 
     449  type(ctrl_out),save :: o_sed_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     450  'sed_ss','Sedmet. COSS','', (/ ('', i=1, 10) /))
     451  type(ctrl_out),save :: o_sed_dust     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     452  'sed_dust','Sedmet. CODU','', (/ ('', i=1, 10) /))
     453  type(ctrl_out),save :: o_sed_dustsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     454  'sed_dustsco','Sedmet. SCDU','', (/ ('', i=1, 10) /))
     455  type(ctrl_out),save :: o_g2p_gas     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     456  'g2p_gas','Gas2particle gas sink','', (/ ('', i=1, 10) /))
     457  type(ctrl_out),save :: o_g2p_aer     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     458  'g2p_aer','Gas2particle tr2 src','', (/ ('', i=1, 10) /))
    479459!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    480460! histrac
    481461
    482   type(ctrl_out),save :: o_fluxbb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    483   'fluxbb','Flux BB','mg/m2/s', (/ ('', i=1, 9) /))
    484 
    485   type(ctrl_out),save :: o_fluxff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    486   'fluxff','Flux FF','mg/m2/s', (/ ('', i=1, 9) /))
    487 
    488   type(ctrl_out),save :: o_fluxbcbb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    489   'fluxbcbb','Flux BC-BB','mg/m2/s', (/ ('', i=1, 9) /))
    490 
    491   type(ctrl_out),save :: o_fluxbcff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    492   'fluxbcff','Flux BC-FF','mg/m2/s', (/ ('', i=1, 9) /))
    493 
    494   type(ctrl_out),save :: o_fluxbcnff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    495   'fluxbcnff','Flux BC-NFF','mg/m2/s', (/ ('', i=1, 9) /))
    496 
    497   type(ctrl_out),save :: o_fluxbcba     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    498   'fluxbcba','Flux BC-BA','mg/m2/s', (/ ('', i=1, 9) /))
    499 
    500   type(ctrl_out),save :: o_fluxbc     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    501   'fluxbc','Flux BC','mg/m2/s', (/ ('', i=1, 9) /))
    502 
    503   type(ctrl_out),save :: o_fluxombb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    504   'fluxombb','Flux OM-BB','mg/m2/s', (/ ('', i=1, 9) /))
    505 
    506   type(ctrl_out),save :: o_fluxomff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    507   'fluxomff','Flux OM-FF','mg/m2/s', (/ ('', i=1, 9) /))
    508 
    509   type(ctrl_out),save :: o_fluxomnff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    510   'fluxomnff','Flux OM-NFF','mg/m2/s', (/ ('', i=1, 9) /))
    511 
    512   type(ctrl_out),save :: o_fluxomba     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    513   'fluxomba','Flux OM-BA','mg/m2/s', (/ ('', i=1, 9) /))
    514 
    515   type(ctrl_out),save :: o_fluxomnat     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    516   'fluxomnat','Flux OM-NT','mg/m2/s', (/ ('', i=1, 9) /))
    517 
    518   type(ctrl_out),save :: o_fluxom     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    519   'fluxom','Flux OM','mg/m2/s', (/ ('', i=1, 9) /))
    520 
    521   type(ctrl_out),save :: o_fluxh2sff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    522   'fluxh2sff','Flux H2S FF','mgS/m2/s', (/ ('', i=1, 9) /))
    523 
    524   type(ctrl_out),save :: o_fluxh2snff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    525   'fluxh2snff','Flux H2S non-FF','mgS/m2/s', (/ ('', i=1, 9) /))
    526 
    527   type(ctrl_out),save :: o_fluxso2ff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    528   'fluxso2ff','Flux SO2 FF','mgS/m2/s', (/ ('', i=1, 9) /))
    529 
    530   type(ctrl_out),save :: o_fluxso2nff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    531   'fluxso2nff','Flux SO2 non-FF','mgS/m2/s', (/ ('', i=1, 9) /))
    532 
    533   type(ctrl_out),save :: o_fluxso2bb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    534   'fluxso2bb','Flux SO2 BB','mgS/m2/s', (/ ('', i=1, 9) /))
    535 
    536   type(ctrl_out),save :: o_fluxso2vol     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    537   'fluxso2vol','Flux SO2 Vol','mgS/m2/s', (/ ('', i=1, 9) /))
    538 
    539   type(ctrl_out),save :: o_fluxso2ba     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    540   'fluxso2ba','Flux SO2 Ba','mgS/m2/s', (/ ('', i=1, 9) /))
    541 
    542   type(ctrl_out),save :: o_fluxso2     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    543   'fluxso2','Flux SO2','mgS/m2/s', (/ ('', i=1, 9) /))
    544 
    545   type(ctrl_out),save :: o_fluxso4ff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    546   'fluxso4ff','Flux SO4 FF','mgS/m2/s', (/ ('', i=1, 9) /))
    547 
    548   type(ctrl_out),save :: o_fluxso4nff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    549   'fluxso4nff','Flux SO4 non-FF','mgS/m2/s', (/ ('', i=1, 9) /))
    550 
    551   type(ctrl_out),save :: o_fluxso4bb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    552   'fluxso4bb','Flux SO4 BB','mgS/m2/s', (/ ('', i=1, 9) /))
    553 
    554   type(ctrl_out),save :: o_fluxso4ba     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    555   'fluxso4ba','Flux SO4 Ba','mgS/m2/s', (/ ('', i=1, 9) /))
    556 
    557   type(ctrl_out),save :: o_fluxso4     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    558   'fluxso4','Flux SO4','mgS/m2/s', (/ ('', i=1, 9) /))
    559 
    560   type(ctrl_out),save :: o_fluxdms     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    561   'fluxdms','Flux DMS','mgS/m2/s', (/ ('', i=1, 9) /))
    562 
    563   type(ctrl_out),save :: o_fluxh2sbio     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    564   'fluxh2sbio','Flux H2S Bio','mgS/m2/s', (/ ('', i=1, 9) /))
    565 
    566   type(ctrl_out),save :: o_fluxdustec     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    567   'fluxdustec','Flux Dust EC','mg/m2/s', (/ ('', i=1, 9) /))
    568 
    569   type(ctrl_out),save :: o_fluxddfine     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    570   'fluxddfine','DD Fine Mode','mg/m2/s', (/ ('', i=1, 9) /))
    571 
    572   type(ctrl_out),save :: o_fluxddcoa     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    573   'fluxddcoa','DD Coarse Mode','mg/m2/s', (/ ('', i=1, 9) /))
    574 
    575   type(ctrl_out),save :: o_fluxddsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    576   'fluxddsco','DD SCoarse Mode','mg/m2/s', (/ ('', i=1, 9) /))
    577 
    578   type(ctrl_out),save :: o_fluxdd     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    579   'fluxdd','Flux DD','mg/m2/s', (/ ('', i=1, 9) /))
    580 
    581   type(ctrl_out),save :: o_fluxssfine     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    582   'fluxssfine','SS Fine Mode','mg/m2/s', (/ ('', i=1, 9) /))
    583 
    584   type(ctrl_out),save :: o_fluxsscoa     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    585   'fluxsscoa','SS Coarse Mode','mg/m2/s', (/ ('', i=1, 9) /))
    586 
    587   type(ctrl_out),save :: o_fluxss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    588   'fluxss','Flux SS','mg/m2/s', (/ ('', i=1, 9) /))
    589 
    590   type(ctrl_out),save :: o_flux_sparam_ind     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    591   'flux_sparam_ind','Ind emiss','mg/m2/s', (/ ('', i=1, 9) /))
    592 
    593   type(ctrl_out),save :: o_flux_sparam_bb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    594   'flux_sparam_bb','BB emiss','mg/m2/s', (/ ('', i=1, 9) /))
    595 
    596   type(ctrl_out),save :: o_flux_sparam_ff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    597   'flux_sparam_ff','FF emiss','mg/m2/s', (/ ('', i=1, 9) /))
    598 
    599   type(ctrl_out),save :: o_flux_sparam_ddfine     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    600   'flux_sparam_ddfine','DD fine emiss','mg/m2/s', (/ ('', i=1, 9) /))
    601 
    602   type(ctrl_out),save :: o_flux_sparam_ddcoa     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    603   'flux_sparam_ddcoa','DD coarse emiss','mg/m2/s', (/ ('', i=1, 9) /))
    604 
    605   type(ctrl_out),save :: o_flux_sparam_ddsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    606   'flux_sparam_ddsco','DD Scoarse emiss','mg/m2/s', (/ ('', i=1, 9) /))
    607 
    608   type(ctrl_out),save :: o_flux_sparam_ssfine     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    609   'flux_sparam_ssfine','SS fine emiss','mg/m2/s', (/ ('', i=1, 9) /))
    610 
    611 
    612   type(ctrl_out),save :: o_flux_sparam_sscoa     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    613   'flux_sparam_sscoa','SS coarse emiss','mg/m2/s', (/ ('', i=1, 9) /))
    614 
    615   type(ctrl_out),save :: o_u10m_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    616   'u10m_ss','Zonal wind at 10 m SS','m/s', (/ ('', i=1, 9) /))
    617 
    618   type(ctrl_out),save :: o_v10m_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    619   'v10m_ss','Meridional wind at 10 m SS','m/s', (/ ('', i=1, 9) /))
    620 
    621 
    622 !  type(ctrl_out),save :: o_     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    623 !  '','','', (/ ('', i=1, 9) /))
    624 
    625 
    626 
    627 !example  TYPE(ctrl_out), SAVE :: o_psbg = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5 /), &
     462  type(ctrl_out),save :: o_fluxbb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     463  'fluxbb','Flux BB','mg/m2/s', (/ ('', i=1, 10) /))
     464
     465  type(ctrl_out),save :: o_fluxff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     466  'fluxff','Flux FF','mg/m2/s', (/ ('', i=1, 10) /))
     467
     468  type(ctrl_out),save :: o_fluxbcbb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     469  'fluxbcbb','Flux BC-BB','mg/m2/s', (/ ('', i=1, 10) /))
     470
     471  type(ctrl_out),save :: o_fluxbcff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     472  'fluxbcff','Flux BC-FF','mg/m2/s', (/ ('', i=1, 10) /))
     473
     474  type(ctrl_out),save :: o_fluxbcnff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     475  'fluxbcnff','Flux BC-NFF','mg/m2/s', (/ ('', i=1, 10) /))
     476
     477  type(ctrl_out),save :: o_fluxbcba     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     478  'fluxbcba','Flux BC-BA','mg/m2/s', (/ ('', i=1, 10) /))
     479
     480  type(ctrl_out),save :: o_fluxbc     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     481  'fluxbc','Flux BC','mg/m2/s', (/ ('', i=1, 10) /))
     482
     483  type(ctrl_out),save :: o_fluxombb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     484  'fluxombb','Flux OM-BB','mg/m2/s', (/ ('', i=1, 10) /))
     485
     486  type(ctrl_out),save :: o_fluxomff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     487  'fluxomff','Flux OM-FF','mg/m2/s', (/ ('', i=1, 10) /))
     488
     489  type(ctrl_out),save :: o_fluxomnff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     490  'fluxomnff','Flux OM-NFF','mg/m2/s', (/ ('', i=1, 10) /))
     491
     492  type(ctrl_out),save :: o_fluxomba     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     493  'fluxomba','Flux OM-BA','mg/m2/s', (/ ('', i=1, 10) /))
     494
     495  type(ctrl_out),save :: o_fluxomnat     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     496  'fluxomnat','Flux OM-NT','mg/m2/s', (/ ('', i=1, 10) /))
     497
     498  type(ctrl_out),save :: o_fluxom     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     499  'fluxom','Flux OM','mg/m2/s', (/ ('', i=1, 10) /))
     500
     501  type(ctrl_out),save :: o_fluxh2sff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     502  'fluxh2sff','Flux H2S FF','mgS/m2/s', (/ ('', i=1, 10) /))
     503
     504  type(ctrl_out),save :: o_fluxh2snff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     505  'fluxh2snff','Flux H2S non-FF','mgS/m2/s', (/ ('', i=1, 10) /))
     506
     507  type(ctrl_out),save :: o_fluxso2ff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     508  'fluxso2ff','Flux SO2 FF','mgS/m2/s', (/ ('', i=1, 10) /))
     509
     510  type(ctrl_out),save :: o_fluxso2nff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     511  'fluxso2nff','Flux SO2 non-FF','mgS/m2/s', (/ ('', i=1, 10) /))
     512
     513  type(ctrl_out),save :: o_fluxso2bb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     514  'fluxso2bb','Flux SO2 BB','mgS/m2/s', (/ ('', i=1, 10) /))
     515
     516  type(ctrl_out),save :: o_fluxso2vol     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     517  'fluxso2vol','Flux SO2 Vol','mgS/m2/s', (/ ('', i=1, 10) /))
     518
     519  type(ctrl_out),save :: o_fluxso2ba     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     520  'fluxso2ba','Flux SO2 Ba','mgS/m2/s', (/ ('', i=1, 10) /))
     521
     522  type(ctrl_out),save :: o_fluxso2     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     523  'fluxso2','Flux SO2','mgS/m2/s', (/ ('', i=1, 10) /))
     524
     525  type(ctrl_out),save :: o_fluxso4ff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     526  'fluxso4ff','Flux SO4 FF','mgS/m2/s', (/ ('', i=1, 10) /))
     527
     528  type(ctrl_out),save :: o_fluxso4nff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     529  'fluxso4nff','Flux SO4 non-FF','mgS/m2/s', (/ ('', i=1, 10) /))
     530
     531  type(ctrl_out),save :: o_fluxso4bb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     532  'fluxso4bb','Flux SO4 BB','mgS/m2/s', (/ ('', i=1, 10) /))
     533
     534  type(ctrl_out),save :: o_fluxso4ba     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     535  'fluxso4ba','Flux SO4 Ba','mgS/m2/s', (/ ('', i=1, 10) /))
     536
     537  type(ctrl_out),save :: o_fluxso4     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     538  'fluxso4','Flux SO4','mgS/m2/s', (/ ('', i=1, 10) /))
     539
     540  type(ctrl_out),save :: o_fluxdms     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     541  'fluxdms','Flux DMS','mgS/m2/s', (/ ('', i=1, 10) /))
     542
     543  type(ctrl_out),save :: o_fluxh2sbio     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     544  'fluxh2sbio','Flux H2S Bio','mgS/m2/s', (/ ('', i=1, 10) /))
     545
     546  type(ctrl_out),save :: o_fluxdustec     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     547  'fluxdustec','Flux Dust EC','mg/m2/s', (/ ('', i=1, 10) /))
     548
     549  type(ctrl_out),save :: o_fluxddfine     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     550  'fluxddfine','DD Fine Mode','mg/m2/s', (/ ('', i=1, 10) /))
     551
     552  type(ctrl_out),save :: o_fluxddcoa     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     553  'fluxddcoa','DD Coarse Mode','mg/m2/s', (/ ('', i=1, 10) /))
     554
     555  type(ctrl_out),save :: o_fluxddsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     556  'fluxddsco','DD SCoarse Mode','mg/m2/s', (/ ('', i=1, 10) /))
     557
     558  type(ctrl_out),save :: o_fluxdd     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     559  'fluxdd','Flux DD','mg/m2/s', (/ ('', i=1, 10) /))
     560
     561  type(ctrl_out),save :: o_fluxssfine     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     562  'fluxssfine','SS Fine Mode','mg/m2/s', (/ ('', i=1, 10) /))
     563
     564  type(ctrl_out),save :: o_fluxsscoa     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     565  'fluxsscoa','SS Coarse Mode','mg/m2/s', (/ ('', i=1, 10) /))
     566
     567  type(ctrl_out),save :: o_fluxss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     568  'fluxss','Flux SS','mg/m2/s', (/ ('', i=1, 10) /))
     569
     570  type(ctrl_out),save :: o_flux_sparam_ind     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     571  'flux_sparam_ind','Ind emiss','mg/m2/s', (/ ('', i=1, 10) /))
     572
     573  type(ctrl_out),save :: o_flux_sparam_bb     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     574  'flux_sparam_bb','BB emiss','mg/m2/s', (/ ('', i=1, 10) /))
     575
     576  type(ctrl_out),save :: o_flux_sparam_ff     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     577  'flux_sparam_ff','FF emiss','mg/m2/s', (/ ('', i=1, 10) /))
     578
     579  type(ctrl_out),save :: o_flux_sparam_ddfine     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     580  'flux_sparam_ddfine','DD fine emiss','mg/m2/s', (/ ('', i=1, 10) /))
     581
     582  type(ctrl_out),save :: o_flux_sparam_ddcoa     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     583  'flux_sparam_ddcoa','DD coarse emiss','mg/m2/s', (/ ('', i=1, 10) /))
     584
     585  type(ctrl_out),save :: o_flux_sparam_ddsco     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     586  'flux_sparam_ddsco','DD Scoarse emiss','mg/m2/s', (/ ('', i=1, 10) /))
     587
     588  type(ctrl_out),save :: o_flux_sparam_ssfine     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     589  'flux_sparam_ssfine','SS fine emiss','mg/m2/s', (/ ('', i=1, 10) /))
     590
     591  type(ctrl_out),save :: o_flux_sparam_sscoa     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     592  'flux_sparam_sscoa','SS coarse emiss','mg/m2/s', (/ ('', i=1, 10) /))
     593
     594  type(ctrl_out),save :: o_u10m_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     595  'u10m_ss','Zonal wind at 10 m SS','m/s', (/ ('', i=1, 10) /))
     596
     597  type(ctrl_out),save :: o_v10m_ss     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     598  'v10m_ss','Meridional wind at 10 m SS','m/s', (/ ('', i=1, 10) /))
     599
     600!  type(ctrl_out),save :: o_     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10, 10 /), &
     601!  '','','', (/ ('', i=1, 10) /))
     602
     603!example  TYPE(ctrl_out), SAVE :: o_psbg = ctrl_out((/ 11, 11, 11, 11, 11, 11, 5, 5, 5, 11 /), &
    628604!       'psbg', 'Pressure sfce below ground', '%', (/ "inst(X)", "inst(X)", "inst(X)", &
    629 !       "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)"  /))
    630 
    631 
    632 
     605!       "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)"  /))
  • LMDZ5/trunk/libf/phylmd/Dust/spla_output_write.h

    r2630 r2752  
    66! traceur
    77
    8       CALL histwrite_phy( o_taue550          , diff_aod550_tot          )
    9       CALL histwrite_phy( o_taue670          , diag_aod670_tot   )
     8      CALL histwrite_phy( o_taue550          , diff_aod550_tot  )
     9      CALL histwrite_phy( o_taue670          , diag_aod670_tot  )
    1010      CALL histwrite_phy( o_taue865          , diag_aod865_tot  )
    1111      IF(id_fine>0) CALL histwrite_phy( o_taue550_tr2      , diff_aod550_tr2  )
    12       IF(id_fine>0) CALL histwrite_phy( o_taue670_tr2      , diag_aod670_tr2   )
     12      IF(id_fine>0) CALL histwrite_phy( o_taue670_tr2      , diag_aod670_tr2  )
    1313      IF(id_fine>0) CALL histwrite_phy( o_taue865_tr2      , diag_aod865_tr2  )
    14       IF(id_coss>0) CALL histwrite_phy( o_taue550_ss       , diag_aod550_ss       )
     14      IF(id_coss>0) CALL histwrite_phy( o_taue550_ss       , diag_aod550_ss   )
    1515      IF(id_coss>0) CALL histwrite_phy( o_taue670_ss       , diag_aod670_ss   )
    16       IF(id_coss>0) CALL histwrite_phy( o_taue865_ss       , diag_aod865_ss    )
    17       IF(id_codu>0) CALL histwrite_phy( o_taue550_dust     , diag_aod550_dust   )
     16      IF(id_coss>0) CALL histwrite_phy( o_taue865_ss       , diag_aod865_ss   )
     17      IF(id_codu>0) CALL histwrite_phy( o_taue550_dust     , diag_aod550_dust )
    1818      IF(id_codu>0) CALL histwrite_phy( o_taue670_dust     , diag_aod670_dust )
    1919      IF(id_codu>0) CALL histwrite_phy( o_taue865_dust     , diag_aod865_dust )
     
    2222      IF(id_scdu>0) CALL histwrite_phy( o_taue865_dustsco  , diag_aod865_dustsco )
    2323      CALL histwrite_phy( o_taue550_aqua     , aod550_aqua  )
    24       CALL histwrite_phy( o_taue550_terra    , aod550_terra    )
    25       CALL histwrite_phy( o_taue670_aqua     , aod670_aqua     )
    26       CALL histwrite_phy( o_taue670_terra    , aod670_terra    )
    27       CALL histwrite_phy( o_taue865_aqua     , aod865_aqua    )
    28       CALL histwrite_phy( o_taue865_terra    , aod865_terra    )
     24      CALL histwrite_phy( o_taue550_terra    , aod550_terra )
     25      CALL histwrite_phy( o_taue670_aqua     , aod670_aqua  )
     26      CALL histwrite_phy( o_taue670_terra    , aod670_terra )
     27      CALL histwrite_phy( o_taue865_aqua     , aod865_aqua  )
     28      CALL histwrite_phy( o_taue865_terra    , aod865_terra )
    2929
    3030      IF(id_fine>0) CALL histwrite_phy( o_taue550_fine_aqua     ,aod550_tr2_aqua   )
    3131      IF(id_fine>0) CALL histwrite_phy( o_taue670_fine_aqua     ,aod670_tr2_aqua   )
    3232      IF(id_fine>0) CALL histwrite_phy( o_taue865_fine_aqua     ,aod865_tr2_aqua   )
    33       IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_aqua      ,aod550_ss_aqua    )
    34       IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_aqua      ,aod670_ss_aqua    )
    35       IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_aqua      ,aod865_ss_aqua    )
    36       IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_aqua    ,aod550_dust_aqua    )
    37       IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_aqua    ,aod670_dust_aqua    )
    38       IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_aqua    ,aod865_dust_aqua    )
    39       IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_aqua ,aod670_dustsco_aqua    )
    40       IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_aqua ,aod550_dustsco_aqua    )
    41       IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_aqua ,aod865_dustsco_aqua    )
     33      IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_aqua      ,aod550_ss_aqua   )
     34      IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_aqua      ,aod670_ss_aqua   )
     35      IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_aqua      ,aod865_ss_aqua   )
     36      IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_aqua    ,aod550_dust_aqua   )
     37      IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_aqua    ,aod670_dust_aqua   )
     38      IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_aqua    ,aod865_dust_aqua   )
     39      IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_aqua ,aod670_dustsco_aqua   )
     40      IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_aqua ,aod550_dustsco_aqua   )
     41      IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_aqua ,aod865_dustsco_aqua   )
    4242     
    4343      IF(id_fine>0) CALL histwrite_phy( o_taue550_fine_terra     ,aod550_tr2_terra   )
    4444      IF(id_fine>0) CALL histwrite_phy( o_taue670_fine_terra     ,aod670_tr2_terra   )
    4545      IF(id_fine>0) CALL histwrite_phy( o_taue865_fine_terra     ,aod865_tr2_terra   )
    46       IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_terra      ,aod550_ss_terra    )
    47       IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_terra      ,aod670_ss_terra    )
    48       IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_terra      ,aod865_ss_terra    )
    49       IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_terra    ,aod550_dust_terra    )
    50       IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_terra    ,aod670_dust_terra    )
    51       IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_terra    ,aod865_dust_terra    )
    52       IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_terra ,aod670_dustsco_terra    )
    53       IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_terra ,aod550_dustsco_terra    )
    54       IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_terra ,aod865_dustsco_terra    )
     46      IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_terra      ,aod550_ss_terra   )
     47      IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_terra      ,aod670_ss_terra   )
     48      IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_terra      ,aod865_ss_terra   )
     49      IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_terra    ,aod550_dust_terra   )
     50      IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_terra    ,aod670_dust_terra   )
     51      IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_terra    ,aod865_dust_terra   )
     52      IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_terra ,aod670_dustsco_terra   )
     53      IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_terra ,aod550_dustsco_terra   )
     54      IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_terra ,aod865_dustsco_terra   )
    5555
    5656
     
    133133      IF(id_codu>0) CALL histwrite_phy( o_d_tr_uscav04 , d_tr_uscav04 )
    134134      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_uscav05 , d_tr_uscav05 )
    135       IF(id_prec>0) CALL histwrite_phy( o_d_tr_insc01 , d_tr_insc01 ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     135      IF(id_prec>0) CALL histwrite_phy( o_d_tr_insc01 , d_tr_insc01 )
    136136      IF(id_fine>0) CALL histwrite_phy( o_d_tr_insc02 , d_tr_insc02 )
    137137      IF(id_coss>0) CALL histwrite_phy( o_d_tr_insc03 , d_tr_insc03 )
     
    152152      IF(id_coss>0) CALL histwrite_phy( o_d_tr_ls03 , d_tr_ls03 )
    153153      IF(id_codu>0) CALL histwrite_phy( o_d_tr_ls04 , d_tr_ls04 )
    154       IF(id_scdu>0) CALL histwrite_phy( o_d_tr_ls05 , d_tr_ls05 )  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    155      
     154      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_ls05 , d_tr_ls05 )
    156155
    157156      IF(id_prec>0) CALL histwrite_phy( o_d_tr_dyn01 , d_tr_dyn01 )
     
    159158      IF(id_coss>0) CALL histwrite_phy( o_d_tr_dyn03 , d_tr_dyn03 )
    160159      IF(id_codu>0) CALL histwrite_phy( o_d_tr_dyn04 , d_tr_dyn04 )
    161       IF(id_scdu>0) CALL histwrite_phy( o_d_tr_dyn05 , d_tr_dyn05 )  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     160      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_dyn05 , d_tr_dyn05 )
    162161     
    163162      IF(id_prec>0) CALL histwrite_phy( o_d_tr_cl01 , d_tr_cl01 )
     
    165164      IF(id_coss>0) CALL histwrite_phy( o_d_tr_cl03 , d_tr_cl03 )
    166165      IF(id_codu>0) CALL histwrite_phy( o_d_tr_cl04 , d_tr_cl04 )
    167       IF(id_scdu>0) CALL histwrite_phy( o_d_tr_cl05 , d_tr_cl05 )  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     166      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_cl05 , d_tr_cl05 )
    168167      IF(id_prec>0) CALL histwrite_phy( o_d_tr_th01 , d_tr_th01 )
    169168      IF(id_fine>0) CALL histwrite_phy( o_d_tr_th02 , d_tr_th02 )
    170169      IF(id_coss>0) CALL histwrite_phy( o_d_tr_th03 , d_tr_th03 )
    171170      IF(id_codu>0) CALL histwrite_phy( o_d_tr_th04 , d_tr_th04 )
    172       IF(id_scdu>0) CALL histwrite_phy( o_d_tr_th05 , d_tr_th05 )  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    173 
     171      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_th05 , d_tr_th05 )
    174172
    175173      IF(id_coss>0) CALL histwrite_phy( o_sed_ss , sed_ss)
     
    232230      CALL histwrite_phy( o_u10m_ss             , u10m_ss              )
    233231      CALL histwrite_phy( o_v10m_ss             , v10m_ss              )
    234 
    235 !      CALL histwrite_phy(  ,)
    236 !
    237 
    238 
    239 
    240 
Note: See TracChangeset for help on using the changeset viewer.