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

    r633 r634  
    2424c     CALL histwrite(nid_hf,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d)
    2525C
    26       CALL gr_fi_ecrit(1, klon,iim,jjmp1, paire_ter, zx_tmp_2d)
    27       CALL histwrite(nid_hf,"aireTER",itau_w,
    28      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     26cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, paire_ter, zx_tmp_2d)
     27      CALL histwrite_phy(nid_hf,"aireTER",itau_w,paire_ter)
    2928c
    3029      DO i=1, klon
    … …  
    3231      ENDDO
    3332c
    34       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    35       CALL histwrite(nid_hf,"contfracATM",itau_w,
    36      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     33cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     34      CALL histwrite_phy(nid_hf,"contfracATM",itau_w,zx_tmp_fi2d)
    3735c
    38       CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d)
    39       CALL histwrite(nid_hf,"contfracOR",itau_w,
    40      $      zx_tmp_2d,iim*jjmp1,ndex2d)
     36cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d)
     37      CALL histwrite_phy(nid_hf,"contfracOR",itau_w,
     38     .                   pctsrf_new(:,is_ter))
    4139c
    42       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
    43       CALL histwrite(nid_hf,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,
    44      .               ndex2d)
     40cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
     41      CALL histwrite_phy(nid_hf,"t2m",itau_w,zt2m)
    4542c
    46       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
    47       CALL histwrite(nid_hf,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,
    48      .               ndex2d)
     43cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
     44      CALL histwrite_phy(nid_hf,"q2m",itau_w,zq2m)
    4945c
    5046      DO i = 1, klon
    5147         zx_tmp_fi2d(i) = paprs(i,1)
    5248      ENDDO
    53       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    54       CALL histwrite(nid_hf,"psol",itau_w,zx_tmp_2d,iim*jjmp1,
    55      .               ndex2d)
     49cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     50      CALL histwrite_phy(nid_hf,"psol",itau_w,zx_tmp_fi2d)
    5651c
    5752      DO i = 1, klon
    5853         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
    5954      ENDDO
    60       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    61       CALL histwrite(nid_hf,"rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     55cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     56      CALL histwrite_phy(nid_hf,"rain",itau_w,zx_tmp_fi2d)
    6257c
    6358c ENSEMBLES BEG
    64       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    65       CALL histwrite(nid_hf,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     59cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
     60      CALL histwrite_phy(nid_hf,"tsol",itau_w,zxtsol)
    6661c
    67       CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
    68       CALL histwrite(nid_hf,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     62cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     63      CALL histwrite_phy(nid_hf,"slp",itau_w,slp)
    6964c
    70       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
    71       CALL histwrite(nid_hf,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     65cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
     66      CALL histwrite_phy(nid_hf,"u10m",itau_w,zu10m)
    7267c
    73       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
    74       CALL histwrite(nid_hf,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     68cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
     69      CALL histwrite_phy(nid_hf,"v10m",itau_w,zv10m)
    7570c
    7671      DO i=1, klon
    7772       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
    7873      ENDDO
    79       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    80       CALL histwrite(nid_hf,"wind10m",itau_w,zx_tmp_2d,
    81      $     iim*jjmp1,ndex2d)
     74cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     75      CALL histwrite_phy(nid_hf,"wind10m",itau_w,zx_tmp_fi2d)
    8276c
    8377      DO k=1, nlevENS
    8478      IF(clev(k).EQ."500") THEN
    85       CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
    86       CALL histwrite(nid_hf,"phi"//clev(k),itau_w,zx_tmp_2d,
    87      $       iim*jjmp1,ndex2d)
     79cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
     80      CALL histwrite_phy(nid_hf,"phi"//clev(k),itau_w,philev(:,k))
    8881      ENDIF !clev(k).EQ."500"
    8982      ENDDO
    … …  
    9386      IF(lev_histhf.GE.2) THEN
    9487c
    95       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    96       CALL histwrite(nid_hf,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     88cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
     89      CALL histwrite_phy(nid_hf,"cldt",itau_w,cldt)
    9790c
    9891      zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
    99       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    100       CALL histwrite(nid_hf,"SWnetOR",itau_w,
    101      $     zx_tmp_2d,iim*jjmp1,ndex2d)
     92      CALL histwrite_phy(nid_hf,"SWnetOR",itau_w,zx_tmp_fi2d)
    10293c
    10394      zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol(1:klon))
    104       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    105       CALL histwrite(nid_hf,"SWdownOR",itau_w,
    106      $     zx_tmp_2d,iim*jjmp1,ndex2d)
     95      CALL histwrite(nid_hf,"SWdownOR",itau_w,zx_tmp_fi2d)
    10796c
    108       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    109       CALL histwrite(nid_hf,"LWdownOR",itau_w,zx_tmp_2d,iim*jjmp1,
    110      $     ndex2d)
     97      CALL histwrite(nid_hf,"LWdownOR",itau_w,zx_tmp_fi2d)
    11198c
    11299c
    … …  
    117104      DO k=1, nlevENS
    118105c
    119       CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlev(:,k),zx_tmp_2d)
    120       CALL histwrite(nid_hf,"t"//clev(k),itau_w,zx_tmp_2d,
    121      $       iim*jjmp1,ndex2d)
     106cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlev(:,k),zx_tmp_2d)
     107      CALL histwrite_phy(nid_hf,"t"//clev(k),itau_w,tlev(:,k))
    122108c
    123109      IF(clev(k).NE."500") THEN !clev(k).NE."500"
    124       CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
    125       CALL histwrite(nid_hf,"phi"//clev(k),itau_w,zx_tmp_2d,
    126      $       iim*jjmp1,ndex2d)
     110cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
     111      CALL histwrite_phy(nid_hf,"phi"//clev(k),itau_w,philev(:,k))
    127112      ENDIF !clev(k).NE."500"
    128113c
    129       CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlev(:,k),zx_tmp_2d)
    130       CALL histwrite(nid_hf,"q"//clev(k),itau_w,zx_tmp_2d,
    131      $       iim*jjmp1,ndex2d)
     114cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlev(:,k),zx_tmp_2d)
     115      CALL histwrite_phy(nid_hf,"q"//clev(k),itau_w,qlev(:,k))
    132116c
    133117      IF(1.EQ.0) THEN
    134       CALL gr_fi_ecrit(1, klon,iim,jjmp1, rhlev(:,k),zx_tmp_2d)
    135       CALL histwrite(nid_hf,"rh"//clev(k),itau_w,zx_tmp_2d,
    136      $       iim*jjmp1,ndex2d)
     118cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, rhlev(:,k),zx_tmp_2d)
     119      CALL histwrite_phy(nid_hf,"rh"//clev(k),itau_w,rhlev(:,k))
    137120      ENDIF !1.EQ.0
    138121c
    139       CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
    140       CALL histwrite(nid_hf,"u"//clev(k),itau_w,zx_tmp_2d,
    141      $       iim*jjmp1,ndex2d)
     122cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
     123      CALL histwrite_phy(nid_hf,"u"//clev(k),itau_w,ulev(:,k))
    142124c
    143       CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
    144       CALL histwrite(nid_hf,"v"//clev(k),itau_w,zx_tmp_2d,
    145      $       iim*jjmp1,ndex2d)
     125cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
     126      CALL histwrite_phy(nid_hf,"v"//clev(k),itau_w,vlev(:,k))
    146127c
    147128      ENDDO !nlevENS
    148129c
    149130      IF(1.EQ.0) THEN
    150       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    151       CALL histwrite(nid_hf,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     131cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
     132      CALL histwrite_phy(nid_hf,"cdrm",itau_w,cdragm)
    152133c
    153       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    154       CALL histwrite(nid_hf,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     134cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
     135      CALL histwrite_phy(nid_hf,"cdrh",itau_w,cdragh)
    155136      ENDIF !(1.EQ.0) THEN
    156137c
Note: See TracChangeset for help on using the changeset viewer.