158 |
! humidite de l'air au contact de la surface |
! humidite de l'air au contact de la surface |
159 |
|
|
160 |
REAL, save:: qsol(klon) ! column-density of water in soil, in kg m-2 |
REAL, save:: qsol(klon) ! column-density of water in soil, in kg m-2 |
161 |
REAL, save:: fsnow(klon, nbsrf) ! \'epaisseur neigeuse |
|
162 |
|
REAL, save:: fsnow(klon, nbsrf) |
163 |
|
! column-density of mass of snow at the surface, in kg m-2 |
164 |
|
|
165 |
REAL, save:: falbe(klon, nbsrf) ! albedo visible par type de surface |
REAL, save:: falbe(klon, nbsrf) ! albedo visible par type de surface |
166 |
|
|
167 |
! Param\`etres de l'orographie \`a l'\'echelle sous-maille (OESM) : |
! Param\`etres de l'orographie \`a l'\'echelle sous-maille (OESM) : |
220 |
REAL sens(klon) ! flux de chaleur sensible au sol |
REAL sens(klon) ! flux de chaleur sensible au sol |
221 |
real dflux_t(klon) ! derivee du flux de chaleur sensible au sol |
real dflux_t(klon) ! derivee du flux de chaleur sensible au sol |
222 |
REAL, save:: dlw(klon) ! derivative of infra-red flux |
REAL, save:: dlw(klon) ! derivative of infra-red flux |
|
REAL bils(klon) ! bilan de chaleur au sol |
|
223 |
REAL fder(klon) ! d\'erive de flux (sensible et latente) |
REAL fder(klon) ! d\'erive de flux (sensible et latente) |
224 |
REAL ve(klon) ! integr. verticale du transport meri. de l'energie |
REAL ve(klon) ! integr. verticale du transport meri. de l'energie |
225 |
REAL vq(klon) ! integr. verticale du transport meri. de l'eau |
REAL vq(klon) ! integr. verticale du transport meri. de l'eau |
846 |
ENDDO |
ENDDO |
847 |
ENDDO |
ENDDO |
848 |
|
|
|
! Calculer le bilan du sol et la d\'erive de temp\'erature (couplage) |
|
|
DO i = 1, klon |
|
|
bils(i) = radsol(i) + sens(i) + zxfluxlat(i) |
|
|
ENDDO |
|
|
|
|
849 |
! Param\'etrisation de l'orographie \`a l'\'echelle sous-maille : |
! Param\'etrisation de l'orographie \`a l'\'echelle sous-maille : |
850 |
|
|
851 |
IF (ok_orodr) THEN |
IF (ok_orodr) THEN |
981 |
CALL histwrite_phy("sols", solsw) |
CALL histwrite_phy("sols", solsw) |
982 |
CALL histwrite_phy("rls", sollw) |
CALL histwrite_phy("rls", sollw) |
983 |
CALL histwrite_phy("solldown", sollwdown) |
CALL histwrite_phy("solldown", sollwdown) |
984 |
CALL histwrite_phy("bils", bils) |
CALL histwrite_phy("bils", radsol + sens + zxfluxlat) |
985 |
CALL histwrite_phy("sens", sens) |
CALL histwrite_phy("sens", sens) |
986 |
CALL histwrite_phy("fder", fder) |
CALL histwrite_phy("fder", fder) |
987 |
CALL histwrite_phy("zxfqcalving", sum(fqcalving * pctsrf, dim = 2)) |
CALL histwrite_phy("zxfqcalving", sum(fqcalving * pctsrf, dim = 2)) |