Ignore:
Timestamp:
Aug 17, 2006, 5:41:51 PM (18 years ago)
Author:
Laurent Fairhead
Message:

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/V3_test/libf/phylmd/ini_histmth.h

    r684 r704  
    1313         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
    1414c
    15          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
    16          DO i = 1, iim
    17             zx_lon(i,1) = rlon(i+1)
    18             zx_lon(i,jjmp1) = rlon(i+1)
    19          ENDDO
    20 c
    21 cIM cf. AM 081204 BEG region
    22 c         imin_histmth=1
    23 c         imax_histmth=iim
    24 c         jmin_histmth=1
    25 c         jmax_histmth=jjmp1
    26 c         do i=1,iim-1
    27 c            if(zx_lon(i,1).lt.lonmin_histmth) imin_histmth=i
    28 c            if(zx_lon(i,1).le.lonmax_histmth) imax_histmth=i+1
    29 c         enddo
    30 cIM 220704 BAD??       do j=1,jjmp1-1
    31 c         do j=1,jjmp1
    32 c            if(zx_lat(1,j).ge.latmin_histmth) jmax_histmth=j
    33 c            if(zx_lat(1,j).gt.latmax_histmth) jmin_histmth=j
    34 c         enddo
    35 c
    36 c         print*,'On stoke le fichier instantanne sur, ',
    37 c     s   imin_histmth,imax_histmth,jmin_histmth,jmax_histmth
    38 c         print*,'On stoke le fichier instantanne sur, ',
    39 c     s   zx_lon(imin_histmth,1),zx_lon(imax_histmth,1),
    40 c     s   zx_lat(1,jmin_histmth),zx_lat(1,jmax_histmth)
    41 cIM cf. AM 081204 END region
    42 c
     15cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
     16cym         DO i = 1, iim
     17cym            zx_lon(i,1) = rlon(i+1)
     18cym            zx_lon(i,jjmp1) = rlon(i+1)
     19cym         ENDDO
    4320         DO ll=1,klev
    4421            znivsig(ll)=float(ll)
    4522         ENDDO
    46          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
    47          CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
    48      .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
     23cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
     24cym         CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
     25cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
     26cym     .                 nhori, nid_mth)
     27
     28           CALL histbeg_phy("histmth",itau_phy, zjulian, dtime,
    4929     .                 nhori, nid_mth)
    50 c
    51 cIM cf. AM 081204 BEG region
    52 c         CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
    53 c     .                 imin_histmth,imax_histmth-imin_histmth+1
    54 c     .                 ,jmin_histmth,jmax_histmth-jmin_histmth+1
    55 c     .                 , itau_phy, zjulian, dtime,
    56 c     .                 nhori, nid_mth)
    57 cIM 081204 END
    58 c
     30     
     31       write(*,*)'Mensuel ', itau_phy, zjulian
    5932         CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
    6033     .                 klev, presnivs/100., nvert)
    … …  
    6639         CALL histdef(nid_mth, "phis",
    6740     .                "Surface geop. height", "m2/s2",
    68      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     41     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    6942     .                "once",  zstophy,zout)
    7043c
    7144         CALL histdef(nid_mth, "aire", "Grid area", "-",
    72      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     45     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    7346     .                "once",  zstophy,zout)
    7447c
    7548         CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
    7649     .         "% "//clnsurf(is_ter), "%", 
    77      .         iim,jjmp1,nhori, 1,1,1, -99, 32,
     50     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    7851     .         "once", zstophy,zout)
    7952c
    8053         CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
    8154     .         "Fraction "//clnsurf(is_ter), "1", 
    82      .         iim,jjmp1,nhori, 1,1,1, -99, 32,
     55     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    8356     .         "once", zstophy,zout)
    8457c
    … …  
    8659c
    8760         CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
    88      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     61     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    8962     .                "ave(X)", zstophy,zout)
    9063c
    9164         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
    92      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     65     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    9366     .                "ave(X)", zstophy,zout)
    9467c
    9568         CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
    96      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     69     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    9770     .                "ave(X)", zstophy,zout)
    9871c
    … …  
    10073         CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
    10174     .                "K",
    102      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     75     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10376     .                t2mincels, zstophy,zout)
    10477c
    10578         CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
    10679     .                "K",
    107      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     80     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10881     .                t2maxcels, zstophy,zout)
    10982c
    11083c        CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
    11184c    .                "K",
    112 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     85c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11386c    .                "ave(X)", zstophy,zout)
    11487c
    11588         CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
    116      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     89     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11790     .                "ave(X)", zstophy,zout)
    11891c
    11992         CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
    120      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     93     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12194     .                "ave(X)", zstophy,zout)
    12295c
    12396c ENSEMBLES END
    12497         CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "kg/kg",
    125      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     98     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12699     .                "ave(X)", zstophy,zout)
    127100c
    128101         CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
    129      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     102     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    130103     .                "ave(X)", zstophy,zout)
    131104c
    132105         CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
    133      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     106     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    134107     .                "ave(X)", zstophy,zout)
    135108c
    136109         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
    137      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     110     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    138111     .                "ave(X)", zstophy,zout)
    139112c
    140113         CALL histdef(nid_mth, "qsurf", "Surface Air humidity", "kg/kg",
    141      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     114     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    142115     .                "ave(X)", zstophy,zout)
    143116c
    144117      if (.not. ok_veget) then
    145118         CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
    146      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     119     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    147120     .                "ave(X)", zstophy,zout)
    148121      endif
    … …  
    150123         CALL histdef(nid_mth, "ndayrain",
    151124     .                "Number of day with rain (liq+sol)", "-",
    152      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     125     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    153126     .                "inst(X)", zstomth,zout)
    154127c
    155128         CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
    156129     .                "kg/(s*m2)",
    157      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     130     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    158131     .                "ave(X)", zstophy,zout)
    159132c
    160133         CALL histdef(nid_mth, "plul", "Large-scale Precip.",
    161134     .   "kg/(s*m2)",
    162      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     135     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    163136     .                "ave(X)", zstophy,zout)
    164137c
    165138         CALL histdef(nid_mth, "pluc", "Convective Precip.",
    166139     .   "kg/(s*m2)",
    167      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     140     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    168141     .                "ave(X)", zstophy,zout)
    169142c
    170143         CALL histdef(nid_mth, "snow", "Snow fall", "kg/(s*m2)",
    171      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     144     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    172145     .                "ave(X)", zstophy,zout)
    173146c
    174147c        CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
    175 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     148c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    176149c    .                "ave(X)", zstophy,zout)
    177150c
    178151         CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)",
    179      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     152     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    180153     .                "ave(X)", zstophy,zout)
    181154c
    182155         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
    183      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     156     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    184157     .                "ave(X)", zstophy,zout)
    185158c
    186159         CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
    187      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     160     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    188161     .                "ave(X)", zstophy,zout)
    189162c
    190163         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
    191      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     164     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    192165     .                "ave(X)", zstophy,zout)
    193166c
    194167         CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
    195      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     168     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    196169     .                "ave(X)", zstophy,zout)
    197170c
    198171         CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
    199      .                iim,jjmp1,nhori, 1,1,1,-99,
     172     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    200173     .                32, "ave(X)", zstophy,zout)
    201174c
    202175         CALL histdef(nid_mth, "SWupTOAclr",
    203176     .                "SWup clear sky at TOA","W/m2",
    204      .                iim,jjmp1,nhori, 1,1,1,-99,
     177     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    205178     .                32, "ave(X)", zstophy,zout)
    206179c
    207180         CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
    208      .                iim,jjmp1,nhori, 1,1,1,-99,
     181     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    209182     .                32, "ave(X)", zstophy,zout)
    210183c
    211184         CALL histdef(nid_mth, "SWdnTOAclr",
    212185     .                "SWdn clear sky at TOA","W/m2",
    213      .                iim,jjmp1,nhori, 1,1,1,-99,
     186     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    214187     .                32, "ave(X)", zstophy,zout)
    215188c
    216189         CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
    217      .                iim,jjmp1,nhori, 1,1,1,-99,
     190     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    218191     .                32, "ave(X)", zstophy,zout)
    219192c
    220193         CALL histdef(nid_mth, "SWup200clr",
    221194     .                "SWup clear sky at 200mb","W/m2",
    222      .                iim,jjmp1,nhori, 1,1,1,-99,
     195     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    223196     .                32, "ave(X)", zstophy,zout)
    224197c
    225198         CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
    226      .                iim,jjmp1,nhori, 1,1,1,-99,
     199     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    227200     .                32, "ave(X)", zstophy,zout)
    228201c
    229202         CALL histdef(nid_mth, "SWdn200clr",
    230203     .                "SWdn clear sky at 200mb","W/m2",
    231      .                iim,jjmp1,nhori, 1,1,1,-99,
     204     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    232205     .                32, "ave(X)", zstophy,zout)
    233206c
    234207         CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
    235      .                iim,jjmp1,nhori, 1,1,1,-99,
     208     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    236209     .                32, "ave(X)", zstophy,zout)
    237210c
    238211         CALL histdef(nid_mth, "LWup200clr",
    239212     .                "LWup clear sky at 200mb","W/m2",
    240      .                iim,jjmp1,nhori, 1,1,1,-99,
     213     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    241214     .                32, "ave(X)", zstophy,zout)
    242215c
    243216         CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
    244      .                iim,jjmp1,nhori, 1,1,1,-99,
     217     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    245218     .                32, "ave(X)", zstophy,zout)
    246219c
    247220         CALL histdef(nid_mth, "LWdn200clr",
    248221     .                "LWdn clear sky at 200mb","W/m2",
    249      .                iim,jjmp1,nhori, 1,1,1,-99,
     222     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    250223     .                32, "ave(X)", zstophy,zout)
    251224c
    252225         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
    253      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     226     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    254227     .                "ave(X)", zstophy,zout)
    255228c
    256229         CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
    257      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     230     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    258231     .                "ave(X)", zstophy,zout)
    259232c
    260233         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
    261      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     234     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    262235     .                "ave(X)", zstophy,zout)
    263236c
    264237         CALL histdef(nid_mth, "radsol", "Rayonnement au sol", "W/m2",
    265      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     238     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    266239     .                "ave(X)", zstophy,zout)
    267240c
    268241         CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
    269      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     242     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    270243     .                "ave(X)", zstophy,zout)
    271244c
    272245         CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
    273      .                iim,jjmp1,nhori, 1,1,1,-99,
     246     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    274247     .                32, "ave(X)", zstophy,zout)
    275248c
    276249         CALL histdef(nid_mth, "SWupSFCclr",
    277250     .                "SWup clear sky at surface","W/m2",
    278      .                iim,jjmp1,nhori, 1,1,1,-99,
     251     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    279252     .                32, "ave(X)", zstophy,zout)
    280253c
    281254         CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
    282      .                iim,jjmp1,nhori, 1,1,1,-99,
     255     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    283256     .                32, "ave(X)", zstophy,zout)
    284257c
    285258         CALL histdef(nid_mth, "SWdnSFCclr",
    286259     .                "SWdn clear sky at surface","W/m2",
    287      .                iim,jjmp1,nhori, 1,1,1,-99,
     260     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    288261     .                32, "ave(X)", zstophy,zout)
    289262c
    290263         CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface",
    291      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     264     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    292265     .                "ave(X)", zstophy,zout)
    293266c
    294267         CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface",
    295      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     268     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    296269     .                "ave(X)", zstophy,zout)
    297270c
    298271         CALL histdef(nid_mth,"LWupSFCclr",
    299272     .                "CS Upwd. IR rad. at surface",
    300      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     273     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    301274     .                "ave(X)", zstophy,zout)
    302275c
    303276         CALL histdef(nid_mth,"LWdnSFCclr",
    304277     .                "Down. CS IR rad. at surface",
    305      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     278     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    306279     .                "ave(X)", zstophy,zout)
    307280c
    308281         CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
    309      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     282     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    310283     .                "ave(X)", zstophy,zout)
    311284c
    … …  
    315288c
    316289         CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
    317      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     290     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    318291     .                "ave(X)", zstophy,zout)
    319292c
    320293c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
    321 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     294c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    322295c    .                "ave(X)", zstophy,zout)
    323296c
    324297c        CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
    325 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     298c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    326299c    .                "ave(X)", zstophy,zout)
    327300c
    328301         CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
    329      .                "W/m2",iim,jjmp1,nhori, 1,1,1, -99, 32,
     302     .                "W/m2",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    330303     .                "ave(X)", zstophy,zout)
    331304c
    332305         CALL histdef(nid_mth, "fqcalving","Ice Calving",
    333      .                "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
    334      .                "ave(X)", zstophy,zout)
    335 c
     306     .                "kg/m2/s",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
     307     .                "ave(X)", zstophy,zout)
     308c
     309         CALL histdef(nid_mth, "fqfonte","Land ice melt",
     310     .                "kg/m2/s",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
     311     .                "ave(X)", zsto,zout)
     312
    336313         DO nsrf = 1, nbsrf
    337314           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
    338315     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
    339      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     316     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    340317     $         "ave(X)", zstophy,zout)
    341318
    342319           call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
    343320     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
    344      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     321     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    345322     $         "ave(X)", zstophy,zout)
    346323         ENDDO
    … …  
    351328           call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
    352329     $         "% "//clnsurf(nsrf), "%", 
    353      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     330     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    354331     $         "ave(X)", zstophy,zout)
    355332c
    356333           call histdef(nid_mth, "fract_"//clnsurf(nsrf),
    357334     $         "Fraction "//clnsurf(nsrf), "1", 
    358      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     335     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    359336     $         "ave(X)", zstophy,zout)
    360337          ENDIF !nsrf.GT.2
    … …  
    362339           call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
    363340     $         "Temperature "//clnsurf(nsrf), "K", 
    364      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     341     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    365342     $         "ave(X)", zstophy,zout)
    366343C
    367344           call histdef(nid_mth, "sens_"//clnsurf(nsrf),
    368345     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
    369      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     346     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    370347     $         "ave(X)", zstophy,zout)
    371348c
    372349           call histdef(nid_mth, "lat_"//clnsurf(nsrf),
    373350     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
    374      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     351     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    375352     $         "ave(X)", zstophy,zout)
    376353C
    377354           call histdef(nid_mth, "flw_"//clnsurf(nsrf),
    378355     $         "LW "//clnsurf(nsrf), "W/m2",
    379      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     356     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    380357     $         "ave(X)", zstophy,zout)
    381358c
    382359           call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
    383360     $         "SW "//clnsurf(nsrf), "W/m2",
    384      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     361     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    385362     $         "ave(X)", zstophy,zout)
    386363C
    387364           call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
    388365     $         "Bilan sol "//clnsurf(nsrf), "W/m2",
    389      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     366     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    390367     $         "ave(X)", zstophy,zout)
    391368C
    392369           call histdef(nid_mth, "wbilo_"//clnsurf(nsrf),
    393370     $         "Bilan eau "//clnsurf(nsrf), "kg/(m2*s)",
    394      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     371     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    395372     $         "ave(X)", zstophy,zout)
    396373C
    … …  
    400377           call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
    401378     $         "Evap*pourc "//clnsurf(nsrf), "kg/(s*m2)",
    402      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     379     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    403380     $         "ave(X)", zstophy,zout)
    404381C
    405382           call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
    406383     $         "Precip*pourc "//clnsurf(nsrf), "kg/(s*m2)",
    407      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     384     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    408385     $         "ave(X)", zstophy,zout)
    409386C
    … …  
    414391c
    415392         CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
    416      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     393     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    417394     .                "ave(X)", zstophy,zout)
    418395c
    419396         CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
    420      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     397     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    421398     .                "ave(X)", zstophy,zout)
    422399c
    423400         CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
    424      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     401     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    425402     .                "ave(X)", zstophy,zout)
    426403c
    427404         CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
    428      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     405     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    429406     .                "ave(X)", zstophy,zout)
    430407c
    431408         CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
    432      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     409     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    433410     .                "ave(X)", zstophy,zout)
    434411c
    435412         CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
    436      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     413     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    437414     .                "ave(X)", zstophy,zout)
    438415c
    439416         CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
    440      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     417     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    441418     .                "ave(X)", zstophy,zout)
    442419c
    443420         CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
    444      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     421     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    445422     .                "ave(X)", zstophy,zout)
    446423c
    447424         CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
    448      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     425     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    449426     .                "ave(X)", zstophy,zout)
    450427c
    451428         CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
    452      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     429     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    453430     .                "ave(X)", zstophy,zout)
    454431c
    455432         CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
    456      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     433     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    457434     .                "ave(X)", zstophy,zout)
    458435c
    459436         CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
    460      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     437     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    461438     .                "ave(X)", zstophy,zout)
    462439c
    463440         CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
    464      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     441     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    465442     .                "ave(X)", zstophy,zout)
    466443cKE43
    … …  
    468445c
    469446         CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
    470      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     447     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    471448     .                "ave(X)", zstophy,zout)
    472449c
    473450         CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
    474      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     451     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    475452     .                "ave(X)", zstophy,zout)
    476453c
    477454         CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
    478      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     455     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    479456     .                "ave(X)", zstophy,zout)
    480457c
    481458         CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
    482      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     459     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    483460     .                "ave(X)", zstophy,zout)
    484461c
    485462         CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
    486      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     463     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    487464     .                "ave(X)", zstophy,zout)
    488465c
    … …  
    493470c HBTM2
    494471         CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
    495      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     472     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    496473     .                "ave(X)", zstophy,zout)
    497474         CALL histdef(nid_mth, "s_pblt", "t at Boundary Layer Height",
    498475     .                "K",
    499      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     476     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    500477     .                "ave(X)", zstophy,zout)
    501478         CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
    502      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     479     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    503480     .                "ave(X)", zstophy,zout)
    504481         CALL histdef(nid_mth, "s_capCL", "Conv avlbl pot ener
    505482     .                for ABL", "J/m2",
    506      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     483     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    507484     .                "ave(X)", zstophy,zout)
    508485         CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
    509      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     486     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    510487     .                "ave(X)", zstophy,zout)
    511488         CALL histdef(nid_mth, "s_cteiCL", "Instability criteria
    512489     .                (ABL)", "K",
    513      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     490     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    514491     .                "ave(X)", zstophy,zout)
    515492         CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
    516      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     493     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    517494     .                "ave(X)", zstophy,zout)
    518495         CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
    519      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     496     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    520497     .                "ave(X)", zstophy,zout)
    521498         CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
    522      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     499     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    523500     .                "ave(X)", zstophy,zout)
    524501         CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
    525      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     502     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    526503     .                "ave(X)", zstophy,zout)
    527504c
    … …  
    545522          CALL histdef(nid_mth, "u"//bb2,
    546523     .                 "Zonal wind "//bb2//"mb","m/s",
    547      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     524     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    548525     .                "inst(X)", zout,zout)
    549526c
    550527          CALL histdef(nid_mth, "v"//bb2,
    551528     .                 "Meridional wind "//bb2//"mb","m/s",
    552      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     529     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    553530     .                "inst(X)", zout,zout)
    554531c
    555532          CALL histdef(nid_mth, "w"//bb2,
    556533     .                 "Vertical wind "//bb2//"mb","m/s",
    557      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     534     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    558535     .                "inst(X)", zout,zout)
    559536c
    560537          CALL histdef(nid_mth, "phi"//bb2,
    561538     .                 "Geopotential "//bb2//"mb","m",
    562      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     539     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    563540     .                "inst(X)", zout,zout)
    564541c
    … …  
    571548       CALL histdef(nid_mth, "fluxo",
    572549     .              "Flux turbulents ocean-atmosphere", "W/m2",
    573      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     550     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    574551     .              "ave(X)", zstophy,zout)
    575552c
    576553       CALL histdef(nid_mth, "fluxg",
    577554     .              "Flux turbulents ocean-glace de mer","W/m2",
    578      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     555     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    579556     .              "ave(X)", zstophy,zout)
    580557c
    581558       CALL histdef(nid_mth, "t_oce_sic",
    582559     .              "Temp. mixte oce-sic","K",
    583      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     560     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    584561     .              "ave(X)", zstophy,zout)
    585562c
    … …  
    588565        CALL histdef(nid_mth, "lmt_bils",
    589566     .       "Bilan au sol atmosphere forcee", "W/m2",
    590      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     567     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    591568     .              "ave(X)", zstophy,zout)
    592569c
    … …  
    595572        CALL histdef(nid_mth, "slab_bils",
    596573     .       "Bilan au sol Slab", "W/m2",
    597      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     574     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    598575     .              "ave(X)", zstophy,zout)
    599576c
    600577        CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
    601      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     578     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    602579     .              "ave(X)", zstophy,zout)
    603580c
    604581        CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
    605      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     582     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    606583     .              "ave(X)", zstophy,zout)
    607584c
    608585        CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
    609      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     586     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    610587     .              "ave(X)", zstophy,zout)
    611588c
    … …  
    621598         CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
    622599     .                "kg/kg",
    623      .                iim,jjmp1,nhori, klev,1,klev, nvert, 32,
     600     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    624601     .                "ave(X)", zstophy,zout)
    625602c
    626603         CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
    627      .                iim,jjmp1,nhori, klev,1,klev, nvert, 32,
     604     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    628605     .                "ave(X)", zstophy,zout)
    629606c
    630607         CALL histdef(nid_mth, "temp", "Air temperature", "K",
    631      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     608     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    632609     .                "ave(X)", zstophy,zout)
    633610c
    634611         CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
    635      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     612     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    636613     .                "ave(X)", zstophy,zout)
    637614c
    638615c        CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
    639 c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     616c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    640617c    .                "ave(X)", zstophy,zout)
    641618c
    642619         CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
    643      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     620     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    644621     .                "ave(X)", zstophy,zout)
    645622c
    646623         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
    647      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     624     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    648625     .                "ave(X)", zstophy,zout)
    649626c
    650627         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
    651      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     628     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    652629     .                "ave(X)", zstophy,zout)
    653630c
    654631         CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
    655      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     632     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    656633     .                "ave(X)", zstophy,zout)
    657634c
    658635         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
    659      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     636     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    660637     .                "ave(X)", zstophy,zout)
    661638c
    662639        CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
    663      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     640     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    664641     .                "ave(X)", zstophy,zout)
    665642c
    666643         CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
    667644     .                , "-",
    668      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     645     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    669646     .                "ave(X)", zstophy,zout)
    670647c
    671648         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
    672      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     649     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    673650     .                "ave(X)", zstophy,zout)
    674651c
    675652         CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
    676      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     653     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    677654     .                "ave(X)", zstophy,zout)
    678655c
    679656         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
    680      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     657     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    681658     .                "ave(X)", zstophy,zout)
    682659c
    683660         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
    684      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     661     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    685662     .                "ave(X)", zstophy,zout)
    686663c
    687664         CALL histdef(nid_mth, "dqphy", "Physics dQ", "(kg/kg)/s",
    688      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     665     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    689666     .                "ave(X)", zstophy,zout)
    690667c
    … …  
    697674           call histdef(nid_mth, "albe_"//clnsurf(nsrf),
    698675     $         "Albedo surf. "//clnsurf(nsrf), "-", 
    699      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     676     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    700677     $         "ave(X)", zstophy,zout)
    701678c
    702679           call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
    703680     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
    704      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     681     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    705682     $         "ave(X)", zstophy,zout)
    706683c
    707684         CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
    708      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     685     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    709686     .                "ave(X)", zstophy,zout)
    710687c
    … …  
    712689c
    713690         CALL histdef(nid_mth, "albs", "Surface albedo", "-",
    714      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     691     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    715692     .                "ave(X)", zstophy,zout)
    716693         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
    717      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     694     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    718695     .                "ave(X)", zstophy,zout)
    719696c
    … …  
    725702     .                "Convective Cloud Liquid water content"
    726703     .                , "kg/kg",
    727      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     704     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    728705     .                "ave(X)", zstophy,zout)
    729706c
    730707         CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
    731      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     708     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    732709     .                "ave(X)", zstophy,zout)
    733710c
    734711         CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
    735      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     712     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    736713     .                "ave(X)", zstophy,zout)
    737714c
    738715         CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
    739      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     716     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    740717     .                "ave(X)", zstophy,zout)
    741718c
    742719         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
    743      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     720     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    744721     .                "ave(X)", zstophy,zout)
    745722c
    746723         CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "(kg/kg)/s",
    747      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     724     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    748725     .                "ave(X)", zstophy,zout)
    749726c
    750727         CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
    751      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     728     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    752729     .                "ave(X)", zstophy,zout)
    753730c
    754731c        CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
    755 c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     732c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    756733c    .                "ave(X)", zstophy,zout)
    757734c
    758735         CALL histdef(nid_mth, "dqcon", "Convection dQ", "(kg/kg)/s",
    759      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     736     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    760737     .                "ave(X)", zstophy,zout)
    761738c
    762739         CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
    763      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     740     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    764741     .                "ave(X)", zstophy,zout)
    765742c
    766743         CALL histdef(nid_mth, "dtlschr",
    767      $       "Large-scale condensational heating rate", "K/s",iim,jjmp1
    768      $       ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
     744     $     "Large-scale condensational heating rate", "K/s",iim,jjphy_nb
     745     $     ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
    769746c
    770747         CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "(kg/kg)/s",
    771      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     748     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    772749     .                "ave(X)", zstophy,zout)
    773750c
    774751         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
    775      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     752     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    776753     .                "ave(X)", zstophy,zout)
    777754c
    778755         CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
    779      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     756     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    780757     .                "ave(X)", zstophy,zout)
    781758c
    782759         CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
    783      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     760     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    784761     .                "ave(X)", zstophy,zout)
    785762c
    786763         CALL histdef(nid_mth, "dqeva","Reevaporation dQ","(kg/kg)/s",
    787      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     764     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    788765     .                "ave(X)", zstophy,zout)
    789766
    790767         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
    791      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     768     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    792769     .                "ave(X)", zstophy,zout)
    793770
    794771         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
    795      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     772     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    796773     .                "ave(X)", zstophy,zout)
    797774
    798775c
    799776         CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
    800      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     777     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    801778     .                "ave(X)", zstophy,zout)
    802779
    803780         CALL histdef(nid_mth,"dqajs","Dry adjust. dQ","(kg/kg)/s",
    804      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     781     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    805782     .                "ave(X)", zstophy,zout)
    806783c
    807784         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
    808      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     785     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    809786     .                "ave(X)", zstophy,zout)
    810787c
    811788         CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
    812      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     789     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    813790     .                "ave(X)", zstophy,zout)
    814791c
    815792         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
    816      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     793     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    817794     .                "ave(X)", zstophy,zout)
    818795c
    819796         CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
    820      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     797     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    821798     .                "ave(X)", zstophy,zout)
    822799c
    823800         CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
    824      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     801     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    825802     .                "ave(X)", zstophy,zout)
    826803c
    827804         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
    828      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     805     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    829806     .                "ave(X)", zstophy,zout)
    830807c
    831808         CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
    832      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     809     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    833810     .                "ave(X)", zstophy,zout)
    834811c
    835812         IF (ok_orodr) THEN
    836813         CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
    837      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     814     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    838815     .                "ave(X)", zstophy,zout)
    839816c
    840817         CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
    841      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     818     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    842819     .                "ave(X)", zstophy,zout)
    843820c
    … …  
    846823         IF (ok_orolf) THEN
    847824         CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
    848      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     825     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    849826     .                "ave(X)", zstophy,zout)
    850827c
    851828         CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
    852      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     829     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    853830     .                "ave(X)", zstophy,zout)
    854831         ENDIF
    … …  
    859836         WRITE(str2,'(i2.2)') iq
    860837         CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
    861      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     838     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    862839     .                "ave(X)", zstophy,zout)
    863840         ELSE
    … …  
    875852c
    876853         CALL histdef(nid_mth, "phis", "Surface geop. height", "m",
    877      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     854     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    878855     .                "once",  zstophy,zout)
    879856c
    880857         CALL histdef(nid_mth, "aire", "Grid area", "-",
    881      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     858     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    882859     .                "once",  zstophy,zout)
    883860c
    884861         CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
    885862     .         "% "//clnsurf(is_ter), "%", 
    886      .         iim,jjmp1,nhori, 1,1,1, -99, 32,
     863     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    887864     .         "once", zstophy,zout)
    888865c
    889866         CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
    890867     .         "Fraction "//clnsurf(is_ter), "1", 
    891      .         iim,jjmp1,nhori, 1,1,1, -99, 32,
     868     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    892869     .         "once", zstophy,zout)
    893870c
    894871         CALL histdef(nid_mth, "pourc_"//clnsurf(is_lic),
    895872     .         "% "//clnsurf(is_lic), "%", 
    896      .         iim,jjmp1,nhori, 1,1,1, -99, 32,
     873     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    897874     .         "once", zstophy,zout)
    898875c
    899876         CALL histdef(nid_mth, "fract_"//clnsurf(is_lic),
    900877     .         "Fraction "//clnsurf(is_lic), "1", 
    901      .         iim,jjmp1,nhori, 1,1,1, -99, 32,
     878     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    902879     .         "once", zstophy,zout)
    903880c
    904881         CALL histdef(nid_mth, "pourc_"//clnsurf(is_sic),
    905882     .         "% "//clnsurf(is_sic), "%",
    906      .         iim,jjmp1,nhori, 1,1,1, -99, 32,
     883     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    907884     .                "ave(X)", zstophy,zout)
    908885c
    909886         CALL histdef(nid_mth, "fract_"//clnsurf(is_sic),
    910887     .         "Fraction "//clnsurf(is_sic), "1",
    911      .         iim,jjmp1,nhori, 1,1,1, -99, 32,
     888     .         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    912889     .                "ave(X)", zstophy,zout)
    913890c
    914891         CALL histdef(nid_mth, "flat", "Latent heat flux", "W/m2",
    915      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     892     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    916893     .                "ave(X)", zstophy,zout)
    917894c
    918895         CALL histdef(nid_mth, "bilTOA",
    919896     .                "Net radiation at TOA", "W/m2",
    920      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     897     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    921898     .                "ave(X)", zstophy,zout)
    922899c
    … …  
    925902c    .                "TOA reflected clear-sky shortwave
    926903c    .                radiation (method 1)", "W/m2",
    927 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     904c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    928905c    .                "ave(X)", zstophy,zout)
    929906c
    930907          CALL histdef(nid_mth, "dtdia","dtdia","K/s",
    931      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     908     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    932909     .                "ave(X)", zstophy,zout)
    933910c
    934911          CALL histdef(nid_mth, "dtsw","dtsw","K/s",
    935      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     912     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    936913     .                "ave(X)", zstophy,zout)
    937914c
    938915          CALL histdef(nid_mth, "dtlw","dtlw","K/s",
    939      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     916     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    940917     .                "ave(X)", zstophy,zout)
    941918c
    942919          CALL histdef(nid_mth, "dtcon","dtcon","K/s",
    943      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     920     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    944921     .                "ave(X)", zstophy,zout)
    945922c
    946923          CALL histdef(nid_mth, "dtajs","dtajs","K/s",
    947      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     924     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    948925     .                "ave(X)", zstophy,zout)
    949926c
    950927          CALL histdef(nid_mth, "dtlsc","dtlsc","K/s",
    951      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     928     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    952929     .                "ave(X)", zstophy,zout)
    953930c
    954931          CALL histdef(nid_mth, "dqdia","dqdia","(kg/kg)/s",
    955      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     932     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    956933     .                "ave(X)", zstophy,zout)
    957934c
    958935          CALL histdef(nid_mth, "duoli",
    959936     .         "duoli","m/s2",
    960      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     937     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    961938     .                "ave(X)", zstophy,zout)
    962939c
    963940          CALL histdef(nid_mth, "dvoli",
    964941     .         "dvoli","m/s2",
    965      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     942     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    966943     .                "ave(X)", zstophy,zout)
    967944c
    968945          CALL histdef(nid_mth, "ducon",
    969946     .         "ducon","m/s2",
    970      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     947     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    971948     .                "ave(X)", zstophy,zout)
    972949c
    973950          CALL histdef(nid_mth, "dvcon",
    974951     .         "dvcon","m/s2",
    975      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     952     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    976953     .                "ave(X)", zstophy,zout)
    977954c
    978955          CALL histdef(nid_mth, "dutot",
    979956     .         "dutot","m/s2",
    980      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     957     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    981958     .                "ave(X)", zstophy,zout)
    982959c
    983960          CALL histdef(nid_mth, "dvtot",
    984961     .         "dvtot","m/s2",
    985      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     962     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    986963     .                "ave(X)", zstophy,zout)
    987964c
    988965c         CALL histdef(nid_mth, "cldtrue",
    989966c    .         "TRUE Cloud fraction ","%",
    990 c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     967c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    991968c    .         "inst(X)", zout,zout)
    992969c
    993970c         CALL histdef(nid_mth, "cldTOA",
    994971c    .         "Cloud amount (satellite view)","%",
    995 c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     972c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    996973c    .         "ave(X)", zstophy,zout)
    997974c
    998975c         CALL histdef(nid_mth, "cldSFC",
    999976c    .         "Cloud amount (surface view)","%",
    1000 c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     977c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    1001978c    .         "ave(X)", zstophy,zout)
    1002979c
    1003980          CALL histdef(nid_mth, "GWDtaux",
    1004981     .         "GWD induced Ewd wind stress","Pa",
    1005      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     982     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    1006983     .         "ave(X)", zstophy,zout)
    1007984c
    1008985          CALL histdef(nid_mth, "GWDtauy",
    1009986     .         "GWD induced Nwd wind stress","Pa",
    1010      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     987     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    1011988     .         "ave(X)", zstophy,zout)
    1012989c
    … …  
    1014991c
    1015992         CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
    1016      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     993     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    1017994     .                "ave(X)", zstorad,zout)
    1018995c
    1019996         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
    1020      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     997     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    1021998     .                "ave(X)", zstophy,zout)
    1022999c
    10231000         CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
    1024      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1001     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10251002     .                "ave(X)", zstorad,zout)
    10261003c
    10271004         CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
    10281005     .                "K",
    1029      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1006     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10301007     .                t2mincels, zstophy,zout)
    10311008c
    10321009         CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
    10331010     .                "K",
    1034      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1011     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10351012     .                t2maxcels, zstophy,zout)
    10361013c
    10371014c        CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
    10381015c    .                "K",
    1039 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1016c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10401017c    .                "ave(X)", zstophy,zout)
    10411018c
    10421019         CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
    1043      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1020     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10441021     .                "ave(X)", zstorad,zout)
    10451022c
    10461023         CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
    1047      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1024     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10481025     .                "ave(X)", zstophy,zout)
    10491026
    10501027c
    10511028         CALL histdef(nid_mth, "q2m","Specific humidity 2m","kg/kg",
    1052      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1029     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10531030     .                "ave(X)", zstorad,zout)
    10541031c
    10551032         CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
    1056      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1033     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10571034     .                "ave(X)", zstorad,zout)
    10581035c
    10591036         CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
    1060      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1037     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10611038     .                "ave(X)", zstorad,zout)
    10621039c
    10631040         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
    1064      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1041     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10651042     .                "ave(X)", zstorad,zout)
    10661043c
    10671044         CALL histdef(nid_mth, "qsurf","Surface Air humidity","kg/kg",
    1068      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1045     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10691046     .                "ave(X)", zstophy,zout)
    10701047c
    10711048      if (.not. ok_veget) then
    10721049         CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
    1073      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1050     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10741051     .                "ave(X)", zstophy,zout)
    10751052      endif
    … …  
    10771054         CALL histdef(nid_mth, "ndayrain",
    10781055     .                "Number of day with rain (liq+sol)", "-",
    1079      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1056     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10801057     .                "inst(X)", zstomth,zout)
    10811058c
    10821059         CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
    10831060     .                "kg/(s*m2)",
    1084      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1061     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10851062     .                "ave(X)", zstophy,zout)
    10861063c
    10871064         CALL histdef(nid_mth, "plul", "Large-scale Precip.",
    10881065     .   "kg/(s*m2)",
    1089      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1066     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10901067     .                "ave(X)", zstophy,zout)
    10911068c
    10921069         CALL histdef(nid_mth, "pluc", "Convective Precip.",
    10931070     .   "kg/(s*m2)",
    1094      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1071     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10951072     .                "ave(X)", zstophy,zout)
    10961073c
    10971074         CALL histdef(nid_mth, "snowf", "Snow fall", "kg/(s*m2)",
    1098      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1075     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    10991076     .                "ave(X)", zstophy,zout)
    11001077c
    11011078c        CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
    1102 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1079c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11031080c    .                "ave(X)", zstophy,zout)
    11041081c
    11051082         CALL histdef(nid_mth, "evap",
    11061083     .                "Evaporation plus sublimation rate", "kg/(s*m2)",
    1107      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1084     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11081085     .                "ave(X)", zstophy,zout)
    11091086c
    11101087         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
    1111      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1088     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11121089     .                "ave(X)", zstophy,zout)
    11131090c
    11141091         CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
    1115      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1092     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11161093     .                "ave(X)", zstophy,zout)
    11171094c
    11181095         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
    1119      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1096     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11201097     .                "ave(X)", zstophy,zout)
    11211098c
    11221099         CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
    1123      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1100     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11241101     .                "ave(X)", zstophy,zout)
    11251102c
    11261103         CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
    1127      .                iim,jjmp1,nhori, 1,1,1,-99,
     1104     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11281105     .                32, "ave(X)", zstophy,zout)
    11291106c
    11301107         CALL histdef(nid_mth, "SWupTOAclr",
    11311108     .                "SWup clear sky at TOA","W/m2",
    1132      .                iim,jjmp1,nhori, 1,1,1,-99,
     1109     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11331110     .                32, "ave(X)", zstophy,zout)
    11341111c
    11351112         CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
    1136      .                iim,jjmp1,nhori, 1,1,1,-99,
     1113     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11371114     .                32, "ave(X)", zstophy,zout)
    11381115c
    11391116         CALL histdef(nid_mth, "SWdnTOAclr",
    11401117     .                "SWdn clear sky at TOA","W/m2",
    1141      .                iim,jjmp1,nhori, 1,1,1,-99,
     1118     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11421119     .                32, "ave(X)", zstophy,zout)
    11431120c
    11441121         CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
    1145      .                iim,jjmp1,nhori, 1,1,1,-99,
     1122     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11461123     .                32, "ave(X)", zstophy,zout)
    11471124c
    11481125         CALL histdef(nid_mth, "SWup200clr",
    11491126     .                "SWup clear sky at 200mb","W/m2",
    1150      .                iim,jjmp1,nhori, 1,1,1,-99,
     1127     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11511128     .                32, "ave(X)", zstophy,zout)
    11521129c
    11531130         CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
    1154      .                iim,jjmp1,nhori, 1,1,1,-99,
     1131     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11551132     .                32, "ave(X)", zstophy,zout)
    11561133c
    11571134         CALL histdef(nid_mth, "SWdn200clr",
    11581135     .                "SWdn clear sky at 200mb","W/m2",
    1159      .                iim,jjmp1,nhori, 1,1,1,-99,
     1136     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11601137     .                32, "ave(X)", zstophy,zout)
    11611138c
    11621139         CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
    1163      .                iim,jjmp1,nhori, 1,1,1,-99,
     1140     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11641141     .                32, "ave(X)", zstophy,zout)
    11651142c
    11661143         CALL histdef(nid_mth, "LWup200clr",
    11671144     .                "LWup clear sky at 200mb","W/m2",
    1168      .                iim,jjmp1,nhori, 1,1,1,-99,
     1145     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11691146     .                32, "ave(X)", zstophy,zout)
    11701147c
    11711148         CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
    1172      .                iim,jjmp1,nhori, 1,1,1,-99,
     1149     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11731150     .                32, "ave(X)", zstophy,zout)
    11741151c
    11751152         CALL histdef(nid_mth, "LWdn200clr",
    11761153     .                "LWdn clear sky at 200mb","W/m2",
    1177      .                iim,jjmp1,nhori, 1,1,1,-99,
     1154     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    11781155     .                32, "ave(X)", zstophy,zout)
    11791156c
    11801157         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
    1181      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1158     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11821159     .                "ave(X)", zstophy,zout)
    11831160c
    11841161         CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
    1185      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1162     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11861163     .                "ave(X)", zstophy,zout)
    11871164c
    11881165         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
    1189      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1166     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11901167     .                "ave(X)", zstophy,zout)
    11911168c
    11921169         CALL histdef(nid_mth, "radsol", "Rayonnement au sol", "W/m2",
    1193      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1170     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11941171     .                "ave(X)", zstophy,zout)
    11951172c
    11961173         CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
    1197      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1174     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    11981175     .                "ave(X)", zstophy,zout)
    11991176c
    12001177         CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
    1201      .                iim,jjmp1,nhori, 1,1,1,-99,
     1178     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    12021179     .                32, "ave(X)", zstophy,zout)
    12031180c
    12041181         CALL histdef(nid_mth, "SWupSFCclr",
    12051182     .                "SWup clear sky at surface","W/m2",
    1206      .                iim,jjmp1,nhori, 1,1,1,-99,
     1183     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    12071184     .                32, "ave(X)", zstophy,zout)
    12081185c
    12091186         CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
    1210      .                iim,jjmp1,nhori, 1,1,1,-99,
     1187     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    12111188     .                32, "ave(X)", zstophy,zout)
    12121189c
    12131190         CALL histdef(nid_mth, "SWdnSFCclr",
    12141191     .                "SWdn clear sky at surface","W/m2",
    1215      .                iim,jjmp1,nhori, 1,1,1,-99,
     1192     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    12161193     .                32, "ave(X)", zstophy,zout)
    12171194c
    12181195         CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface",
    1219      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     1196     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12201197     .                "ave(X)", zstophy,zout)
    12211198c
    12221199         CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface",
    1223      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     1200     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12241201     .                "ave(X)", zstophy,zout)
    12251202c
    12261203         CALL histdef(nid_mth,"LWupSFCclr",
    12271204     .                "CS Upwd. IR rad. at surface",
    1228      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     1205     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12291206     .                "ave(X)", zstophy,zout)
    12301207c
    12311208         CALL histdef(nid_mth,"LWdnSFCclr",
    12321209     .                "Down. CS IR rad. at surface",
    1233      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     1210     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12341211     .                "ave(X)", zstophy,zout)
    12351212c
    12361213         CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
    1237      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1214     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12381215     .                "ave(X)", zstophy,zout)
    12391216c
    12401217         CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
    1241      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1218     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12421219     .                "ave(X)", zstophy,zout)
    12431220c
    12441221         CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
    1245      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1222     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12461223     .                "ave(X)", zstophy,zout)
    12471224c
    12481225c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
    1249 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1226c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12501227c    .                "ave(X)", zstophy,zout)
    12511228c
    12521229c        CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
    1253 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1230c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12541231c    .                "ave(X)", zstophy,zout)
    12551232c
    12561233         CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
    1257      .                "W/m2",iim,jjmp1,nhori, 1,1,1, -99, 32,
     1234     .                "W/m2",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12581235     .                "ave(X)", zstophy,zout)
    12591236c
    12601237         CALL histdef(nid_mth, "fqcalving","Ice Calving",
    1261      .                "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
     1238     .                "kg/m2/s",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12621239     .                "ave(X)", zstophy,zout)
    12631240c
    … …  
    12651242           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
    12661243     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
    1267      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1244     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12681245     $         "ave(X)", zstophy,zout)
    12691246
    12701247           call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
    12711248     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
    1272      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1249     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12731250     $         "ave(X)", zstophy,zout)
    12741251         ENDDO
    … …  
    12761253           call histdef(nid_mth, "pourc_"//clnsurf(is_oce),
    12771254     $         "% "//clnsurf(is_oce), "%", 
    1278      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1255     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12791256     $         "ave(X)", zstophy,zout)
    12801257c
    12811258           call histdef(nid_mth, "fract_"//clnsurf(is_oce),
    12821259     $         "Fraction "//clnsurf(is_oce), "1", 
    1283      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1260     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12841261     $         "ave(X)", zstophy,zout)
    12851262C
    … …  
    12881265           call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
    12891266     $         "Temperature "//clnsurf(nsrf), "K", 
    1290      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1267     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12911268     $         "ave(X)", zstophy,zout)
    12921269C
    12931270           call histdef(nid_mth, "sens_"//clnsurf(nsrf),
    12941271     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
    1295      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1272     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    12961273     $         "ave(X)", zstophy,zout)
    12971274c
    12981275           call histdef(nid_mth, "lat_"//clnsurf(nsrf),
    12991276     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
    1300      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1277     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13011278     $         "ave(X)", zstophy,zout)
    13021279C
    13031280           call histdef(nid_mth, "flw_"//clnsurf(nsrf),
    13041281     $         "LW "//clnsurf(nsrf), "W/m2",
    1305      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1282     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13061283     $         "ave(X)", zstophy,zout)
    13071284c
    13081285           call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
    13091286     $         "SW "//clnsurf(nsrf), "W/m2",
    1310      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1287     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13111288     $         "ave(X)", zstophy,zout)
    13121289C
    13131290           call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
    13141291     $         "Bilan sol "//clnsurf(nsrf), "W/m2",
    1315      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1292     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13161293     $         "ave(X)", zstophy,zout)
    13171294C
    13181295           call histdef(nid_mth, "wbilo_"//clnsurf(nsrf),
    13191296     $         "Bilan eau "//clnsurf(nsrf), "kg/(m2*s)",
    1320      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1297     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13211298     $         "ave(X)", zstophy,zout)
    13221299C
    … …  
    13261303           call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
    13271304     $         "Evap*pourc "//clnsurf(nsrf), "kg/(s*m2)",
    1328      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1305     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13291306     $         "ave(X)", zstophy,zout)
    13301307C
    13311308           call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
    13321309     $         "Precip*pourc "//clnsurf(nsrf), "kg/(s*m2)",
    1333      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1310     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13341311     $         "ave(X)", zstophy,zout)
    13351312C
    … …  
    13401317c
    13411318         CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
    1342      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1319     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13431320     .                "ave(X)", zstophy,zout)
    13441321c
    13451322         CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
    1346      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1323     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13471324     .                "ave(X)", zstophy,zout)
    13481325c
    13491326         CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
    1350      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1327     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13511328     .                "ave(X)", zstophy,zout)
    13521329c
    13531330         CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
    1354      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1331     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13551332     .                "ave(X)", zstophy,zout)
    13561333c
    13571334         CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
    1358      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1335     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13591336     .                "ave(X)", zstophy,zout)
    13601337c
    13611338         CALL histdef(nid_mth, "cldt", "Total cloudiness", "%",
    1362      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1339     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13631340     .                "ave(X)", zstophy,zout)
    13641341c
    13651342         CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
    1366      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1343     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13671344     .                "ave(X)", zstophy,zout)
    13681345c
    13691346         CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
    1370      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1347     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13711348     .                "ave(X)", zstophy,zout)
    13721349c
    13731350         CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
    1374      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1351     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13751352     .                "ave(X)", zstophy,zout)
    13761353c
    13771354         CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
    1378      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1355     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13791356     .                "ave(X)", zstophy,zout)
    13801357c
    13811358         CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
    1382      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1359     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13831360     .                "ave(X)", zstophy,zout)
    13841361c
    13851362         CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
    1386      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1363     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13871364     .                "ave(X)", zstophy,zout)
    13881365c
    13891366         CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
    1390      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1367     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13911368     .                "ave(X)", zstophy,zout)
    13921369cKE43
    … …  
    13941371c
    13951372         CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
    1396      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1373     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    13971374     .                "ave(X)", zstophy,zout)
    13981375c
    13991376         CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
    1400      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1377     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14011378     .                "ave(X)", zstophy,zout)
    14021379c
    14031380         CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
    1404      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1381     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14051382     .                "ave(X)", zstophy,zout)
    14061383c
    14071384         CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
    1408      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1385     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14091386     .                "ave(X)", zstophy,zout)
    14101387c
    14111388         CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
    1412      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1389     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14131390     .                "ave(X)", zstophy,zout)
    14141391c
    14151392         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
    1416      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1393     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    14171394     .                "ave(X)", zstophy,zout)
    14181395c
    … …  
    14231400c HBTM2
    14241401         CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
    1425      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1402     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14261403     .                "ave(X)", zstophy,zout)
    14271404         CALL histdef(nid_mth, "s_pblt", "T at Boundary Layer Height",
    14281405     .                "K",
    1429      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1406     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14301407     .                "ave(X)", zstophy,zout)
    14311408         CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
    1432      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1409     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14331410     .                "ave(X)", zstophy,zout)
    14341411         CALL histdef(nid_mth,"s_capCL","Conv avlbl pot ener for ABL",
    14351412     .                "J/m2",
    1436      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1413     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14371414     .                "ave(X)", zstophy,zout)
    14381415         CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
    1439      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1416     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14401417     .                "ave(X)", zstophy,zout)
    14411418         CALL histdef(nid_mth,"s_cteiCL","Instability criteria (ABL)",
    14421419     .                "K",
    1443      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1420     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14441421     .                "ave(X)", zstophy,zout)
    14451422         CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
    1446      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1423     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14471424     .                "ave(X)", zstophy,zout)
    14481425         CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
    1449      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1426     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14501427     .                "ave(X)", zstophy,zout)
    14511428         CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
    1452      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1429     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14531430     .                "ave(X)", zstophy,zout)
    14541431         CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
    1455      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1432     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14561433     .                "ave(X)", zstophy,zout)
    14571434c
    … …  
    14751452          CALL histdef(nid_mth, "u"//bb2,
    14761453     .                 "Zonal wind "//bb2//"mb","m/s",
    1477      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1454     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14781455     .                "inst(X)", zout,zout)
    14791456c
    14801457          CALL histdef(nid_mth, "v"//bb2,
    14811458     .                 "Meridional wind "//bb2//"mb","m/s",
    1482      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1459     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14831460     .                "inst(X)", zout,zout)
    14841461c
    14851462          CALL histdef(nid_mth, "w"//bb2,
    14861463     .                 "Vertical wind "//bb2//"mb","m/s",
    1487      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1464     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14881465     .                "inst(X)", zout,zout)
    14891466c
    14901467          CALL histdef(nid_mth, "phi"//bb2,
    14911468     .                 "Geopotential "//bb2//"mb","m",
    1492      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1469     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    14931470     .                "inst(X)", zout,zout)
    14941471c
    … …  
    15021479       CALL histdef(nid_mth, "fluxo",
    15031480     .              "Flux turbulents ocean-atmosphere", "W/m2",
    1504      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     1481     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    15051482     .              "ave(X)", zstophy,zout)
    15061483c
    15071484       CALL histdef(nid_mth, "fluxg",
    15081485     .              "Flux turbulents ocean-glace de mer","W/m2",
    1509      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     1486     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    15101487     .              "ave(X)", zstophy,zout)
    15111488c
    15121489       CALL histdef(nid_mth, "t_oce_sic",
    15131490     .              "Temp. mixte oce-sic","K",
    1514      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     1491     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    15151492     .              "ave(X)", zstophy,zout)
    15161493c
    … …  
    15191496        CALL histdef(nid_mth, "lmt_bils",
    15201497     .       "Bilan au sol atmosphere forcee", "W/m2",
    1521      .       iim,jjmp1,nhori, 1,1,1, -99, 32,
     1498     .       iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    15221499     .       "ave(X)", zstophy,zout)
    15231500c
    … …  
    15261503        CALL histdef(nid_mth, "slab_bils",
    15271504     .              "Bilan au sol Slab", "W/m2",
    1528      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     1505     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    15291506     .              "ave(X)", zstophy,zout)
    15301507c
    15311508        CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
    1532      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     1509     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    15331510     .              "ave(X)", zstophy,zout)
    15341511c
    15351512        CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
    1536      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     1513     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    15371514     .              "ave(X)", zstophy,zout)
    15381515c
    15391516        CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
    1540      .              iim,jjmp1,nhori, 1,1,1, -99, 32,
     1517     .              iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    15411518     .              "ave(X)", zstophy,zout)
    15421519c
    … …  
    15521529         CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
    15531530     .                "kg/kg",
    1554      .                iim,jjmp1,nhori, klev,1,klev, nvert, 32,
     1531     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    15551532     .                "ave(X)", zstophy,zout)
    15561533c
    15571534         CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
    1558      .                iim,jjmp1,nhori, klev,1,klev, nvert, 32,
     1535     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    15591536     .                "ave(X)", zstophy,zout)
    15601537c
    15611538         CALL histdef(nid_mth, "temp", "Air temperature", "K",
    1562      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1539     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    15631540     .                "ave(X)", zstophy,zout)
    15641541c
    15651542         CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
    1566      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1543     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    15671544     .                "ave(X)", zstophy,zout)
    15681545c
    15691546c        CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
    1570 c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1547c    .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    15711548c    .                "ave(X)", zstophy,zout)
    15721549c
    15731550         CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
    1574      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1551     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    15751552     .                "ave(X)", zstophy,zout)
    15761553c
    15771554         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
    1578      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1555     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    15791556     .                "ave(X)", zstophy,zout)
    15801557c
    15811558         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
    1582      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1559     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    15831560     .                "ave(X)", zstophy,zout)
    15841561c
    15851562         CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
    1586      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1563     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    15871564     .                "ave(X)", zstophy,zout)
    15881565c
    15891566         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
    1590      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1567     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    15911568     .                "ave(X)", zstophy,zout)
    15921569c
    15931570        CALL histdef(nid_mth, "rneb", "Cloud fraction", "%",
    1594      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1571     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    15951572     .                "ave(X)", zstophy,zout)
    15961573c
    15971574         CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
    15981575     .                , "-",
    1599      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1576     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16001577     .                "ave(X)", zstophy,zout)
    16011578c
    16021579         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
    1603      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1580     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16041581     .                "ave(X)", zstophy,zout)
    16051582c
    16061583         CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
    1607      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1584     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16081585     .                "ave(X)", zstophy,zout)
    16091586c
    16101587         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
    1611      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1588     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16121589     .                "ave(X)", zstophy,zout)
    16131590c
    16141591         CALL histdef(nid_mth, "dqphy","Physics dQ","(kg/kg)/s",
    1615      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1592     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16161593     .                "ave(X)", zstophy,zout)
    16171594c
    … …  
    16241601           call histdef(nid_mth, "albe_"//clnsurf(nsrf),
    16251602     $         "Albedo surf. "//clnsurf(nsrf), "-", 
    1626      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1603     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    16271604     $         "ave(X)", zstophy,zout)
    16281605c
    16291606           call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
    16301607     $         "Rugosite "//clnsurf(nsrf), "m", 
    1631      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     1608     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    16321609     $         "ave(X)", zstophy,zout)
    16331610c
    16341611         CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
    1635      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1612     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    16361613     .                "ave(X)", zstophy,zout)
    16371614c
    … …  
    16391616c
    16401617         CALL histdef(nid_mth, "albs", "Surface albedo", "-",
    1641      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1618     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    16421619     .                "ave(X)", zstophy,zout)
    16431620         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
    1644      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1621     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    16451622     .                "ave(X)", zstophy,zout)
    16461623c
    … …  
    16521629     .                "Convective Cloud Liquid water content"
    16531630     .                , "kg/kg",
    1654      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1631     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16551632     .                "ave(X)", zstophy,zout)
    16561633c
    16571634         CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
    1658      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1635     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16591636     .                "ave(X)", zstophy,zout)
    16601637c
    16611638         CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
    1662      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1639     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16631640     .                "ave(X)", zstophy,zout)
    16641641c
    16651642         CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
    1666      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1643     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16671644     .                "ave(X)", zstophy,zout)
    16681645c
    16691646         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
    1670      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1647     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16711648     .                "ave(X)", zstophy,zout)
    16721649c
    16731650         CALL histdef(nid_mth, "dqdyn","Dynamics dQ","(kg/kg)/s",
    1674      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1651     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16751652     .                "ave(X)", zstophy,zout)
    16761653c
    16771654         CALL histdef(nid_mth, "dqcon","Convection dQ","(kg/kg)/s",
    1678      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1655     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16791656     .                "ave(X)", zstophy,zout)
    16801657c
    16811658         CALL histdef(nid_mth, "dtlschr",
    1682      $       "Large-scale condensational heating rate", "K/s",iim,jjmp1
     1659     $       "Large-scale condensational heating rate", "K/s",
     1660     $       iim,jjphy_nb
    16831661     $       ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
    16841662c
    16851663         CALL histdef(nid_mth, "dqlsc","Condensation dQ","(kg/kg)/s",
    1686      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1664     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16871665     .                "ave(X)", zstophy,zout)
    16881666c
    16891667         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
    1690      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1668     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16911669     .                "ave(X)", zstophy,zout)
    16921670c
    16931671         CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
    1694      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1672     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16951673     .                "ave(X)", zstophy,zout)
    16961674c
    16971675         CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
    1698      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1676     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    16991677     .                "ave(X)", zstophy,zout)
    17001678c
    17011679         CALL histdef(nid_mth,"dqeva","Reevaporation dQ","(kg/kg)/s",
    1702      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1680     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17031681     .                "ave(X)", zstophy,zout)
    17041682
    17051683         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
    1706      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1684     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17071685     .                "ave(X)", zstophy,zout)
    17081686
    17091687         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
    1710      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1688     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17111689     .                "ave(X)", zstophy,zout)
    17121690
    17131691c
    17141692         CALL histdef(nid_mth, "dqajs","Dry adjust. dQ","(kg/kg)/s",
    1715      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1693     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17161694     .                "ave(X)", zstophy,zout)
    17171695c
    17181696         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
    1719      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1697     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17201698     .                "ave(X)", zstophy,zout)
    17211699c
    17221700         CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
    1723      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1701     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17241702     .                "ave(X)", zstophy,zout)
    17251703c
    17261704         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
    1727      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1705     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17281706     .                "ave(X)", zstophy,zout)
    17291707c
    17301708         CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
    1731      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1709     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17321710     .                "ave(X)", zstophy,zout)
    17331711c
    17341712         CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
    1735      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1713     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17361714     .                "ave(X)", zstophy,zout)
    17371715c
    17381716         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
    1739      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1717     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17401718     .                "ave(X)", zstophy,zout)
    17411719c
    17421720         CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
    1743      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1721     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17441722     .                "ave(X)", zstophy,zout)
    17451723c
    17461724         IF (ok_orodr) THEN
    17471725         CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
    1748      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1726     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17491727     .                "ave(X)", zstophy,zout)
    17501728c
    17511729         CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
    1752      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1730     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17531731     .                "ave(X)", zstophy,zout)
    17541732c
    … …  
    17571735         IF (ok_orolf) THEN
    17581736         CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
    1759      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1737     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17601738     .                "ave(X)", zstophy,zout)
    17611739c
    17621740         CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
    1763      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1741     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17641742     .                "ave(X)", zstophy,zout)
    17651743         ENDIF
    … …  
    17691747         iiq=niadv(iq)
    17701748         CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
    1771      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1749     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    17721750     .                "ave(X)", zstophy,zout)
    17731751         ENDDO
Note: See TracChangeset for help on using the changeset viewer.