Ignore:
Timestamp:
Mar 24, 2004, 9:43:54 AM (20 years ago)
Author:
lmdzadmin
Message:

IM: ajout lev_hist(hf,day,mth) niveaux output definis dans physiq.def et

variables ENSEMBLES + rearrangement

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histmth.h

    r486 r506  
    33      ndex2d = 0
    44      ndex3d = 0
    5 c
    6 c Champs 2D:
    7 c
     5
    86         zsto = dtime
    97         zout = dtime * ecrit_mth
    … …  
    119
    1210      i = NINT(zout/zsto)
     11c
     12      IF(lev_histmth.GE.1) THEN
     13c
     14c Champs 2D:
     15c
    1316      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
    1417      CALL histwrite(nid_mth,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    … …  
    1720      CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
    1821      CALL histwrite(nid_mth,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    19 
     22c
     23      zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_ter)
     24      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     25      CALL histwrite(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
     26     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     27c
    2028      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
    2129      CALL histwrite(nid_mth,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    … …  
    2735      CALL histwrite(nid_mth,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    2836c
     37c ENSEMBLES BEG
     38      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     39      CALL histwrite(nid_mth,"t2m_min",itau_w,zx_tmp_2d,
     40     .               iim*jjmp1,ndex2d)
     41c
     42      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     43      CALL histwrite(nid_mth,"t2m_max",itau_w,zx_tmp_2d,
     44     .               iim*jjmp1,ndex2d)
     45c
     46c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
     47c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
     48c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
     49c     CALL histwrite(nid_mth,"tsoil",itau_w,zx_tmp_2d,
     50c    .               iim*jjmp1,ndex2d)
     51c
     52      DO i=1, klon
     53       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
     54      ENDDO
     55      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     56      CALL histwrite(nid_mth,"wind10m",itau_w,zx_tmp_2d,
     57     $     iim*jjmp1,ndex2d)
     58c
     59      DO i = 1, klon
     60         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
     61      ENDDO
     62      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     63      CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     64c
     65c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
     66c     CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
     67c    .                                   iim*jjmp1*klev,ndex3d)
     68c
     69c ENSEMBLES END
     70c
    2971      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
    3072      CALL histwrite(nid_mth,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    … …  
    5799      CALL histwrite(nid_mth,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    58100c
     101      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
     102      CALL histwrite(nid_mth,"ndayrain",itau_w,zx_tmp_2d,
     103     $               iim*jjmp1,ndex2d)
     104c
    59105      DO i = 1, klon
    60106         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
    … …  
    77123c
    78124cIM: 071003
    79       zx_tmp_fi2d(1:klon)=evap(1:klon)*86400.
    80       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     125c     zx_tmp_fi2d(1:klon)=evap(1:klon)*86400.
     126c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     127      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
    81128      CALL histwrite(nid_mth,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    82129c
    … …  
    84131      CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    85132c
     133      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
     134      CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     135c
    86136      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    87137      CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    88138c
     139      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
     140      CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     141c
     142c     zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
     143cIM 080304   zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 2 )
     144      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
     145      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     146      CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
     147     .                               iim*jjmp1,ndex2d)
     148c
     149c
     150c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
     151cIM 080304   zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 2 )
     152      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
     153      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     154      CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
     155     .                               iim*jjmp1,ndex2d)
     156c
     157c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
     158cIM 080304   zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 2 )
     159      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
     160      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     161      CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
     162     .                               iim*jjmp1,ndex2d)
     163c
     164c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
     165cIM 080304   zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 2 )
     166      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
     167      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     168      CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
     169     .                               iim*jjmp1,ndex2d)
     170c
     171      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
     172      CALL histwrite(nid_mth,"SWup200",itau_w,zx_tmp_2d,
     173     $     iim*jjmp1,ndex2d)
     174c
     175      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
     176      CALL histwrite(nid_mth,"SWup200clr",itau_w,zx_tmp_2d,
     177     $     iim*jjmp1,ndex2d)
     178c
     179      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
     180      CALL histwrite(nid_mth,"SWdn200",itau_w,zx_tmp_2d,
     181     $     iim*jjmp1,ndex2d)
     182c
     183      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
     184      CALL histwrite(nid_mth,"SWdn200clr",itau_w,zx_tmp_2d,
     185     $     iim*jjmp1,ndex2d)
     186c
     187      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
     188      CALL histwrite(nid_mth,"LWup200",itau_w,zx_tmp_2d,
     189     $     iim*jjmp1,ndex2d)
     190c
     191      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
     192      CALL histwrite(nid_mth,"LWup200clr",itau_w,zx_tmp_2d,
     193     $     iim*jjmp1,ndex2d)
     194c
     195      zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
     196      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     197      CALL histwrite(nid_mth,"LWdn200",itau_w,zx_tmp_2d,
     198     $     iim*jjmp1,ndex2d)
     199c
     200      zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
     201      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     202      CALL histwrite(nid_mth,"LWdn200clr",itau_w,zx_tmp_2d,
     203     $     iim*jjmp1,ndex2d)
     204c
     205c200 END
    89206      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
    90207      CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    91208c
     209      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
     210      CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     211c
    92212      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
    93213      CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    94214c
    95       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    96       CALL histwrite(nid_mth,"LWdnSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    97      .               ndex2d)
    98 cIM: 071003
     215      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
     216      CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     217c
     218c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
     219      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
     220      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     221      CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
     222     .                               iim*jjmp1,ndex2d)
     223c
     224c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
     225      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
     226      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     227      CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
     228     .                               iim*jjmp1,ndex2d)
     229c
     230c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
     231      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
     232      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     233      CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
     234     .                               iim*jjmp1,ndex2d)
     235c
     236c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
     237      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
     238      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     239      CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
     240     .                               iim*jjmp1,ndex2d)
     241c
    99242      zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
    100243      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    101244      CALL histwrite(nid_mth,"LWupSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    102245     .               ndex2d)
    103 cLWupSFC
    104       CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
    105       CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    106 c
    107       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
    108       CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    109 c
    110       CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
    111       CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    112 c
    113       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
    114       CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     246c
     247      zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
     248      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     249      CALL histwrite(nid_mth,"LWupSFCclr",itau_w,zx_tmp_2d,iim*jjmp1,
     250     .               ndex2d)
     251c
     252      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
     253      CALL histwrite(nid_mth,"LWdnSFC",itau_w,zx_tmp_2d,iim*jjmp1,
     254     .               ndex2d)
     255c
     256      DO i=1, klon
     257       sollwdownclr(i)=lwdn0(i,1)
     258      ENDDO !i=1, klon
     259c
     260      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
     261      CALL histwrite(nid_mth,"LWdnSFCclr",itau_w,zx_tmp_2d,
     262     $               iim*jjmp1,ndex2d)
    115263c
    116264      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
    117265      CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    118266c
    119       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
     267      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
     268c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
     269      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    120270      CALL histwrite(nid_mth,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    121271c
    … …  
    144294       CALL histwrite(nid_mth,"fqcalving",itau_w,zx_tmp_2d,iim*jjmp1,
    145295     $                ndex2d)
    146 c
    147 
    148       DO nsrf = 1, nbsrf
    149 CYYY
    150         zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
    151         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    152         CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
    153      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    154 C
    155         zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
    156         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    157         CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
    158      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    159 C
    160         zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
    161         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    162         CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
    163      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    164 C
    165         zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
    166         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    167         CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
    168      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    169 C
    170         zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
    171         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    172         CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
    173      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    174 c
    175         zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
    176         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    177         CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
    178      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    179 C
    180         zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
    181         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    182         CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
    183      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    184 C
     296cIM: 171003
     297       DO nsrf = 1, nbsrf
    185298        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
    186299        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    … …  
    192305        CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
    193306     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     307       ENDDO
     308cIM: 171003
     309c
     310cIM      if ( lev_histday.gt.1 ) then
     311      DO nsrf = 1, nbsrf
     312CYYY
     313       IF(nsrf.GT.1) THEN
     314        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
     315        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     316        CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
     317     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     318       ENDIF !nsrf.GT.1
    194319C
     320        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
     321        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     322        CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
     323     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     324C
     325        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
     326        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     327        CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
     328     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     329C
     330        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
     331        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     332        CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
     333     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     334C
     335        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
     336        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     337        CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
     338     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     339c
     340        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
     341        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     342        CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
     343     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     344C
     345        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
     346        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     347        CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
     348     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     349C
     350
     351      END DO 
     352cIM      endif !lev_histday
     353
     354c
     355      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
     356      CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     357c
     358      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
     359      CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     360c
     361      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
     362      CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     363c
     364      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
     365      CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     366c
     367      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
     368      CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     369c
     370      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
     371      CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     372c
     373      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
     374      CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     375c
     376      zx_tmp_fi2d(1:klon) = flwp(1:klon)
     377      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     378      CALL histwrite(nid_mth,"lwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     379c
     380      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
     381      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     382      CALL histwrite(nid_mth,"iwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     383c
     384      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
     385      CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     386c
     387      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
     388      CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     389c
     390      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
     391      CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     392c
     393      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
     394      CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     395cKE43
     396      IF (iflag_con.GE.3) THEN ! sb
     397c
     398      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
     399      CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     400c
     401      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
     402      CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     403c
     404      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
     405      CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     406c
     407      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
     408      CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     409c
     410c34EK
     411c
     412      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
     413      CALL histwrite(nid_mth,"prw",itau_w,zx_tmp_2d,
     414     .               iim*jjmp1,ndex2d)
     415c
     416      ENDIF !iflag_con.GE.3
     417c
     418c
     419c
     420c  Champs interpolles sur des niveaux de pression
     421      DO k=1, nlevENS
     422c
     423        CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
     424        CALL histwrite(nid_mth,"u"//clev(k),itau_w,zx_tmp_2d,
     425     $       iim*jjmp1,ndex2d)
     426c
     427        CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
     428        CALL histwrite(nid_mth,"v"//clev(k),itau_w,zx_tmp_2d,
     429     $       iim*jjmp1,ndex2d)
     430c
     431      CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlev(:,k),zx_tmp_2d)
     432      CALL histwrite(nid_mth,"w"//clev(k),itau_w,zx_tmp_2d,
     433     $     iim*jjmp1,ndex2d)
     434c
     435      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
     436      CALL histwrite(nid_mth,"phi"//clev(k),itau_w,zx_tmp_2d,
     437     $     iim*jjmp1,ndex2d)
     438c
     439      ENDDO !nlevENS
     440      ENDIF !lev_histmth.GE.1
     441c
     442      IF(lev_histmth.GE.2) THEN
     443c
     444c Champs 3D:
     445C
     446c
     447      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
     448      CALL histwrite(nid_mth,"lwcon",itau_w,zx_tmp_3d,
     449     .               iim*jjmp1*klev,ndex3d)
     450c
     451      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
     452      CALL histwrite(nid_mth,"iwcon",itau_w,zx_tmp_3d,
     453     .               iim*jjmp1*klev,ndex3d)
     454c
     455      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
     456      CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
     457     .                                   iim*jjmp1*klev,ndex3d)
     458c
     459c#ifdef histmthNMC
     460c#include "write_histmthNMC.h"
     461c#endif
     462c
     463      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
     464      CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
     465     .                                   iim*jjmp1*klev,ndex3d)
     466c
     467cIM: 071003
     468c     zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
     469c    .                         (1-qx(1:klon,1:klev,ivap))
     470c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
     471c     CALL histwrite(nid_mth,"wvap",itau_w,zx_tmp_3d,
     472c    .                                   iim*jjmp1*klev,ndex3d)
     473c
     474      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
     475      CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
     476     .                                   iim*jjmp1*klev,ndex3d)
     477c
     478      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
     479      CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
     480     .                                   iim*jjmp1*klev,ndex3d)
     481c
     482      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
     483      CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
     484     .                                   iim*jjmp1*klev,ndex3d)
     485c
     486      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
     487      CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
     488     .                                   iim*jjmp1*klev,ndex3d)
     489c
     490      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
     491      CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
     492     .                                   iim*jjmp1*klev,ndex3d)
     493c
     494      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
     495      CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
     496     .                                   iim*jjmp1*klev,ndex3d)
     497c
     498      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
     499      CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
     500     .                                   iim*jjmp1*klev,ndex3d)
     501c
     502      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
     503      CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
     504     .                                   iim*jjmp1*klev,ndex3d)
     505c
     506      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
     507      CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
     508     .                                   iim*jjmp1*klev,ndex3d)
     509c
     510cIM 22.03.04 BEG
     511      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
     512      CALL histwrite(nid_mth,"dtphy",itau_w,zx_tmp_3d,
     513     .                                   iim*jjmp1*klev,ndex3d)
     514c
     515      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),
     516     .                 zx_tmp_3d)
     517      CALL histwrite(nid_mth,"dqphy",itau_w,zx_tmp_3d,
     518     .                                   iim*jjmp1*klev,ndex3d)
     519cIM 22.03.04 END
     520c
     521      ENDIF !lev_histmth.GE.2
     522c
     523      IF(lev_histmth.GE.3) THEN
     524c
     525       DO nsrf=1, nbsrf
     526c
    195527        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
    196528        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    … …  
    207539      CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w
    208540     $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
    209 
    210       END DO 
    211 cXXX      DO i = 1, klon
    212 cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
    213 cXXX      ENDDO
    214 cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    215 cXXX      CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     541c
     542      ENDDO !nsrf=1, nbsrf
    216543c
    217544      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
    218545      CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     546c
    219547      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
    220548      CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    221549c
    222       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    223       CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    224 c
    225       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    226       CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    227 c
    228       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
    229       CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    230 c
    231       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
    232       CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    233 c
    234       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
    235       CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    236 c
    237       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    238       CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    239 c
    240       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
    241       CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    242 c
    243       zx_tmp_fi2d(1:klon) = flwp(1:klon)
    244       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    245       CALL histwrite(nid_mth,"lwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    246 c
    247       zx_tmp_fi2d(1:klon) = fiwp(1:klon)
    248       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    249       CALL histwrite(nid_mth,"iwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    250 c
    251       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
    252       CALL histwrite(nid_mth,"lwcon",itau_w,zx_tmp_3d,
    253      .               iim*jjmp1*klev,ndex3d)
    254 c
    255       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
    256       CALL histwrite(nid_mth,"iwcon",itau_w,zx_tmp_3d,
    257      .               iim*jjmp1*klev,ndex3d)
    258 c
    259       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
    260       CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    261 c
    262       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
    263       CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    264 c
    265       CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
    266       CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    267 c
    268       CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
    269       CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    270 cKE43
    271       IF (iflag_con .GE. 3) THEN ! sb
    272 c
    273       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
    274       CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    275 c
    276       CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
    277       CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    278 c
    279       CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
    280       CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    281 c
    282       CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
    283       CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     550      ENDIF !lev_histmth.GE.3
     551c
     552c
     553      IF(lev_histmth.GE.4) THEN
     554c
     555      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
     556      CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
     557     .                                   iim*jjmp1*klev,ndex3d)
     558c
     559      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
     560      CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
     561     .                                   iim*jjmp1*klev,ndex3d)
    284562c
    285563      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
    … …  
    295573     .                                   iim*jjmp1*klev,ndex3d)
    296574c
    297       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
    298       CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
    299      .                                   iim*jjmp1*klev,ndex3d)
    300 c
    301       CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
    302       CALL histwrite(nid_mth,"prw",itau_w,zx_tmp_2d,
    303      .               iim*jjmp1,ndex2d)
    304 c
    305       ENDIF
    306 c34EK
    307 c
    308 c Champs 3D:
    309 C
    310       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    311       CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
    312      .                                   iim*jjmp1*klev,ndex3d)
    313 c
    314       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
    315       CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
    316      .                                   iim*jjmp1*klev,ndex3d)
    317 cIM: 071003
    318       zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
    319      .                         (1-qx(1:klon,1:klev,ivap))
    320       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
    321       CALL histwrite(nid_mth,"wvap",itau_w,zx_tmp_3d,
    322      .                                   iim*jjmp1*klev,ndex3d)
    323 c
    324       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    325       CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
    326      .                                   iim*jjmp1*klev,ndex3d)
    327 c
    328       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    329       CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
    330      .                                   iim*jjmp1*klev,ndex3d)
    331 c
    332       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    333       CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
    334      .                                   iim*jjmp1*klev,ndex3d)
    335 c
    336       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    337       CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
    338      .                                   iim*jjmp1*klev,ndex3d)
    339 c
    340       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    341       CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
    342      .                                   iim*jjmp1*klev,ndex3d)
    343 c
    344       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
    345       CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
    346      .                                   iim*jjmp1*klev,ndex3d)
    347 c
    348       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
    349       CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
    350      .                                   iim*jjmp1*klev,ndex3d)
    351 c
    352       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
    353       CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
    354      .                                   iim*jjmp1*klev,ndex3d)
    355 c
    356       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
    357       CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
    358      .                                   iim*jjmp1*klev,ndex3d)
    359 c
    360       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)
    361       CALL histwrite(nid_mth,"oliq",itau_w,zx_tmp_3d,
    362      .                                   iim*jjmp1*klev,ndex3d)
    363 c
    364       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
    365       CALL histwrite(nid_mth,"dtphy",itau_w,zx_tmp_3d,
    366      .                                   iim*jjmp1*klev,ndex3d)
    367575c
    368576      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
    … …  
    504712      ENDIF
    505713C
    506       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
    507       CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
    508      .                                   iim*jjmp1*klev,ndex3d)
    509 c
    510714      IF (nqmax.GE.3) THEN
    511715      DO iq=1,nqmax-2
    … …  
    522726      ENDIF
    523727c
    524 c     zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
    525       zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 2 )
    526       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    527       CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
    528      .                               iim*jjmp1,ndex2d)
    529 c
    530 c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
    531       zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
    532       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    533       CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
    534      .                               iim*jjmp1,ndex2d)
    535 c
    536 c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
    537       zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 2 )
    538       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    539       CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
    540      .                               iim*jjmp1,ndex2d)
    541 c
    542 c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
    543       zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
    544       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    545       CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
    546      .                               iim*jjmp1,ndex2d)
    547 c
    548 c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
    549       zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 2 )
    550       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    551       CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
    552      .                               iim*jjmp1,ndex2d)
    553 c
    554 c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
    555       zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
    556       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    557       CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
    558      .                               iim*jjmp1,ndex2d)
    559 c
    560 c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
    561       zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 2 )
    562       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    563       CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
    564      .                               iim*jjmp1,ndex2d)
    565 c
    566 c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
    567       zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
    568       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    569       CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
    570      .                               iim*jjmp1,ndex2d)
    571 c
    572 c  Champs interpolles sur des niveaux de pression
    573 
    574       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
    575       CALL histwrite(nid_mth,"u850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    576 
    577       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
    578       CALL histwrite(nid_mth,"v850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    579 
    580       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
    581       CALL histwrite(nid_mth,"u500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    582 
    583       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
    584       CALL histwrite(nid_mth,"v500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    585 
    586       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
    587       CALL histwrite(nid_mth,"u200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    588 
    589       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
    590       CALL histwrite(nid_mth,"v200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    591 
    592       CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
    593       CALL histwrite(nid_mth,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    594 
    595       CALL gr_fi_ecrit(1, klon,iim,jjmp1, w500,zx_tmp_2d)
    596       CALL histwrite(nid_mth,"w500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     728      ENDIF !lev_histmth.GE.4
    597729c
    598730      if (ok_sync) then
Note: See TracChangeset for help on using the changeset viewer.