Ignore:
Timestamp:
Oct 19, 2001, 12:32:41 PM (23 years ago)
Author:
lmdzadmin
Message:

Passage des 2 albedos de surface vis et nir
mise en commentaire de la nouvelle surcouche IOIPSL
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

    r271 r282  
    244244      REAL falbe(klon,nbsrf)
    245245      SAVE falbe                  ! albedo par type de surface
     246      REAL falblw(klon,nbsrf)
     247      SAVE falblw                 ! albedo par type de surface
     248
    246249c
    247250c
     
    366369      REAL albsol(klon)
    367370      SAVE albsol                 ! albedo du sol total
     371      REAL albsollw(klon)
     372      SAVE albsollw                 ! albedo du sol total
     373
    368374      REAL wo(klon,klev)
    369375      SAVE wo                     ! ozone
     
    770776         zout = dtime * FLOAT(ecrit_day)
    771777C Essai writephys
    772          nom_fichier = 'histday1'
    773          call writephy_ini(fid_day,nom_fichier,klon,iim,jjmp1,klev,
    774      .                     rlon,rlat, presnivs,
    775      .                     zjulian, dtime)
    776          call writephy_def(prof2d_on, fid_day, "once", zsto, zout, 0)
    777          call writephy_def(prof3d_on, fid_day, "once", zsto, zout,
    778      .                                                         klev)
    779          call writephy_def(prof2d_av, fid_day, "ave(X)", zsto, zout, 0)
    780          call writephy_def(prof3d_av, fid_day, "ave(X)", zsto, zout,
    781      .                                                         klev)
     778c        nom_fichier = 'histday1'
     779c        call writephy_ini(fid_day,nom_fichier,klon,iim,jjmp1,klev,
     780c    .                     rlon,rlat, presnivs,
     781c    .                     zjulian, dtime)
     782c        call writephy_def(prof2d_on, fid_day, "once", zsto, zout, 0)
     783c        call writephy_def(prof3d_on, fid_day, "once", zsto, zout,
     784c    .                                                         klev)
     785c        call writephy_def(prof2d_av, fid_day, "ave(X)", zsto, zout, 0)
     786c        call writephy_def(prof3d_av, fid_day, "ave(X)", zsto, zout,
     787c    .                                                         klev)
    782788 
    783789c
     
    11621168     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    11631169     .                "ave(X)", zsto,zout)
     1170         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
     1171     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1172     .                "ave(X)", zsto,zout)
    11641173c
    11651174         CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
     
    15891598     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    15901599     .                "inst(X)", zsto,zout)
     1600         CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-",
     1601     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     1602     .                "inst(X)", zsto,zout)
    15911603c
    15921604c
     
    17951807     e            ok_veget, ocean, npas, nexca, ftsol,
    17961808     $            soil_model,ftsoil,
    1797      $            paprs,pplay,radsol, fsnow,fqsol,fevap,falbe,fluxlat,
     1809     $            paprs,pplay,radsol, fsnow,fqsol,fevap,falbe,falblw,
     1810     $            fluxlat,
    17981811     e            rain_fall, snow_fall, solsw, sollw, sollwdown, fder,
    17991812     e            rlon, rlat, cufi, cvfi, frugs,
     
    21732186     .             + falbe(i,is_ter) * pctsrf(i,is_ter)
    21742187     .             + falbe(i,is_sic) * pctsrf(i,is_sic)
     2188         albsollw(i) = falblw(i,is_oce) * pctsrf(i,is_oce)
     2189     .               + falblw(i,is_lic) * pctsrf(i,is_lic)
     2190     .               + falblw(i,is_ter) * pctsrf(i,is_ter)
     2191     .               + falblw(i,is_sic) * pctsrf(i,is_sic)
    21752192      ENDDO
    21762193      CALL radlwsw ! nouveau rayonnement (compatible Arpege-IFS)
    21772194     e            (dist, rmu0, fract, co2_ppm, solaire,
    2178      e             paprs, pplay,zxtsol,albsol, t_seri,q_seri,wo,
     2195     e             paprs, pplay,zxtsol,albsol, albsollw, t_seri,q_seri,
     2196     e             wo,
    21792197     e             cldfra, cldemi, cldtau,
    21802198     s             heat,heat0,cool,cool0,radsol,albpla,
     
    23742392         vartitle= 'Surface geop. height'
    23752393         varunits= '-'
    2376          call writephy(fid_day,prof2d_on,varname,pphis,vartitle,
    2377      .                                                    varunits)
     2394c        call writephy(fid_day,prof2d_on,varname,pphis,vartitle,
     2395c    .                                                    varunits)
    23782396c
    23792397         i = NINT(zout/zsto)
     
    23832401         vartitle= 'Grid area'
    23842402         varunits= '-'
    2385          call writephy(fid_day,prof2d_on,varname,paire,vartitle,
    2386      .                                                    varunits)
     2403c        call writephy(fid_day,prof2d_on,varname,paire,vartitle,
     2404c    .                                                    varunits)
    23872405C
    23882406      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    23892407      CALL histwrite(nid_day,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2390       call writephy(fid_day,prof2d_av,'tsol',zxtsol,
    2391      .              'Surface Temperature','K')
     2408c     call writephy(fid_day,prof2d_av,'tsol',zxtsol,
     2409c    .              'Surface Temperature','K')
    23922410c
    23932411C
     
    23952413      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
    23962414      CALL histwrite(nid_day,"tter",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2397       call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
    2398      .              'Surface Temperature','K')
     2415c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
     2416c    .              'Surface Temperature','K')
    23992417C
    24002418      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
    24012419      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    24022420      CALL histwrite(nid_day,"tlic",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2403       call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
    2404      .              'Surface Temperature','K')
     2421c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
     2422c    .              'Surface Temperature','K')
    24052423C
    24062424      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
    24072425      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    24082426      CALL histwrite(nid_day,"toce",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2409       call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
    2410      .              'Surface Temperature','K')
     2427c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
     2428c    .              'Surface Temperature','K')
    24112429C
    24122430      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
    24132431      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
    24142432      CALL histwrite(nid_day,"tsic",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2415       call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
    2416      .              'Surface Temperature','K')
     2433c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
     2434c    .              'Surface Temperature','K')
    24172435C
    24182436      DO i = 1, klon
     
    24252443      vartitle= 'pression au sol'
    24262444      varunits= 'hPa'
    2427       call writephy(fid_day,prof2d_av,varname,zx_tmp_fi2d,vartitle,
    2428      .                                                    varunits)
     2445c     call writephy(fid_day,prof2d_av,varname,zx_tmp_fi2d,vartitle,
     2446c    .                                                    varunits)
    24292447c
    24302448      DO i = 1, klon
     
    24332451      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    24342452      CALL histwrite(nid_day,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2435       call writephy(fid_day,prof2d_av,'rain',zx_tmp_fi2d,
    2436      .              'Precipitation','mm/day')
     2453c     call writephy(fid_day,prof2d_av,'rain',zx_tmp_fi2d,
     2454c    .              'Precipitation','mm/day')
    24372455
    24382456
     
    24402458      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
    24412459      CALL histwrite(nid_day,"snow",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2442       call writephy(fid_day,prof2d_av,'snow',snow_fall,
    2443      .              'Snow','mm/day')
     2460c     call writephy(fid_day,prof2d_av,'snow',snow_fall,
     2461c    .              'Snow','mm/day')
    24442462c
    24452463      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
    24462464      CALL histwrite(nid_day,"snow_cov",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2447       call writephy(fid_day,prof2d_av,'snow_cov',zxsnow,
    2448      .              'Snow cover','mm')
     2465c     call writephy(fid_day,prof2d_av,'snow_cov',zxsnow,
     2466c    .              'Snow cover','mm')
    24492467c
    24502468      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
    24512469      CALL histwrite(nid_day,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2452       call writephy(fid_day,prof2d_av,'evap',evap,
    2453      .              'Evaporation','mm/day')
     2470c     call writephy(fid_day,prof2d_av,'evap',evap,
     2471c    .              'Evaporation','mm/day')
    24542472c
    24552473      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
    24562474      CALL histwrite(nid_day,"tops",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2457       call writephy(fid_day,prof2d_av,'tops',topsw,
    2458      .              'Solar rad. at TOA','W/m2')
     2475c     call writephy(fid_day,prof2d_av,'tops',topsw,
     2476c    .              'Solar rad. at TOA','W/m2')
    24592477c
    24602478      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    24612479      CALL histwrite(nid_day,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2462       call writephy(fid_day,prof2d_av,'topl',toplw,
    2463      .              'IR rad. at TOA','W/m2')
     2480c     call writephy(fid_day,prof2d_av,'topl',toplw,
     2481c    .              'IR rad. at TOA','W/m2')
    24642482c
    24652483      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
    24662484      CALL histwrite(nid_day,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2467       call writephy(fid_day,prof2d_av,'sols',solsw,
    2468      .              'Solar rad. at surf.','W/m2')
     2485c     call writephy(fid_day,prof2d_av,'sols',solsw,
     2486c    .              'Solar rad. at surf.','W/m2')
    24692487c
    24702488      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
    24712489      CALL histwrite(nid_day,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2472       call writephy(fid_day,prof2d_av,'soll',sollw,
    2473      .              'IR rad. at surface','W/m2')
     2490c     call writephy(fid_day,prof2d_av,'soll',sollw,
     2491c    .              'IR rad. at surface','W/m2')
    24742492c
    24752493      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
    24762494      CALL histwrite(nid_day,"solldown",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2477       call writephy(fid_day,prof2d_av,'solldown',sollwdown,
    2478      .              'Down. IR rad. at surface','W/m2')
     2495c     call writephy(fid_day,prof2d_av,'solldown',sollwdown,
     2496c    .              'Down. IR rad. at surface','W/m2')
    24792497c
    24802498      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
    24812499      CALL histwrite(nid_day,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2482       call writephy(fid_day,prof2d_av,'bils',bils,
    2483      .              'Surf. total heat flux','W/m2')
     2500c     call writephy(fid_day,prof2d_av,'bils',bils,
     2501c    .              'Surf. total heat flux','W/m2')
    24842502c
    24852503      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
    24862504      CALL histwrite(nid_day,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2487       call writephy(fid_day,prof2d_av,'sens',sens,
    2488      .              'Sensible heat flux','W/m2')
     2505c     call writephy(fid_day,prof2d_av,'sens',sens,
     2506c    .              'Sensible heat flux','W/m2')
    24892507c
    24902508      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
    24912509      CALL histwrite(nid_day,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2492       call writephy(fid_day,prof2d_av,'fder',fder,
    2493      .              'Heat flux derivation','W/m2')
     2510c     call writephy(fid_day,prof2d_av,'fder',fder,
     2511c    .              'Heat flux derivation','W/m2')
    24942512c
    24952513c
     
    25002518        CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itap,
    25012519     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    2502         call writephy(fid_day,prof2d_av,'pourc_'//clnsurf(nsrf),
    2503      .                pctsrf( 1 : klon, nsrf),
    2504      .                'Fraction'//clnsurf(nsrf),'-')
     2520c       call writephy(fid_day,prof2d_av,'pourc_'//clnsurf(nsrf),
     2521c    .                pctsrf( 1 : klon, nsrf),
     2522c    .                'Fraction'//clnsurf(nsrf),'-')
    25052523C
    25062524        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
     
    25082526        CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itap,
    25092527     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    2510         call writephy(fid_day,prof2d_av,'tsol_'//clnsurf(nsrf),
    2511      .                ftsol( 1 : klon, nsrf),
    2512      .                'Surf. Temp'//clnsurf(nsrf),'K')
     2528c       call writephy(fid_day,prof2d_av,'tsol_'//clnsurf(nsrf),
     2529c    .                ftsol( 1 : klon, nsrf),
     2530c    .                'Surf. Temp'//clnsurf(nsrf),'K')
    25132531C
    25142532        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
     
    25162534        CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itap,
    25172535     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    2518         call writephy(fid_day,prof2d_av,'sens_'//clnsurf(nsrf),
    2519      .                fluxt( 1 : klon, 1, nsrf),
    2520      .                'Sensible heat flux '//clnsurf(nsrf),'W/m2')
     2536c       call writephy(fid_day,prof2d_av,'sens_'//clnsurf(nsrf),
     2537c    .                fluxt( 1 : klon, 1, nsrf),
     2538c    .                'Sensible heat flux '//clnsurf(nsrf),'W/m2')
    25212539
    25222540        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
     
    25242542        CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itap,
    25252543     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    2526         call writephy(fid_day,prof2d_av,'lat_'//clnsurf(nsrf),
    2527      .                fluxlat( 1 : klon, nsrf),
    2528      .                'Latent heat flux '//clnsurf(nsrf),'W/m2')
     2544c       call writephy(fid_day,prof2d_av,'lat_'//clnsurf(nsrf),
     2545c    .                fluxlat( 1 : klon, nsrf),
     2546c    .                'Latent heat flux '//clnsurf(nsrf),'W/m2')
    25292547C
    25302548        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
     
    25322550        CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itap,
    25332551     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    2534         call writephy(fid_day,prof2d_av,'taux_'//clnsurf(nsrf),
    2535      .                fluxu( 1 : klon, 1, nsrf),
    2536      .                'Zonal wind stress '//clnsurf(nsrf),'Pa')
     2552c       call writephy(fid_day,prof2d_av,'taux_'//clnsurf(nsrf),
     2553c    .                fluxu( 1 : klon, 1, nsrf),
     2554c    .                'Zonal wind stress '//clnsurf(nsrf),'Pa')
    25372555C     
    25382556        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
     
    25402558        CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itap,
    25412559     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    2542         call writephy(fid_day,prof2d_av,'tauy_'//clnsurf(nsrf),
    2543      .                fluxv( 1 : klon, 1, nsrf),
    2544      .                'Meridional wind stress '//clnsurf(nsrf),'Pa')
     2560c       call writephy(fid_day,prof2d_av,'tauy_'//clnsurf(nsrf),
     2561c    .                fluxv( 1 : klon, 1, nsrf),
     2562c    .                'Meridional wind stress '//clnsurf(nsrf),'Pa')
    25452563C
    25462564        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
     
    25482566        CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itap,
    25492567     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    2550         call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf),
    2551      .                falbe( 1 : klon, nsrf),
    2552      .                'Albedo surf.'//clnsurf(nsrf),'-')
     2568c       call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf),
     2569c    .                falbe( 1 : klon, nsrf),
     2570c    .                'Albedo surf. SW'//clnsurf(nsrf),'-')
     2571c       call writephy(fid_day,prof2d_av,'alblw_'//clnsurf(nsrf),
     2572c    .                falblw( 1 : klon, nsrf),
     2573c    .                'Albedo surf. LW'//clnsurf(nsrf),'-')
    25532574C
    25542575        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
     
    25562577        CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itap,
    25572578     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    2558         call writephy(fid_day,prof2d_av,'rugs_'//clnsurf(nsrf),
    2559      .                frugs( 1 : klon, nsrf),
    2560      .                'Rugosity '//clnsurf(nsrf),' - ')
     2579c       call writephy(fid_day,prof2d_av,'rugs_'//clnsurf(nsrf),
     2580c    .                frugs( 1 : klon, nsrf),
     2581c    .                'Rugosity '//clnsurf(nsrf),' - ')
    25612582C
    25622583      END DO 
     
    25702591      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
    25712592      CALL histwrite(nid_day,"cldl",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2572       call writephy(fid_day,prof2d_av,'cldl',cldl,
    2573      .              'Low-level cloudiness','-')
     2593c     call writephy(fid_day,prof2d_av,'cldl',cldl,
     2594c    .              'Low-level cloudiness','-')
    25742595c
    25752596      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
    25762597      CALL histwrite(nid_day,"cldm",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2577       call writephy(fid_day,prof2d_av,'cldm',cldm,
    2578      .              'Mid-level cloudiness','-')
     2598c     call writephy(fid_day,prof2d_av,'cldm',cldm,
     2599c    .              'Mid-level cloudiness','-')
    25792600c
    25802601      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
    25812602      CALL histwrite(nid_day,"cldh",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2582       call writephy(fid_day,prof2d_av,'cldh',cldh,
    2583      .              'High-level cloudiness','-')
     2603c     call writephy(fid_day,prof2d_av,'cldh',cldh,
     2604c    .              'High-level cloudiness','-')
    25842605c
    25852606      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    25862607      CALL histwrite(nid_day,"cldt",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2587       call writephy(fid_day,prof2d_av,'cldt',cldt,
    2588      .              'Total cloudiness','-')
     2608c     call writephy(fid_day,prof2d_av,'cldt',cldt,
     2609c    .              'Total cloudiness','-')
    25892610c
    25902611      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
    25912612      CALL histwrite(nid_day,"cldq",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    2592       call writephy(fid_day,prof2d_av,'cldq',cldq,
    2593      .              'Cloud liquid water path','-')
     2613c     call writephy(fid_day,prof2d_av,'cldq',cldq,
     2614c    .              'Cloud liquid water path','-')
    25942615c
    25952616c Champs 3D:
     
    26022623      vartitle= 'temperature 3D'
    26032624      varunits= 'K'
    2604       call writephy(fid_day,prof3d_av,varname,t_seri,vartitle,varunits)
     2625c     call writephy(fid_day,prof3d_av,varname,t_seri,vartitle,varunits)
    26052626c
    26062627      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
    26072628      CALL histwrite(nid_day,"ovap",itap,zx_tmp_3d,
    26082629     .                                   iim*jjmp1*klev,ndex3d)
    2609       call writephy(fid_day,prof3d_av,'ovap',qx(1,1,ivap),
    2610      .              'Specific humidity','Kg/Kg')
     2630c     call writephy(fid_day,prof3d_av,'ovap',qx(1,1,ivap),
     2631c    .              'Specific humidity','Kg/Kg')
    26112632c
    26122633      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    26132634      CALL histwrite(nid_day,"geop",itap,zx_tmp_3d,
    26142635     .                                   iim*jjmp1*klev,ndex3d)
    2615       call writephy(fid_day,prof3d_av,'geop',zphi,
    2616      .              'Geopotential height','m')
     2636c     call writephy(fid_day,prof3d_av,'geop',zphi,
     2637c    .              'Geopotential height','m')
    26172638c
    26182639      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    26192640      CALL histwrite(nid_day,"vitu",itap,zx_tmp_3d,
    26202641     .                                   iim*jjmp1*klev,ndex3d)
    2621       call writephy(fid_day,prof3d_av,'vitu',u_seri,
    2622      .              'Zonal wind','m/s')
     2642c     call writephy(fid_day,prof3d_av,'vitu',u_seri,
     2643c    .              'Zonal wind','m/s')
    26232644c
    26242645      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    26252646      CALL histwrite(nid_day,"vitv",itap,zx_tmp_3d,
    26262647     .                                   iim*jjmp1*klev,ndex3d)
    2627       call writephy(fid_day,prof3d_av,'vitv',v_seri,
    2628      .              'Meridional wind','m/s')
     2648c     call writephy(fid_day,prof3d_av,'vitv',v_seri,
     2649c    .              'Meridional wind','m/s')
    26292650c
    26302651      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    26312652      CALL histwrite(nid_day,"vitw",itap,zx_tmp_3d,
    26322653     .                                   iim*jjmp1*klev,ndex3d)
    2633       call writephy(fid_day,prof3d_av,'vitw',omega,
    2634      .              'Vertical wind','m/s')
     2654c     call writephy(fid_day,prof3d_av,'vitw',omega,
     2655c    .              'Vertical wind','m/s')
    26352656c
    26362657      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    26372658      CALL histwrite(nid_day,"pres",itap,zx_tmp_3d,
    26382659     .                                   iim*jjmp1*klev,ndex3d)
    2639       call writephy(fid_day,prof3d_av,'pres',pplay,
    2640      .              'Air pressure','Pa')
     2660c     call writephy(fid_day,prof3d_av,'pres',pplay,
     2661c    .              'Air pressure','Pa')
    26412662
    26422663c
    26432664      if (ok_sync) then
    2644         call writephy_sync(fid_day)
     2665c       call writephy_sync(fid_day)
    26452666        call histsync(nid_day)
    26462667      endif
     
    28092830      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
    28102831      CALL histwrite(nid_mth,"albs",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     2832      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
     2833      CALL histwrite(nid_mth,"albslw",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    28112834c
    28122835      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
     
    31793202      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
    31803203      CALL histwrite(nid_ins,"albs",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     3204      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
     3205      CALL histwrite(nid_ins,"albslw",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
    31813206c
    31823207      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
Note: See TracChangeset for help on using the changeset viewer.