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_histday.h

    r633 r634  
    1414
    1515         i = NINT(zout/zsto)
    16          CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
     16cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
    1717c
    1818      IF(lev_histday.GE.1) THEN
    1919c
    20        CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     20       CALL histwrite_phy(nid_day,"phis",itau_w,pphis)
    2121         varname = 'phis'
    2222         vartitle= 'Surface geop. height'
    … …  
    2424c
    2525         i = NINT(zout/zsto)
    26          CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
    27        CALL histwrite(nid_day,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     26cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
     27       CALL histwrite_phy(nid_day,"aire",itau_w,airephy)
    2828         varname = 'aire'
    2929         vartitle= 'Grid area'
    … …  
    3434      ENDDO
    3535c
    36       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    37       CALL histwrite(nid_day,"contfracATM",itau_w,
    38      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    39 c
    40       CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d)
    41       CALL histwrite(nid_day,"contfracOR",itau_w,
    42      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    43 c
    44       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    45       CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    46 C
    47       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    48       CALL histwrite(nid_day,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     36cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     37      CALL histwrite_phy(nid_day,"contfracATM",itau_w,zx_tmp_fi2d)
     38c
     39cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d)
     40      CALL histwrite_phy(nid_day,"contfracOR",itau_w,
     41     &                   pctsrf_new(:,is_ter))
     42c
     43cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     44      CALL histwrite_phy(nid_day,"tsol",itau_w,zxtsol)
     45C
     46cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     47      CALL histwrite_phy(nid_day,"t2m",itau_w,zt2m)
    4948c
    5049c   En attendant un eventuel debugage.
    5150
    52       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    53       CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
    54      .               iim*jjmp1,ndex2d)
    55 c
    56       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    57       CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
    58      .               iim*jjmp1,ndex2d)
     51cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     52      CALL histwrite_phy(nid_day,"t2m_min",itau_w,zt2m)
     53c
     54cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     55      CALL histwrite_phy(nid_day,"t2m_max",itau_w,zt2m)
    5956c
    6057      DO i = 1, klon
    6158         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
    6259      ENDDO
    63       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    64       CALL histwrite(nid_day,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     60cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     61      CALL histwrite_phy(nid_day,"plul",itau_w,zx_tmp_fi2d)
    6562c
    6663      DO i = 1, klon
    6764         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
    6865      ENDDO
    69       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    70       CALL histwrite(nid_day,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    71 c
    72       CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_lsc,zx_tmp_2d)
    73       CALL histwrite(nid_day,"snowl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     66cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     67      CALL histwrite_phy(nid_day,"pluc",itau_w,zx_tmp_fi2d)
     68c
     69cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_lsc,zx_tmp_2d)
     70      CALL histwrite_phy(nid_day,"snowl",itau_w,snow_lsc)
    7471c
    7572      IF(1.EQ.0) THEN
    76       CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_con,zx_tmp_2d)
    77       CALL histwrite(nid_day,"snowc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     73cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_con,zx_tmp_2d)
     74      CALL histwrite_phy(nid_day,"snowc",itau_w,snow_con)
    7875      ENDIF
    7976c
    80       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
    81       CALL histwrite(nid_day,"flat",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     77cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
     78      CALL histwrite_phy(nid_day,"flat",itau_w,zxfluxlat)
    8279c
    8380      DO i = 1, klon
    … …  
    8582      ENDDO
    8683c
    87       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    88       CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    89 c
    90       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zq2m,zx_tmp_2d)
    91       CALL histwrite(nid_day,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    92 c
    93       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zu10m,zx_tmp_2d)
    94       CALL histwrite(nid_day,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,
    95      .               ndex2d)
    96 c
    97       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zv10m,zx_tmp_2d)
    98       CALL histwrite(nid_day,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,
    99      .               ndex2d)
     84cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     85      CALL histwrite_phy(nid_day,"sicf",itau_w,zx_tmp_fi2d)
     86c
     87cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zq2m,zx_tmp_2d)
     88      CALL histwrite_phy(nid_day,"q2m",itau_w,zq2m)
     89c
     90cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zu10m,zx_tmp_2d)
     91      CALL histwrite_phy(nid_day,"u10m",itau_w,zu10m)
     92c
     93cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zv10m,zx_tmp_2d)
     94      CALL histwrite_phy(nid_day,"v10m",itau_w,zv10m)
    10095c
    10196      DO i=1, klon
    10297       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
    10398      ENDDO
    104       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    105       CALL histwrite(nid_day,"wind10m",itau_w,zx_tmp_2d,
    106      $     iim*jjmp1,ndex2d)
    107 c
    108       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    109       CALL histwrite(nid_day,"wind10max",itau_w,zx_tmp_2d,
    110      .               iim*jjmp1,ndex2d)
     99cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)     
     100      CALL histwrite_phy(nid_day,"wind10m",itau_w,zx_tmp_fi2d)
     101c
     102cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     103      CALL histwrite_phy(nid_day,"wind10max",itau_w,zx_tmp_fi2d)
    111104c
    112105      DO i=1, klon
    113106         zx_tmp_fi2d(i) = paprs(i,1)
    114107      ENDDO
    115       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    116       CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     108cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     109      CALL histwrite_phy(nid_day,"psol",itau_w,zx_tmp_fi2d)
    117110c
    118111      DO i=1, klon
    119112         zx_tmp_fi2d(i) = (rain_fall(i) + snow_fall(i))
    120113      ENDDO
    121       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    122       CALL histwrite(nid_day,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    123 c
    124       CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
    125       CALL histwrite(nid_day,"snowf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    126 c
    127       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
    128       CALL histwrite(nid_day,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
    129      .               ndex2d)
    130 c
    131       CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
    132       CALL histwrite(nid_day,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    133 c
    134       CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
    135       CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    136 c
    137       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    138       CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    139 c
    140       CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
    141       CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    142 c
    143       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
    144       CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    145 c
    146       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
    147       CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    148 c
    149       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
    150       CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    151 c
    152       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
    153       CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    154 c
    155       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    156       CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    157 c
    158       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
    159       CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    160 c
    161       CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
    162       CALL histwrite(nid_day,"prw",itau_w,zx_tmp_2d,
    163      .               iim*jjmp1,ndex2d)
     114cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     115      CALL histwrite_phy(nid_day,"precip",itau_w,zx_tmp_fi2d)
     116c
     117cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
     118      CALL histwrite_phy(nid_day,"snowf",itau_w,snow_fall)
     119c
     120cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
     121      CALL histwrite_phy(nid_day,"snow_mass",itau_w,zxsnow)
     122c
     123cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
     124      CALL histwrite_phy(nid_day,"evap",itau_w,evap)
     125c
     126cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
     127      CALL histwrite_phy(nid_day,"tops",itau_w,topsw)
     128c
     129cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
     130      CALL histwrite_phy(nid_day,"topl",itau_w,toplw)
     131c
     132cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
     133      CALL histwrite_phy(nid_day,"sols",itau_w,solsw)
     134c
     135cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
     136      CALL histwrite_phy(nid_day,"soll",itau_w,sollw)
     137c
     138cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
     139      CALL histwrite_phy(nid_day,"cldl",itau_w,cldl)
     140c
     141cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
     142      CALL histwrite_phy(nid_day,"cldm",itau_w,cldm)
     143c
     144cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
     145      CALL histwrite_phy(nid_day,"cldh",itau_w,cldh)
     146c
     147cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
     148      CALL histwrite_phy(nid_day,"cldt",itau_w,cldt)
     149c
     150cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
     151      CALL histwrite_phy(nid_day,"cldq",itau_w,cldq)
     152c
     153cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
     154      CALL histwrite_phy(nid_day,"prw",itau_w,prw)
    164155
    165156c   Ecriture de champs dynamiques sur des niveaux de pression
    166157      DO k=1, nlevENS
    167158c
    168         CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
    169         CALL histwrite(nid_day,"u"//clev(k),itau_w,zx_tmp_2d,
    170      $       iim*jjmp1,ndex2d)
    171 c
    172         CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
    173         CALL histwrite(nid_day,"v"//clev(k),itau_w,zx_tmp_2d,
    174      $       iim*jjmp1,ndex2d)
     159cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
     160        CALL histwrite_phy(nid_day,"u"//clev(k),itau_w,ulev(:,k))
     161c
     162cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
     163        CALL histwrite_phy(nid_day,"v"//clev(k),itau_w,vlev(:,k))
    175164c
    176165      ENDDO !nlevENS
    … …  
    179168      IF(clev(k).EQ."500") THEN
    180169c
    181       CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlev(:,indENS(k)),zx_tmp_2d)
    182       CALL histwrite(nid_day,"w500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     170cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlev(:,indENS(k)),zx_tmp_2d)
     171cym      CALL histwrite_phy(nid_day,"w500",itau_w,wlev(:,indENS(k)))
     172         CALL histwrite_phy(nid_day,"w500",itau_w,wlev(:,k))
     173
    183174c
    184175c     DO i=1, klon
    … …  
    188179c
    189180cBAD  CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,indENS(k)),zx_tmp_2d)
    190       CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
    191       CALL histwrite(nid_day,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     181cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
     182      CALL histwrite_phy(nid_day,"phi500",itau_w,philev(:,k))
    192183c
    193184      ENDIF !clev(k).EQ.500
    194185      ENDDO !k=1, nlevENS
    195186c
    196       CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
    197       CALL histwrite(nid_day,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    198 
    199       CALL gr_fi_ecrit(1,klon,iim,jjmp1,cape,zx_tmp_2d)
    200       CALL histwrite(nid_day,"cape_max",itau_w,zx_tmp_2d,
    201      .               iim*jjmp1,ndex2d)
    202 c
    203       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    204       CALL histwrite(nid_day,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
    205      .               ndex2d)
     187cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     188      CALL histwrite_phy(nid_day,"slp",itau_w,slp)
     189
     190cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,cape,zx_tmp_2d)
     191      CALL histwrite_phy(nid_day,"cape_max",itau_w,cape)
     192c
     193cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
     194      CALL histwrite_phy(nid_day,"solldown",itau_w,sollwdown)
    206195c
    207196      DO i=1, klon
    208197        zx_tmp_fi2d(i)=-1*sens(i)
    209198      ENDDO
    210       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    211       CALL histwrite(nid_day,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     199cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     200      CALL histwrite_phy(nid_day,"sens",itau_w,zx_tmp_fi2d)
    212201c
    213202      zx_tmp_fi2d(1 : klon) = swdn(1 : klon, 1)
    214203
    215       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    216       CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
    217      .                               iim*jjmp1,ndex2d)
     204cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     205      CALL histwrite_phy(nid_day, "SWdnSFC",itau_w,zx_tmp_fi2d)
    218206c
    219207      ENDIF !lev_histday.GE.1
    … …  
    221209      IF(lev_histday.GE.2) THEN
    222210c
    223       CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
    224       CALL histwrite(nid_day,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    225 c
    226       CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
    227       CALL histwrite(nid_day,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     211cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
     212      CALL histwrite_phy(nid_day,"bils",itau_w, bils)
     213c
     214cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
     215      CALL histwrite_phy(nid_day,"fder",itau_w,fder)
    228216c
    229217      ENDIF !lev_histday.GE.2
    … …  
    237225c Champs 3D:
    238226c
    239       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    240       CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
    241      .                                   iim*jjmp1*klev,ndex3d)
    242 c
    243       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
    244       CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d,
    245      .                                   iim*jjmp1*klev,ndex3d)
    246 c
    247       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    248       CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
    249      .                                   iim*jjmp1*klev,ndex3d)
    250 c
    251       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    252       CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
    253      .                                   iim*jjmp1*klev,ndex3d)
    254 c
    255       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    256       CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
    257      .                                   iim*jjmp1*klev,ndex3d)
    258 c
    259       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    260       CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
    261      .                                   iim*jjmp1*klev,ndex3d)
    262 c
    263       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    264       CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
    265      .                                   iim*jjmp1*klev,ndex3d)
     227cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
     228      CALL histwrite_phy(nid_day,"temp",itau_w,t_seri)
     229c
     230cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
     231      CALL histwrite_phy(nid_day,"ovap",itau_w,qx(:,:,ivap))
     232c
     233cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
     234      CALL histwrite_phy(nid_day,"geop",itau_w,zphi)
     235c
     236cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
     237      CALL histwrite_phy(nid_day,"vitu",itau_w,u_seri)
     238c
     239cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
     240      CALL histwrite_phy(nid_day,"vitv",itau_w,v_seri)
     241c
     242cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
     243      CALL histwrite_phy(nid_day,"vitw",itau_w,omega)
     244c
     245cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
     246      CALL histwrite_phy(nid_day,"pres",itau_w,pplay)
    266247c
    267248      ENDIF !lev_histday.GE.3
    … …  
    276257c=================================================================
    277258      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, klevp1 )
    278       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    279       CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
    280      .                               iim*jjmp1,ndex2d)
     259cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     260      CALL histwrite_phy(nid_day, "SWupTOA",itau_w,zx_tmp_fi2d)
    281261c
    282262      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, 1 )
    283       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    284       CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
    285      .                               iim*jjmp1,ndex2d)
     263cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     264      CALL histwrite_phy(nid_day, "SWupSFC",itau_w,zx_tmp_fi2d)
    286265c
    287266      zx_tmp_fi2d(1 : klon) = swdn( 1 : klon, klevp1 )
    288       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    289       CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
    290      .                               iim*jjmp1,ndex2d)
     267cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     268      CALL histwrite_phy(nid_day, "SWdnTOA",itau_w,zx_tmp_fi2d)
    291269c
    292270      zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, klevp1 )
    293       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    294       CALL histwrite(nid_day, "SWupTOAclr",itau_w,zx_tmp_2d,
    295      .                               iim*jjmp1,ndex2d)
     271cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     272      CALL histwrite_phy(nid_day, "SWupTOAclr",itau_w,zx_tmp_fi2d)
    296273c
    297274      zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, 1 )
    298       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    299       CALL histwrite(nid_day, "SWupSFCclr",itau_w,zx_tmp_2d,
    300      .                               iim*jjmp1,ndex2d)
     275cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     276      CALL histwrite_phy(nid_day, "SWupSFCclr",itau_w,zx_tmp_fi2d)
    301277c
    302278      zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, klevp1 )
    303       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    304       CALL histwrite(nid_day, "SWdnTOAclr",itau_w,zx_tmp_2d,
    305      .                               iim*jjmp1,ndex2d)
     279cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     280      CALL histwrite_phy(nid_day, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
    306281c
    307282      zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, 1 )
    308       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    309       CALL histwrite(nid_day, "SWdnSFCclr",itau_w,zx_tmp_2d,
    310      .                               iim*jjmp1,ndex2d)
     283cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     284      CALL histwrite_phy(nid_day, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
    311285c
    312286      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
    313       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
    314       CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     287cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
     288      CALL histwrite_phy(nid_day,"tter",itau_w,zx_tmp_fi2d)
    315289c
    316290      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
    317       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    318       CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     291cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     292      CALL histwrite_phy(nid_day,"tlic",itau_w,zx_tmp_fi2d)
    319293c
    320294      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
    321       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    322       CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     295cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     296      CALL histwrite_phy(nid_day,"toce",itau_w,zx_tmp_fi2d)
    323297c
    324298      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
    325       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    326       CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     299cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     300      CALL histwrite_phy(nid_day,"tsic",itau_w,zx_tmp_fi2d)
    327301c
    328302      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
    329       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    330       CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
    331      .               ndex2d)
     303cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     304      CALL histwrite_phy(nid_day,"t2mter",itau_w,zx_tmp_fi2d)
    332305c
    333306      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
    334       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    335       CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
    336      .               iim*jjmp1,ndex2d)
     307cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     308      CALL histwrite_phy(nid_day,"t2mter_min",itau_w,zx_tmp_fi2d)
    337309c
    338310      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
    339       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    340       CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
    341      .               iim*jjmp1,ndex2d)
     311cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     312      CALL histwrite_phy(nid_day,"t2mter_max",itau_w,zx_tmp_fi2d)
    342313c
    343314      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
    344       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    345       CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
    346      .               ndex2d)
     315cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     316      CALL histwrite_phy(nid_day,"t2mlic",itau_w,zx_tmp_fi2d)
    347317c
    348318      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
    349       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    350       CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
    351      .               ndex2d)
     319cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     320      CALL histwrite_phy(nid_day,"t2moce",itau_w,zx_tmp_fi2d)
    352321c
    353322      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
    354       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    355       CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
    356      .               ndex2d)
     323cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     324      CALL histwrite_phy(nid_day,"t2msic",itau_w,zx_tmp_fi2d)
    357325c
    358326      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
    359       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    360       CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
    361      .               ndex2d)
     327cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     328      CALL histwrite_phy(nid_day,"u10mter",itau_w,zx_tmp_fi2d)
    362329c
    363330      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
    364       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    365       CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
    366      .               ndex2d)
     331cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     332      CALL histwrite_phy(nid_day,"v10mter",itau_w,zx_tmp_fi2d)
    367333c
    368334      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
    369       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    370       CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
    371      .               ndex2d)
     335cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     336      CALL histwrite_phy(nid_day,"u10mlic",itau_w,zx_tmp_fi2d)
    372337c
    373338      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
    374       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    375       CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
    376      .               ndex2d)
     339cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     340      CALL histwrite_phy(nid_day,"v10mlic",itau_w,zx_tmp_fi2d)
    377341c
    378342      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
    379       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    380       CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
    381      .               ndex2d)
     343cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     344      CALL histwrite_phy(nid_day,"u10moce",itau_w,zx_tmp_fi2d)
    382345c
    383346      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
    384       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    385       CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
    386      .               ndex2d)
     347cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     348      CALL histwrite_phy(nid_day,"v10moce",itau_w,zx_tmp_fi2d)
    387349c
    388350      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
    389       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    390       CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
    391      .               ndex2d)
     351cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     352      CALL histwrite_phy(nid_day,"u10msic",itau_w,zx_tmp_fi2d)
    392353C
    393354      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
    394       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    395       CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
    396      .               ndex2d)
     355cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     356      CALL histwrite_phy(nid_day,"v10msic",itau_w,zx_tmp_fi2d)
    397357C
    398358      DO nsrf = 1, nbsrf
    399359C
    400360        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
    401         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    402         CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itau_w,
    403      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     361cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     362        CALL histwrite_phy(nid_day,"pourc_"//clnsurf(nsrf),itau_w,
     363     $                     zx_tmp_fi2d)
    404364C
    405365        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
    406         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    407         CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itau_w,
    408      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     366cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     367        CALL histwrite_phy(nid_day,"tsol_"//clnsurf(nsrf),itau_w,
     368     $                     zx_tmp_fi2d)
    409369C
    410370        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
    411         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    412         CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itau_w,
    413      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     371cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     372        CALL histwrite_phy(nid_day,"sens_"//clnsurf(nsrf),itau_w,
     373     $                     zx_tmp_fi2d)
    414374C 
    415375        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
    416         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    417         CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itau_w,
    418      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     376cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     377        CALL histwrite_phy(nid_day,"lat_"//clnsurf(nsrf),itau_w,
     378     $                     zx_tmp_fi2d)
    419379C
    420380        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
    421         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    422         CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itau_w,
    423      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     381cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     382        CALL histwrite_phy(nid_day,"taux_"//clnsurf(nsrf),itau_w,
     383     $                     zx_tmp_fi2d)
    424384C     
    425385        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
    426         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    427         CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itau_w,
    428      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     386cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     387        CALL histwrite_phy(nid_day,"tauy_"//clnsurf(nsrf),itau_w,
     388     $                     zx_tmp_fi2d)
    429389C
    430390        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
    431         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    432         CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itau_w,
    433      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     391cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     392        CALL histwrite_phy(nid_day,"albe_"//clnsurf(nsrf),itau_w,
     393     $                     zx_tmp_fi2d)
    434394C
    435395        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
    436         CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    437         CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itau_w,
    438      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     396cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     397        CALL histwrite_phy(nid_day,"rugs_"//clnsurf(nsrf),itau_w,
     398     $                     zx_tmp_fi2d)
    439399C
    440400      END DO 
Note: See TracChangeset for help on using the changeset viewer.