Ignore:
Timestamp:
Jun 4, 2007, 4:34:47 PM (17 years ago)
Author:
Laurent Fairhead
Message:

Merge entre la version V3_conv et le HEAD
YM, JG, LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/write_histmth.h

    r684 r766  
    1616c Champs 2D:
    1717c
    18       CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
    19       CALL histwrite(nid_mth,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     18cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
     19      CALL histwrite_phy(nid_mth,"phis",itau_w,pphis)
    2020C
    2121cIM 310804   CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
    22       CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
    23       CALL histwrite(nid_mth,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     22cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
     23      CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
    2424c
    2525      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100.
    26       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    27       CALL histwrite(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
    28      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     26cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     27      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
     28     $                   zx_tmp_fi2d)
    2929c
    3030      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)
    31       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    32       CALL histwrite(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
    33      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    34 c
    35       CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
    36       CALL histwrite(nid_mth,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     31cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     32      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
     33     $                   zx_tmp_fi2d)
     34c
     35cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     36      CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
    3737
    38       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    39       CALL histwrite(nid_mth,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    40 c
    41       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
    42       CALL histwrite(nid_mth,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     38cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     39      CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
     40c
     41cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
     42      CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
    4343c
    4444c ENSEMBLES BEG
    45       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    46       CALL histwrite(nid_mth,"t2m_min",itau_w,zx_tmp_2d,
    47      .               iim*jjmp1,ndex2d)
    48 c
    49       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    50       CALL histwrite(nid_mth,"t2m_max",itau_w,zx_tmp_2d,
    51      .               iim*jjmp1,ndex2d)
     45cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     46      CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
     47c
     48cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     49      CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
    5250c
    5351c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
    5452c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
    5553c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
    56 c     CALL histwrite(nid_mth,"tsoil",itau_w,zx_tmp_2d,
     54c     CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
    5755c    .               iim*jjmp1,ndex2d)
    5856c
     
    6058       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
    6159      ENDDO
    62       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    63       CALL histwrite(nid_mth,"wind10m",itau_w,zx_tmp_2d,
    64      $     iim*jjmp1,ndex2d)
     60cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     61      CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
    6562c
    6663      DO i = 1, klon
    6764         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
    6865      ENDDO
    69       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    70       CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     66cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     67      CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
    7168c
    7269c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    73 c     CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
     70c     CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
    7471c    .                                   iim*jjmp1*klev,ndex3d)
    7572c
    7673c ENSEMBLES END
    7774c
    78       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
    79       CALL histwrite(nid_mth,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    80 c
    81       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
    82       CALL histwrite(nid_mth,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    83 C
    84       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
    85       CALL histwrite(nid_mth,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     75cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
     76      CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
     77c
     78cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
     79      CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
     80C
     81cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
     82      CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
    8683C
    8784      DO i = 1, klon
    8885         zx_tmp_fi2d(i) = paprs(i,1)
    8986      ENDDO
    90       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    91       CALL histwrite(nid_mth,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    92 c
    93       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
    94       CALL histwrite(nid_mth,"qsurf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     87cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     88      CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
     89c
     90cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
     91      CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
    9592c
    9693      if (.not. ok_veget) then
    97           CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
    98           CALL histwrite(nid_mth,"qsol",itau_w,zx_tmp_2d,iim*jjmp1
    99      $        ,ndex2d)
     94cym          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
     95          CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
    10096      endif
    10197c
     
    10399         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
    104100      ENDDO
    105       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    106       CALL histwrite(nid_mth,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    107 c
    108       CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
    109       CALL histwrite(nid_mth,"ndayrain",itau_w,zx_tmp_2d,
    110      $               iim*jjmp1,ndex2d)
     101cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     102      CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
     103c
     104cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
     105      CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
    111106c
    112107      DO i = 1, klon
    113108         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
    114109      ENDDO
    115       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    116       CALL histwrite(nid_mth,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     110cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     111      CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
    117112c
    118113      DO i = 1, klon
    119114         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
    120115      ENDDO
    121       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    122       CALL histwrite(nid_mth,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    123 c
    124       CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
    125       CALL histwrite(nid_mth,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     116cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     117      CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
     118c
     119cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
     120      CALL histwrite_phy(nid_mth,"snow",itau_w,snow_fall)
    126121c
    127122c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
     
    129124c    .               ndex2d)
    130125c
    131       CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
    132       CALL histwrite(nid_mth,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    133 c
    134       CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
    135       CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    136 c
    137       CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
    138       CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    139 c
    140       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    141       CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    142 c
    143       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
    144       CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     126cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
     127      CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
     128c
     129cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
     130      CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
     131c
     132cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
     133      CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
     134c
     135cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
     136      CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
     137c
     138cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
     139      CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
    145140c
    146141      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
    147       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    148       CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
    149      .                               iim*jjmp1,ndex2d)
     142cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     143      CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d)
     144c
    150145c
    151146      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
    152       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    153       CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
    154      .                               iim*jjmp1,ndex2d)
     147cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     148      CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
    155149c
    156150      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
    157       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    158       CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
    159      .                               iim*jjmp1,ndex2d)
     151cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     152      CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
    160153c
    161154      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
    162       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    163       CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
    164      .                               iim*jjmp1,ndex2d)
    165 c
    166       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
    167       CALL histwrite(nid_mth,"SWup200",itau_w,zx_tmp_2d,
    168      $     iim*jjmp1,ndex2d)
    169 c
    170       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
    171       CALL histwrite(nid_mth,"SWup200clr",itau_w,zx_tmp_2d,
    172      $     iim*jjmp1,ndex2d)
    173 c
    174       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
    175       CALL histwrite(nid_mth,"SWdn200",itau_w,zx_tmp_2d,
    176      $     iim*jjmp1,ndex2d)
    177 c
    178       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
    179       CALL histwrite(nid_mth,"SWdn200clr",itau_w,zx_tmp_2d,
    180      $     iim*jjmp1,ndex2d)
    181 c
    182       CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
    183       CALL histwrite(nid_mth,"LWup200",itau_w,zx_tmp_2d,
    184      $     iim*jjmp1,ndex2d)
    185 c
    186       CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
    187       CALL histwrite(nid_mth,"LWup200clr",itau_w,zx_tmp_2d,
    188      $     iim*jjmp1,ndex2d)
     155cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     156      CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
     157c
     158cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
     159      CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
     160c
     161cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
     162      CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
     163c
     164cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
     165      CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
     166c
     167cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
     168      CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
     169c
     170cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
     171      CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
     172c
     173cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
     174      CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
    189175c
    190176      zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
    191       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    192       CALL histwrite(nid_mth,"LWdn200",itau_w,zx_tmp_2d,
    193      $     iim*jjmp1,ndex2d)
     177cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     178      CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
    194179c
    195180      zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
    196       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    197       CALL histwrite(nid_mth,"LWdn200clr",itau_w,zx_tmp_2d,
    198      $     iim*jjmp1,ndex2d)
     181cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     182      CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
    199183c
    200184c200 END
    201       CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
    202       CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    203 c
    204       CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
    205       CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    206 c
    207       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
    208       CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    209 c
    210       CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
    211       CALL histwrite(nid_mth,"radsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    212 c
    213       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
    214       CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     185cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
     186      CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
     187c
     188cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
     189      CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
     190c
     191cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
     192      CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
     193c
     194cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
     195      CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol)
     196c
     197cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
     198      CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
    215199c
    216200c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
    217201      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
    218       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    219       CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
    220      .                               iim*jjmp1,ndex2d)
     202cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     203      CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
    221204c
    222205c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
    223206      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
    224       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    225       CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
    226      .                               iim*jjmp1,ndex2d)
     207cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     208      CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
    227209c
    228210c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
    229211      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
    230       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    231       CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
    232      .                               iim*jjmp1,ndex2d)
     212cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     213      CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
    233214c
    234215c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
    235216      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
    236       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    237       CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
    238      .                               iim*jjmp1,ndex2d)
     217cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     218      CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
    239219c
    240220cIM 080904   zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
    241221cIM 130904   zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
    242222      zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
    243       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    244       CALL histwrite(nid_mth,"LWupSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    245      .               ndex2d)
     223cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     224      CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
    246225c
    247226      DO i=1, klon
     
    252231cIM 130904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon)
    253232      zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
    254       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    255       CALL histwrite(nid_mth,"LWupSFCclr",itau_w,zx_tmp_2d,iim*jjmp1,
    256      .               ndex2d)
    257 c
    258       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    259       CALL histwrite(nid_mth,"LWdnSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    260      .               ndex2d)
     233cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     234      CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
     235c
     236cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
     237      CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
    261238c
    262239c     DO i=1, klon
     
    265242c     ENDDO !i=1, klon
    266243c
    267       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
    268       CALL histwrite(nid_mth,"LWdnSFCclr",itau_w,zx_tmp_2d,
    269      $               iim*jjmp1,ndex2d)
    270 c
    271       CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
    272       CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     244cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
     245      CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
     246c
     247cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
     248      CALL histwrite_phy(nid_mth,"bils",itau_w,bils)
    273249c
    274250      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
    275251c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
    276       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    277       CALL histwrite(nid_mth,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    278 c
    279       CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
    280       CALL histwrite(nid_mth,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     252cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     253      CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d)
     254c
     255cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
     256      CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
    281257c
    282258c
     
    294270c
    295271c
    296        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
    297        CALL histwrite(nid_mth,"ffonte",itau_w,zx_tmp_2d,iim*jjmp1,
    298      $                ndex2d)
    299 c
    300        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
    301        CALL histwrite(nid_mth,"fqcalving",itau_w,zx_tmp_2d,iim*jjmp1,
    302      $                ndex2d)
    303 c
     272cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
     273       CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
     274c
     275cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
     276       CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
     277
     278cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqfonte,zx_tmp_2d)
     279       CALL histwrite_phy(nid_mth,"fqfonte",itau_w,zxfqfonte)
     280
    304281       DO nsrf = 1, nbsrf
    305282        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
    306         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    307         CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
    308      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     283cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     284        CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
     285     $      zx_tmp_fi2d)
    309286C     
    310287        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
    311         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    312         CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
    313      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     288cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     289        CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
     290     $      zx_tmp_fi2d)
    314291       ENDDO
    315292c
     
    318295       IF(nsrf.GT.2) THEN
    319296        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
    320         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    321         CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
    322      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     297cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     298        CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
     299     &                     zx_tmp_fi2d)
    323300c
    324301        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
    325         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    326         CALL histwrite(nid_mth,"fract_"//clnsurf(nsrf),itau_w,
    327      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     302cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     303        CALL histwrite_phy(nid_mth,"fract_"//clnsurf(nsrf),itau_w,
     304     $      zx_tmp_fi2d)
    328305       ENDIF !nsrf.GT.2
    329306C
    330307        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
    331         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    332         CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
    333      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     308cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     309        CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
     310     $      zx_tmp_fi2d)
    334311C
    335312        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
    336         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    337         CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
    338      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     313cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     314        CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
     315     $      zx_tmp_fi2d)
    339316C
    340317        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
    341         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    342         CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
    343      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     318cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     319        CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
     320     $      zx_tmp_fi2d)
    344321C
    345322        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
    346         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    347         CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
    348      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     323cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     324        CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
     325     $      zx_tmp_fi2d)
    349326c
    350327        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
    351         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    352         CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
    353      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     328cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     329        CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
     330     $      zx_tmp_fi2d)
    354331C
    355332        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
    356         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    357         CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
    358      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     333cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     334        CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
     335     $      zx_tmp_fi2d)
    359336C
    360337        zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
    361         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    362         CALL histwrite(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
    363      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     338cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     339        CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
     340     $      zx_tmp_fi2d)
    364341C
    365342       IF(1.EQ.0) THEN
    366343C
    367344        zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf)
    368         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    369         CALL histwrite(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
    370      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     345cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     346        CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
     347     $      zx_tmp_fi2d)
    371348C
    372349        zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon))
    373350     $                        *pctsrf(1:klon,nsrf)
    374         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    375         CALL histwrite(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
    376      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     351cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     352        CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
     353     $      zx_tmp_fi2d)
    377354C
    378355       ENDIF !(1.EQ.0) THEN
     
    380357      END DO 
    381358c
    382       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    383       CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    384 c
    385       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    386       CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    387 c
    388       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
    389       CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    390 c
    391       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
    392       CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    393 c
    394       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
    395       CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    396 c
    397       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    398       CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    399 c
    400       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
    401       CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     359cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
     360      CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
     361c
     362cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
     363      CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
     364c
     365cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
     366      CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
     367c
     368cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
     369      CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
     370c
     371cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
     372      CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
     373c
     374cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
     375      CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt)
     376c
     377cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
     378      CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
    402379c
    403380      zx_tmp_fi2d(1:klon) = flwp(1:klon)
    404       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    405       CALL histwrite(nid_mth,"lwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     381cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     382      CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
    406383c
    407384      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
    408       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    409       CALL histwrite(nid_mth,"iwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    410 c
    411       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
    412       CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    413 c
    414       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
    415       CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    416 c
    417       CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
    418       CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    419 c
    420       CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
    421       CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     385cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     386      CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
     387c
     388cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
     389      CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
     390c
     391cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
     392      CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
     393c
     394cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
     395      CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
     396c
     397cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
     398      CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
    422399cKE43
    423400      IF (iflag_con.GE.3) THEN ! sb
    424401c
    425       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
    426       CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    427 c
    428       CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
    429       CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    430 c
    431       CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
    432       CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    433 c
    434       CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
    435       CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     402cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
     403      CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
     404c
     405cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
     406      CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
     407c
     408cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
     409      CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
     410c
     411cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
     412      CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
    436413c
    437414c34EK
    438415c
    439       CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
    440       CALL histwrite(nid_mth,"prw",itau_w,zx_tmp_2d,
    441      .               iim*jjmp1,ndex2d)
     416cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
     417      CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
    442418c
    443419      ENDIF !iflag_con.GE.3
     
    448424c
    449425c
    450       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
    451       CALL histwrite(nid_mth,"s_pblh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    452 c
    453       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
    454       CALL histwrite(nid_mth,"s_pblt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    455 c
    456       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
    457       CALL histwrite(nid_mth,"s_lcl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    458 c
    459       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
    460       CALL histwrite(nid_mth,"s_capCL",itau_w,zx_tmp_2d,iim*jjmp1,
    461      .               ndex2d)
    462 c
    463       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
    464       CALL histwrite(nid_mth,"s_oliqCL",itau_w,zx_tmp_2d,iim*jjmp1,
    465      .               ndex2d)
    466 c
    467       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
    468       CALL histwrite(nid_mth,"s_cteiCL",itau_w,zx_tmp_2d,iim*jjmp1,
    469      .               ndex2d)
    470 c
    471       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
    472       CALL histwrite(nid_mth,"s_therm",itau_w,zx_tmp_2d,iim*jjmp1,
    473      .               ndex2d)
    474 c
    475       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
    476       CALL histwrite(nid_mth,"s_trmb1",itau_w,zx_tmp_2d,iim*jjmp1,
    477      .               ndex2d)
    478 c
    479       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
    480       CALL histwrite(nid_mth,"s_trmb2",itau_w,zx_tmp_2d,iim*jjmp1,
    481      .               ndex2d)
    482 c
    483       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
    484       CALL histwrite(nid_mth,"s_trmb3",itau_w,zx_tmp_2d,iim*jjmp1,
    485      .               ndex2d)
     426cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
     427      CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh)
     428c
     429cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
     430      CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt)
     431c
     432cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
     433      CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl)
     434c
     435cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
     436      CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL)
     437c
     438cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
     439      CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL)
     440c
     441cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
     442      CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL)
     443c
     444cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
     445      CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm)
     446c
     447cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
     448      CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1)
     449c
     450cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
     451      CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2)
     452c
     453cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
     454      CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3)
    486455c
    487456cIM cf. AM 081204 END
     
    502471     $   bb2.EQ."500".OR.bb2.EQ."200") THEN
    503472c
    504         CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
    505         CALL histwrite(nid_mth,"u"//bb2,itau_w,zx_tmp_2d,
    506      $       iim*jjmp1,ndex2d)
    507 c
    508         CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
    509         CALL histwrite(nid_mth,"v"//bb2,itau_w,zx_tmp_2d,
    510      $       iim*jjmp1,ndex2d)
    511 c
    512       CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
    513       CALL histwrite(nid_mth,"w"//bb2,itau_w,zx_tmp_2d,
    514      $     iim*jjmp1,ndex2d)
    515 c
    516       CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
    517       CALL histwrite(nid_mth,"phi"//bb2,itau_w,zx_tmp_2d,
    518      $     iim*jjmp1,ndex2d)
     473cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
     474        CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2))
     475c
     476cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
     477        CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2))
     478c
     479cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
     480      CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2))
     481c
     482cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
     483      CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2))
    519484c
    520485      ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
     
    529494      ENDDO
    530495c
    531       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    532       CALL histwrite(nid_mth,"fluxo",itau_w,zx_tmp_2d,
    533      .     iim*jjmp1,ndex2d)
     496cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     497      CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d)
    534498c
    535499      DO i=1, klon
     
    541505      ENDDO
    542506c
    543       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    544       CALL histwrite(nid_mth,"fluxg",itau_w,zx_tmp_2d,
    545      .     iim*jjmp1,ndex2d)
     507cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     508      CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d)
    546509c
    547510      DO i=1, klon
     
    556519      ENDDO
    557520c
    558       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    559       CALL histwrite(nid_mth,"t_oce_sic",itau_w,zx_tmp_2d,
    560      .     iim*jjmp1,ndex2d)
     521cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     522      CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d)
    561523c
    562524      IF (OCEAN.EQ.'force ') THEN
     
    571533      ENDIF
    572534      ENDDO
    573       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    574       CALL histwrite(nid_mth,"lmt_bils",itau_w,
    575      .     zx_tmp_2d,iim*jjmp1,ndex2d)
     535cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     536      CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,zx_tmp_fi2d)
    576537c
    577538      ELSE IF (OCEAN.EQ.'slab  ') THEN
     
    586547      ENDIF
    587548      ENDDO
    588       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    589       CALL histwrite(nid_mth,"slab_bils",itau_w,
    590      .     zx_tmp_2d,iim*jjmp1,ndex2d)
     549cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     550      CALL histwrite_phy(nid_mth,"slab_bils",itau_w,zx_tmp_fi2d)
    591551c
    592552      DO i=1, klon
     
    598558       ENDIF
    599559      ENDDO !i=1, klon
    600       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    601       CALL histwrite(nid_mth,"tslab",itau_w,zx_tmp_2d,
    602      .     iim*jjmp1,ndex2d)
    603 c
    604       CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
    605       CALL histwrite(nid_mth,"seaice",itau_w,zx_tmp_2d,
    606      .     iim*jjmp1,ndex2d)
    607 c
    608       CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
    609       CALL histwrite(nid_mth,"siceh",itau_w,zx_tmp_2d,
    610      .     iim*jjmp1,ndex2d)
     560cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     561      CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d)
     562c
     563cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
     564      CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice)
     565c
     566cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
     567      CALL histwrite_phy(nid_mth,"siceh",itau_w, seaice/1000.)
    611568c
    612569      ENDIF !OCEAN.EQ.force/slab
     
    619576C
    620577c
    621       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
    622       CALL histwrite(nid_mth,"lwcon",itau_w,zx_tmp_3d,
    623      .               iim*jjmp1*klev,ndex3d)
    624 c
    625       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
    626       CALL histwrite(nid_mth,"iwcon",itau_w,zx_tmp_3d,
    627      .               iim*jjmp1*klev,ndex3d)
    628 c
    629       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    630       CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
    631      .                                   iim*jjmp1*klev,ndex3d)
     578cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
     579      CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
     580c
     581cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
     582      CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
     583c
     584cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
     585      CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
    632586c
    633587c#ifdef histmthNMC
     
    635589c#endif
    636590c
    637       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
    638       CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
    639      .                                   iim*jjmp1*klev,ndex3d)
     591cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
     592      CALL histwrite_phy(nid_mth,"ovap",itau_w,qx(:,:,ivap))
    640593c
    641594c     zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
     
    645598c    .                                   iim*jjmp1*klev,ndex3d)
    646599c
    647       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    648       CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
    649      .                                   iim*jjmp1*klev,ndex3d)
    650 c
    651       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    652       CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
    653      .                                   iim*jjmp1*klev,ndex3d)
    654 c
    655       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    656       CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
    657      .                                   iim*jjmp1*klev,ndex3d)
    658 c
    659       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    660       CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
    661      .                                   iim*jjmp1*klev,ndex3d)
    662 c
    663       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    664       CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
    665      .                                   iim*jjmp1*klev,ndex3d)
    666 c
    667       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
    668       CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
    669      .                                   iim*jjmp1*klev,ndex3d)
    670 c
    671       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
    672       CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
    673      .                                   iim*jjmp1*klev,ndex3d)
    674 c
    675       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
    676       CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
    677      .                                   iim*jjmp1*klev,ndex3d)
    678 c
     600cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
     601      CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
     602c
     603cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
     604      CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
     605c
     606cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
     607      CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
     608c
     609cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
     610      CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
     611c
     612cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
     613      CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
     614c
     615cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
     616      CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra)
     617c
     618cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
     619      CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
     620c
     621cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
     622      CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
     623
    679624      DO k=1, klev
    680625       DO i=1, klon
     
    691636c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
    692637c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d)
    693       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    694       CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
    695      .                                   iim*jjmp1*klev,ndex3d)
    696 c
    697       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
    698       CALL histwrite(nid_mth,"dtphy",itau_w,zx_tmp_3d,
    699      .                                   iim*jjmp1*klev,ndex3d)
    700 c
    701       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),
    702      .                 zx_tmp_3d)
    703       CALL histwrite(nid_mth,"dqphy",itau_w,zx_tmp_3d,
    704      .                                   iim*jjmp1*klev,ndex3d)
     638cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     639      CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d)
     640c
     641cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
     642      CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
     643c
     644cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),
     645cym     .                 zx_tmp_3d)
     646      CALL histwrite_phy(nid_mth,"dqphy",itau_w, d_qx(:,:,ivap))
     647cIM 22.03.04 END
    705648c
    706649      ENDIF !lev_histmth.GE.2
     
    711654c
    712655        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
    713         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    714         CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
    715      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     656cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     657        CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
     658     $      zx_tmp_fi2d)
    716659C
    717660        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
    718         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    719         CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
    720      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     661cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     662        CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
     663     $      zx_tmp_fi2d)
    721664c
    722665      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
    723       CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
    724       CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w
    725      $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
     666cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
     667      CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
     668     $    ,zx_tmp_fi2d)
    726669c
    727670      ENDDO !nsrf=1, nbsrf
    728671c
    729       CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
    730       CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    731 c
    732       CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
    733       CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     672cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
     673      CALL histwrite_phy(nid_mth,"albs",itau_w,albsol)
     674c
     675cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
     676      CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw)
    734677c
    735678      ENDIF !lev_histmth.GE.3
     
    738681      IF(lev_histmth.GE.4) THEN
    739682c
    740       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
    741       CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
    742      .                                   iim*jjmp1*klev,ndex3d)
    743 c
     683cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
     684      CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
     685
    744686      IF (iflag_con.GE.3) THEN
    745687c
    746       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
    747       CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
    748      .                                   iim*jjmp1*klev,ndex3d)
    749 c
    750       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
    751       CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
    752      .                                   iim*jjmp1*klev,ndex3d)
    753 c
    754       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
    755       CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
    756      .                                   iim*jjmp1*klev,ndex3d)
    757 c
    758       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
    759       CALL histwrite(nid_mth,"dnwd0",itau_w,zx_tmp_3d,
    760      .                                   iim*jjmp1*klev,ndex3d)
     688cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
     689      CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
     690c
     691cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
     692      CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
     693c
     694cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
     695      CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
     696c
     697cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
     698      CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
    761699c
    762700      ENDIF !iflag_con.GE.3
    763701c
    764       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
    765       CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
    766      .                                   iim*jjmp1*klev,ndex3d)
    767 c
    768       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
    769       CALL histwrite(nid_mth,"dqdyn",itau_w,zx_tmp_3d,
    770      .                                   iim*jjmp1*klev,ndex3d)
     702cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
     703      CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
     704c
     705cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
     706      CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
    771707c
    772708cIM: 101003 : K/30min ==> K/s
    773709      zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
    774       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    775       CALL histwrite(nid_mth,"dtcon",itau_w,zx_tmp_3d,
    776      .                                   iim*jjmp1*klev,ndex3d)
     710cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     711      CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
    777712c
    778713      zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
    779       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    780       CALL histwrite(nid_mth,"dqcon",itau_w,zx_tmp_3d,
    781      .                                   iim*jjmp1*klev,ndex3d)
     714cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     715      CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
    782716c
    783717cIM: 101003 : K/30min ==> K/s
    784718      zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
    785       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    786       CALL histwrite(nid_mth,"dtlsc",itau_w,zx_tmp_3d,
    787      .                                   iim*jjmp1*klev,ndex3d)
     719cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     720      CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
     721cIM: 071003
    788722cIM: 101003 : K/30min ==> K/s
    789723      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
    790724     .                             d_t_eva(1:klon,1:klev))/pdtphys
    791       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
    792       CALL histwrite(nid_mth,"dtlschr",itau_w,zx_tmp_3d,
    793      .                                   iim*jjmp1*klev,ndex3d)
     725cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
     726      CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
    794727c
    795728      zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
    796       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    797       CALL histwrite(nid_mth,"dqlsc",itau_w,zx_tmp_3d,
    798      .                                   iim*jjmp1*klev,ndex3d)
     729cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     730      CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d)
    799731c
    800732cIM: 101003 : K/30min ==> K/s
    801733      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
    802       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    803       CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
    804      .                                   iim*jjmp1*klev,ndex3d)
     734cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     735      CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
    805736c
    806737      zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
    807       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    808       CALL histwrite(nid_mth,"dqvdf",itau_w,zx_tmp_3d,
    809      .                                   iim*jjmp1*klev,ndex3d)
     738cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     739      CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
    810740c
    811741cIM: 101003 : K/30min ==> K/s
    812742      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
    813       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    814       CALL histwrite(nid_mth,"dteva",itau_w,zx_tmp_3d,
    815      .                                   iim*jjmp1*klev,ndex3d)
     743cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     744      CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
    816745c
    817746      zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
    818       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    819       CALL histwrite(nid_mth,"dqeva",itau_w,zx_tmp_3d,
    820      .                                   iim*jjmp1*klev,ndex3d)
     747cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     748      CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d)
    821749c
    822750      zpt_conv = 0.
    823751      where (ptconv) zpt_conv = 1.
    824       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
    825       CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d,
    826      .                                   iim*(jjmp1)*klev,ndex3d)
    827 c
    828       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
    829       CALL histwrite(nid_mth,"ratqs",itau_w,zx_tmp_3d,
    830      .                                   iim*(jjmp1)*klev,ndex3d)
     752cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
     753      CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
     754c
     755cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
     756      CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
    831757c
    832758cIM: 101003 : K/30min ==> K/s
    833759      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
    834       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    835       CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
    836      .                                   iim*jjmp1*klev,ndex3d)
     760cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     761      CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
    837762c
    838763      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
    839       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    840       CALL histwrite(nid_mth,"dqajs",itau_w,zx_tmp_3d,
    841      .                                   iim*jjmp1*klev,ndex3d)
     764cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     765      CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d)
    842766c
    843767cIM: 101003 : K/day ==> K/s
    844768      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
    845       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    846       CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
    847      .                                   iim*jjmp1*klev,ndex3d)
     769cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     770      CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
    848771c
    849772cIM: 101003 : K/day ==> K/s     
    850773      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
    851       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    852       CALL histwrite(nid_mth,"dtsw0",itau_w,zx_tmp_3d,
    853      .                                   iim*jjmp1*klev,ndex3d)
     774cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     775      CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
    854776c
    855777cIM: 101003 : K/day ==> K/s     
    856778      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
    857       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    858       CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
    859      .                                   iim*jjmp1*klev,ndex3d)
     779cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     780      CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
    860781c
    861782cIM: 101003 : K/day ==> K/s     
    862783      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
    863       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    864       CALL histwrite(nid_mth,"dtlw0",itau_w,zx_tmp_3d,
    865      .                                   iim*jjmp1*klev,ndex3d)
     784cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     785      CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
    866786c
    867787cIM: 101003 : deja en K/s     
    868788      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
    869       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    870       CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
    871      .                                   iim*jjmp1*klev,ndex3d)
     789cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     790      CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
    872791c
    873792      zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
    874       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    875       CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
    876      .                                   iim*jjmp1*klev,ndex3d)
     793cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     794      CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
    877795c
    878796      zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
    879       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    880       CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
    881      .                                   iim*jjmp1*klev,ndex3d)
     797cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     798      CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
    882799c
    883800      IF (ok_orodr) THEN
    884801      zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
    885       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    886       CALL histwrite(nid_mth,"duoro",itau_w,zx_tmp_3d,
    887      .                                   iim*jjmp1*klev,ndex3d)
     802cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     803      CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
    888804c
    889805      zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
    890       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    891       CALL histwrite(nid_mth,"dvoro",itau_w,zx_tmp_3d,
    892      .                                   iim*jjmp1*klev,ndex3d)
     806cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     807      CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
    893808c
    894809      ENDIF
     
    896811      IF (ok_orolf) THEN
    897812      zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
    898       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    899       CALL histwrite(nid_mth,"dulif",itau_w,zx_tmp_3d,
    900      .                                   iim*jjmp1*klev,ndex3d)
     813cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     814      CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d)
    901815c
    902816      zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
    903       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    904       CALL histwrite(nid_mth,"dvlif",itau_w,zx_tmp_3d,
    905      .                                   iim*jjmp1*klev,ndex3d)
     817cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     818      CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
    906819      ENDIF
    907820C
    908821      IF (nqmax.GE.3) THEN
    909822      DO iq=3,nqmax
    910          CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
    911          CALL histwrite(nid_mth,tnom(iq),itau_w,zx_tmp_3d,
    912      .                                   iim*jjmp1*klev,ndex3d)
     823cym         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
     824         CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
    913825      ENDDO
    914826      ENDIF
     
    922834c Champs 2D:
    923835c
    924       CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis/RG,zx_tmp_2d)
    925       CALL histwrite(nid_mth,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     836cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis/RG,zx_tmp_2d)
     837      CALL histwrite_phy(nid_mth,"phis",itau_w,pphis/RG)
    926838C
    927839cIM 310804   CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
    928       CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
    929       CALL histwrite(nid_mth,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     840cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
     841      CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
    930842c
    931843      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100.
    932       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
    933       CALL histwrite(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
    934      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     844cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
     845      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
     846     $      zx_tmp_fi2d)
    935847c
    936848      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)
    937       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
    938       CALL histwrite(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
    939      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     849cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
     850      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
     851     $      zx_tmp_fi2d)
    940852c
    941853      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)* 100.
    942       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
    943       CALL histwrite(nid_mth,"pourc_"//clnsurf(is_lic),itau_w,
    944      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     854cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
     855      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_lic),itau_w,
     856     $      zx_tmp_fi2d)
    945857c
    946858      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)
    947       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
    948       CALL histwrite(nid_mth,"fract_"//clnsurf(is_lic),itau_w,
    949      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     859cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
     860      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_lic),itau_w,
     861     $      zx_tmp_fi2d)
    950862c
    951863      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)* 100.
    952       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    953       CALL histwrite(nid_mth,"pourc_"//clnsurf(is_sic),itau_w,
    954      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     864cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     865      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_sic),itau_w,
     866     $      zx_tmp_fi2d)
    955867c
    956868      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)
    957       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    958       CALL histwrite(nid_mth,"fract_"//clnsurf(is_sic),itau_w,
    959      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    960 c
    961       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
    962       CALL histwrite(nid_mth,"flat",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     869cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     870      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_sic),itau_w,
     871     $      zx_tmp_fi2d)
     872c
     873cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
     874      CALL histwrite_phy(nid_mth,"flat",itau_w,zxfluxlat)
    963875c
    964876      zx_tmp_fi2d(1 : klon) = topsw(1:klon)-toplw(1:klon)
    965       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    966       CALL histwrite(nid_mth,"bilTOA",itau_w,
    967      .               zx_tmp_2d,iim*jjmp1,ndex2d)
     877cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     878      CALL histwrite_phy(nid_mth,"bilTOA",itau_w,zx_tmp_fi2d)
    968879c
    969880c     IF(itap.EQ.ecrit_mth) THEN
     
    976887c       ENDDO
    977888c       CALL gr_fi_ecrit(1, klon,iim,jjmp1, tops01,zx_tmp_2d)
    978 c       CALL histwrite(nid_mth,"tops01",itau_w,zx_tmp_2d,
     889c       CALL histwrite_phy(nid_mth,"tops01",itau_w,zx_tmp_2d,
    979890c    .       iim*jjmp1,ndex2d)
    980891c     ENDIF !itap.EQ.ecrit_mth
     
    989900       ENDDO !l
    990901c
    991       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t,zx_tmp_3d)
    992       CALL histwrite(nid_mth,"dtdia",itau_w,zx_tmp_3d,
    993      $       iim*jjmp1*klev,ndex3d)
     902cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t,zx_tmp_3d)
     903      CALL histwrite_phy(nid_mth,"dtdia",itau_w,d_t)
    994904c
    995905c temperature tendency due to SW radiation
     
    1000910       ENDDO !l
    1001911c
    1002       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
    1003       CALL histwrite(nid_mth,"dtsw",itau_w,zx_tmp_3d,
    1004      $       iim*jjmp1*klev,ndex3d)
     912cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
     913      CALL histwrite_phy(nid_mth,"dtsw",itau_w,zx_tmp_fi3d)
    1005914c
    1006915c temperature tendency due to LW radiation
     
    1011920       ENDDO !l
    1012921c
    1013       CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
    1014       CALL histwrite(nid_mth,"dtlw",itau_w,zx_tmp_3d,
    1015      $       iim*jjmp1*klev,ndex3d)
     922cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
     923      CALL histwrite_phy(nid_mth,"dtlw",itau_w,zx_tmp_fi3d)
    1016924c
    1017925c  temperature tendency due to moist convective processes
     
    1022930       ENDDO !l
    1023931c
    1024       CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
    1025       CALL histwrite(nid_mth,"dtcon",itau_w,zx_tmp_3d,
    1026      $       iim*jjmp1*klev,ndex3d)
     932cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
     933      CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
    1027934c
    1028935c temperature tendency due to dry convective processes
     
    1033940       ENDDO !l
    1034941c
    1035       CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
    1036       CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
    1037      $       iim*jjmp1*klev,ndex3d)
     942cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
     943      CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
    1038944c
    1039945c  temperature tendency due to large scale precipitation
     
    1044950       ENDDO !l
    1045951c
    1046       CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
    1047       CALL histwrite(nid_mth,"dtlsc",itau_w,zx_tmp_3d,
    1048      $       iim*jjmp1*klev,ndex3d)
     952cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
     953      CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
    1049954c
    1050955c  total moisture tendency due to diabatic processes
     
    1056961       ENDDO !l
    1057962c
    1058       CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
    1059       CALL histwrite(nid_mth,"dqdia",itau_w,zx_tmp_3d,
    1060      $       iim*jjmp1*klev,ndex3d)
     963cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
     964      CALL histwrite_phy(nid_mth,"dqdia",itau_w,zx_tmp_fi3d)
    1061965c
    1062966cBEG 210404 BEG
     
    1071975      ENDIF !(ok_orodr.AND.ok_orolf) THEN
    1072976c
    1073       CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u_oli,zx_tmp_3d)
    1074       CALL histwrite(nid_mth,"duoli",itau_w,zx_tmp_3d,
    1075      $       iim*jjmp1*klev,ndex3d)
    1076 c
    1077       CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v_oli,zx_tmp_3d)
    1078       CALL histwrite(nid_mth,"dvoli",itau_w,zx_tmp_3d,
    1079      $       iim*jjmp1*klev,ndex3d)
     977cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u_oli,zx_tmp_3d)
     978      CALL histwrite_phy(nid_mth,"duoli",itau_w,d_u_oli)
     979c
     980cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v_oli,zx_tmp_3d)
     981      CALL histwrite_phy(nid_mth,"dvoli",itau_w,d_v_oli)
    1080982c
    1081983      zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
    1082       CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1083       CALL histwrite(nid_mth,"ducon",itau_w,zx_tmp_3d,
    1084      $       iim*jjmp1*klev,ndex3d)
     984cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     985      CALL histwrite_phy(nid_mth,"ducon",itau_w,zx_tmp_fi3d)
    1085986c
    1086987      zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
    1087       CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1088       CALL histwrite(nid_mth,"dvcon",itau_w,zx_tmp_3d,
    1089      $       iim*jjmp1*klev,ndex3d)
    1090 c
    1091       CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u,zx_tmp_3d)
    1092       CALL histwrite(nid_mth,"dutot",itau_w,zx_tmp_3d,
    1093      $       iim*jjmp1*klev,ndex3d)
    1094 c
    1095       CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v,zx_tmp_3d)
    1096       CALL histwrite(nid_mth,"dvtot",itau_w,zx_tmp_3d,
    1097      $       iim*jjmp1*klev,ndex3d)
     988cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     989      CALL histwrite_phy(nid_mth,"dvcon",itau_w,zx_tmp_fi3d)
     990c
     991cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u,zx_tmp_3d)
     992      CALL histwrite_phy(nid_mth,"dutot",itau_w,d_u)
     993c
     994cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v,zx_tmp_3d)
     995      CALL histwrite_phy(nid_mth,"dvtot",itau_w,d_v)
    1098996
    1099997cEND 210404 END
     
    1101999c cldtrue
    11021000c     CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldtrue,zx_tmp_3d)
    1103 c     CALL histwrite(nid_mth,"cldtrue",itau_w,zx_tmp_3d,
     1001c     CALL histwrite_phy(nid_mth,"cldtrue",itau_w,zx_tmp_3d,
    11041002c    $       iim*jjmp1*klev,ndex3d)
    11051003c
    11061004c     CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldTOA,zx_tmp_3d)
    1107 c     CALL histwrite(nid_mth,"cldTOA",itau_w,zx_tmp_3d,
     1005c     CALL histwrite_phy(nid_mth,"cldTOA",itau_w,zx_tmp_3d,
    11081006c    $       iim*jjmp1*klev,ndex3d)
    11091007c
    11101008c     CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldSFC,zx_tmp_3d)
    1111 c     CALL histwrite(nid_mth,"cldSFC",itau_w,zx_tmp_3d,
     1009c     CALL histwrite_phy(nid_mth,"cldSFC",itau_w,zx_tmp_3d,
    11121010c    $       iim*jjmp1*klev,ndex3d)
    11131011c
     
    11161014       zx_tmp_fi2d(i)=zustrdr(i)+zustrli(i)
    11171015      ENDDO
    1118       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1119       CALL histwrite(nid_mth,"GWDtaux",
    1120      $               itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1016cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1017      CALL histwrite_phy(nid_mth,"GWDtaux",itau_w,zx_tmp_fi2d)
    11211018c
    11221019      DO i=1, klon
     
    11241021       zx_tmp_fi2d(i)=zvstrdr(i)+zvstrli(i)
    11251022      ENDDO
    1126       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1127       CALL histwrite(nid_mth,"GWDtauy",
    1128      $               itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1023cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1024      CALL histwrite_phy(nid_mth,"GWDtauy",itau_w,zx_tmp_fi2d)
    11291025c
    11301026c220404 END
    11311027cAMIP2 END
    11321028c
    1133       CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
    1134       CALL histwrite(nid_mth,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1029cym     CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     1030      CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
    11351031
    1136       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    1137       CALL histwrite(nid_mth,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1138 c
    1139       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
    1140       CALL histwrite(nid_mth,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1032cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     1033      CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
     1034c
     1035cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
     1036      CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
    11411037c
    11421038c ENSEMBLES BEG
    1143       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    1144       CALL histwrite(nid_mth,"t2m_min",itau_w,zx_tmp_2d,
    1145      .               iim*jjmp1,ndex2d)
    1146 c
    1147       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    1148       CALL histwrite(nid_mth,"t2m_max",itau_w,zx_tmp_2d,
    1149      .               iim*jjmp1,ndex2d)
     1039cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     1040      CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
     1041c
     1042cym     CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     1043      CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
    11501044c
    11511045c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
    11521046c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
    11531047c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
    1154 c     CALL histwrite(nid_mth,"tsoil",itau_w,zx_tmp_2d,
     1048c     CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
    11551049c    .               iim*jjmp1,ndex2d)
    11561050c
     
    11581052       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
    11591053      ENDDO
    1160       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1161       CALL histwrite(nid_mth,"wind10m",itau_w,zx_tmp_2d,
    1162      $     iim*jjmp1,ndex2d)
     1054cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1055      CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
    11631056c
    11641057      DO i = 1, klon
    11651058         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
    11661059      ENDDO
    1167       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1168       CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1060cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1061      CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
    11691062c
    11701063c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    1171 c     CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
     1064c     CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
    11721065c    .                                   iim*jjmp1*klev,ndex3d)
    11731066c
    11741067c ENSEMBLES END
    11751068c
    1176       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
    1177       CALL histwrite(nid_mth,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1178 c
    1179       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
    1180       CALL histwrite(nid_mth,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1181 C
    1182       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
    1183       CALL histwrite(nid_mth,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1069cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
     1070      CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
     1071c
     1072cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
     1073      CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
     1074C
     1075cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
     1076      CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
    11841077C
    11851078      DO i = 1, klon
    11861079         zx_tmp_fi2d(i) = paprs(i,1)
    11871080      ENDDO
    1188       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1189       CALL histwrite(nid_mth,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1190 c
    1191       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
    1192       CALL histwrite(nid_mth,"qsurf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1081cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1082      CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
     1083c
     1084cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
     1085      CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
    11931086c
    11941087      if (.not. ok_veget) then
    1195           CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
    1196           CALL histwrite(nid_mth,"qsol",itau_w,zx_tmp_2d,iim*jjmp1
    1197      $        ,ndex2d)
     1088cym          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
     1089          CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
    11981090      endif
    11991091c
     
    12011093         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
    12021094      ENDDO
    1203       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1204       CALL histwrite(nid_mth,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1205 c
    1206       CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
    1207       CALL histwrite(nid_mth,"ndayrain",itau_w,zx_tmp_2d,
    1208      $               iim*jjmp1,ndex2d)
     1095cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1096      CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
     1097c
     1098cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
     1099      CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
    12091100c
    12101101      DO i = 1, klon
    12111102         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
    12121103      ENDDO
    1213       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1214       CALL histwrite(nid_mth,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1104cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1105      CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
    12151106c
    12161107      DO i = 1, klon
    12171108         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
    12181109      ENDDO
    1219       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1220       CALL histwrite(nid_mth,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1221 c
    1222       CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
    1223       CALL histwrite(nid_mth,"snowf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1110cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1111      CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
     1112c
     1113cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
     1114      CALL histwrite_phy(nid_mth,"snowf",itau_w,snow_fall)
    12241115c
    12251116c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
    1226 c     CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
     1117c     CALL histwrite_phy(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
    12271118c    .               ndex2d)
    12281119c
    1229       CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
    1230       CALL histwrite(nid_mth,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1231 c
    1232       CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
    1233       CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1234 c
    1235       CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
    1236       CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1237 c
    1238       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    1239       CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1240 c
    1241       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
    1242       CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1120cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
     1121      CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
     1122c
     1123cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
     1124      CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
     1125c
     1126cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
     1127      CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
     1128c
     1129cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
     1130      CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
     1131c
     1132cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
     1133      CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
    12431134c
    12441135      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
    1245       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1246       CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
    1247      .                               iim*jjmp1,ndex2d)
     1136cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1137      CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d)
    12481138c
    12491139c
    12501140      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
    1251       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1252       CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
    1253      .                               iim*jjmp1,ndex2d)
     1141cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1142      CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
    12541143c
    12551144      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
    1256       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1257       CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
    1258      .                               iim*jjmp1,ndex2d)
     1145cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1146      CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
    12591147c
    12601148      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
    1261       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1262       CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
    1263      .                               iim*jjmp1,ndex2d)
    1264 c
    1265       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
    1266       CALL histwrite(nid_mth,"SWup200",itau_w,zx_tmp_2d,
    1267      $     iim*jjmp1,ndex2d)
    1268 c
    1269       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
    1270       CALL histwrite(nid_mth,"SWup200clr",itau_w,zx_tmp_2d,
    1271      $     iim*jjmp1,ndex2d)
    1272 c
    1273       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
    1274       CALL histwrite(nid_mth,"SWdn200",itau_w,zx_tmp_2d,
    1275      $     iim*jjmp1,ndex2d)
    1276 c
    1277       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
    1278       CALL histwrite(nid_mth,"SWdn200clr",itau_w,zx_tmp_2d,
    1279      $     iim*jjmp1,ndex2d)
    1280 c
    1281       CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
    1282       CALL histwrite(nid_mth,"LWup200",itau_w,zx_tmp_2d,
    1283      $     iim*jjmp1,ndex2d)
    1284 c
    1285       CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
    1286       CALL histwrite(nid_mth,"LWup200clr",itau_w,zx_tmp_2d,
    1287      $     iim*jjmp1,ndex2d)
     1149cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1150      CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
     1151c
     1152cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
     1153      CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
     1154c
     1155cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
     1156      CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
     1157c
     1158cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
     1159      CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
     1160c
     1161cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
     1162      CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
     1163c
     1164cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
     1165      CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
     1166c
     1167cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
     1168      CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
    12881169c
    12891170      zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
    1290       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1291       CALL histwrite(nid_mth,"LWdn200",itau_w,zx_tmp_2d,
    1292      $     iim*jjmp1,ndex2d)
     1171cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1172      CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
    12931173c
    12941174      zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
    1295       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1296       CALL histwrite(nid_mth,"LWdn200clr",itau_w,zx_tmp_2d,
    1297      $     iim*jjmp1,ndex2d)
     1175cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1176      CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
    12981177c
    12991178c200 END
    1300       CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
    1301       CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1302 c
    1303       CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
    1304       CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1305 c
    1306       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
    1307       CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1308 c
    1309       CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
    1310       CALL histwrite(nid_mth,"radsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1311 c
    1312       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
    1313       CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1179cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
     1180      CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
     1181c
     1182cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
     1183      CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
     1184c
     1185cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
     1186      CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
     1187c
     1188cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
     1189      CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol)
     1190c
     1191cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
     1192      CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
    13141193c
    13151194c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
    13161195      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
    1317       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1318       CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
    1319      .                               iim*jjmp1,ndex2d)
     1196cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1197      CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
    13201198c
    13211199c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
    13221200      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
    1323       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1324       CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
    1325      .                               iim*jjmp1,ndex2d)
     1201cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1202      CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
    13261203c
    13271204c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
    13281205      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
    1329       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1330       CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
    1331      .                               iim*jjmp1,ndex2d)
     1206cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1207      CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
    13321208c
    13331209c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
    13341210      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
    1335       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1336       CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
    1337      .                               iim*jjmp1,ndex2d)
     1211cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1212      CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
    13381213c
    13391214cIM 080904   zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
    13401215cIM 130904   zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
    13411216      zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
    1342       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    1343       CALL histwrite(nid_mth,"LWupSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    1344      .               ndex2d)
     1217cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     1218      CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
    13451219c
    13461220      DO i=1, klon
     
    13511225cIM 130904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon)
    13521226      zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
    1353       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    1354       CALL histwrite(nid_mth,"LWupSFCclr",itau_w,zx_tmp_2d,iim*jjmp1,
    1355      .               ndex2d)
    1356 c
    1357       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    1358       CALL histwrite(nid_mth,"LWdnSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    1359      .               ndex2d)
     1227cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     1228      CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
     1229c
     1230cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
     1231      CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
    13601232c
    13611233c     DO i=1, klon
     
    13641236c     ENDDO !i=1, klon
    13651237c
    1366       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
    1367       CALL histwrite(nid_mth,"LWdnSFCclr",itau_w,zx_tmp_2d,
    1368      $               iim*jjmp1,ndex2d)
    1369 c
    1370       CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
    1371       CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1238cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
     1239      CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
     1240c
     1241cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
     1242      CALL histwrite_phy(nid_mth,"bils",itau_w,bils)
    13721243c
    13731244      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
    13741245c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
    1375       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1376       CALL histwrite(nid_mth,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1377 c
    1378       CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
    1379       CALL histwrite(nid_mth,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1246cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1247      CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d)
     1248c
     1249cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
     1250      CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
    13801251c
    13811252c
     
    13841255c      ENDDO
    13851256c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1386 c      CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1257c      CALL histwrite_phy(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    13871258c
    13881259c      DO i = 1, klon
     
    13901261c      ENDDO
    13911262c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1392 c      CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1393 c
    1394 c
    1395        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
    1396        CALL histwrite(nid_mth,"ffonte",itau_w,zx_tmp_2d,iim*jjmp1,
    1397      $                ndex2d)
    1398 c
    1399        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
    1400        CALL histwrite(nid_mth,"fqcalving",itau_w,zx_tmp_2d,iim*jjmp1,
    1401      $                ndex2d)
     1263c      CALL histwrite_phy(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1264c
     1265c
     1266cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
     1267       CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
     1268c
     1269cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
     1270       CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
    14021271c
    14031272       DO nsrf = 1, nbsrf
    14041273        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
    1405         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1406         CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
    1407      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1274cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1275        CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
     1276     $      zx_tmp_fi2d)
    14081277C     
    14091278        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
    1410         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1411         CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
    1412      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1279cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1280        CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
     1281     $      zx_tmp_fi2d)
    14131282       ENDDO
    14141283c
    14151284        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)*100.
    1416         CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1417         CALL histwrite(nid_mth,"pourc_"//clnsurf(is_oce),itau_w,
    1418      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1285cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1286        CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_oce),itau_w,
     1287     $      zx_tmp_fi2d)
    14191288c
    14201289        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)
    1421         CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    1422         CALL histwrite(nid_mth,"fract_"//clnsurf(is_oce),itau_w,
    1423      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1290cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     1291        CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_oce),itau_w,
     1292     $      zx_tmp_fi2d)
    14241293c
    14251294      DO nsrf = 1, nbsrf
     
    14271296C
    14281297        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
    1429         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1430         CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
    1431      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1298cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1299        CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
     1300     $      zx_tmp_fi2d)
    14321301C
    14331302        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
    1434         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1435         CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
    1436      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1303cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1304        CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
     1305     $      zx_tmp_fi2d)
    14371306C
    14381307        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
    1439         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1440         CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
    1441      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1308cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1309        CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
     1310     $      zx_tmp_fi2d)
    14421311C
    14431312        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
    1444         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1445         CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
    1446      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1313cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1314        CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
     1315     $      zx_tmp_fi2d)
    14471316c
    14481317        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
    1449         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1450         CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
    1451      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1318cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1319        CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
     1320     $      zx_tmp_fi2d)
    14521321C
    14531322        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
    1454         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1455         CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
    1456      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1323cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1324        CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
     1325     $      zx_tmp_fi2d)
    14571326C
    14581327        zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
    1459         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1460         CALL histwrite(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
    1461      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1328cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1329        CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
     1330     $      zx_tmp_fi2d)
    14621331C
    14631332       IF(1.EQ.0) THEN
    14641333C
    14651334        zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf)
    1466         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1467         CALL histwrite(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
    1468      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1335cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1336        CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
     1337     $      zx_tmp_fi2d)
    14691338C
    14701339        zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon))
    14711340     $                        *pctsrf(1:klon,nsrf)
    1472         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1473         CALL histwrite(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
    1474      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1341cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1342        CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
     1343     $      zx_tmp_fi2d)
    14751344C
    14761345       ENDIF !(1.EQ.0) THEN
     
    14781347      END DO 
    14791348c
    1480       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    1481       CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1482 c
    1483       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    1484       CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1485 c
    1486       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
    1487       CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1488 c
    1489       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
    1490       CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1491 c
    1492       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
    1493       CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1494 c
    1495       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt*100.,zx_tmp_2d)
    1496       CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1497 c
    1498       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
    1499       CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1349cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
     1350      CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
     1351c
     1352cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
     1353      CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
     1354c
     1355cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
     1356      CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
     1357c
     1358cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
     1359      CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
     1360c
     1361cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
     1362      CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
     1363c
     1364cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt*100.,zx_tmp_2d)
     1365      CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt*100.)
     1366c
     1367cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
     1368      CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
    15001369c
    15011370      zx_tmp_fi2d(1:klon) = flwp(1:klon)
    1502       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1503       CALL histwrite(nid_mth,"lwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1371cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1372      CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
    15041373c
    15051374      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
    1506       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1507       CALL histwrite(nid_mth,"iwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1508 c
    1509       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
    1510       CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1511 c
    1512       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
    1513       CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1514 c
    1515       CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
    1516       CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1517 c
    1518       CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
    1519       CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1375cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1376      CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
     1377c
     1378cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
     1379      CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
     1380c
     1381cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
     1382      CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
     1383c
     1384cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
     1385      CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
     1386c
     1387cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
     1388      CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
    15201389cKE43
    15211390      IF (iflag_con.GE.3) THEN ! sb
    15221391c
    1523       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
    1524       CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1525 c
    1526       CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
    1527       CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1528 c
    1529       CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
    1530       CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1531 c
    1532       CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
    1533       CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1392cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
     1393      CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
     1394c
     1395cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
     1396      CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
     1397c
     1398cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
     1399      CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
     1400c
     1401cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
     1402      CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
    15341403c
    15351404c34EK
    15361405c
    1537       CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
    1538       CALL histwrite(nid_mth,"prw",itau_w,zx_tmp_2d,
    1539      .               iim*jjmp1,ndex2d)
     1406cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
     1407      CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
    15401408c
    15411409      ENDIF !iflag_con.GE.3
     
    15451413c
    15461414c
    1547       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
    1548       CALL histwrite(nid_mth,"s_pblh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1549 c
    1550       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
    1551       CALL histwrite(nid_mth,"s_pblt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1552 c
    1553       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
    1554       CALL histwrite(nid_mth,"s_lcl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1555 c
    1556       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
    1557       CALL histwrite(nid_mth,"s_capCL",itau_w,zx_tmp_2d,iim*jjmp1,
    1558      .               ndex2d)
    1559 c
    1560       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
    1561       CALL histwrite(nid_mth,"s_oliqCL",itau_w,zx_tmp_2d,iim*jjmp1,
    1562      .               ndex2d)
    1563 c
    1564       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
    1565       CALL histwrite(nid_mth,"s_cteiCL",itau_w,zx_tmp_2d,iim*jjmp1,
    1566      .               ndex2d)
    1567 c
    1568       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
    1569       CALL histwrite(nid_mth,"s_therm",itau_w,zx_tmp_2d,iim*jjmp1,
    1570      .               ndex2d)
    1571 c
    1572       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
    1573       CALL histwrite(nid_mth,"s_trmb1",itau_w,zx_tmp_2d,iim*jjmp1,
    1574      .               ndex2d)
    1575 c
    1576       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
    1577       CALL histwrite(nid_mth,"s_trmb2",itau_w,zx_tmp_2d,iim*jjmp1,
    1578      .               ndex2d)
    1579 c
    1580       CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
    1581       CALL histwrite(nid_mth,"s_trmb3",itau_w,zx_tmp_2d,iim*jjmp1,
    1582      .               ndex2d)
     1415cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
     1416      CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh)
     1417c
     1418cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
     1419      CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt)
     1420c
     1421cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
     1422      CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl)
     1423c
     1424cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
     1425      CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL)
     1426c
     1427cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
     1428      CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL)
     1429c
     1430cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
     1431      CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL)
     1432c
     1433cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
     1434      CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm)
     1435c
     1436cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
     1437      CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1)
     1438c
     1439cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
     1440      CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2)
     1441c
     1442cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
     1443      CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3)
    15831444c
    15841445cIM cf. AM 081204 END
     
    16001461     $   .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
    16011462c
    1602         CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
    1603         CALL histwrite(nid_mth,"u"//bb2,itau_w,zx_tmp_2d,
    1604      $       iim*jjmp1,ndex2d)
    1605 c
    1606         CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
    1607         CALL histwrite(nid_mth,"v"//bb2,itau_w,zx_tmp_2d,
    1608      $       iim*jjmp1,ndex2d)
    1609 c
    1610       CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
    1611       CALL histwrite(nid_mth,"w"//bb2,itau_w,zx_tmp_2d,
    1612      $     iim*jjmp1,ndex2d)
    1613 c
    1614       CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
    1615       CALL histwrite(nid_mth,"phi"//bb2,itau_w,zx_tmp_2d,
    1616      $     iim*jjmp1,ndex2d)
     1463cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
     1464        CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2))
     1465c
     1466cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
     1467        CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2))
     1468c
     1469cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
     1470      CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2))
     1471c
     1472cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
     1473      CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2))
    16171474c
    16181475      ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
     
    16271484      ENDDO
    16281485c
    1629       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1630       CALL histwrite(nid_mth,"fluxo",itau_w,zx_tmp_2d,
    1631      .     iim*jjmp1,ndex2d)
     1486cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1487      CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d)
    16321488c
    16331489      DO i=1, klon
     
    16391495      ENDDO
    16401496c
    1641       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1642       CALL histwrite(nid_mth,"fluxg",itau_w,zx_tmp_2d,
    1643      .     iim*jjmp1,ndex2d)
     1497cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1498      CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d)
    16441499c
    16451500      DO i=1, klon
     
    16541509      ENDDO
    16551510c
    1656       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1657       CALL histwrite(nid_mth,"t_oce_sic",itau_w,zx_tmp_2d,
    1658      .     iim*jjmp1,ndex2d)
     1511cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1512      CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d)
    16591513c
    16601514      IF (OCEAN.EQ.'force ') THEN
     
    16691523      ENDIF
    16701524      ENDDO
    1671       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1672       CALL histwrite(nid_mth,"lmt_bils",itau_w,
    1673      .     zx_tmp_2d,iim*jjmp1,ndex2d)
     1525cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1526      CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,
     1527     .     zx_tmp_fi2d)
    16741528c
    16751529      ELSE IF (OCEAN.EQ.'slab  ') THEN
     
    16851539      ENDIF
    16861540      ENDDO
    1687       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1688       CALL histwrite(nid_mth,"slab_bils",itau_w,
    1689      .     zx_tmp_2d,iim*jjmp1,ndex2d)
     1541cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1542      CALL histwrite_phy(nid_mth,"slab_bils",itau_w,
     1543     .     zx_tmp_fi2d)
    16901544c
    16911545      DO i=1, klon
     
    16971551       ENDIF
    16981552      ENDDO !i=1, klon
    1699       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    1700       CALL histwrite(nid_mth,"tslab",itau_w,zx_tmp_2d,
    1701      .     iim*jjmp1,ndex2d)
    1702 c
    1703       CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
    1704       CALL histwrite(nid_mth,"seaice",itau_w,zx_tmp_2d,
    1705      .     iim*jjmp1,ndex2d)
    1706 c
    1707       CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
    1708       CALL histwrite(nid_mth,"siceh",itau_w,zx_tmp_2d,
    1709      .     iim*jjmp1,ndex2d)
     1553cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     1554      CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d)
     1555c
     1556cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
     1557      CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice)
     1558c
     1559cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
     1560      CALL histwrite_phy(nid_mth,"siceh",itau_w,seaice/1000.)
    17101561c
    17111562      ENDIF !OCEAN.EQ.force/slab
     
    17181569C
    17191570c
    1720       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
    1721       CALL histwrite(nid_mth,"lwcon",itau_w,zx_tmp_3d,
    1722      .               iim*jjmp1*klev,ndex3d)
    1723 c
    1724       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
    1725       CALL histwrite(nid_mth,"iwcon",itau_w,zx_tmp_3d,
    1726      .               iim*jjmp1*klev,ndex3d)
    1727 c
    1728       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    1729       CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
    1730      .                                   iim*jjmp1*klev,ndex3d)
     1571cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
     1572      CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
     1573c
     1574cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
     1575      CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
     1576c
     1577cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
     1578      CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
    17311579c
    17321580c#ifdef histmthNMC
     
    17341582c#endif
    17351583c
    1736       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
    1737       CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
    1738      .                                   iim*jjmp1*klev,ndex3d)
     1584cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
     1585      CALL histwrite_phy(nid_mth,"ovap",itau_w, qx(:,:,ivap))
    17391586c
    17401587c     zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
    17411588c    .                         (1-qx(1:klon,1:klev,ivap))
    17421589c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
    1743 c     CALL histwrite(nid_mth,"wvap",itau_w,zx_tmp_3d,
     1590c     CALL histwrite_phy(nid_mth,"wvap",itau_w,zx_tmp_3d,
    17441591c    .                                   iim*jjmp1*klev,ndex3d)
    17451592c
    1746       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    1747       CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
    1748      .                                   iim*jjmp1*klev,ndex3d)
    1749 c
    1750       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    1751       CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
    1752      .                                   iim*jjmp1*klev,ndex3d)
    1753 c
    1754       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    1755       CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
    1756      .                                   iim*jjmp1*klev,ndex3d)
    1757 c
    1758       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    1759       CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
    1760      .                                   iim*jjmp1*klev,ndex3d)
    1761 c
    1762       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    1763       CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
    1764      .                                   iim*jjmp1*klev,ndex3d)
    1765 c
    1766       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra*100., zx_tmp_3d)
    1767       CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
    1768      .                                   iim*jjmp1*klev,ndex3d)
    1769 c
    1770       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
    1771       CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
    1772      .                                   iim*jjmp1*klev,ndex3d)
    1773 c
    1774       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
    1775       CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
    1776      .                                   iim*jjmp1*klev,ndex3d)
     1593cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
     1594      CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
     1595c
     1596cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
     1597      CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
     1598c
     1599cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
     1600      CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
     1601c
     1602cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
     1603      CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
     1604c
     1605cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
     1606      CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
     1607c
     1608cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra*100., zx_tmp_3d)
     1609      CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra*100.)
     1610c
     1611cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
     1612      CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
     1613c
     1614cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
     1615      CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
    17771616c
    17781617      DO k=1, klev
     
    17901629c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
    17911630c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d)
    1792       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1793       CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
    1794      .                                   iim*jjmp1*klev,ndex3d)
    1795 c
    1796       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
    1797       CALL histwrite(nid_mth,"dtphy",itau_w,zx_tmp_3d,
    1798      .                                   iim*jjmp1*klev,ndex3d)
    1799 c
    1800       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),
    1801      .                 zx_tmp_3d)
    1802       CALL histwrite(nid_mth,"dqphy",itau_w,zx_tmp_3d,
    1803      .                                   iim*jjmp1*klev,ndex3d)
     1631cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1632      CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d)
     1633c
     1634cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
     1635      CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
     1636c
     1637cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),zx_tmp_3d)
     1638      CALL histwrite_phy(nid_mth,"dqphy",itau_w,d_qx(:,:,ivap))
    18041639c
    18051640      ENDIF !lev_histmth.GE.2
     
    18101645c
    18111646        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
    1812         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1813         CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
    1814      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1647cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1648        CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
     1649     $      zx_tmp_fi2d)
    18151650C
    18161651        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
    1817         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    1818         CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
    1819      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     1652cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     1653        CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
     1654     $      zx_tmp_fi2d)
    18201655c
    18211656      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
    1822       CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
    1823       CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w
    1824      $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
     1657cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
     1658      CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
     1659     $    ,agesno)
    18251660c
    18261661      ENDDO !nsrf=1, nbsrf
    18271662c
    1828       CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
    1829       CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1830 c
    1831       CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
    1832       CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     1663cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
     1664      CALL histwrite_phy(nid_mth,"albs",itau_w,albsol)
     1665c
     1666cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
     1667      CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw)
    18331668c
    18341669      ENDIF !lev_histmth.GE.3
     
    18371672      IF(lev_histmth.GE.4) THEN
    18381673c
    1839       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
    1840       CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
    1841      .                                   iim*jjmp1*klev,ndex3d)
     1674cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
     1675      CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
    18421676c
    18431677      IF (iflag_con.GE.3) THEN
    18441678c
    1845       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
    1846       CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
    1847      .                                   iim*jjmp1*klev,ndex3d)
    1848 c
    1849       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
    1850       CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
    1851      .                                   iim*jjmp1*klev,ndex3d)
    1852 c
    1853       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
    1854       CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
    1855      .                                   iim*jjmp1*klev,ndex3d)
    1856 c
    1857       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
    1858       CALL histwrite(nid_mth,"dnwd0",itau_w,zx_tmp_3d,
    1859      .                                   iim*jjmp1*klev,ndex3d)
     1679cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
     1680      CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
     1681c
     1682cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
     1683      CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
     1684c
     1685cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
     1686      CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
     1687c
     1688cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
     1689      CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
    18601690c
    18611691      ENDIF !iflag_con.GE.3
    18621692c
    1863       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
    1864       CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
    1865      .                                   iim*jjmp1*klev,ndex3d)
    1866 c
    1867       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
    1868       CALL histwrite(nid_mth,"dqdyn",itau_w,zx_tmp_3d,
    1869      .                                   iim*jjmp1*klev,ndex3d)
     1693cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
     1694      CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
     1695c
     1696cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
     1697      CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
    18701698c
    18711699cIM: 101003 : K/30min ==> K/s
    18721700c
    18731701      zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
    1874       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1875       CALL histwrite(nid_mth,"dqcon",itau_w,zx_tmp_3d,
    1876      .                                   iim*jjmp1*klev,ndex3d)
     1702cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1703      CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
    18771704c
    18781705cIM: 101003 : K/30min ==> K/s
    18791706      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
    18801707     .                             d_t_eva(1:klon,1:klev))/pdtphys
    1881       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
    1882       CALL histwrite(nid_mth,"dtlschr",itau_w,zx_tmp_3d,
    1883      .                                   iim*jjmp1*klev,ndex3d)
     1708cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
     1709      CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
    18841710c
    18851711      zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
    1886       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1887       CALL histwrite(nid_mth,"dqlsc",itau_w,zx_tmp_3d,
    1888      .                                   iim*jjmp1*klev,ndex3d)
     1712cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1713      CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d)
    18891714c
    18901715cIM: 101003 : K/30min ==> K/s
    18911716      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
    1892       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1893       CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
    1894      .                                   iim*jjmp1*klev,ndex3d)
     1717cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1718      CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
    18951719c
    18961720      zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
    1897       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1898       CALL histwrite(nid_mth,"dqvdf",itau_w,zx_tmp_3d,
    1899      .                                   iim*jjmp1*klev,ndex3d)
     1721cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1722      CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
    19001723c
    19011724cIM: 101003 : K/30min ==> K/s
    19021725      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
    1903       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1904       CALL histwrite(nid_mth,"dteva",itau_w,zx_tmp_3d,
    1905      .                                   iim*jjmp1*klev,ndex3d)
     1726cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1727      CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
    19061728c
    19071729      zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
    1908       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1909       CALL histwrite(nid_mth,"dqeva",itau_w,zx_tmp_3d,
    1910      .                                   iim*jjmp1*klev,ndex3d)
     1730cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1731      CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d)
    19111732c
    19121733      zpt_conv = 0.
    19131734      where (ptconv) zpt_conv = 1.
    1914       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
    1915       CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d,
    1916      .                                   iim*(jjmp1)*klev,ndex3d)
    1917 c
    1918       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
    1919       CALL histwrite(nid_mth,"ratqs",itau_w,zx_tmp_3d,
    1920      .                                   iim*(jjmp1)*klev,ndex3d)
     1735cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
     1736      CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
     1737c
     1738cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
     1739      CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
    19211740c
    19221741      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
    1923       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1924       CALL histwrite(nid_mth,"dqajs",itau_w,zx_tmp_3d,
    1925      .                                   iim*jjmp1*klev,ndex3d)
     1742cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1743      CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d)
    19261744c
    19271745cIM: 101003 : K/day ==> K/s
    19281746      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
    1929       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1930       CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
    1931      .                                   iim*jjmp1*klev,ndex3d)
     1747cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1748      CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
    19321749c
    19331750cIM: 101003 : K/day ==> K/s     
    19341751      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
    1935       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1936       CALL histwrite(nid_mth,"dtsw0",itau_w,zx_tmp_3d,
    1937      .                                   iim*jjmp1*klev,ndex3d)
     1752cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1753      CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
    19381754c
    19391755cIM: 101003 : K/day ==> K/s     
    19401756      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
    1941       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1942       CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
    1943      .                                   iim*jjmp1*klev,ndex3d)
     1757cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1758      CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
    19441759c
    19451760cIM: 101003 : K/day ==> K/s     
    19461761      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
    1947       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1948       CALL histwrite(nid_mth,"dtlw0",itau_w,zx_tmp_3d,
    1949      .                                   iim*jjmp1*klev,ndex3d)
     1762cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1763      CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
    19501764c
    19511765cIM: 101003 : deja en K/s     
    19521766      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
    1953       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1954       CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
    1955      .                                   iim*jjmp1*klev,ndex3d)
     1767cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1768      CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
    19561769c
    19571770      zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
    1958       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1959       CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
    1960      .                                   iim*jjmp1*klev,ndex3d)
     1771cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1772      CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
    19611773c
    19621774      zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
    1963       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1964       CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
    1965      .                                   iim*jjmp1*klev,ndex3d)
     1775cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1776      CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
    19661777c
    19671778      IF (ok_orodr) THEN
    19681779      zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
    1969       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1970       CALL histwrite(nid_mth,"duoro",itau_w,zx_tmp_3d,
    1971      .                                   iim*jjmp1*klev,ndex3d)
     1780cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1781      CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
    19721782c
    19731783      zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
    1974       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1975       CALL histwrite(nid_mth,"dvoro",itau_w,zx_tmp_3d,
    1976      .                                   iim*jjmp1*klev,ndex3d)
     1784cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1785      CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
    19771786c
    19781787      ENDIF
     
    19801789      IF (ok_orolf) THEN
    19811790      zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
    1982       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1983       CALL histwrite(nid_mth,"dulif",itau_w,zx_tmp_3d,
    1984      .                                   iim*jjmp1*klev,ndex3d)
     1791cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1792      CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d)
    19851793c
    19861794      zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
    1987       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    1988       CALL histwrite(nid_mth,"dvlif",itau_w,zx_tmp_3d,
    1989      .                                   iim*jjmp1*klev,ndex3d)
     1795cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     1796      CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
    19901797      ENDIF
    19911798C
    19921799      IF (nqmax.GE.3) THEN
    19931800      DO iq=3,nqmax
    1994          CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
    1995          CALL histwrite(nid_mth,tnom(iq),itau_w,zx_tmp_3d,
    1996      .                                   iim*jjmp1*klev,ndex3d)
     1801cym         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
     1802         CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
    19971803      ENDDO
    19981804      ENDIF
     
    20031809c
    20041810      if (ok_sync) then
     1811c$OMP MASTER
    20051812        call histsync(nid_mth)
     1813c$OMP END MASTER
    20061814      endif
    20071815c
Note: See TracChangeset for help on using the changeset viewer.