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

    r467 r506  
    1212         i = NINT(zout/zsto)
    1313         CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
     14c
     15      IF(lev_histday.GE.1) THEN
     16c
    1417       CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    1518         varname = 'phis'
    1619         vartitle= 'Surface geop. height'
    1720         varunits= '-'
    18 c        call writephy(fid_day,prof2d_on,varname,pphis,vartitle,
    19 c    .                                                    varunits)
    2021c
    2122         i = NINT(zout/zsto)
    … …  
    2526         vartitle= 'Grid area'
    2627         varunits= '-'
    27 c        call writephy(fid_day,prof2d_on,varname,paire,vartitle,
    28 c    .                                                    varunits)
    29 C
     28c
     29      DO i=1, klon
     30       zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
     31      ENDDO
     32c
     33      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     34      CALL histwrite(nid_day,"contfracATM",itau_w,
     35     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     36c
     37      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d)
     38      CALL histwrite(nid_day,"contfracOR",itau_w,
     39     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     40c
    3041      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    3142      CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    32 c     call writephy(fid_day,prof2d_av,'tsol',zxtsol,
    33 c    .              'Surface Temperature','K')
    34 c
    35 C
    3643C
    3744      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    3845      CALL histwrite(nid_day,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    3946c
     47c   En attendant un eventuel debugage.
     48
     49      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     50      CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
     51     .               iim*jjmp1,ndex2d)
     52c
     53      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     54      CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
     55     .               iim*jjmp1,ndex2d)
     56c
     57      DO i = 1, klon
     58         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
     59      ENDDO
     60      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     61      CALL histwrite(nid_day,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     62c
     63      DO i = 1, klon
     64         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
     65      ENDDO
     66      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     67      CALL histwrite(nid_day,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     68c
     69      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_lsc,zx_tmp_2d)
     70      CALL histwrite(nid_day,"snowl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     71c
     72      IF(1.EQ.0) THEN
     73      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_con,zx_tmp_2d)
     74      CALL histwrite(nid_day,"snowc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     75      ENDIF
     76c
     77      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
     78      CALL histwrite(nid_day,"flat",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     79c
     80      DO i = 1, klon
     81         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
     82      ENDDO
     83c
     84      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     85      CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     86c
    4087      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zq2m,zx_tmp_2d)
    4188      CALL histwrite(nid_day,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    … …  
    4996     .               ndex2d)
    5097c
    51       DO i = 1, klon
     98      DO i=1, klon
     99       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
     100      ENDDO
     101      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     102      CALL histwrite(nid_day,"wind10m",itau_w,zx_tmp_2d,
     103     $     iim*jjmp1,ndex2d)
     104c
     105      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     106      CALL histwrite(nid_day,"wind10max",itau_w,zx_tmp_2d,
     107     .               iim*jjmp1,ndex2d)
     108c
     109      DO i=1, klon
    52110         zx_tmp_fi2d(i) = paprs(i,1)
    53111      ENDDO
    54112      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    55113      CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    56 c Essai writephys
    57       varname = 'psol'
    58       vartitle= 'pression au sol'
    59       varunits= 'hPa'
    60 c     call writephy(fid_day,prof2d_av,varname,zx_tmp_fi2d,vartitle,
    61 c    .                                                    varunits)
    62 c
    63       DO i = 1, klon
     114c
     115      DO i=1, klon
    64116         zx_tmp_fi2d(i) = (rain_fall(i) + snow_fall(i))
    65117      ENDDO
    66118      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    67119      CALL histwrite(nid_day,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    68 c     call writephy(fid_day,prof2d_av,'rain',zx_tmp_fi2d,
    69 c    .              'Precipitation','mm/day')
    70 
    71 
    72120c
    73121      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
    74       CALL histwrite(nid_day,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    75 c     call writephy(fid_day,prof2d_av,'snow',snow_fall,
    76 c    .              'Snow','mm/day')
     122      CALL histwrite(nid_day,"snowf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    77123c
    78124      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
    79125      CALL histwrite(nid_day,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
    80126     .               ndex2d)
    81 c     call writephy(fid_day,prof2d_av,'snow_mass',zxsnow,
    82 c    .              'Snow cover','mm')
    83127c
    84128      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
    85129      CALL histwrite(nid_day,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    86 c     call writephy(fid_day,prof2d_av,'evap',evap,
    87 c    .              'Evaporation','mm/day')
    88130c
    89131      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
    90132      CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    91 c     call writephy(fid_day,prof2d_av,'tops',topsw,
    92 c    .              'Solar rad. at TOA','W/m2')
    93133c
    94134      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    95135      CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    96 c     call writephy(fid_day,prof2d_av,'topl',toplw,
    97 c    .              'IR rad. at TOA','W/m2')
    98136c
    99137      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
    100138      CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    101 c     call writephy(fid_day,prof2d_av,'sols',solsw,
    102 c    .              'Solar rad. at surf.','W/m2')
    103139c
    104140      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
    105141      CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    106 c     call writephy(fid_day,prof2d_av,'soll',sollw,
    107 c    .              'IR rad. at surface','W/m2')
     142c
     143      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
     144      CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     145c
     146      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
     147      CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     148c
     149      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
     150      CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     151c
     152      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
     153      CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     154c
     155      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
     156      CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     157c
     158      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
     159      CALL histwrite(nid_day,"prw",itau_w,zx_tmp_2d,
     160     .               iim*jjmp1,ndex2d)
     161
     162c   Ecriture de champs dynamiques sur des niveaux de pression
     163      DO k=1, nlevENS
     164c
     165        CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
     166        CALL histwrite(nid_day,"u"//clev(k),itau_w,zx_tmp_2d,
     167     $       iim*jjmp1,ndex2d)
     168c
     169        CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
     170        CALL histwrite(nid_day,"v"//clev(k),itau_w,zx_tmp_2d,
     171     $       iim*jjmp1,ndex2d)
     172c
     173      ENDDO !nlevENS
     174c
     175      DO k=1, nlevENS
     176      IF(clev(k).EQ."500") THEN
     177c
     178      CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlev(:,indENS(k)),zx_tmp_2d)
     179      CALL histwrite(nid_day,"w500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     180c
     181c     DO i=1, klon
     182c     print*,' ind, phi 500hPa',i,indENS(k),philev(:,indENS(k))
     183c     print*,' ind, phi 500hPa',i,indENS(k)
     184c     ENDDO
     185c
     186cBAD  CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,indENS(k)),zx_tmp_2d)
     187      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
     188      CALL histwrite(nid_day,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     189c
     190      ENDIF !clev(k).EQ.500
     191      ENDDO !k=1, nlevENS
     192c
     193      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     194      CALL histwrite(nid_day,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     195
     196      CALL gr_fi_ecrit(1,klon,iim,jjmp1,cape,zx_tmp_2d)
     197      CALL histwrite(nid_day,"cape_max",itau_w,zx_tmp_2d,
     198     .               iim*jjmp1,ndex2d)
    108199c
    109200      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    110201      CALL histwrite(nid_day,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
    111202     .               ndex2d)
    112 c     call writephy(fid_day,prof2d_av,'solldown',sollwdown,
    113 c    .              'Down. IR rad. at surface','W/m2')
     203c
     204      DO i=1, klon
     205        zx_tmp_fi2d(i)=-1*sens(i)
     206      ENDDO
     207      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     208      CALL histwrite(nid_day,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     209c
     210      zx_tmp_fi2d(1 : klon) = swdn(1 : klon, 1)
     211
     212      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     213      CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
     214     .                               iim*jjmp1,ndex2d)
     215c
     216      ENDIF !lev_histday.GE.1
     217c
     218      IF(lev_histday.GE.2) THEN
    114219c
    115220      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
    116221      CALL histwrite(nid_day,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    117 c     call writephy(fid_day,prof2d_av,'bils',bils,
    118 c    .              'Surf. total heat flux','W/m2')
    119 c
    120       CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
    121       CALL histwrite(nid_day,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    122 c     call writephy(fid_day,prof2d_av,'sens',sens,
    123 c    .              'Sensible heat flux','W/m2')
    124222c
    125223      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
    126224      CALL histwrite(nid_day,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    127 c     call writephy(fid_day,prof2d_av,'fder',fder,
    128 c    .              'Heat flux derivation','W/m2')
    129 c
    130 c
    131 cXXX      DO i = 1, klon
    132 cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
    133 cXXX      ENDDO
    134 cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    135 cXXX      CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    136 c
    137       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
    138       CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    139 c     call writephy(fid_day,prof2d_av,'cldl',cldl,
    140 c    .              'Low-level cloudiness','-')
    141 c
    142       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
    143       CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    144 c     call writephy(fid_day,prof2d_av,'cldm',cldm,
    145 c    .              'Mid-level cloudiness','-')
    146 c
    147       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
    148       CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    149 c     call writephy(fid_day,prof2d_av,'cldh',cldh,
    150 c    .              'High-level cloudiness','-')
    151 c
    152       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    153       CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    154 c     call writephy(fid_day,prof2d_av,'cldt',cldt,
    155 c    .              'Total cloudiness','-')
    156 c
    157       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
    158       CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    159 c     call writephy(fid_day,prof2d_av,'cldq',cldq,
    160 c    .              'Cloud liquid water path','-')
    161 c
    162 c     zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
    163       zx_tmp_fi2d(1 : klon) = swup( 1 : klon, 2 )
    164 
     225c
     226      ENDIF !lev_histday.GE.2
     227c
     228      IF(lev_histday.GE.3) THEN
     229c=================================================================
     230c   ECRITURE DES CHAMPS 3D
     231c=================================================================
     232
     233c
     234c Champs 3D:
     235c
     236      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
     237      CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
     238     .                                   iim*jjmp1*klev,ndex3d)
     239c
     240      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
     241      CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d,
     242     .                                   iim*jjmp1*klev,ndex3d)
     243c
     244      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
     245      CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
     246     .                                   iim*jjmp1*klev,ndex3d)
     247c
     248      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
     249      CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
     250     .                                   iim*jjmp1*klev,ndex3d)
     251c
     252      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
     253      CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
     254     .                                   iim*jjmp1*klev,ndex3d)
     255c
     256      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
     257      CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
     258     .                                   iim*jjmp1*klev,ndex3d)
     259c
     260      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
     261      CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
     262     .                                   iim*jjmp1*klev,ndex3d)
     263c
     264      ENDIF !lev_histday.GE.3
     265c=================================================================
     266c   FIN ECRITURE DES CHAMPS 3D
     267c=================================================================
     268      IF(lev_histday.ge.4) THEN
     269c=================================================================
     270c
     271c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
     272c
     273c=================================================================
     274      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, klevp1 )
    165275      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    166276      CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
    167277     .                               iim*jjmp1,ndex2d)
    168278c
    169 c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
    170279      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, 1 )
    171 
    172280      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    173281      CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
    174282     .                               iim*jjmp1,ndex2d)
    175283c
    176 c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
    177       zx_tmp_fi2d(1 : klon) = swdn( 1 : klon, 2 )
    178 
     284      zx_tmp_fi2d(1 : klon) = swdn( 1 : klon, klevp1 )
    179285      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    180286      CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
    181287     .                               iim*jjmp1,ndex2d)
    182288c
    183 c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
    184       zx_tmp_fi2d(1 : klon) = swdn( 1 : klon, 1 )
    185 
    186       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    187       CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
    188      .                               iim*jjmp1,ndex2d)
    189 
    190 c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
    191       zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, 2 )
    192 
     289      zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, klevp1 )
    193290      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    194291      CALL histwrite(nid_day, "SWupTOAclr",itau_w,zx_tmp_2d,
    195292     .                               iim*jjmp1,ndex2d)
    196 
    197 c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
     293c
    198294      zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, 1 )
    199 
    200295      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    201296      CALL histwrite(nid_day, "SWupSFCclr",itau_w,zx_tmp_2d,
    202297     .                               iim*jjmp1,ndex2d)
    203  
    204 c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
    205       zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, 2 )
    206 
     298c
     299      zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, klevp1 )
    207300      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    208301      CALL histwrite(nid_day, "SWdnTOAclr",itau_w,zx_tmp_2d,
    209302     .                               iim*jjmp1,ndex2d)
    210  
    211 c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
     303c
    212304      zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, 1 )
    213 
    214305      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    215306      CALL histwrite(nid_day, "SWdnSFCclr",itau_w,zx_tmp_2d,
    216307     .                               iim*jjmp1,ndex2d)
    217 cIM
    218       CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
    219       CALL histwrite(nid_day,"prw",itau_w,zx_tmp_2d,
    220      .               iim*jjmp1,ndex2d)
    221 
    222 c   Ecriture de champs dynamiques sur des niveaux de pression
    223 
    224       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
    225       CALL histwrite(nid_day,"u850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    226 
    227       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
    228       CALL histwrite(nid_day,"v850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    229 
    230       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
    231       CALL histwrite(nid_day,"u500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    232 
    233       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
    234       CALL histwrite(nid_day,"v500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    235 
    236       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
    237       CALL histwrite(nid_day,"u200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    238 
    239       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
    240       CALL histwrite(nid_day,"v200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    241 
    242       CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
    243       CALL histwrite(nid_day,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    244 
    245       CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
    246       CALL histwrite(nid_day,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    247 
    248       CALL gr_fi_ecrit(1,klon,iim,jjmp1,cape,zx_tmp_2d)
    249       CALL histwrite(nid_day,"cape_max",itau_w,zx_tmp_2d,
    250      .               iim*jjmp1,ndex2d)
    251 
    252 
    253 c
    254 c=================================================================
    255 c   CI DESSOUS, SORTIES A UTILISER EN DEBUG
    256 c=================================================================
    257       if ( lev_histday.gt.1 ) then
    258 
    259 c   En attendant un eventuel debugage.
    260 
    261       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    262       CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
    263      .               iim*jjmp1,ndex2d)
    264 c
    265       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    266       CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
    267      .               iim*jjmp1,ndex2d)
    268 c
    269 c=================================================================
    270 c=================================================================
    271 c=================================================================
    272 c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
    273 c=================================================================
    274 
     308c
    275309      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
    276310      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
    277311      CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    278 c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
    279 c    .              'Surface Temperature','K')
    280 C
     312c
    281313      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
    282314      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    283315      CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    284 c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
    285 c    .              'Surface Temperature','K')
    286 C
     316c
    287317      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
    288318      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    289319      CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    290 c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
    291 c    .              'Surface Temperature','K')
    292 C
     320c
    293321      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
    294322      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    295323      CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    296 c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
    297 c    .              'Surface Temperature','K')
    298 C
    299 cccIM
     324c
    300325      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    301326      CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
    … …  
    368393C
    369394      DO nsrf = 1, nbsrf
    370 CXXX
     395C
    371396        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
    372397        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
    373398        CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itau_w,
    374399     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    375 c       call writephy(fid_day,prof2d_av,'pourc_'//clnsurf(nsrf),
    376 c    .                pctsrf( 1 : klon, nsrf),
    377 c    .                'Fraction'//clnsurf(nsrf),'-')
    378400C
    379401        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
    … …  
    381403        CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itau_w,
    382404     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    383 c       call writephy(fid_day,prof2d_av,'tsol_'//clnsurf(nsrf),
    384 c    .                ftsol( 1 : klon, nsrf),
    385 c    .                'Surf. Temp'//clnsurf(nsrf),'K')
    386405C
    387406        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
    … …  
    389408        CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itau_w,
    390409     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    391 c       call writephy(fid_day,prof2d_av,'sens_'//clnsurf(nsrf),
    392 c    .                fluxt( 1 : klon, 1, nsrf),
    393 c    .                'Sensible heat flux '//clnsurf(nsrf),'W/m2')
    394410C 
    395411        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
    … …  
    397413        CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itau_w,
    398414     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    399 c       call writephy(fid_day,prof2d_av,'lat_'//clnsurf(nsrf),
    400 c    .                fluxlat( 1 : klon, nsrf),
    401 c    .                'Latent heat flux '//clnsurf(nsrf),'W/m2')
    402415C
    403416        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
    … …  
    405418        CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itau_w,
    406419     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    407 c       call writephy(fid_day,prof2d_av,'taux_'//clnsurf(nsrf),
    408 c    .                fluxu( 1 : klon, 1, nsrf),
    409 c    .                'Zonal wind stress '//clnsurf(nsrf),'Pa')
    410420C     
    411421        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
    … …  
    413423        CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itau_w,
    414424     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    415 c       call writephy(fid_day,prof2d_av,'tauy_'//clnsurf(nsrf),
    416 c    .                fluxv( 1 : klon, 1, nsrf),
    417 c    .                'Meridional wind stress '//clnsurf(nsrf),'Pa')
    418425C
    419426        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
    … …  
    421428        CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itau_w,
    422429     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    423 c       call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf),
    424 c    .                falbe( 1 : klon, nsrf),
    425 c    .                'Albedo surf. SW'//clnsurf(nsrf),'-')
    426 c       call writephy(fid_day,prof2d_av,'alblw_'//clnsurf(nsrf),
    427 c    .                falblw( 1 : klon, nsrf),
    428 c    .                'Albedo surf. LW'//clnsurf(nsrf),'-')
    429430C
    430431        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
    … …  
    432433        CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itau_w,
    433434     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    434 c       call writephy(fid_day,prof2d_av,'rugs_'//clnsurf(nsrf),
    435 c    .                frugs( 1 : klon, nsrf),
    436 c    .                'Rugosity '//clnsurf(nsrf),' - ')
    437435C
    438436      END DO 
    … …  
    440438c   FIN ECRITURE DES CHAMPS SUR LES SOUS SURFACES
    441439c=================================================================
    442 
    443 c=================================================================
    444 c   ECRITURE DES CHAMPS 3D
    445 c=================================================================
    446 
    447 c
    448 c Champs 3D:
    449 c
    450       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    451       CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
    452      .                                   iim*jjmp1*klev,ndex3d)
    453 c Essai writephys
    454       varname = 'temp'
    455       vartitle= 'temperature 3D'
    456       varunits= 'K'
    457 c     call writephy(fid_day,prof3d_av,varname,t_seri,vartitle,varunits)
    458 c
    459       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
    460       CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d,
    461      .                                   iim*jjmp1*klev,ndex3d)
    462 c     call writephy(fid_day,prof3d_av,'ovap',qx(1,1,ivap),
    463 c    .              'Specific humidity','Kg/Kg')
    464 c
    465       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    466       CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
    467      .                                   iim*jjmp1*klev,ndex3d)
    468 c     call writephy(fid_day,prof3d_av,'geop',zphi,
    469 c    .              'Geopotential height','m')
    470 c
    471       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    472       CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
    473      .                                   iim*jjmp1*klev,ndex3d)
    474 c     call writephy(fid_day,prof3d_av,'vitu',u_seri,
    475 c    .              'Zonal wind','m/s')
    476 c
    477       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    478       CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
    479      .                                   iim*jjmp1*klev,ndex3d)
    480 c     call writephy(fid_day,prof3d_av,'vitv',v_seri,
    481 c    .              'Meridional wind','m/s')
    482 c
    483       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    484       CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
    485      .                                   iim*jjmp1*klev,ndex3d)
    486 c     call writephy(fid_day,prof3d_av,'vitw',omega,
    487 c    .              'Vertical wind','m/s')
    488 c
    489       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    490       CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
    491      .                                   iim*jjmp1*klev,ndex3d)
    492 c     call writephy(fid_day,prof3d_av,'pres',pplay,
    493 c    .              'Air pressure','Pa')
    494 cccIM
    495 
    496 c=================================================================
    497 c   FIN ECRITURE DES CHAMPS 3D
    498 c=================================================================
    499 
    500 c=================================================================
    501 c=================================================================
    502       endif !  lev_histday.gt.1
    503 c=================================================================
    504 
    505 
     440      ENDIF !lev_histday.GE.4
     441c
     442c=================================================================
     443c=================================================================
     444c=================================================================
     445c
    506446      if (ok_sync) then
    507 c       call writephy_sync(fid_day)
    508447        call histsync(nid_day)
    509448      endif
Note: See TracChangeset for help on using the changeset viewer.