/[lmdze]/trunk/phylmd/physiq.f90
ViewVC logotype

Diff of /trunk/phylmd/physiq.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 212 by guez, Thu Jan 12 12:31:31 2017 UTC revision 215 by guez, Tue Mar 28 12:46:28 2017 UTC
# Line 151  contains Line 151  contains
151      ! soil temperature of surface fraction      ! soil temperature of surface fraction
152    
153      REAL, save:: fevap(klon, nbsrf) ! evaporation      REAL, save:: fevap(klon, nbsrf) ! evaporation
154      REAL, save:: fluxlat(klon, nbsrf)      REAL fluxlat(klon, nbsrf)
155    
156      REAL, save:: fqsurf(klon, nbsrf)      REAL, save:: fqsurf(klon, nbsrf)
157      ! humidite de l'air au contact de la surface      ! humidite de l'air au contact de la surface
158    
159      REAL, save:: qsol(klon)      REAL, save:: qsol(klon) ! column-density of water in soil, in kg m-2
160      ! column-density of water in soil, in kg m-2      REAL, save:: fsnow(klon, nbsrf) ! \'epaisseur neigeuse
   
     REAL, save:: fsnow(klon, nbsrf) ! epaisseur neigeuse  
161      REAL, save:: falbe(klon, nbsrf) ! albedo visible par type de surface      REAL, save:: falbe(klon, nbsrf) ! albedo visible par type de surface
162    
163      ! Param\`etres de l'orographie \`a l'\'echelle sous-maille (OESM) :      ! Param\`etres de l'orographie \`a l'\'echelle sous-maille (OESM) :
# Line 206  contains Line 204  contains
204      REAL, save:: pfrac_1nucl(klon, llm)      REAL, save:: pfrac_1nucl(klon, llm)
205      ! Produits des coefs lessi nucl (alpha = 1)      ! Produits des coefs lessi nucl (alpha = 1)
206    
207      REAL frac_impa(klon, llm) ! fractions d'aerosols lessivees (impaction)      REAL frac_impa(klon, llm) ! fraction d'a\'erosols lessiv\'es (impaction)
208      REAL frac_nucl(klon, llm) ! idem (nucleation)      REAL frac_nucl(klon, llm) ! idem (nucleation)
209    
210      REAL, save:: rain_fall(klon)      REAL, save:: rain_fall(klon)
# Line 276  contains Line 274  contains
274      REAL cldl(klon), cldm(klon), cldh(klon) ! nuages bas, moyen et haut      REAL cldl(klon), cldm(klon), cldh(klon) ! nuages bas, moyen et haut
275      REAL cldt(klon), cldq(klon) ! nuage total, eau liquide integree      REAL cldt(klon), cldq(klon) ! nuage total, eau liquide integree
276    
277      REAL zxqsurf(klon), zxsnow(klon), zxfluxlat(klon)      REAL zxqsurf(klon), zxfluxlat(klon)
278    
279      REAL dist, mu0(klon), fract(klon)      REAL dist, mu0(klon), fract(klon)
280      real longi      real longi
# Line 387  contains Line 385  contains
385      REAL ztsol(klon)      REAL ztsol(klon)
386    
387      REAL d_t_ec(klon, llm)      REAL d_t_ec(klon, llm)
388      ! tendance due \`a la conversion Ec en énergie thermique      ! tendance due \`a la conversion d'\'energie cin\'etique en
389        ! énergie thermique
     REAL ZRCPD  
390    
391      REAL, save:: t2m(klon, nbsrf), q2m(klon, nbsrf)      REAL, save:: t2m(klon, nbsrf), q2m(klon, nbsrf)
392      ! temperature and humidity at 2 m      ! temperature and humidity at 2 m
# Line 944  contains Line 941  contains
941    
942      ! Calculer l'hydrologie de la surface      ! Calculer l'hydrologie de la surface
943      zxqsurf = sum(fqsurf * pctsrf, dim = 2)      zxqsurf = sum(fqsurf * pctsrf, dim = 2)
     zxsnow = sum(fsnow * pctsrf, dim = 2)  
944    
945      ! Calculer le bilan du sol et la d\'erive de temp\'erature (couplage)      ! Calculer le bilan du sol et la d\'erive de temp\'erature (couplage)
946      DO i = 1, klon      DO i = 1, klon
# Line 1044  contains Line 1040  contains
1040      ! conversion Ec en énergie thermique      ! conversion Ec en énergie thermique
1041      DO k = 1, llm      DO k = 1, llm
1042         DO i = 1, klon         DO i = 1, klon
1043            ZRCPD = RCPD * (1. + RVTMP2 * q_seri(i, k))            d_t_ec(i, k) = 0.5 / (RCPD * (1. + RVTMP2 * q_seri(i, k))) &
           d_t_ec(i, k) = 0.5 / ZRCPD &  
1044                 * (u(i, k)**2 + v(i, k)**2 - u_seri(i, k)**2 - v_seri(i, k)**2)                 * (u(i, k)**2 + v(i, k)**2 - u_seri(i, k)**2 - v_seri(i, k)**2)
1045            t_seri(i, k) = t_seri(i, k) + d_t_ec(i, k)            t_seri(i, k) = t_seri(i, k) + d_t_ec(i, k)
1046            d_t_ec(i, k) = d_t_ec(i, k) / dtphys            d_t_ec(i, k) = d_t_ec(i, k) / dtphys
# Line 1156  contains Line 1151  contains
1151      CALL histwrite_phy("dtvdf", d_t_vdf)      CALL histwrite_phy("dtvdf", d_t_vdf)
1152      CALL histwrite_phy("dqvdf", d_q_vdf)      CALL histwrite_phy("dqvdf", d_q_vdf)
1153      CALL histwrite_phy("rhum", zx_rh)      CALL histwrite_phy("rhum", zx_rh)
1154        CALL histwrite_phy("d_t_ec", d_t_ec)
1155        CALL histwrite_phy("dtsw0", heat0 / 86400.)
1156        CALL histwrite_phy("dtlw0", - cool0 / 86400.)
1157        CALL histwrite_phy("msnow", sum(fsnow * pctsrf, dim = 2))
1158    
1159      if (ok_instan) call histsync(nid_ins)      if (ok_instan) call histsync(nid_ins)
1160    
Legend:
Removed from v.212  
changed lines
  Added in v.215
  ViewVC Help
Powered by ViewVC 1.1.21