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

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

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r644 r704  
    1414      IF(lev_histhf.GE.1) THEN
    1515c
    16       CALL gr_fi_ecrit(1, klon,iim,jjmp1, paire_ter, zx_tmp_2d)
    17       CALL histwrite(nid_hf,"aireTER",itau_w,
    18      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     16cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, paire_ter, zx_tmp_2d)
     17      CALL histwrite_phy(nid_hf,"aireTER",itau_w,paire_ter)
    1918c
    2019      DO i=1, klon
    … …  
    2221      ENDDO
    2322c
    24       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    25       CALL histwrite(nid_hf,"contfracATM",itau_w,
    26      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    27 c
    28       CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d)
    29       CALL histwrite(nid_hf,"contfracOR",itau_w,
    30      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    31 c
    32       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
    33       CALL histwrite(nid_hf,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,
    34      .               ndex2d)
    35 c
    36       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
    37       CALL histwrite(nid_hf,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,
    38      .               ndex2d)
     23cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     24      CALL histwrite_phy(nid_hf,"contfracATM",itau_w,zx_tmp_fi2d)
     25c
     26cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d)
     27      CALL histwrite_phy(nid_hf,"contfracOR",itau_w,
     28     .                   pctsrf_new(:,is_ter))
     29c
     30cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
     31      CALL histwrite_phy(nid_hf,"t2m",itau_w,zt2m)
     32c
     33cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
     34      CALL histwrite_phy(nid_hf,"q2m",itau_w,zq2m)
    3935c
    4036      DO i = 1, klon
    4137         zx_tmp_fi2d(i) = paprs(i,1)
    4238      ENDDO
    43       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    44       CALL histwrite(nid_hf,"psol",itau_w,zx_tmp_2d,iim*jjmp1,
    45      .               ndex2d)
     39cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     40      CALL histwrite_phy(nid_hf,"psol",itau_w,zx_tmp_fi2d)
    4641c
    4742      DO i = 1, klon
    4843         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
    4944      ENDDO
    50       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    51       CALL histwrite(nid_hf,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     45cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     46      CALL histwrite_phy(nid_hf,"precip",itau_w,zx_tmp_fi2d)
    5247c
    5348c ENSEMBLES BEG
    54       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    55       CALL histwrite(nid_hf,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    56 c
    57       CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
    58       CALL histwrite(nid_hf,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    59 c
    60 c
    61       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
    62       CALL histwrite(nid_hf,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    63 c
    64       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
    65       CALL histwrite(nid_hf,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     49cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     50      CALL histwrite_phy(nid_hf,"tsol",itau_w,zxtsol)
     51c
     52cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     53      CALL histwrite_phy(nid_hf,"slp",itau_w,slp)
     54c
     55cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
     56      CALL histwrite_phy(nid_hf,"u10m",itau_w,zu10m)
     57c
     58cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
     59      CALL histwrite_phy(nid_hf,"v10m",itau_w,zv10m)
    6660c
    6761      DO i=1, klon
    6862       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
    6963      ENDDO
    70       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    71       CALL histwrite(nid_hf,"wind10m",itau_w,zx_tmp_2d,
    72      $     iim*jjmp1,ndex2d)
     64cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     65      CALL histwrite_phy(nid_hf,"wind10m",itau_w,zx_tmp_fi2d)
    7366c
    7467      DO k=1, nlevSTD
    … …  
    8578      IF(bb2.EQ."500") THEN
    8679c
    87       CALL gr_fi_ecrit(1, klon,iim,jjmp1, philevSTD(:,k),zx_tmp_2d)
    88       CALL histwrite(nid_hf,"phi"//bb2,itau_w,zx_tmp_2d,
    89      $       iim*jjmp1,ndex2d)
     80cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philevSTD(:,k),zx_tmp_2d)
     81      CALL histwrite_phy(nid_hf,"phi"//bb2,itau_w, philevSTD(:,k))
    9082      ENDIF !bb2.EQ."500"
    9183      ENDDO
    … …  
    10193cIM 140904 END
    10294cIM 140904   CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    103       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    104       CALL histwrite(nid_hf,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     95cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     96      CALL histwrite_phy(nid_hf,"cldt",itau_w,zx_tmp_fi2d)
    10597c
    10698      zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
    107       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    108       CALL histwrite(nid_hf,"SWnetOR",itau_w,
    109      $     zx_tmp_2d,iim*jjmp1,ndex2d)
     99cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     100      CALL histwrite_phy(nid_hf,"SWnetOR",itau_w, zx_tmp_fi2d)
    110101c
    111102      zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol(1:klon))
    112       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    113       CALL histwrite(nid_hf,"SWdownOR",itau_w,
    114      $     zx_tmp_2d,iim*jjmp1,ndex2d)
    115 c
    116       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    117       CALL histwrite(nid_hf,"LWdownOR",itau_w,zx_tmp_2d,iim*jjmp1,
    118      $     ndex2d)
     103cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     104      CALL histwrite_phy(nid_hf,"SWdownOR",itau_w, zx_tmp_fi2d)
     105c
     106cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
     107      CALL histwrite_phy(nid_hf,"LWdownOR",itau_w,sollwdown)
    119108c
    120109c
    … …  
    137126     $   bb2.EQ."500".OR.bb2.EQ."200") THEN
    138127c
    139       CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
    140       CALL histwrite(nid_hf,"t"//bb2,itau_w,zx_tmp_2d,
    141      $       iim*jjmp1,ndex2d)
     128cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
     129      CALL histwrite_phy(nid_hf,"t"//bb2,itau_w, tlevSTD(:,k))
    142130c
    143131      IF(bb2.NE."500") THEN !bb2.NE."500"
    144132c
    145       CALL gr_fi_ecrit(1, klon,iim,jjmp1,philevSTD(:,k),zx_tmp_2d)
    146       CALL histwrite(nid_hf,"phi"//bb2,itau_w,zx_tmp_2d,
    147      $       iim*jjmp1,ndex2d)
     133cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,philevSTD(:,k),zx_tmp_2d)
     134      CALL histwrite_phy(nid_hf,"phi"//bb2,itau_w,philevSTD(:,k))
    148135      ENDIF !bb2.NE."500"
    149136c
    150       CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d)
    151       CALL histwrite(nid_hf,"q"//bb2,itau_w,zx_tmp_2d,
    152      $       iim*jjmp1,ndex2d)
    153 c
    154       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
    155       CALL histwrite(nid_hf,"u"//bb2,itau_w,zx_tmp_2d,
    156      $       iim*jjmp1,ndex2d)
    157 c
    158       CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d)
    159       CALL histwrite(nid_hf,"v"//bb2,itau_w,zx_tmp_2d,
    160      $       iim*jjmp1,ndex2d)
     137cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d)
     138      CALL histwrite_phy(nid_hf,"q"//bb2,itau_w, qlevSTD(:,k))
     139c
     140cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
     141      CALL histwrite_phy(nid_hf,"u"//bb2,itau_w,ulevSTD(:,k))
     142c
     143cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d)
     144      CALL histwrite_phy(nid_hf,"v"//bb2,itau_w, vlevSTD(:,k))
    161145c
    162146      ENDIF !bb2.EQ."850".OR.bb2.EQ."700"
    … …  
    164148c
    165149      IF(1.EQ.0) THEN
    166       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    167       CALL histwrite(nid_hf,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    168 c
    169       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    170       CALL histwrite(nid_hf,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     150cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
     151      CALL histwrite_phy(nid_hf,"cdrm",itau_w,cdragm)
     152c
     153cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
     154      CALL histwrite_phy(nid_hf,"cdrh",itau_w,cdragh)
    171155      ENDIF !(1.EQ.0) THEN
    172156c
    … …  
    199183        IF(bb2.EQ."850".or.bb2.EQ."200") THEN
    200184c
    201         CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
    202         CALL histwrite(nid_hf,"u"//bb2,itau_w,zx_tmp_2d,
    203      $       iim*jjmp1,ndex2d)
    204 c
    205         CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d)
    206         CALL histwrite(nid_hf,"v"//bb2,itau_w,zx_tmp_2d,
    207      $       iim*jjmp1,ndex2d)
     185cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
     186        CALL histwrite_phy(nid_hf,"u"//bb2,itau_w,ulevSTD(:,k))
     187c
     188cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d)
     189        CALL histwrite_phy(nid_hf,"v"//bb2,itau_w,vlevSTD(:,k))
    208190c
    209191        ENDIF !bb2.EQ."850".or.bb2.EQ."200"
    … …  
    211193      ENDDO !nlevSTD
    212194c
    213       CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    214       CALL histwrite(nid_hf,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     195cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
     196      CALL histwrite_phy(nid_hf,"topl",itau_w,toplw)
    215197c
    216198      DO i = 1, klon
    217199         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
    218200      ENDDO
    219       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    220       CALL histwrite(nid_hf,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    221 c
    222       CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
    223       CALL histwrite(nid_hf,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     201cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     202      CALL histwrite_phy(nid_hf,"precip",itau_w,zx_tmp_fi2d)
     203c
     204cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     205      CALL histwrite_phy(nid_hf,"slp",itau_w,slp)
    224206c
    225207cIM 280504 BEG
    … …  
    241223        IF(bb2.EQ."850".or.bb2.EQ."500") THEN
    242224c
    243          CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
    244          CALL histwrite(nid_hf,"t"//bb2,itau_w,zx_tmp_2d,
    245      $       iim*jjmp1,ndex2d)
    246 c
    247          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d)
    248          CALL histwrite(nid_hf,"q"//bb2,itau_w,zx_tmp_2d,
    249      $       iim*jjmp1,ndex2d)
     225cym         CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
     226         CALL histwrite_phy(nid_hf,"t"//bb2,itau_w,tlevSTD(:,k))
     227c
     228cym         CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d)
     229         CALL histwrite_phy(nid_hf,"q"//bb2,itau_w, qlevSTD(:,k))
    250230c
    251231        ENDIF !bb2.EQ."850".or.bb2.EQ."500"
    … …  
    253233        IF(bb2.EQ."500") THEN
    254234c
    255          CALL gr_fi_ecrit(1, klon,iim,jjmp1, philevSTD(:,k),zx_tmp_2d)
    256          CALL histwrite(nid_hf,"phi"//bb2,itau_w,zx_tmp_2d,
    257      $       iim*jjmp1,ndex2d)
    258 c
    259         CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlevSTD(:,k),zx_tmp_2d)
    260         CALL histwrite(nid_hf,"w"//bb2,itau_w,zx_tmp_2d,
    261      $       iim*jjmp1,ndex2d)
     235cym         CALL gr_fi_ecrit(1, klon,iim,jjmp1, philevSTD(:,k),zx_tmp_2d)
     236         CALL histwrite_phy(nid_hf,"phi"//bb2,itau_w, philevSTD(:,k))
     237c
     238cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlevSTD(:,k),zx_tmp_2d)
     239        CALL histwrite_phy(nid_hf,"w"//bb2,itau_w, wlevSTD(:,k))
    262240c
    263241        ENDIF !bb2.EQ."500"
    … …  
    265243        IF(bb3.EQ."50") THEN
    266244c
    267          CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
    268          CALL histwrite(nid_hf,"t"//bb3,itau_w,zx_tmp_2d,
    269      $       iim*jjmp1,ndex2d)
    270 c
    271         CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
    272         CALL histwrite(nid_hf,"u"//bb3,itau_w,zx_tmp_2d,
    273      $       iim*jjmp1,ndex2d)
    274 c
    275         CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d)
    276         CALL histwrite(nid_hf,"v"//bb3,itau_w,zx_tmp_2d,
    277      $       iim*jjmp1,ndex2d)
     245cym         CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
     246         CALL histwrite_phy(nid_hf,"t"//bb3,itau_w, tlevSTD(:,k))
     247c
     248cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
     249        CALL histwrite_phy(nid_hf,"u"//bb3,itau_w, ulevSTD(:,k))
     250c
     251cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d)
     252        CALL histwrite_phy(nid_hf,"v"//bb3,itau_w, vlevSTD(:,k))
    278253c
    279254        ENDIF !bb3.EQ."50"
    … …  
    281256       ENDDO !k=1, nlevSTD
    282257c
    283       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
    284       CALL histwrite(nid_hf,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    285 c
    286       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
    287       CALL histwrite(nid_hf,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    288 c
    289       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
    290       CALL histwrite(nid_hf,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,
    291      $               ndex2d)
    292 c
    293       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
    294       CALL histwrite(nid_hf,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,
    295      $               ndex2d)
     258cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
     259      CALL histwrite_phy(nid_hf,"u10m",itau_w,zu10m)
     260c
     261cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
     262      CALL histwrite_phy(nid_hf,"v10m",itau_w,zv10m)
     263c
     264cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
     265      CALL histwrite_phy(nid_hf,"t2m",itau_w, zt2m)
     266c
     267cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
     268      CALL histwrite_phy(nid_hf,"q2m",itau_w,zq2m)
    296269c
    297270      DO i = 1, klon
    298271         zx_tmp_fi2d(i) = paprs(i,1)
    299272      ENDDO
    300       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    301       CALL histwrite(nid_hf,"psol",itau_w,zx_tmp_2d,iim*jjmp1,
    302      $               ndex2d)
     273cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     274      CALL histwrite_phy(nid_hf,"psol",itau_w,zx_tmp_fi2d)
    303275c
    304276cIM 140904 BEG
    … …  
    308280cIM 140904 END
    309281cIM 140904   CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    310       CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    311       CALL histwrite(nid_hf,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     282cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     283      CALL histwrite_phy(nid_hf,"cldt",itau_w,zx_tmp_fi2d)
    312284c
    313285      zx_tmp_fi2d(1 : klon) = zxfluxu(1 : klon, 1)
    314       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    315       CALL histwrite(nid_hf,"taux",itau_w,
    316      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     286cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     287      CALL histwrite_phy(nid_hf,"taux",itau_w,zx_tmp_fi2d)
    317288C
    318289      zx_tmp_fi2d(1 : klon) = zxfluxv(1 : klon, 1)
    319       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    320       CALL histwrite(nid_hf,"tauy",itau_w,
    321      $      zx_tmp_2d,iim*jjmp1,ndex2d)
    322 c
    323       CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
    324       CALL histwrite(nid_hf,"prw",itau_w,zx_tmp_2d,
    325      $               iim*jjmp1,ndex2d)
     290cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     291      CALL histwrite_phy(nid_hf,"tauy",itau_w,zx_tmp_fi2d)
     292c
     293cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
     294      CALL histwrite_phy(nid_hf,"prw",itau_w,prw)
    326295c
    327296      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
    328       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    329       CALL histwrite(nid_hf,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    330 c
    331       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
    332       CALL histwrite(nid_hf,"flat",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     297cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     298      CALL histwrite_phy(nid_hf,"sens",itau_w,zx_tmp_fi2d)
     299c
     300cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
     301      CALL histwrite_phy(nid_hf,"flat",itau_w,zxfluxlat)
    333302c
    334303c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
    335 c     CALL histwrite(nid_hf,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
    336 c    $               ndex2d)
     304c     CALL histwrite_phy(nid_hf,"snow_mass",itau_w,zsnow_mass)
    337305c
    338306      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
    339       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    340       CALL histwrite(nid_hf, "SWdnSFC",itau_w,zx_tmp_2d,
    341      $                               iim*jjmp1,ndex2d)
     307cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     308      CALL histwrite_phy(nid_hf, "SWdnSFC",itau_w,zx_tmp_fi2d)
    342309c
    343310      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
    344       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    345       CALL histwrite(nid_hf, "SWupSFC",itau_w,zx_tmp_2d,
    346      $                               iim*jjmp1,ndex2d)
    347 c
    348       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    349       CALL histwrite(nid_hf,"LWdnSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    350      $               ndex2d)
     311cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     312      CALL histwrite_phy(nid_hf, "SWupSFC",itau_w,zx_tmp_fi2d)
     313c
     314cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
     315      CALL histwrite_phy(nid_hf,"LWdnSFC",itau_w,sollwdown)
    351316c
    352317cIM 090904   zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
    353318cIM 130904   zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
    354319      zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
    355       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    356       CALL histwrite(nid_hf,"LWupSFC",itau_w,zx_tmp_2d,iim*jjmp1,
    357      $               ndex2d)                   
     320cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     321      CALL histwrite_phy(nid_hf,"LWupSFC",itau_w,zx_tmp_fi2d)                   
    358322c
    359323      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
    360       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    361       CALL histwrite(nid_hf, "SWdnTOA",itau_w,zx_tmp_2d,
    362      $                               iim*jjmp1,ndex2d)
     324cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     325      CALL histwrite_phy(nid_hf, "SWdnTOA",itau_w,zx_tmp_fi2d)
    363326c
    364327      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
    365       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    366       CALL histwrite(nid_hf, "SWupTOA",itau_w,zx_tmp_2d,
    367      $                               iim*jjmp1,ndex2d)
    368 c
    369       CALL gr_fi_ecrit(nbteta,klon,iim,jjmp1,PVteta,zx_tmp_3dte)
     328cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     329      CALL histwrite_phy(nid_hf, "SWupTOA",itau_w,zx_tmp_fi2d)
     330c
     331      IF (monocpu) THEN
     332cym      CALL gr_fi_ecrit(nbteta,klon,iim,jjmp1,PVteta,zx_tmp_3dte)
    370333      DO k=1, nbteta
    371        CALL histwrite(nid_hf,"PV"//ctetaSTD(k),
    372      .      itau_w,zx_tmp_3dte(:,:,k),iim*jjmp1,ndex2d)
     334       CALL histwrite_phy(nid_hf,"PV"//ctetaSTD(k),
     335     .      itau_w,PVteta)
    373336      ENDDO !k=1, nbteta
     337      ENDIF
    374338c
    375339c
Note: See TracChangeset for help on using the changeset viewer.