Changeset 433 for LMDZ.3.3


Ignore:
Timestamp:
Dec 19, 2002, 5:46:39 PM (21 years ago)
Author:
lmdzadmin
Message:

Convergence avec la version de Ionela dec 2002

YOMCST.? : suppression RI0 (IM)
albedo.F : facteur 1.2 sur le nouveau calcul (IM)
clesphys.h : rajout de différentes ctes (concentration des gaz) (IM)
clmain.F : separation des flux LW, SW (JLD)

remplace qsurf par yqsol (IM)

conf_phys.F90 : rajout de différentes ctes (gaz + orbite) (IM)
convect3.F : DPINV+SIGD*0.5*(EVAP(1)+EVAP(2)) (SBL)
cv3_routines.F:
cvparam3.h : compatibilite avec conema3 TEMPORAIRE (FH)
phyetat0.F : lecture de co2_ppm et solaire pour tests de coherence
phyredem.F : co2_ppm et solaire passé en common
physiq.F : separation flux LW, SW

rajout diagnostiques (slp, w500)
suppression iflag_con = 4
clwcon0=qcondc (FH)
position dU "ENDIF ! ok_cvl"

radlwsw.F : passage des concentrations gaz dans un common (IM)

PEMIS(i) = 1.0 (JLD pour cohérence ORCHIDEE)

stdlevvar.F90 :
suphec.F : suppression init. des ctes orbitales (IM)

nouvelles E/S (ini_hist..., write_hist...)

Location:
LMDZ.3.3/branches/rel-LF/libf/phylmd
Files:
19 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/YOMCST.h

    r98 r433  
    88      REAL RA,RG,R1SA
    99! A1.3 Radiation
    10       REAL RSIGMA,RI0
     10!     REAL RSIGMA,RI0
     11      REAL RSIGMA
    1112! A1.4 Thermodynamic gas phase
    1213      REAL R,RMD,RMV,RD,RV,RCPD,RCPV,RCVD,RCVV
    … …  
    2021      REAL RALPD,RBETD,RGAMD
    2122!
     23!    S      ,RSIGMA,RI0
    2224      COMMON/YOMCST/RPI   ,RCLUM ,RHPLA ,RKBOL ,RNAVO
    2325     S      ,RDAY  ,REA   ,REPSM ,RSIYEA,RSIDAY,ROMEGA
    2426     s      ,R_ecc, R_peri, R_incl
    2527     S      ,RA    ,RG    ,R1SA
    26      S      ,RSIGMA,RI0
     28     S      ,RSIGMA
    2729     S      ,R     ,RMD   ,RMV   ,RD    ,RV    ,RCPD
    2830     S      ,RCPV  ,RCVD  ,RCVV  ,RKAPPA,RETV

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/YOMCST.inc

    r230 r433  
    88      REAL :: RA,RG,R1SA
    99! A1.3 Radiation
    10       REAL :: RSIGMA,RI0
     10!     REAL :: RSIGMA,RI0
     11      REAL :: RSIGMA
    1112! A1.4 Thermodynamic gas phase
    1213      REAL :: R,RMD,RMV,RD,RV,RCPD,RCPV,RCVD,RCVV
    … …  
    2021      REAL :: RALPD,RBETD,RGAMD
    2122!
     23!IM  & ,R1SA ,RSIGMA,RI0,R ,RMD   ,RMV   ,RD    ,RV    ,RCPD ,RCPV,RCVD &
    2224      COMMON/YOMCST/RPI ,RCLUM, RHPLA, RKBOL, RNAVO ,RDAY  ,REA &
    2325     & ,REPSM ,RSIYEA,RSIDAY,ROMEGA ,R_ecc, R_peri, R_incl, RA    ,RG &
    24      & ,R1SA ,RSIGMA,RI0,R ,RMD   ,RMV   ,RD    ,RV    ,RCPD ,RCPV,RCVD &
     26     & ,R1SA ,RSIGMA,R ,RMD   ,RMV   ,RD    ,RV    ,RCPD ,RCPV,RCVD &
    2527     & ,RCVV  ,RKAPPA,RETV ,RCW   ,RCS ,RLVTT ,RLSTT ,RLMLT ,RTT ,RATM &
    2628     & ,RESTT ,RALPW ,RBETW ,RGAMW ,RALPS ,RBETS ,RGAMS ,RALPD ,RBETD ,RGAMD

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/albedo.F

    r406 r433  
     1c
     2c $Header$
     3c
    14      SUBROUTINE alboc(rjour,rlat,albedo)
    25      IMPLICIT none
    … …  
    1922ccc      PARAMETER (fmagic=0.7)
    2023cccIM => a remplacer 
    21 cccIM   PARAMETER (fmagic=1.1)
     24        PARAMETER (fmagic=1.1)
    2225c       PARAMETER (fmagic=1.0)
    23         PARAMETER (fmagic=0.7)
     26c       PARAMETER (fmagic=0.7)
    2427      INTEGER npts ! il controle la precision de l'integration
    2528      PARAMETER (npts=120) ! 120 correspond a l'interval 6 minutes
    … …  
    9295      rmu = aa + bb * COS(0.0)
    9396      rmu = MAX(0.0, rmu)
    94       alb = 0.058/(rmu + 0.30)
     97cIM cf. PB  alb = 0.058/(rmu + 0.30)
     98c     alb = 0.058/(rmu + 0.30) * 1.5
     99      alb = 0.058/(rmu + 0.30) * 1.2
    95100      srmu = rmu
    96101      salb = alb * rmu
    … …  
    101106         rmu = aa + bb * COS(FLOAT(k)/FLOAT(npts)*zpi)
    102107         rmu = MAX(0.0, rmu)
    103          alb = 0.058/(rmu + 0.30)
     108cIM cf. PB      alb = 0.058/(rmu + 0.30)
     109c        alb = 0.058/(rmu + 0.30) * 1.5
     110         alb = 0.058/(rmu + 0.30) * 1.2
    104111         srmu = srmu + rmu * 2.0
    105112         salb = salb + alb*rmu * 2.0
    … …  
    136143ccc      PARAMETER (fmagic=0.7)
    137144cccIM => a remplacer 
    138 cccIM   PARAMETER (fmagic=1.1)
     145        PARAMETER (fmagic=1.1)
    139146c       PARAMETER (fmagic=1.0)
    140         PARAMETER (fmagic=0.7)
     147c       PARAMETER (fmagic=0.7)
    141148c
    142149      REAL fauxo
    … …  
    164171      DO i = 1, klon
    165172         rmu0(i) = MAX(rmu0(i),0.0)
    166          albedo(i) = 0.058/(rmu0(i) + 0.30)
     173cIM cf. PB      albedo(i) = 0.058/(rmu0(i) + 0.30)
     174c        albedo(i) = 0.058/(rmu0(i) + 0.30) * 1.5
     175         albedo(i) = 0.058/(rmu0(i) + 0.30) * 1.2
    167176         albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
    168177      ENDDO

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/clesphys.h

    r2 r433  
    44       LOGICAL ok_limitvrai
    55       INTEGER nbapp_rad, iflag_con
     6       REAL co2_ppm, solaire
     7       REAL*8 RCO2, RCH4, RN2O, RCFC11, RCFC12 
    68
    7        COMMON/clesphys/cycle_diurne, soil_model, new_oliq ,
     9       COMMON/clesphys/cycle_diurne, soil_model, new_oliq,
    810     ,     ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad, iflag_con
    9 c
     11     ,     , co2_ppm, solaire, RCO2, RCH4, RN2O, RCFC11, RCFC12

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F

    r419 r433  
    119119      real rain_f(klon), snow_f(klon)
    120120      REAL fder(klon)
    121       REAL sollw(klon), solsw(klon), sollwdown(klon)
     121cIM cf. JLD   REAL sollw(klon), solsw(klon), sollwdown(klon)
     122      REAL sollw(klon,nbsrf), solsw(klon,nbsrf), sollwdown(klon)
    122123      REAL rugos(klon,nbsrf)
    123124C la nouvelle repartition des surfaces sortie de l'interface
    … …  
    407408        ytaux(j) = flux_u(i,1,nsrf)
    408409        ytauy(j) = flux_v(i,1,nsrf)
    409 c$$$        ysolsw(j) = solsw(i)
    410        ysolsw(j) = (1 - albe(i,nsrf))
    411      $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter)
    412      $      - pctsrf(i, is_lic) *albe(i,is_lic)
    413      $      - pctsrf(i, is_oce) *albe(i,is_oce)
    414      $      - pctsrf(i, is_sic) *albe(i,is_sic)
    415      $      ) * solsw(i)
    416         ysollw(j) = sollw(i)
     410c$$$    ysolsw(j) = solsw(i)
     411cIM cf. JLD
     412        ysolsw(j) = solsw(i,nsrf)
     413c       ysolsw(j) = (1 - albe(i,nsrf))
     414c    $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter)
     415c    $      - pctsrf(i, is_lic) *albe(i,is_lic)
     416c    $      - pctsrf(i, is_oce) *albe(i,is_oce)
     417c    $      - pctsrf(i, is_sic) *albe(i,is_sic)
     418c    $      ) * solsw(i)
     419cIM cf. JLD     ysollw(j) = sollw(i)
     420        ysollw(j) = sollw(i,nsrf)
    417421        ysollwdown(j) = sollwdown(i)
    418422        yrugos(j) = rugos(i,nsrf)
    … …  
    421425        yv1(j) = v1lay(i)
    422426c$$$        yrads(j) = totalflu(i)
    423         yrads(j) = (1 - albe(i,nsrf))
    424      $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter)
    425      $      - pctsrf(i, is_lic) *albe(i,is_lic)
    426      $      - pctsrf(i, is_oce) *albe(i,is_oce)
    427      $      - pctsrf(i, is_sic) *albe(i,is_sic)
    428      $      ) * solsw(i) + sollw(i)
     427cIM cf. JLD
     428        yrads(j) =  ysolsw(j)+ ysollw(j)
     429c       yrads(j) = (1 - albe(i,nsrf))
     430c    $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter)
     431c    $      - pctsrf(i, is_lic) *albe(i,is_lic)
     432c    $      - pctsrf(i, is_oce) *albe(i,is_oce)
     433c    $      - pctsrf(i, is_sic) *albe(i,is_sic)
     434cIM cf. JLD  $      ) * solsw(i) + sollw(i)
     435c    $      ) * solsw(i) + ysollw(j)
    429436        ypaprs(j,klev+1) = paprs(i,klev+1)
    430437      END DO
    … …  
    453460      CALL coefkz(nsrf, knon, ypaprs, ypplay,
    454461     .            yts, yrugos, yu, yv, yt, yq,
     462cIM remplace qsurf par yqsol
     463     .            yqsol,
    455464     .            ycoefm, ycoefh)
    456465      CALL coefkz2(nsrf, knon, ypaprs, ypplay,yt,
    … …  
    595604c
    596605#undef T2m     
     606#define T2m     
    597607#ifdef T2m
    598608ccc diagnostic t,q a 2m et u, v a 10m
    … …  
    11361146     .                  ts, rugos,
    11371147     .                  u,v,t,q,
     1148cIM remplace qsurf par yqsol
     1149     .                  qsol,
    11381150     .                  pcfm, pcfh)
    11391151      IMPLICIT none
    … …  
    12261238      REAL gamt(2:klev)
    12271239c essai qsurf
    1228       real qsurf(klon)
     1240cIM   real qsurf(klon)
     1241      real qsol(klon)
    12291242c
    12301243      LOGICAL appel1er
    … …  
    12641277      ENDDO
    12651278
     1279cIM remplace qsurf par qsol
     1280      IF(nsrf.NE.1) THEN
    12661281      do i = 1, knon
    1267         qsurf(i) = qsatl(ts(i))/paprs(i,1)
     1282cIM     qsurf(i) = qsatl(ts(i))/paprs(i,1)
     1283        qsol(i) = qsatl(ts(i))/paprs(i,1)
    12681284      enddo
     1285      ENDIF
    12691286
    12701287c
    … …  
    12981315c Calculer le frottement au sol (Cdrag)
    12991316c
    1300 c     DO i = 1, knon
    1301       DO i = 1, klon
     1317      DO i = 1, knon
    13021318       u1(i) = u(i,1)
    13031319       v1(i) = v(i,1)
    … …  
    13091325      CALL clcdrag(klon, knon, nsrf, zxli,
    13101326     $             u1, v1, t1, q1, z1,
    1311      $             ts, qsurf, rugos,
     1327     $             ts, qsol, rugos,
    13121328     $             pcfm1, pcfh1)
     1329cIM  $             ts, qsurf, rugos,
    13131330C
    13141331      DO i = 1, knon

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/conf_phys.F90

    r395 r433  
    1313#include "fisrtilp.inc"
    1414#include "nuage.h"
     15#include "YOMCST.inc"
     16!IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
     17include "clesphys.inc"
    1518!
    1619! Configuration de la "physique" de LMDZ a l'aide de la fonction
    … …  
    9598  call getin('if_ebil', if_ebil)
    9699!!
     100!! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
     101!!
     102!Config Key  = R_ecc
     103!Config Desc = Excentricite
     104!Config Def  = 0.016715
     105!Config Help =
     106!               
     107!valeur AMIP II
     108  R_ecc = 0.016715
     109  call getin('R_ecc', R_ecc)
     110!!
     111!Config Key  = R_peri
     112!Config Desc = Equinoxe
     113!Config Def  =
     114!Config Help =
     115!               
     116!
     117!valeur AMIP II
     118  R_peri = 102.7
     119  call getin('R_peri', R_peri)
     120!!
     121!Config Key  = R_incl
     122!Config Desc = Inclinaison
     123!Config Def  =
     124!Config Help =
     125!               
     126!
     127!valeur AMIP II
     128  R_incl = 23.441
     129  call getin('R_incl', R_incl)
     130!!
     131!Config Key  = solaire
     132!Config Desc = Constante solaire en W/m2
     133!Config Def  = 1365.
     134!Config Help =
     135!               
     136!
     137!valeur AMIP II
     138  solaire = 1365.
     139  call getin('solaire', solaire)
     140!!
     141!Config Key  = co2_ppm
     142!Config Desc = concentration du gaz carbonique en ppmv
     143!Config Def  = 348.
     144!Config Help =
     145!               
     146!
     147!valeur AMIP II
     148  co2_ppm = 348.
     149  call getin('co2_ppm', co2_ppm)
     150!!
     151!Config Key  = RCO2
     152!Config Desc = Concentration du CO2
     153!Config Def  = co2_ppm * 1.0e-06  * 44.011/28.97
     154!Config Def  = 348. * 1.0e-06  * 44.011/28.97
     155!Config Help =
     156!               
     157! RCO2 = 5.286789092164308E-04
     158!ancienne valeur
     159  RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97 ! pour co2_ppm=348.
     160
     161  call getin('RCO2', RCO2)
     162!!
     163!Config Key  = RCH4
     164!Config Desc = Concentration du CH4
     165!Config Def  = 1.65E-06* 16.043/28.97
     166!Config Help =
     167!               
     168!
     169!valeur AMIP II
     170  RCH4 = 1.65E-06* 16.043/28.97
     171! RCH4 = 9.137366240938903E-07
     172!
     173!ancienne valeur
     174! RCH4 = 1.72E-06* 16.043/28.97
     175  call getin('RCH4', RCH4)
     176!!
     177!Config Key  = RN2O
     178!Config Desc = Concentration du N2O
     179!Config Def  = 306.E-09* 44.013/28.97
     180!Config Help =
     181!               
     182!
     183!valeur AMIP II
     184  RN2O = 306.E-09* 44.013/28.97
     185! RN2O = 4.648939592682085E-07
     186!
     187!ancienne valeur
     188! RN2O = 310.E-09* 44.013/28.97
     189  call getin('RN2O', RN2O)
     190!!
     191!Config Key  = RCFC11
     192!Config Desc = Concentration du CFC11
     193!Config Def  = 280.E-12* 137.3686/28.97
     194!Config Help =
     195!               
     196!
     197  RCFC11 = 280.E-12* 137.3686/28.97
     198! RCFC11 = 1.327690990680013E-09
     199  call getin('RCFC11', RCFC11)
     200!!
     201!Config Key  = RCFC12
     202!Config Desc = Concentration du CFC12
     203!Config Def  = 484.E-12* 120.9140/28.97
     204!Config Help =
     205!               
     206!
     207  RCFC12 = 484.E-12* 120.9140/28.97
     208! RCFC12 = 2.020102726958923E-09
     209  call getin('RCFC12', RCFC12)
     210!!
     211!! Constante solaire & Parametres orbitaux & taux gaz effet de serre END
     212!!
    97213!! KE
    98214!
    … …  
    292408  write(numout,*)' Sortie instantanee = ', ok_instan
    293409  write(numout,*)' Sortie bilan d''energie, if_ebil =', if_ebil
     410!IM constantes physiques BEG
     411  write(numout,*)' Excentricite = ',R_ecc
     412  write(numout,*)' Equinoxe = ',R_peri
     413  write(numout,*)' Inclinaison =',R_incl
     414  write(numout,*)' Constante solaire =',solaire
     415  write(numout,*)' co2_ppm =',co2_ppm
     416  write(numout,*)' RCO2 = ',RCO2
     417  write(numout,*)' RCH4 = ',RCH4
     418  write(numout,*)' RN2O =  ',RN2O
     419  write(numout,*)' RCFC11 =  ',RCFC11
     420  write(numout,*)' RCFC12 =  ',RCFC12
     421!IM constantes physiques END
    294422  write(numout,*)' epmax = ', epmax
    295423  write(numout,*)' ok_adj_ema = ', ok_adj_ema

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/convect3.F

    r408 r433  
    10421042        FR(1)=0.1*MP(2)*(RP(2)-RR(1))*
    10431043Ccorrection bug conservation eau
    1044 C     1    DPINV+SIGD*0.5*(EVAP(1)+EVAP(2))
    1045      1    DPINV+SIGD*EVAP(1)
     1044C    1    DPINV+SIGD*0.5*(EVAP(1)+EVAP(2))
     1045     1    DPINV+SIGD*0.5*(EVAP(1)+EVAP(2))
     1046cIM cf. SBL
     1047C    1    DPINV+SIGD*EVAP(1)
    10461048        FR(1)=FR(1)+0.1*AM*(RR(2)-RR(1))*DPINV
    10471049        FU(1)=FU(1)+0.1*DPINV*(MP(2)*(UP(2)-U(1))+AM*(U(2)-U(1)))

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/cv3_routines.F

    r418 r433  
    2626
    2727#include "cvparam3.h"
     28#include "conema3.h"
    2829
    2930      integer nd
    … …  
    4748      pbcrit = 150.0
    4849      ptcrit = 500.0
    49       epmax  = 0.993
     50cIM cf. FH     epmax  = 0.993
    5051
    5152      omtrain = 45.0 ! used also for snow (no disctinction rain/snow)
    … …  
    864865#include "cvthermo.h"
    865866#include "cvparam3.h"
     867#include "conema3.h"
    866868
    867869c inputs:
    … …  
    21282130#include "cvparam3.h"
    21292131#include "cvflag.h"
     2132#include "conema3.h"
    21302133
    21312134c inputs:
    … …  
    29362939     :         *rrd*tvp(il,i)/p(il,i)/100./delta            ! cld
    29372940          siga(il,i) = min(siga(il,i),1.0)                  ! cld
     2941cIM cf. FH
     2942         if (iflag_clw.eq.0) then
    29382943          qcondc(il,i)=siga(il,i)*clw(il,i)*(1.-ep(il,i))   ! cld
    29392944     :           + (1.-siga(il,i))*qcond(il,i)              ! cld
     2945         else if (iflag_clw.eq.1) then
     2946          qcondc(il,i)=qcond(il,i)              ! cld
     2947         endif
     2948
    29402949        enddo                                               ! cld
    29412950       enddo                                                ! cld

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/cvparam3.h

    r418 r433  
    99      integer noff, minorig, nl, nlp, nlm
    1010      real sigd, spfac
    11       real pbcrit, ptcrit, epmax
     11cIM cf. FH : pour compatibilite avec conema3 TEMPORAIRE   real pbcrit, ptcrit, epmax
     12      real pbcrit, ptcrit
    1213      real omtrain
    1314      real dtovsh, dpbase, dttrig
    … …  
    1819      COMMON /cvparam3/  noff, minorig, nl, nlp, nlm
    1920     :                ,  sigd, spfac
    20      :                ,pbcrit, ptcrit, epmax
     21cIM cf. FH : pour compatibilite avec conema3 TEMPORAIRE  :                ,pbcrit, ptcrit, epmax
     22     :                ,pbcrit, ptcrit
    2123     :                ,omtrain
    2224     :                ,dtovsh, dpbase, dttrig

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histday.h

    r418 r433  
    5353     .                "ave(X)", zsto,zout)
    5454c
    55          CALL histdef(nid_day, "tter", "Surface Temperature", "K",
    56      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    57      .                "ave(X)", zsto,zout)
    58 c
    59          CALL histdef(nid_day, "tlic", "Surface Temperature", "K",
    60      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    61      .                "ave(X)", zsto,zout)
    62 c
    63          CALL histdef(nid_day, "toce", "Surface Temperature", "K",
    64      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    65      .                "ave(X)", zsto,zout)
    66 c
    67          CALL histdef(nid_day, "tsic", "Surface Temperature", "K",
    68      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    69      .                "ave(X)", zsto,zout)
    70 c
    71 cccIM
    7255c
    7356         CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
    … …  
    7558     .                "ave(X)", zsto,zout)
    7659c
    77          CALL histdef(nid_day, "t2mter", "Temp.terre 2m", "K",
    78      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    79      .                "ave(X)", zsto,zout)
    80 c
    81          CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
    82      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    83      .                "ave(X)", zsto,zout)
    84 c
    85          CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
    86      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    87      .                "ave(X)", zsto,zout)
    88 c
    89          CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
    90      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    91      .                "ave(X)", zsto,zout)
    92 c
     60         CALL histdef(nid_day, "q2m", "Specific humidity", "Kg/Kg",
     61     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     62     .                "ave(X)", zsto,zout)
     63c
     64         CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
     65     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     66     .                "ave(X)", zsto,zout)
     67c
     68         CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
     69     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     70     .                "ave(X)", zsto,zout)
     71c
     72         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
     73     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     74     .                "ave(X)", zsto,zout)
     75c
     76         CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
     77     .                , "kg/s",
     78     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     79     .                "ave(X)", zsto,zout)
     80c
     81         CALL histdef(nid_day, "snow", "Snow fall", "kg/s",
     82     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     83     .                "ave(X)", zsto,zout)
     84c
     85         CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
     86     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     87     .                "ave(X)", zsto,zout)
     88c
     89         CALL histdef(nid_day, "evap", "Evaporation", "kg/s",
     90     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     91     .                "ave(X)", zsto,zout)
     92c
     93         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
     94     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     95     .                "ave(X)", zsto,zout)
     96c
     97         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
     98     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     99     .                "ave(X)", zsto,zout)
     100c
     101         CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
     102     .                "W/m2",
     103     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     104     .                "ave(X)", zsto,zout)
     105c
     106         CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
     107     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     108     .                "ave(X)", zsto,zout)
     109c
     110         CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
     111     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     112     .                "ave(X)", zsto,zout)
     113c
     114         CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
     115     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     116     .                "ave(X)", zsto,zout)
     117c
     118         CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
     119     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     120     .                "ave(X)", zsto,zout)
     121c
     122         CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2",
     123     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     124     .                "ave(X)", zsto,zout)
     125c
     126c        CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa",
     127c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     128c    .                "ave(X)", zsto,zout)
     129c
     130c        CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa",
     131c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     132c    .                "ave(X)", zsto,zout)
     133c
     134c        CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
     135c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     136c    .                "ave(X)", zsto,zout)
     137c
     138         CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
     139     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     140     .                "ave(X)", zsto,zout)
     141c
     142         CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
     143     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     144     .                "ave(X)", zsto,zout)
     145c
     146         CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
     147     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     148     .                "ave(X)", zsto,zout)
     149c
     150         CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
     151     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     152     .                "ave(X)", zsto,zout)
     153c
     154         CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-",
     155     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     156     .                "ave(X)", zsto,zout)
     157c
     158         CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
     159     .                iim,jjmp1,nhori, 1,1,1,-99,
     160     .                32, "ave(X)", zsto,zout)
     161c
     162         CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
     163     .                iim,jjmp1,nhori, 1,1,1,-99,
     164     .                32, "ave(X)", zsto,zout)
     165c
     166         CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
     167     .                iim,jjmp1,nhori, 1,1,1,-99,
     168     .                32, "ave(X)", zsto,zout)
     169c
     170         CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
     171     .                iim,jjmp1,nhori, 1,1,1,-99,
     172     .                32, "ave(X)", zsto,zout)
     173c
     174c  Champs dynamiques sur niveaux de pression
     175
     176         CALL histdef(nid_day, "u850", "Zonal wind 850mb", "m/s",
     177     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     178     .                "ave(X)", zsto,zout)
     179
     180         CALL histdef(nid_day, "v850", "Meridional wind 850mb", "m/s",
     181     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     182     .                "ave(X)", zsto,zout)
     183c
     184         CALL histdef(nid_day, "u500", "Zonal wind 500mb", "m/s",
     185     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     186     .                "ave(X)", zsto,zout)
     187
     188         CALL histdef(nid_day, "v500", "Meridional wind 500mb", "m/s",
     189     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     190     .                "ave(X)", zsto,zout)
     191
     192         CALL histdef(nid_day, "u200", "Zonal wind 200mb", "m/s",
     193     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     194     .                "ave(X)", zsto,zout)
     195
     196         CALL histdef(nid_day, "v200", "Meridional wind 200mb", "m/s",
     197     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     198     .                "ave(X)", zsto,zout)
     199
     200         CALL histdef(nid_day,"phi500", "Geopotentiel à 500mb", "m2/s2",
     201     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     202     .                "ave(X)", zsto,zout)
     203
     204         CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
     205     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     206     .                "ave(X)", zsto,zout)
     207
     208
     209c=================================================================
     210c   CI DESSOUS, SORTIES A UTILISER EN DEBUG
     211c=================================================================
     212      if ( lev_histday.gt.1 ) then
     213
     214c Champs retires momentannéement en attendant un hypothetique
     215c debugage
     216
    93217         CALL histdef(nid_day, "t2m_min", "Temp. 2m min.",
    94218     .                "K",
    … …  
    101225     .                t2maxcels, zsto,zout)
    102226c
     227c=================================================================
     228c=================================================================
     229c=================================================================
     230c   INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
     231c=================================================================
     232
     233         CALL histdef(nid_day, "tter", "Surface Temperature", "K",
     234     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     235     .                "ave(X)", zsto,zout)
     236c
     237         CALL histdef(nid_day, "tlic", "Surface Temperature", "K",
     238     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     239     .                "ave(X)", zsto,zout)
     240c
     241         CALL histdef(nid_day, "toce", "Surface Temperature", "K",
     242     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     243     .                "ave(X)", zsto,zout)
     244c
     245         CALL histdef(nid_day, "tsic", "Surface Temperature", "K",
     246     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     247     .                "ave(X)", zsto,zout)
     248c
     249cccIM
     250c
     251         CALL histdef(nid_day, "t2mter", "Temp.terre 2m", "K",
     252     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     253     .                "ave(X)", zsto,zout)
     254c
     255         CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
     256     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     257     .                "ave(X)", zsto,zout)
     258c
     259         CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
     260     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     261     .                "ave(X)", zsto,zout)
     262c
     263         CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
     264     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     265     .                "ave(X)", zsto,zout)
     266c
    103267         CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
    104268     .                "K",
    … …  
    111275     .                t2maxcels, zsto,zout)
    112276c
    113          CALL histdef(nid_day, "q2m", "Specific humidity", "Kg/Kg",
    114      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    115      .                "ave(X)", zsto,zout)
    116 c
    117          CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
    118      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    119      .                "ave(X)", zsto,zout)
    120 c
    121277         CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
    122278     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    … …  
    135291     .                "ave(X)", zsto,zout)
    136292c
    137          CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
    138      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    139      .                "ave(X)", zsto,zout)
    140 c
    141293         CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m",
    142294     .                "m/s", iim,jjmp1,nhori, 1,1,1, -99, 32,
    … …  
    155307     .                "ave(X)", zsto,zout)
    156308c
    157          CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
    158      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    159      .                "ave(X)", zsto,zout)
    160 c
    161          CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
    162      .                , "kg/s",
    163      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    164      .                "ave(X)", zsto,zout)
    165 c
    166          CALL histdef(nid_day, "snow", "Snow fall", "kg/s",
    167      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    168      .                "ave(X)", zsto,zout)
    169 c
    170          CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
    171      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    172      .                "ave(X)", zsto,zout)
    173 c
    174          CALL histdef(nid_day, "evap", "Evaporation", "kg/s",
    175      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    176      .                "ave(X)", zsto,zout)
    177 c
    178          CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
    179      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    180      .                "ave(X)", zsto,zout)
    181 c
    182          CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
    183      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    184      .                "ave(X)", zsto,zout)
    185 c
    186          CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
    187      .                "W/m2",
    188      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    189      .                "ave(X)", zsto,zout)
    190 c
    191          CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
    192      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    193      .                "ave(X)", zsto,zout)
    194 c
    195          CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
    196      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
    197      .                "ave(X)", zsto,zout)
    198 c
    199          CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
    200      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    201      .                "ave(X)", zsto,zout)
    202 c
    203          CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
    204      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    205      .                "ave(X)", zsto,zout)
    206 c
    207          CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2",
    208      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    209      .                "ave(X)", zsto,zout)
    210 c
    211          CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa",
    212      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    213      .                "ave(X)", zsto,zout)
    214 c
    215          CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa",
    216      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    217      .                "ave(X)", zsto,zout)
    218309c
    219310CXXX PB flux pour chaque sous surface
    … …  
    264355         END DO
    265356           
    266          CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
    267      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    268      .                "ave(X)", zsto,zout)
    269 c
    270          CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
    271      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    272      .                "ave(X)", zsto,zout)
    273 c
    274          CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
    275      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    276      .                "ave(X)", zsto,zout)
    277 c
    278          CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
    279      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    280      .                "ave(X)", zsto,zout)
    281 c
    282          CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
    283      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    284      .                "ave(X)", zsto,zout)
    285 c
    286          CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-",
    287      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    288      .                "ave(X)", zsto,zout)
    289 c
     357c=================================================================
     358c   FIN DES CHAMPS SUR LES SOUS SURFACES
     359c=================================================================
     360c=================================================================
     361c   FIN INITIALISATION DES CHAMPS 3D
     362c=================================================================
     363
    290364c Champs 3D:
    291365c
    … …  
    337411         CALL histend(nid_day)
    338412c
     413c=================================================================
     414c   FIN INITIALISATION DES CHAMPS 3D
     415c=================================================================
     416
     417
     418c=================================================================
     419c=================================================================
     420      endif !  lev_histday.gt.1
     421c=================================================================
     422
     423         CALL histend(nid_day)
     424c
    339425         ndex2d = 0
    340426         ndex3d = 0
    341427c
     428c=================================================================
    342429      ENDIF ! fin de test sur ok_journe

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h

    r418 r433  
    167167     $         "ave(X)", zsto,zout)
    168168C
     169           call histdef(nid_mth, "flw_"//clnsurf(nsrf),
     170     $         "LW "//clnsurf(nsrf), "W/m2",
     171     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     172     $         "ave(X)", zsto,zout)
     173cIM cf. JLD
     174           call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
     175     $         "SW "//clnsurf(nsrf), "W/m2",
     176     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     177     $         "ave(X)", zsto,zout)
     178C
     179           call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
     180     $         "Bilan sol "//clnsurf(nsrf), "W/m2",
     181     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     182     $         "ave(X)", zsto,zout)
     183C
    169184           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
    170185     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
    … …  
    301316c
    302317         CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
     318     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     319     .                "ave(X)", zsto,zout)
     320cIM cf. FH
     321         CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
     322     .                , "-",
    303323     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    304324     .                "ave(X)", zsto,zout)
    … …  
    479499     .                iim,jjmp1,nhori, 1,1,1,-99,
    480500     .                32, "ave(X)", zsto,zout)
    481 c
     501cccIM clear sky
     502         CALL histdef(nid_mth, "SWupTOAclr",
     503     .                "SWup clear sky at TOA","W/m2",
     504     .                iim,jjmp1,nhori, 1,1,1,-99,
     505     .                32, "ave(X)", zsto,zout)
     506c
     507         CALL histdef(nid_mth, "SWupSFCclr",
     508     .                "SWup clear sky at surface","W/m2",
     509     .                iim,jjmp1,nhori, 1,1,1,-99,
     510     .                32, "ave(X)", zsto,zout)
     511c
     512         CALL histdef(nid_mth, "SWdnTOAclr",
     513     .                "SWdn clear sky at TOA","W/m2",
     514     .                iim,jjmp1,nhori, 1,1,1,-99,
     515     .                32, "ave(X)", zsto,zout)
     516
     517c Champs interpolles sur des niveaux de pression
     518
     519         CALL histdef(nid_mth, "u850", "Zonal wind 850mb", "m/s",
     520     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     521     .                "ave(X)", zsto,zout)
     522
     523         CALL histdef(nid_mth, "v850", "Meridional wind 850mb", "m/s",
     524     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     525     .                "ave(X)", zsto,zout)
     526c
     527         CALL histdef(nid_mth, "u500", "Zonal wind 500mb", "m/s",
     528     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     529     .                "ave(X)", zsto,zout)
     530
     531         CALL histdef(nid_mth, "v500", "Meridional wind 500mb", "m/s",
     532     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     533     .                "ave(X)", zsto,zout)
     534
     535         CALL histdef(nid_mth, "u200", "Zonal wind 200mb", "m/s",
     536     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     537     .                "ave(X)", zsto,zout)
     538
     539         CALL histdef(nid_mth, "v200", "Meridional wind 200mb", "m/s",
     540     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     541     .                "ave(X)", zsto,zout)
     542
     543         CALL histdef(nid_mth, "phi500","Geopotentiel à 500mb","m2/s2",
     544     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     545     .                "ave(X)", zsto,zout)
     546
     547         CALL histdef(nid_mth, "w500", "Vertical wind 500mb", "Pa/s",
     548     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     549     .                "ave(X)", zsto,zout)
     550
     551
     552c
     553         CALL histdef(nid_mth, "SWdnSFCclr",
     554     .                "SWdn clear sky at surface","W/m2",
     555     .                iim,jjmp1,nhori, 1,1,1,-99,
     556     .                32, "ave(X)", zsto,zout)
    482557      ENDIF
    483558c34EK

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/phyetat0.F

    r373 r433  
    22c $Header$
    33c
    4       SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm,solaire,
     4      SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0,
    55     .            rlat,rlon, pctsrf, tsol,tsoil,deltat,qsol,snow,
    66     .           albe, evap, rain_fall, snow_fall, solsw, sollw,
    … …  
    2525      INTEGER radpas
    2626      REAL rlat(klon), rlon(klon)
    27       REAL co2_ppm
    28       REAL solaire
     27      REAL co2_ppm_etat0
     28      REAL solaire_etat0
    2929      REAL tsol(klon,nbsrf)
    3030      REAL tsoil(klon,nsoilmx,nbsrf)
    … …  
    145145         dtime        = tab_cntrl(1)
    146146         radpas       = tab_cntrl(2)
    147          co2_ppm      = tab_cntrl(3)
    148          solaire      = tab_cntrl(4)
     147         co2_ppm_etat0      = tab_cntrl(3)
     148         solaire_etat0      = tab_cntrl(4)
    149149         iflag_con    = tab_cntrl(5)
    150150         nbapp_rad    = tab_cntrl(6)

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/phyredem.F

    r373 r433  
    11c $Header$
    22c
    3       SUBROUTINE phyredem (fichnom,dtime,radpas,co2_ppm,solaire,
     3      SUBROUTINE phyredem (fichnom,dtime,radpas,
    44     .           rlat,rlon, pctsrf,tsol,tsoil,deltat,qsol,snow,
    55     .           albedo, evap, rain_fall, snow_fall,
    … …  
    2626      INTEGER radpas
    2727      REAL rlat(klon), rlon(klon)
    28       REAL co2_ppm
    29       REAL solaire
    3028      REAL tsol(klon,nbsrf)
    3129      REAL tsoil(klon,nsoilmx,nbsrf)

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

    r422 r433  
    123123      save ok_journe
    124124c      PARAMETER (ok_journe=.true.)
     125
     126      integer lev_histday
     127      save lev_histday
     128      data lev_histday/1/
    125129c
    126130      LOGICAL ok_mensuel ! sortir le fichier mensuel
    … …  
    199203cccIM cf. FH
    200204      real u850(klon),v850(klon),u200(klon),v200(klon)
    201       real u500(klon),v500(klon),phi500(klon)
     205      real u500(klon),v500(klon),phi500(klon),w500(klon)
    202206
    203207      logical ok_hf
    … …  
    248252      SAVE itap                   ! compteur pour la physique
    249253c
    250       REAL co2_ppm
    251       SAVE co2_ppm                ! concentration du CO2
    252 c
    253       REAL solaire
    254       SAVE solaire                ! constante solaire
    255 c
     254      REAL co2_ppm_etat0
     255c
     256      REAL solaire_etat0
     257c
     258      real slp(klon) ! sea level pressure
     259
    256260      REAL ftsol(klon,nbsrf)
    257261      SAVE ftsol                  ! temperature du sol
    … …  
    384388      REAL dlw(klon)    ! derivee infra rouge
    385389      REAL bils(klon) ! bilan de chaleur au sol
     390cIM cf. JLD
     391      REAL wfbils(klon,nbsrf) ! bilan de chaleur au sol, pour chaque
     392C                   type de sous-surface et pondere par la fraction
    386393      REAL fder(klon) ! Derive de flux (sensible et latente)
    387394      save fder
    … …  
    481488      REAL topsw0(klon), toplw0(klon), solsw0(klon), sollw0(klon)
    482489      REAL albpla(klon)
     490cIM cf. JLD
     491      REAL fsollw(klon, nbsrf)   ! bilan flux IR pour chaque sous surface
     492      REAL fsolsw(klon, nbsrf)   ! flux solaire absorb. pour chaque sous surface
    483493c Le rayonnement n'est pas calcule tous les pas, il faut donc
    484494c                      sauvegarder les sorties du rayonnement
    … …  
    529539      REAL cape(klon)           ! CAPE
    530540      SAVE cape
     541cccIM
     542      CHARACTER*40 capemaxcels
     543
    531544      REAL pbase(klon)          ! cloud base pressure
    532545      SAVE pbase
    … …  
    726739     .                  ok_instan, fact_cldcon, facttemps,ok_newmicro,
    727740     .                  iflag_cldcon,ratqsbas,ratqshaut, if_ebil)
     741cIM  .                  , RI0)
    728742
    729743         DO k = 2, nvm          ! pas de vegetation
    … …  
    745759         itaprad = 0
    746760c
    747          CALL phyetat0 ("startphy.nc",dtime,co2_ppm,solaire,
     761         CALL phyetat0 ("startphy.nc",dtime,co2_ppm_etat0,solaire_etat0,
    748762     .       rlat,rlon,pctsrf, ftsol,ftsoil,deltat,fqsol,fsnow,
    749763     .       falbe, fevap, rain_fall,snow_fall,solsw, sollwdown,
    … …  
    799813           ema_workcbmf(i) = 0.
    800814          ENDDO
     815
     816cIM15/11/02 rajout initialisation ibas_con,itop_con cf. SB =>BEG
     817          DO i = 1, klon
     818           ibas_con(i) = 1
     819           itop_con(i) = klev+1
     820          ENDDO
     821cIM15/11/02 rajout initialisation ibas_con,itop_con cf. SB =>END
     822
    801823         ENDIF
     824
    802825c34EK
    803826         IF (ok_orodr) THEN
    … …  
    856879c
    857880cccIM
     881      capemaxcels = 't_max(X)'
    858882      t2mincels = 't_min(X)'
    859883      t2maxcels = 't_max(X)'
    … …  
    10761100      ENDIF
    10771101
     1102      DO nsrf = 1, nbsrf
     1103      DO i = 1, klon
     1104        fsollw(i,nsrf) = sollwdown(i) - RSIGMA*ftsol(i,nsrf)**4
     1105        fsolsw(i,nsrf) = solsw(i)*(1.-falbe(i,nsrf))/(1.-albsol(i))
     1106      ENDDO
     1107      ENDDO
     1108
    10781109      fder = dlw
    10791110
    … …  
    10851116     $            paprs,pplay,radsol, fsnow,fqsol,fevap,falbe,falblw,
    10861117     $            fluxlat,
    1087      e            rain_fall, snow_fall, solsw, sollw, sollwdown, fder,
     1118cIM cf. JLD  e            rain_fall, snow_fall, solsw, sollw, sollwdown, fder,
     1119     e            rain_fall, snow_fall, fsolsw, fsollw, sollwdown, fder,
    10881120     e            rlon, rlat, cufi, cvfi, frugs,
    10891121     e            debut, lafin, agesno,rugoro ,
    … …  
    11661198c        IF (pctsrf(i,nsrf) .GE. EPSFRA) THEN
    11671199            ftsol(i,nsrf) = ftsol(i,nsrf) + d_ts(i,nsrf)
     1200cIM cf. JLD
     1201            wfbils(i,nsrf) = ( fsolsw(i,nsrf) + fsollw(i,nsrf)
     1202     $         + fluxt(i,nsrf) + fluxlat(i,nsrf) ) * pctsrf(i,nsrf)
    11681203            zxtsol(i) = zxtsol(i) + ftsol(i,nsrf)*pctsrf(i,nsrf)
    11691204            zxfluxlat(i) = zxfluxlat(i) + fluxlat(i,nsrf)*pctsrf(i,nsrf)
    … …  
    12711306     .        Ma,cape,tvp,iflagctrl,
    12721307     .        pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,qcondc,wd)
     1308cIM cf. FH
     1309              clwcon0=qcondc
    12731310
    12741311          ELSE ! ok_cvl
    1275 
    1276           if (iflag_con.eq.4) then ! vectorise
    1277           CALL conemav (dtime,paprs,pplay,t_seri,q_seri,
    1278      .        u_seri,v_seri,tr_seri,nbtr,
    1279      .        ema_work1,ema_work2,
    1280      .        d_t_con,d_q_con,d_u_con,d_v_con,d_tr,
    1281      .        rain_con, snow_con, ibas_con, itop_con,
    1282      .        upwd,dnwd,dnwd0,
    1283 c    .        Ma,cape,tvp,(/(nint(rflag(i)),i=1,size(rflag))/),
    1284      .        Ma,cape,tvp,iflagctrl,
    1285      .        pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr)
    1286 
    1287             qcondc=0.0
    1288 
    1289           else
    12901312
    12911313c          print*,'Avant conema OUI'
    … …  
    13031325          print*,'Apres conema3 '
    13041326
     1327          ENDIF ! ok_cvl
     1328
    13051329           IF (.NOT. ok_gust) THEN
    13061330           do i = 1, klon
    … …  
    13091333           ENDIF
    13101334
    1311 c Calculer l'humidite relative pour diagnostique
     1335c =================================================================== c
     1336c Calcul des proprietes des nuages convectifs
    13121337c
    13131338      DO k = 1, klev
    … …  
    13361361     s       (klon,klev,q_seri,zqsat,clwcon0,ptconv,ratqsc,rnebcon0)
    13371362
    1338           endif
    1339 
    1340           ENDIF ! ok_cvl
     1363c =================================================================== c
    13411364
    13421365          DO i = 1, klon
    … …  
    16681691!      albsollw = albsollw1
    16691692      CALL radlwsw ! nouveau rayonnement (compatible Arpege-IFS)
    1670      e            (dist, rmu0, fract, co2_ppm, solaire,
     1693cIM  e            (dist, rmu0, fract, co2_ppm, solaire,
     1694     e            (dist, rmu0, fract,
    16711695     e             paprs, pplay,zxtsol,albsol, albsollw, t_seri,q_seri,
    16721696     e             wo,
    … …  
    16811705      ENDIF
    16821706      itaprad = itaprad + 1
     1707
    16831708c
    16841709c Ajouter la tendance des rayonnements (tous les pas)
    … …  
    19231948      call plevel(klon,klev,.false.,pplay,20000.,v_seri,v200)
    19241949      call plevel(klon,klev,.true. ,pplay,50000.,zphi,phi500)
     1950      call plevel(klon,klev,.true. ,paprs,50000.,omega,w500)
     1951
     1952      slp(:) = paprs(:,1)*exp(pphis(:)/(289.*t_seri(:,1)))
     1953c
     1954
    19251955
    19261956c=============================================================

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/radlwsw.F

    r412 r433  
    1       SUBROUTINE radlwsw(dist, rmu0, fract, co2_ppm, solaire,
     1cIM   SUBROUTINE radlwsw(dist, rmu0, fract, co2_ppm, solaire,
     2      SUBROUTINE radlwsw(dist, rmu0, fract,
    23     .                  paprs, pplay,tsol,albedo, alblw, t,q,wo,
    34     .                  cldfra, cldemi, cldtau,
    … …  
    4344c
    4445      real rmu0(klon), fract(klon), dist
    45       real co2_ppm
    46       real solaire
     46cIM   real co2_ppm
     47cIM   real solaire
     48#include "clesphys.h"
    4749c
    4850      real paprs(klon,klev+1), pplay(klon,klev)
    … …  
    6971      EXTERNAL lw, sw
    7072c
    71       REAL*8 RCO2, RCH4, RN2O, RCFC11, RCFC12
     73cIM ctes ds clesphys.h  REAL*8 RCO2, RCH4, RN2O, RCFC11, RCFC12
    7274      REAL*8 PSCT
    7375c
    … …  
    119121      zdist = dist
    120122c
    121       RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97
    122       RCH4 = 1.72E-06* 16.043/28.97
    123       RN2O = 310.E-09* 44.013/28.97
    124       RCFC11 = 280.E-12* 137.3686/28.97
    125       RCFC12 = 484.E-12* 120.9140/28.97
     123cIM anciennes valeurs
     124c     RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97
     125c
     126cIM : on met RCO2, RCH4, RN2O, RCFC11 et RCFC12 dans clesphys.h /lecture ds conf_phys.F90
     127c     RCH4 = 1.65E-06* 16.043/28.97
     128c     RN2O = 306.E-09* 44.013/28.97
     129c     RCFC11 = 280.E-12* 137.3686/28.97
     130c     RCFC12 = 484.E-12* 120.9140/28.97
     131cIM anciennes valeurs
     132c     RCH4 = 1.72E-06* 16.043/28.97
     133c     RN2O = 310.E-09* 44.013/28.97
     134c
     135c     PRINT*,'IMradlwsw : solaire, co2= ', solaire, co2_ppm
    126136      PSCT = solaire/zdist/zdist
    127137c
    … …  
    138148!         PALBP(i,2) = albedo(iof+i)
    139149         PALBP(i,2) = alblw(iof+i)
    140          PEMIS(i) = 0.96
     150cIM cf. JLD pour etre en accord avec ORCHIDEE il faut mettre PEMIS(i) = 0.96
     151         PEMIS(i) = 1.0
    141152         PVIEW(i) = 1.66
    142153         PPSOL(i) = paprs(iof+i,1)
    … …  
    191202c
    192203c======================================================================
    193       CALL LW(RCO2,RCH4,RN2O,RCFC11,RCFC12,
     204cIM ctes ds clesphys.h   CALL LW(RCO2,RCH4,RN2O,RCFC11,RCFC12,
     205      CALL LW(
    194206     .        PPMB, PDP,
    195207     .        PPSOL,PDT0,PEMIS,
    … …  
    200212     .        ztoplw,zsollw,ztoplw0,zsollw0,
    201213     .        zsollwdown)
    202       CALL SW(PSCT, RCO2, zrmu0, zfract,
     214cIM ctes ds clesphys.h   CALL SW(PSCT, RCO2, zrmu0, zfract,
     215      CALL SW(PSCT, zrmu0, zfract,
    203216     S        PPMB, PDP,
    204217     S        PPSOL, PALBD, PALBP,
    … …  
    243256      RETURN
    244257      END
    245       SUBROUTINE SW(PSCT, RCO2, PRMU0, PFRAC,
     258cIM ctes ds clesphys.h   SUBROUTINE SW(PSCT, RCO2, PRMU0, PFRAC,
     259      SUBROUTINE SW(PSCT, PRMU0, PFRAC,
    246260     S              PPMB, PDP,
    247261     S              PPSOL, PALBD, PALBP,
    … …  
    292306C
    293307      REAL*8 PSCT  ! constante solaire (valeur conseillee: 1370)
    294       REAL*8 RCO2  ! concentration CO2 (IPCC: 353.E-06*44.011/28.97)
     308cIM ctes ds clesphys.h   REAL*8 RCO2  ! concentration CO2 (IPCC: 353.E-06*44.011/28.97)
     309#include "clesphys.h"
    295310C
    296311      REAL*8 PPSOL(KDLON)        ! SURFACE PRESSURE (PA)
    … …  
    373388C
    374389c clear-sky:
    375       CALL SWU(PSCT,RCO2,ZCLDSW0,PPMB,PPSOL,
     390cIM ctes ds clesphys.h  CALL SWU(PSCT,RCO2,ZCLDSW0,PPMB,PPSOL,
     391      CALL SWU(PSCT,ZCLDSW0,PPMB,PPSOL,
    376392     S         PRMU0,PFRAC,PTAVE,PWV,
    377393     S         ZAKI,ZCLD,ZCLEAR,ZDSIG,ZFACT,ZRMU,ZSEC,ZUD)
    … …  
    394410      ENDDO
    395411c cloudy-sky:
    396       CALL SWU(PSCT,RCO2,PCLDSW,PPMB,PPSOL,
     412cIM ctes ds clesphys.h   CALL SWU(PSCT,RCO2,PCLDSW,PPMB,PPSOL,
     413      CALL SWU(PSCT,PCLDSW,PPMB,PPSOL,
    397414     S         PRMU0,PFRAC,PTAVE,PWV,
    398415     S         ZAKI,ZCLD,ZCLEAR,ZDSIG,ZFACT,ZRMU,ZSEC,ZUD)
    … …  
    443460      END
    444461c
    445       SUBROUTINE SWU (PSCT,RCO2,PCLDSW,PPMB,PPSOL,PRMU0,PFRAC,
     462cIM ctes ds clesphys.h   SUBROUTINE SWU (PSCT,RCO2,PCLDSW,PPMB,PPSOL,PRMU0,PFRAC,
     463      SUBROUTINE SWU (PSCT,PCLDSW,PPMB,PPSOL,PRMU0,PFRAC,
    446464     S                PTAVE,PWV,PAKI,PCLD,PCLEAR,PDSIG,PFACT,
    447465     S                PRMU,PSEC,PUD)
    … …  
    457475C
    458476      REAL*8 PSCT
    459       REAL*8 RCO2
     477cIM ctes ds clesphys.h   REAL*8 RCO2
     478#include "clesphys.h"
    460479      REAL*8 PCLDSW(KDLON,KFLEV)
    461480      REAL*8 PPMB(KDLON,KFLEV+1)
    … …  
    24712490      RETURN
    24722491      END
    2473       SUBROUTINE LW(RCO2,RCH4,RN2O,RCFC11,RCFC12,
     2492cIM ctes ds clesphys.h   SUBROUTINE LW(RCO2,RCH4,RN2O,RCFC11,RCFC12,
     2493      SUBROUTINE LW(
    24742494     .              PPMB, PDP,
    24752495     .              PPSOL,PDT0,PEMIS,
    … …  
    25162536C        ORIGINAL : 89-07-14
    25172537C-----------------------------------------------------------------------
    2518       REAL*8 RCO2   ! CO2 CONCENTRATION (IPCC:353.E-06* 44.011/28.97)
    2519       REAL*8 RCH4   ! CH4 CONCENTRATION (IPCC: 1.72E-06* 16.043/28.97)
    2520       REAL*8 RN2O   ! N2O CONCENTRATION (IPCC: 310.E-09* 44.013/28.97)
    2521       REAL*8 RCFC11 ! CFC11 CONCENTRATION (IPCC: 280.E-12* 137.3686/28.97)
    2522       REAL*8 RCFC12 ! CFC12 CONCENTRATION (IPCC: 484.E-12* 120.9140/28.97)
     2538cIM ctes ds clesphys.h
     2539c     REAL*8 RCO2   ! CO2 CONCENTRATION (IPCC:353.E-06* 44.011/28.97)
     2540c     REAL*8 RCH4   ! CH4 CONCENTRATION (IPCC: 1.72E-06* 16.043/28.97)
     2541c     REAL*8 RN2O   ! N2O CONCENTRATION (IPCC: 310.E-09* 44.013/28.97)
     2542c     REAL*8 RCFC11 ! CFC11 CONCENTRATION (IPCC: 280.E-12* 137.3686/28.97)
     2543c     REAL*8 RCFC12 ! CFC12 CONCENTRATION (IPCC: 484.E-12* 120.9140/28.97)
     2544#include "clesphys.h"
    25232545      REAL*8 PCLDLD(KDLON,KFLEV)  ! DOWNWARD EFFECTIVE CLOUD COVER
    25242546      REAL*8 PCLDLU(KDLON,KFLEV)  ! UPWARD EFFECTIVE CLOUD COVER
    … …  
    25822604      ENDDO
    25832605      ENDDO
    2584       CALL LWU(RCO2,RCH4, RN2O, RCFC11, RCFC12,
     2606cIM ctes ds clesphys.h   CALL LWU(RCO2,RCH4, RN2O, RCFC11, RCFC12,
     2607      CALL LWU(
    25852608     S         PAER,PDP,PPMB,PPSOL,ZOZ,PTAVE,PVIEW,PWV,ZABCU)
    25862609      CALL LWBV(ILIM,PDP,PDT0,PEMIS,PPMB,PTL,PTAVE,ZABCU,
    … …  
    26202643      RETURN
    26212644      END
    2622       SUBROUTINE LWU(RCO2, RCH4, RN2O, RCFC11, RCFC12,
     2645cIM ctes ds clesphys.h   SUBROUTINE LWU(RCO2, RCH4, RN2O, RCFC11, RCFC12,
     2646      SUBROUTINE LWU(
    26232647     S               PAER,PDP,PPMB,PPSOL,POZ,PTAVE,PVIEW,PWV,
    26242648     S               PABCU)
    … …  
    26602684C-----------------------------------------------------------------------
    26612685C* ARGUMENTS:
    2662       REAL*8 RCO2
    2663       REAL*8 RCH4, RN2O, RCFC11, RCFC12
     2686cIM ctes ds clesphys.h
     2687c     REAL*8 RCO2
     2688c     REAL*8 RCH4, RN2O, RCFC11, RCFC12
     2689#include "clesphys.h"
    26642690      REAL*8 PAER(KDLON,KFLEV,5)
    26652691      REAL*8 PDP(KDLON,KFLEV)

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/stdlevvar.F90

    r419 r433  
    160160          conv_te(i) = (te_zref(i) - te_zref_p(i))/te_zref_p(i)
    161161          conv_q(i) = (q_zref(i) - q_zref_p(i))/q_zref_p(i)
    162           IF(abs(conv_te(i)).GE.0.0025.AND.abs(conv_q(i)).GE.0.05) THEN
    163             PRINT*,'DIV','i=',i,te_zref_p(i),te_zref(i),conv_te(i), &
    164             q_zref_p(i),q_zref(i),conv_q(i)
    165           ENDIF
     162!         IF(abs(conv_te(i)).GE.0.0025.AND.abs(conv_q(i)).GE.0.05) THEN
     163!           PRINT*,'DIV','i=',i,te_zref_p(i),te_zref(i),conv_te(i), &
     164!           q_zref_p(i),q_zref(i),conv_q(i)
     165!         ENDIF
    166166        ENDDO
    167167!

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/suphec.F

    r393 r433  
    33#include "YOMCST.h"
    44#include "YOETHF.h"
     5cIM cf. JLD
     6       LOGICAL firstcall
     7       SAVE firstcall
     8       DATA firstcall /.TRUE./
     9       IF (firstcall) THEN
     10         PRINT*, 'suphec initialise les constantes du GCM'
     11         firstcall = .FALSE.
     12       ELSE
     13         PRINT*, 'suphec DEJA APPELE '
     14         RETURN
     15       ENDIF
    516C      -----------------------------------------------------------------
    617C
    … …  
    4556c ref      R_peri = 102.04
    4657c ref      R_incl = 23.5
    47       R_ecc = 0.016724
    48       R_peri = 102.04
    49       R_incl = 23.5
     58c
     59cIM 161002 : pour avoir les ctes AMIP II
     60cIM 161002   R_ecc = 0.016724
     61cIM 161002   R_peri = 102.04
     62cIM 161002   R_incl = 23.5
     63cIM on mets R_ecc, R_peri, R_incl dans conf_phys.F90
     64c     R_ecc = 0.016715
     65c     R_peri = 102.7
     66c     R_incl = 23.441
    5067c
    5168      WRITE(UNIT=6,FMT='('' *** Astronomical constants ***'')')
    … …  
    5774      WRITE(UNIT=6,FMT='(''        omega = '',E13.7,'' s-1'')')
    5875     S                  ROMEGA
    59       write(unit=6,fmt='('' excentricite = '',e13.7,''-'')')R_ecc
    60       write(unit=6,fmt='(''     equinoxe = '',e13.7,''-'')')R_peri
    61       write(unit=6,fmt='(''  inclinaison = '',e13.7,''-'')')R_incl
     76c     write(unit=6,fmt='('' excentricite = '',e13.7,''-'')')R_ecc
     77c     write(unit=6,fmt='(''     equinoxe = '',e13.7,''-'')')R_peri
     78c     write(unit=6,fmt='(''  inclinaison = '',e13.7,''-'')')R_incl
    6279C
    6380C     ------------------------------------------------------------------
    … …  
    8299c z.x.li      RSIGMA=2. * RPI**5 * RKBOL**4 /(15.* RCLUM**2 * RHPLA**3)
    83100      rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15.
    84       RI0=1370.
     101cIM init. dans conf_phys.F90   RI0=1365.
    85102      WRITE(UNIT=6,FMT='('' ***        Radiation       ***'')')
    86103      WRITE(UNIT=6,FMT='('' Stefan-Bol.  = '',E13.7,'' W m-2 K-4''
    87104     S )')  RSIGMA
    88       WRITE(UNIT=6,FMT='('' Solar const. = '',E13.7,'' W m-2'')')
    89      S      RI0
     105cIM init. dans conf_phys.F90   WRITE(UNIT=6,FMT='('' Solar const. = '',E13.7,'' W m-2'')')
     106cIM init. dans conf_phys.F90  S      RI0
    90107C
    91108C     -----------------------------------------------------------------

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histday.h

    r418 r433  
    3434c
    3535C
    36       zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
    37       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
    38       CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    39 c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
    40 c    .              'Surface Temperature','K')
    41 C
    42       zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
    43       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    44       CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    45 c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
    46 c    .              'Surface Temperature','K')
    47 C
    48       zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
    49       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    50       CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    51 c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
    52 c    .              'Surface Temperature','K')
    53 C
    54       zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
    55       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    56       CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    57 c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
    58 c    .              'Surface Temperature','K')
    59 C
    60 cccIM
     36C
    6137      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    6238      CALL histwrite(nid_day,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    6339c
    64       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    65       CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
    66      .               iim*jjmp1,ndex2d)
    67 c
    68       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    69       CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
    70      .               iim*jjmp1,ndex2d)
    71 c
    72       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
    73       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    74       CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
    75      .               ndex2d)
    76 c
    77       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
    78       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    79       CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
    80      .               iim*jjmp1,ndex2d)
    81 c
    82       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
    83       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    84       CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
    85      .               iim*jjmp1,ndex2d)
    86 c
    87       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
    88       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    89       CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
    90      .               ndex2d)
    91 c
    92       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
    93       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    94       CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
    95      .               ndex2d)
    96 c
    97       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
    98       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    99       CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
    100      .               ndex2d)
    101 c
    10240      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zq2m,zx_tmp_2d)
    10341      CALL histwrite(nid_day,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    … …  
    11149     .               ndex2d)
    11250c
    113       zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
    114       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    115       CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
    116      .               ndex2d)
    117 c
    118       zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
    119       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    120       CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
    121      .               ndex2d)
    122 c
    123       zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
    124       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    125       CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
    126      .               ndex2d)
    127 c
    128       zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
    129       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    130       CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
    131      .               ndex2d)
    132 c
    133       zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
    134       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    135       CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
    136      .               ndex2d)
    137 c
    138       zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
    139       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    140       CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
    141      .               ndex2d)
    142 c
    143       zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
    144       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    145       CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
    146      .               ndex2d)
    147 C
    148       zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
    149       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    150       CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
    151      .               ndex2d)
    152 C
    15351      DO i = 1, klon
    15452         zx_tmp_fi2d(i) = paprs(i,1)
    … …  
    231129c
    232130c
     131cXXX      DO i = 1, klon
     132cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
     133cXXX      ENDDO
     134cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     135cXXX      CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     136c
     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)
     139c     call writephy(fid_day,prof2d_av,'cldl',cldl,
     140c    .              'Low-level cloudiness','-')
     141c
     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)
     144c     call writephy(fid_day,prof2d_av,'cldm',cldm,
     145c    .              'Mid-level cloudiness','-')
     146c
     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)
     149c     call writephy(fid_day,prof2d_av,'cldh',cldh,
     150c    .              'High-level cloudiness','-')
     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     call writephy(fid_day,prof2d_av,'cldt',cldt,
     155c    .              'Total cloudiness','-')
     156c
     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)
     159c     call writephy(fid_day,prof2d_av,'cldq',cldq,
     160c    .              'Cloud liquid water path','-')
     161c
     162      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
     163      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     164      CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
     165     .                               iim*jjmp1,ndex2d)
     166c
     167      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
     168      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     169      CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
     170     .                               iim*jjmp1,ndex2d)
     171c
     172      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
     173      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     174      CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
     175     .                               iim*jjmp1,ndex2d)
     176c
     177      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
     178      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     179      CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
     180     .                               iim*jjmp1,ndex2d)
     181
     182
     183c   Ecriture de champs dynamiques sur des niveaux de pression
     184
     185      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
     186      CALL histwrite(nid_day,"u850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     187
     188      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
     189      CALL histwrite(nid_day,"v850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     190
     191      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
     192      CALL histwrite(nid_day,"u500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     193
     194      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
     195      CALL histwrite(nid_day,"v500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     196
     197      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
     198      CALL histwrite(nid_day,"u200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     199
     200      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
     201      CALL histwrite(nid_day,"v200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     202
     203      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
     204      CALL histwrite(nid_day,"phi500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     205
     206      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     207      CALL histwrite(nid_day,"slp",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     208
     209
     210c
     211c=================================================================
     212c   CI DESSOUS, SORTIES A UTILISER EN DEBUG
     213c=================================================================
     214      if ( lev_histday.gt.1 ) then
     215
     216c   En attendant un eventuel debugage.
     217
     218      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     219      CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
     220     .               iim*jjmp1,ndex2d)
     221c
     222      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     223      CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
     224     .               iim*jjmp1,ndex2d)
     225c
     226c=================================================================
     227c=================================================================
     228c=================================================================
     229c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
     230c=================================================================
     231
     232      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
     233      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
     234      CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     235c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
     236c    .              'Surface Temperature','K')
     237C
     238      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
     239      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     240      CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     241c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
     242c    .              'Surface Temperature','K')
     243C
     244      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
     245      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     246      CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     247c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
     248c    .              'Surface Temperature','K')
     249C
     250      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
     251      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     252      CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     253c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
     254c    .              'Surface Temperature','K')
     255C
     256cccIM
     257      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     258      CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
     259     .               ndex2d)
     260c
     261      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
     262      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     263      CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
     264     .               iim*jjmp1,ndex2d)
     265c
     266      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
     267      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     268      CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
     269     .               iim*jjmp1,ndex2d)
     270c
     271      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
     272      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     273      CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
     274     .               ndex2d)
     275c
     276      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
     277      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     278      CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
     279     .               ndex2d)
     280c
     281      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
     282      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     283      CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
     284     .               ndex2d)
     285c
     286      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
     287      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     288      CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
     289     .               ndex2d)
     290c
     291      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
     292      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     293      CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
     294     .               ndex2d)
     295c
     296      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
     297      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     298      CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
     299     .               ndex2d)
     300c
     301      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
     302      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     303      CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
     304     .               ndex2d)
     305c
     306      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
     307      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     308      CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
     309     .               ndex2d)
     310c
     311      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
     312      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     313      CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
     314     .               ndex2d)
     315c
     316      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
     317      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     318      CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
     319     .               ndex2d)
     320C
     321      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
     322      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     323      CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
     324     .               ndex2d)
     325C
    233326      DO nsrf = 1, nbsrf
    234327CXXX
    … …  
    301394C
    302395      END DO 
    303 C
    304 cXXX      DO i = 1, klon
    305 cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
    306 cXXX      ENDDO
    307 cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    308 cXXX      CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    309 c
    310       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
    311       CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    312 c     call writephy(fid_day,prof2d_av,'cldl',cldl,
    313 c    .              'Low-level cloudiness','-')
    314 c
    315       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
    316       CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    317 c     call writephy(fid_day,prof2d_av,'cldm',cldm,
    318 c    .              'Mid-level cloudiness','-')
    319 c
    320       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
    321       CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    322 c     call writephy(fid_day,prof2d_av,'cldh',cldh,
    323 c    .              'High-level cloudiness','-')
    324 c
    325       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    326       CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    327 c     call writephy(fid_day,prof2d_av,'cldt',cldt,
    328 c    .              'Total cloudiness','-')
    329 c
    330       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
    331       CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    332 c     call writephy(fid_day,prof2d_av,'cldq',cldq,
    333 c    .              'Cloud liquid water path','-')
     396c=================================================================
     397c   FIN ECRITURE DES CHAMPS SUR LES SOUS SURFACES
     398c=================================================================
     399
     400c=================================================================
     401c   ECRITURE DES CHAMPS 3D
     402c=================================================================
     403
    334404c
    335405c Champs 3D:
    … …  
    380450c    .              'Air pressure','Pa')
    381451cccIM
    382       zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
    383       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    384       CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
    385      .                               iim*jjmp1,ndex2d)
    386 c
    387       zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
    388       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    389       CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
    390      .                               iim*jjmp1,ndex2d)
    391 c
    392       zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
    393       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    394       CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
    395      .                               iim*jjmp1,ndex2d)
    396 c
    397       zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
    398       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    399       CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
    400      .                               iim*jjmp1,ndex2d)
    401 c
     452
     453c=================================================================
     454c   FIN ECRITURE DES CHAMPS 3D
     455c=================================================================
     456
     457c=================================================================
     458c=================================================================
     459      endif !  lev_histday.gt.1
     460c=================================================================
     461
     462
    402463      if (ok_sync) then
    403464c       call writephy_sync(fid_day)
    404465        call histsync(nid_day)
    405466      endif
     467
    406468      ENDIF

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histmth.h

    r418 r433  
    142142        CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
    143143     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     144C
     145        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
     146        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     147        CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
     148     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     149cIM cf. JLD
     150        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
     151        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     152        CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
     153     $      zx_tmp_2d,iim*jjmp1,ndex2d)
     154C
     155        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
     156        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     157        CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
     158     $      zx_tmp_2d,iim*jjmp1,ndex2d)
    144159C
    145160        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
    … …  
    263278      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
    264279      CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
     280     .                                   iim*jjmp1*klev,ndex3d)
     281c
     282      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
     283      CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
    265284     .                                   iim*jjmp1*klev,ndex3d)
    266285c
    … …  
    441460     .                               iim*jjmp1,ndex2d)
    442461c
     462cccIM clear sky
     463      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
     464      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     465      CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
     466     .                               iim*jjmp1,ndex2d)
     467c
     468      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
     469      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     470      CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
     471     .                               iim*jjmp1,ndex2d)
     472c
     473      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
     474      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     475      CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
     476     .                               iim*jjmp1,ndex2d)
     477c
     478      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
     479      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     480      CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
     481     .                               iim*jjmp1,ndex2d)
     482c
     483c  Champs interpolles sur des niveaux de pression
     484
     485      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
     486      CALL histwrite(nid_mth,"u850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     487
     488      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
     489      CALL histwrite(nid_mth,"v850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     490
     491      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
     492      CALL histwrite(nid_mth,"u500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     493
     494      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
     495      CALL histwrite(nid_mth,"v500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     496
     497      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
     498      CALL histwrite(nid_mth,"u200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     499
     500      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
     501      CALL histwrite(nid_mth,"v200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     502
     503      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
     504      CALL histwrite(nid_mth,"phi500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     505
     506      CALL gr_fi_ecrit(1, klon,iim,jjmp1, w500,zx_tmp_2d)
     507      CALL histwrite(nid_mth,"w500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     508
     509
    443510     
    444       ENDIF
     511       ENDIF
    445512
    446513      if (ok_sync) then
Note: See TracChangeset for help on using the changeset viewer.