Ignore:
Timestamp:
May 4, 2005, 5:11:29 PM (19 years ago)
Author:
Laurent Fairhead
Message:

Modifications faites à la physique pour la rendre parallele YM
Une branche de travail LMDZ4_par_0 a été créée provisoirement afin de tester
les modifs pleinement avant leurs inclusions dans le tronc principal
LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4_par_0/libf/phylmd/write_histmth.h

    r633 r634  
    1717c Champs 2D:
    1818c
    19       CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
    20       CALL histwrite(nid_mth,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     19cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
     20      CALL histwrite_phy(nid_mth,"phis",itau_w,pphis)
    2121C
    2222      i = NINT(zout/zsto)
    23       CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
    24       CALL histwrite(nid_mth,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     23cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
     24      CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
    2525c
    2626      zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_ter)
    27       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    28       CALL histwrite(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
    29      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    30 c
    31       CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
    32       CALL histwrite(nid_mth,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     27cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 
     28      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
     29     $      zx_tmp_fi2d)
     30c
     31cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     32      CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
    3333
    34       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    35       CALL histwrite(nid_mth,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    36 c
    37       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
    38       CALL histwrite(nid_mth,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     34cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     35      CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
     36c
     37cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
     38      CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
    3939c
    4040c ENSEMBLES BEG
    41       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    42       CALL histwrite(nid_mth,"t2m_min",itau_w,zx_tmp_2d,
    43      .               iim*jjmp1,ndex2d)
    44 c
    45       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    46       CALL histwrite(nid_mth,"t2m_max",itau_w,zx_tmp_2d,
    47      .               iim*jjmp1,ndex2d)
     41cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     42      CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
     43c
     44cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     45      CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
    4846c
    4947c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
    5048c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
    5149c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
    52 c     CALL histwrite(nid_mth,"tsoil",itau_w,zx_tmp_2d,
     50c     CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
    5351c    .               iim*jjmp1,ndex2d)
    5452c
    … …  
    5654       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
    5755      ENDDO
    58       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    59       CALL histwrite(nid_mth,"wind10m",itau_w,zx_tmp_2d,
    60      $     iim*jjmp1,ndex2d)
     56cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     57      CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
    6158c
    6259      DO i = 1, klon
    6360         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
    6461      ENDDO
    65       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    66       CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     62cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     63      CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
    6764c
    6865c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    69 c     CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
     66c     CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
    7067c    .                                   iim*jjmp1*klev,ndex3d)
    7168c
    7269c ENSEMBLES END
    7370c
    74       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
    75       CALL histwrite(nid_mth,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    76 c
    77       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
    78       CALL histwrite(nid_mth,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    79 C
    80       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
    81       CALL histwrite(nid_mth,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     71cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
     72      CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
     73c
     74cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
     75      CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
     76C
     77cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
     78      CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
    8279C
    8380      DO i = 1, klon
    8481         zx_tmp_fi2d(i) = paprs(i,1)
    8582      ENDDO
    86       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    87       CALL histwrite(nid_mth,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    88 c
    89       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
    90       CALL histwrite(nid_mth,"qsurf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     83cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     84      CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
     85c
     86cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
     87      CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
    9188c
    9289      if (.not. ok_veget) then
    93           CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
    94           CALL histwrite(nid_mth,"qsol",itau_w,zx_tmp_2d,iim*jjmp1
    95      $        ,ndex2d)
     90cym          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
     91          CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
    9692      endif
    9793c
    … …  
    9995         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
    10096      ENDDO
    101       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    102       CALL histwrite(nid_mth,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    103 c
    104       CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
    105       CALL histwrite(nid_mth,"ndayrain",itau_w,zx_tmp_2d,
    106      $               iim*jjmp1,ndex2d)
     97cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     98      CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
     99c
     100cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
     101      CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
    107102c
    108103      DO i = 1, klon
    109104         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
    110105      ENDDO
    111       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    112       CALL histwrite(nid_mth,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     106cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     107      CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
    113108c
    114109      DO i = 1, klon
    115110         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
    116111      ENDDO
    117       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    118       CALL histwrite(nid_mth,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    119 c
    120       CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
    121       CALL histwrite(nid_mth,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    122 c
    123       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
    124       CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
    125      .               ndex2d)
     112cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     113      CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
     114c
     115cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
     116      CALL histwrite_phy(nid_mth,"snow",itau_w,snow_fall)
     117c
     118cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
     119      CALL histwrite_phy(nid_mth,"snow_mass",itau_w,zxsnow)
    126120c
    127121cIM: 071003
    128122c     zx_tmp_fi2d(1:klon)=evap(1:klon)*86400.
    129123c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    130       CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
    131       CALL histwrite(nid_mth,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    132 c
    133       CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
    134       CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    135 c
    136       CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
    137       CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    138 c
    139       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    140       CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    141 c
    142       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
    143       CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     124cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
     125      CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
     126c
     127cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
     128      CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
     129c
     130cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
     131      CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
     132c
     133cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
     134      CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
     135c
     136cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
     137      CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
    144138c
    145139c     zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
    146140cIM 080304   zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 2 )
    147141      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
    148       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    149       CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
    150      .                               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)
    151144c
    152145c
    … …  
    154147cIM 080304   zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 2 )
    155148      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
    156       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    157       CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
    158      .                               iim*jjmp1,ndex2d)
     149cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     150      CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
    159151c
    160152c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
    161153cIM 080304   zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 2 )
    162154      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
    163       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    164       CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
    165      .                               iim*jjmp1,ndex2d)
     155cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     156      CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
    166157c
    167158c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
    168159cIM 080304   zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 2 )
    169160      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
    170       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    171       CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
    172      .                               iim*jjmp1,ndex2d)
    173 c
    174       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
    175       CALL histwrite(nid_mth,"SWup200",itau_w,zx_tmp_2d,
    176      $     iim*jjmp1,ndex2d)
    177 c
    178       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
    179       CALL histwrite(nid_mth,"SWup200clr",itau_w,zx_tmp_2d,
    180      $     iim*jjmp1,ndex2d)
    181 c
    182       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
    183       CALL histwrite(nid_mth,"SWdn200",itau_w,zx_tmp_2d,
    184      $     iim*jjmp1,ndex2d)
    185 c
    186       CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
    187       CALL histwrite(nid_mth,"SWdn200clr",itau_w,zx_tmp_2d,
    188      $     iim*jjmp1,ndex2d)
    189 c
    190       CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
    191       CALL histwrite(nid_mth,"LWup200",itau_w,zx_tmp_2d,
    192      $     iim*jjmp1,ndex2d)
    193 c
    194       CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
    195       CALL histwrite(nid_mth,"LWup200clr",itau_w,zx_tmp_2d,
    196      $     iim*jjmp1,ndex2d)
     161cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     162      CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
     163c
     164cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
     165      CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
     166c
     167cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
     168      CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
     169c
     170cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
     171      CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
     172c
     173cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
     174      CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
     175c
     176cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
     177      CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
     178c
     179cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
     180      CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
    197181c
    198182      zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
    199       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    200       CALL histwrite(nid_mth,"LWdn200",itau_w,zx_tmp_2d,
    201      $     iim*jjmp1,ndex2d)
     183cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     184      CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
    202185c
    203186      zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
    204       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    205       CALL histwrite(nid_mth,"LWdn200clr",itau_w,zx_tmp_2d,
    206      $     iim*jjmp1,ndex2d)
     187cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     188      CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
    207189c
    208190c200 END
    209       CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
    210       CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    211 c
    212       CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
    213       CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    214 c
    215       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
    216       CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    217 c
    218       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
    219       CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     191cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
     192      CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
     193c
     194cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
     195      CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
     196c
     197cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
     198      CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
     199c
     200cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
     201      CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
    220202c
    221203c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
    222204      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
    223       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    224       CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
    225      .                               iim*jjmp1,ndex2d)
     205cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     206      CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
    226207c
    227208c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
    228209      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
    229       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    230       CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
    231      .                               iim*jjmp1,ndex2d)
     210cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     211      CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
    232212c
    233213c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
    234214      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
    235       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    236       CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
    237      .                               iim*jjmp1,ndex2d)
     215cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     216      CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
    238217c
    239218c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
    240219      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
    241       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    242       CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
    243      .                               iim*jjmp1,ndex2d)
     220cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     221      CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
    244222c
    245223      zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
    246       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    247       CALL histwrite(nid_mth,"LWupSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    248      .               ndex2d)
     224cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     225      CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
    249226c
    250227      zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
    251       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    252       CALL histwrite(nid_mth,"LWupSFCclr",itau_w,zx_tmp_2d,iim*jjmp1,
    253      .               ndex2d)
    254 c
    255       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    256       CALL histwrite(nid_mth,"LWdnSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    257      .               ndex2d)
     228cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     229      CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
     230c
     231cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
     232      CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
    258233c
    259234      DO i=1, klon
    … …  
    261236      ENDDO !i=1, klon
    262237c
    263       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
    264       CALL histwrite(nid_mth,"LWdnSFCclr",itau_w,zx_tmp_2d,
    265      $               iim*jjmp1,ndex2d)
     238cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
     239      CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
    266240c
    267241c
    … …  
    271245      zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
    272246c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topswad,zx_tmp_2d)
    273       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
    274       CALL histwrite(nid_mth,"topsad",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     247cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
     248      CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d)
    275249c
    276250      zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
    277251c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solswad,zx_tmp_2d)
    278       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
    279       CALL histwrite(nid_mth,"solsad",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,"solsad",itau_w,zx_tmp_fi2d)
    280254c
    281255      zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
    282256c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topswai,zx_tmp_2d)
    283       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
    284       CALL histwrite(nid_mth,"topsai",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     257cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
     258      CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d)
    285259c
    286260      zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
    287261c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solswai,zx_tmp_2d)
    288       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
    289       CALL histwrite(nid_mth,"solsai",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    290 c     endif
    291 c
    292       CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
    293       CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     262cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
     263      CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d)
     264cym
     265c      endif
     266c
     267cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
     268      CALL histwrite_phy(nid_mth,"bils",itau_w,bils)
    294269c
    295270      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
    296271c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
    297       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    298       CALL histwrite(nid_mth,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    299 c
    300       CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
    301       CALL histwrite(nid_mth,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     272cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     273      CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d)
     274c
     275cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
     276      CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
    302277c
    303278c
    … …  
    315290c
    316291c
    317        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
    318        CALL histwrite(nid_mth,"ffonte",itau_w,zx_tmp_2d,iim*jjmp1,
    319      $                ndex2d)
    320 c
    321        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
    322        CALL histwrite(nid_mth,"fqcalving",itau_w,zx_tmp_2d,iim*jjmp1,
    323      $                ndex2d)
     292cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
     293       CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
     294c
     295cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
     296       CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
    324297cIM: 171003
    325298       DO nsrf = 1, nbsrf
    326299        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
    327         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    328         CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
    329      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     300cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     301        CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
     302     $      zx_tmp_fi2d)
    330303C     
    331304        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
    332         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    333         CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
    334      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     305cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     306        CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
     307     $      zx_tmp_fi2d)
    335308       ENDDO
    336309cIM: 171003
    … …  
    341314       IF(nsrf.GT.1) THEN
    342315        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
    343         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    344         CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
    345      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     316cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     317        CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
     318     $      zx_tmp_fi2d)
    346319       ENDIF !nsrf.GT.1
    347320C
    348321        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
    349         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    350         CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
    351      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     322cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     323        CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
     324     $      zx_tmp_fi2d)
    352325C
    353326        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
    354         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    355         CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
    356      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     327cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     328        CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
     329     $      zx_tmp_fi2d)
    357330C
    358331        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
    359         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    360         CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
    361      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     332cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     333        CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
     334     $      zx_tmp_fi2d)
    362335C
    363336        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
    364         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    365         CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
    366      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     337cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     338        CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
     339     $      zx_tmp_fi2d)
    367340c
    368341        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
    369         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    370         CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
    371      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     342cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     343        CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
     344     $      zx_tmp_fi2d)
    372345C
    373346        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
    374         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    375         CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
    376      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     347cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     348        CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
     349     $      zx_tmp_fi2d)
    377350C
    378351
    … …  
    381354
    382355c
    383       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    384       CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    385 c
    386       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    387       CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    388 c
    389       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
    390       CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    391 c
    392       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
    393       CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    394 c
    395       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
    396       CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    397 c
    398       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    399       CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    400 c
    401       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
    402       CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     356cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
     357      CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
     358c
     359cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
     360      CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
     361c
     362cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
     363      CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
     364c
     365cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
     366      CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
     367c
     368cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
     369      CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
     370c
     371cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
     372      CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt)
     373c
     374cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
     375      CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
    403376c
    404377      zx_tmp_fi2d(1:klon) = flwp(1:klon)
    405       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    406       CALL histwrite(nid_mth,"lwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     378cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     379      CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
    407380c
    408381      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
    409       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    410       CALL histwrite(nid_mth,"iwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    411 c
    412       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
    413       CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    414 c
    415       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
    416       CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    417 c
    418       CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
    419       CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    420 c
    421       CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
    422       CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     382cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     383      CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
     384c
     385cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
     386      CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
     387c
     388cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
     389      CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
     390c
     391cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
     392      CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
     393c
     394cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
     395      CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
    423396cKE43
    424397      IF (iflag_con.GE.3) THEN ! sb
    425398c
    426       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
    427       CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    428 c
    429       CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
    430       CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    431 c
    432       CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
    433       CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    434 c
    435       CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
    436       CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     399cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
     400      CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
     401c
     402cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
     403      CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
     404c
     405cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
     406      CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
     407c
     408cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
     409      CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
    437410c
    438411c34EK
    … …  
    440413      ENDIF !iflag_con.GE.3
    441414c
    442       CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
    443       CALL histwrite(nid_mth,"prw",itau_w,zx_tmp_2d,
    444      .               iim*jjmp1,ndex2d)
     415cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
     416      CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
    445417c
    446418c
    … …  
    450422      DO k=1, nlevENS
    451423c
    452         CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
    453         CALL histwrite(nid_mth,"u"//clev(k),itau_w,zx_tmp_2d,
    454      $       iim*jjmp1,ndex2d)
    455 c
    456         CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
    457         CALL histwrite(nid_mth,"v"//clev(k),itau_w,zx_tmp_2d,
    458      $       iim*jjmp1,ndex2d)
    459 c
    460       CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlev(:,k),zx_tmp_2d)
    461       CALL histwrite(nid_mth,"w"//clev(k),itau_w,zx_tmp_2d,
    462      $     iim*jjmp1,ndex2d)
    463 c
    464       CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
    465       CALL histwrite(nid_mth,"phi"//clev(k),itau_w,zx_tmp_2d,
    466      $     iim*jjmp1,ndex2d)
     424cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
     425        CALL histwrite_phy(nid_mth,"u"//clev(k),itau_w, ulev(:,k))
     426c
     427cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
     428        CALL histwrite_phy(nid_mth,"v"//clev(k),itau_w,vlev(:,k))
     429c
     430cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlev(:,k),zx_tmp_2d)
     431      CALL histwrite_phy(nid_mth,"w"//clev(k),itau_w,wlev(:,k))
     432c
     433cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
     434      CALL histwrite_phy(nid_mth,"phi"//clev(k),itau_w,philev(:,k))
    467435c
    468436      ENDDO !nlevENS
    … …  
    474442C
    475443c
    476       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
    477       CALL histwrite(nid_mth,"lwcon",itau_w,zx_tmp_3d,
    478      .               iim*jjmp1*klev,ndex3d)
    479 c
    480       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
    481       CALL histwrite(nid_mth,"iwcon",itau_w,zx_tmp_3d,
    482      .               iim*jjmp1*klev,ndex3d)
    483 c
    484       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    485       CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
    486      .                                   iim*jjmp1*klev,ndex3d)
     444cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
     445      CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
     446c
     447cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
     448      CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
     449c
     450cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
     451      CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
    487452c
    488453c#ifdef histmthNMC
    … …  
    490455c#endif
    491456c
    492       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
    493       CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
    494      .                                   iim*jjmp1*klev,ndex3d)
     457cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
     458      CALL histwrite_phy(nid_mth,"ovap",itau_w,qx(:,:,ivap))
    495459c
    496460cIM: 071003
    … …  
    501465c    .                                   iim*jjmp1*klev,ndex3d)
    502466c
    503       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    504       CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
    505      .                                   iim*jjmp1*klev,ndex3d)
    506 c
    507       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    508       CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
    509      .                                   iim*jjmp1*klev,ndex3d)
    510 c
    511       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    512       CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
    513      .                                   iim*jjmp1*klev,ndex3d)
    514 c
    515       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    516       CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
    517      .                                   iim*jjmp1*klev,ndex3d)
    518 c
    519       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    520       CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
    521      .                                   iim*jjmp1*klev,ndex3d)
    522 c
    523       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
    524       CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
    525      .                                   iim*jjmp1*klev,ndex3d)
    526 c
    527       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
    528       CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
    529      .                                   iim*jjmp1*klev,ndex3d)
    530 c
    531       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
    532       CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
    533      .                                   iim*jjmp1*klev,ndex3d)
    534 c
    535       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
    536       CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
    537      .                                   iim*jjmp1*klev,ndex3d)
     467cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
     468      CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
     469c
     470cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
     471      CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
     472c
     473cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
     474      CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
     475c
     476cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
     477      CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
     478c
     479cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
     480      CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
     481c
     482cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
     483      CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra)
     484c
     485cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
     486      CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
     487c
     488cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
     489      CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
     490c
     491cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
     492      CALL histwrite_phy(nid_mth,"ozone",itau_w,wo)
    538493c
    539494cIM 22.03.04 BEG
    540       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
    541       CALL histwrite(nid_mth,"dtphy",itau_w,zx_tmp_3d,
    542      .                                   iim*jjmp1*klev,ndex3d)
    543 c
    544       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),
    545      .                 zx_tmp_3d)
    546       CALL histwrite(nid_mth,"dqphy",itau_w,zx_tmp_3d,
    547      .                                   iim*jjmp1*klev,ndex3d)
     495cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
     496      CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
     497c
     498cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),
     499cym     .                 zx_tmp_3d)
     500      CALL histwrite_phy(nid_mth,"dqphy",itau_w, d_qx(:,:,ivap))
    548501cIM 22.03.04 END
    549502c
    … …  
    555508c
    556509        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
    557         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    558         CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
    559      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     510cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     511        CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
     512     $      zx_tmp_fi2d)
    560513C
    561514        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
    562         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    563         CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
    564      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     515cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     516        CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
     517     $      zx_tmp_fi2d)
    565518c
    566519      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
    567       CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
    568       CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w
    569      $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
     520cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
     521      CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
     522     $    ,zx_tmp_fi2d)
    570523c
    571524      ENDDO !nsrf=1, nbsrf
    572525c
    573       CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
    574       CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    575 c
    576       CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
    577       CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     526cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
     527      CALL histwrite_phy(nid_mth,"albs",itau_w,albsol)
     528c
     529cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
     530      CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw)
    578531c
    579532      ENDIF !lev_histmth.GE.3
    … …  
    582535      IF(lev_histmth.GE.4) THEN
    583536c
    584       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
    585       CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
    586      .                                   iim*jjmp1*klev,ndex3d)
     537cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
     538      CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
    587539c
    588540      IF (iflag_con.GE.3) THEN
    589541c
    590       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
    591       CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
    592      .                                   iim*jjmp1*klev,ndex3d)
    593 c
    594       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
    595       CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
    596      .                                   iim*jjmp1*klev,ndex3d)
    597 c
    598       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
    599       CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
    600      .                                   iim*jjmp1*klev,ndex3d)
    601 c
    602       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
    603       CALL histwrite(nid_mth,"dnwd0",itau_w,zx_tmp_3d,
    604      .                                   iim*jjmp1*klev,ndex3d)
    605 c
    606       ENDIF !iflag_con.GE.3
    607 c
    608       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
    609       CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
    610      .                                   iim*jjmp1*klev,ndex3d)
    611 c
    612       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
    613       CALL histwrite(nid_mth,"dqdyn",itau_w,zx_tmp_3d,
    614      .                                   iim*jjmp1*klev,ndex3d)
     542cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
     543      CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
     544c
     545cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
     546      CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
     547c
     548cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
     549      CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
     550c
     551cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
     552      CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
     553c
     554
     555      ENDIF
     556c
     557cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
     558      CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
     559c
     560cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
     561      CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
    615562c
    616563cIM: 101003 : K/30min ==> K/s
    617564      zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
    618       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    619       CALL histwrite(nid_mth,"dtcon",itau_w,zx_tmp_3d,
    620      .                                   iim*jjmp1*klev,ndex3d)
    621 c
    622       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
    623       CALL histwrite(nid_mth,"dqcon",itau_w,zx_tmp_3d,
    624      .                                   iim*jjmp1*klev,ndex3d)
     565cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     566      CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
     567c
     568cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
     569      CALL histwrite_phy(nid_mth,"dqcon",itau_w,d_q_con)
    625570c
    626571cIM: 101003 : K/30min ==> K/s
    627572      zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
    628       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    629       CALL histwrite(nid_mth,"dtlsc",itau_w,zx_tmp_3d,
    630      .                                   iim*jjmp1*klev,ndex3d)
     573cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     574      CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
    631575cIM: 071003
    632576cIM: 101003 : K/30min ==> K/s
    633577      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
    634578     .                             d_t_eva(1:klon,1:klev))/pdtphys
    635       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
    636       CALL histwrite(nid_mth,"dtlschr",itau_w,zx_tmp_3d,
    637      .                                   iim*jjmp1*klev,ndex3d)
    638 c
    639       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
    640       CALL histwrite(nid_mth,"dqlsc",itau_w,zx_tmp_3d,
    641      .                                   iim*jjmp1*klev,ndex3d)
     579cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
     580      CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
     581c
     582cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
     583      CALL histwrite_phy(nid_mth,"dqlsc",itau_w,d_q_lsc)
    642584c
    643585cIM: 101003 : K/30min ==> K/s
    644586      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
    645       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    646       CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
    647      .                                   iim*jjmp1*klev,ndex3d)
    648 c
    649       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
    650       CALL histwrite(nid_mth,"dqvdf",itau_w,zx_tmp_3d,
    651      .                                   iim*jjmp1*klev,ndex3d)
     587cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     588      CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
     589c
     590cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
     591      CALL histwrite_phy(nid_mth,"dqvdf",itau_w,d_q_vdf)
    652592c
    653593cIM: 101003 : K/30min ==> K/s
    654594      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
    655       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    656       CALL histwrite(nid_mth,"dteva",itau_w,zx_tmp_3d,
    657      .                                   iim*jjmp1*klev,ndex3d)
    658 c
    659       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
    660       CALL histwrite(nid_mth,"dqeva",itau_w,zx_tmp_3d,
    661      .                                   iim*jjmp1*klev,ndex3d)
     595cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     596      CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
     597c
     598cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
     599      CALL histwrite_phy(nid_mth,"dqeva",itau_w, d_q_eva)
    662600c
    663601      zpt_conv = 0.
    664602      where (ptconv) zpt_conv = 1.
    665       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
    666       CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d,
    667      .                                   iim*(jjmp1)*klev,ndex3d)
    668 c
    669       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
    670       CALL histwrite(nid_mth,"ratqs",itau_w,zx_tmp_3d,
    671      .                                   iim*(jjmp1)*klev,ndex3d)
     603cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
     604      CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
     605c
     606cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
     607      CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
    672608c
    673609cIM: 101003 : K/30min ==> K/s
    674610      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
    675       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    676       CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
    677      .                                   iim*jjmp1*klev,ndex3d)
    678 c
    679       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
    680       CALL histwrite(nid_mth,"dqajs",itau_w,zx_tmp_3d,
    681      .                                   iim*jjmp1*klev,ndex3d)
     611cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     612      CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
     613c
     614cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
     615      CALL histwrite_phy(nid_mth,"dqajs",itau_w,d_q_ajs)
    682616c
    683617cIM: 101003 : K/day ==> K/s
    684618cIM: LMD_ARMIP3   zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)*pdtphys/RDAY
    685619      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
    686       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    687       CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
    688      .                                   iim*jjmp1*klev,ndex3d)
     620cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     621      CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
    689622c
    690623cIM: 101003 : K/day ==> K/s     
    691624cIM: LMD_ARMIP3   zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)*pdtphys/RDAY
    692625      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
    693       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    694       CALL histwrite(nid_mth,"dtsw0",itau_w,zx_tmp_3d,
    695      .                                   iim*jjmp1*klev,ndex3d)
     626cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     627      CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
    696628c
    697629cIM: 101003 : K/day ==> K/s     
    698630cIM: LMD_ARMIP3     zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)*pdtphys/RDAY
    699631      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
    700       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    701       CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
    702      .                                   iim*jjmp1*klev,ndex3d)
     632cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     633      CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
    703634c
    704635cIM: 101003 : K/day ==> K/s     
    705636cIM: LMD_ARMIP3     zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)*pdtphys/RDAY
    706637      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
    707       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    708       CALL histwrite(nid_mth,"dtlw0",itau_w,zx_tmp_3d,
    709      .                                   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,"dtlw0",itau_w,zx_tmp_fi3d)
    710640c
    711641cIM: 101003 : deja en K/s     
    712642      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
    713       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    714       CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
    715      .                                   iim*jjmp1*klev,ndex3d)
    716 c
    717       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
    718       CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
    719      .                                   iim*jjmp1*klev,ndex3d)
    720 c
    721       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
    722       CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
    723      .                                   iim*jjmp1*klev,ndex3d)
     643cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     644      CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
     645c
     646cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
     647      CALL histwrite_phy(nid_mth,"duvdf",itau_w,d_u_vdf)
     648c
     649cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
     650      CALL histwrite_phy(nid_mth,"dvvdf",itau_w,d_v_vdf)
    724651c
    725652      IF (ok_orodr) THEN
    726       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d)
    727       CALL histwrite(nid_mth,"duoro",itau_w,zx_tmp_3d,
    728      .                                   iim*jjmp1*klev,ndex3d)
    729 c
    730       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d)
    731       CALL histwrite(nid_mth,"dvoro",itau_w,zx_tmp_3d,
    732      .                                   iim*jjmp1*klev,ndex3d)
     653cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d)
     654      CALL histwrite_phy(nid_mth,"duoro",itau_w,d_u_oro)
     655c
     656cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d)
     657      CALL histwrite_phy(nid_mth,"dvoro",itau_w,d_v_oro)
    733658c
    734659      ENDIF
    735660C
    736661      IF (ok_orolf) THEN
    737       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d)
    738       CALL histwrite(nid_mth,"dulif",itau_w,zx_tmp_3d,
    739      .                                   iim*jjmp1*klev,ndex3d)
    740 c
    741       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d)
    742       CALL histwrite(nid_mth,"dvlif",itau_w,zx_tmp_3d,
    743      .                                   iim*jjmp1*klev,ndex3d)
     662cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d)
     663      CALL histwrite_phy(nid_mth,"dulif",itau_w,d_u_lif)
     664c
     665cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d)
     666      CALL histwrite_phy(nid_mth,"dvlif",itau_w, d_v_lif)
    744667      ENDIF
    745668C
    … …  
    748671c
    749672c     IF (ok_ade.OR.ok_aie) THEN
    750       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, re, zx_tmp_3d)
    751       CALL histwrite(nid_mth,"re",itau_w,zx_tmp_3d,
    752      .                                   iim*jjmp1*klev,ndex3d)
    753 c
    754       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fl, zx_tmp_3d)
    755       CALL histwrite(nid_mth,"redenom",itau_w,zx_tmp_3d,
    756      .                                   iim*jjmp1*klev,ndex3d)
    757 c
    758       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtau, zx_tmp_3d)
    759       CALL histwrite(nid_mth,"tau",itau_w,zx_tmp_3d,
    760      .                                   iim*jjmp1*klev,ndex3d)
    761 c
    762       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtaupi, zx_tmp_3d)
    763       CALL histwrite(nid_mth,"taupi",itau_w,zx_tmp_3d,
    764      .                                   iim*jjmp1*klev,ndex3d)
     673cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, re, zx_tmp_3d)
     674      CALL histwrite_phy(nid_mth,"re",itau_w,re)
     675c
     676cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fl, zx_tmp_3d)
     677      CALL histwrite_phy(nid_mth,"redenom",itau_w,fl)
     678c
     679cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtau, zx_tmp_3d)
     680      CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau)
     681c
     682cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtaupi, zx_tmp_3d)
     683      CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi)
    765684c     endif
    766685c
    767686      IF (nqmax.GE.3) THEN
    768687      DO iq=3,nqmax
    769          CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
    770          CALL histwrite(nid_mth,tnom(iq),itau_w,zx_tmp_3d,
    771      .                                   iim*jjmp1*klev,ndex3d)
     688c         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
     689         CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
    772690      ENDDO
    773691      ENDIF
Note: See TracChangeset for help on using the changeset viewer.