/[lmdze]/trunk/phylmd/Interface_surf/calcul_fluxs.f
ViewVC logotype

Diff of /trunk/phylmd/Interface_surf/calcul_fluxs.f

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

trunk/Sources/phylmd/Interface_surf/calcul_fluxs.f revision 207 by guez, Thu Sep 1 10:30:53 2016 UTC trunk/phylmd/Interface_surf/calcul_fluxs.f revision 266 by guez, Thu Apr 19 17:54:55 2018 UTC
# Line 12  contains Line 12  contains
12      ! Cette routine calcule les flux en h et q à l'interface et une      ! Cette routine calcule les flux en h et q à l'interface et une
13      ! température de surface.      ! température de surface.
14    
15      ! L. Fairhead April 2000      ! L. Fairhead, April 2000
16    
17      USE fcttre, ONLY: dqsatl, dqsats, foede, foeew, qsatl, qsats      USE fcttre, ONLY: foede, foeew
18      use nr_util, only: assert_eq      use nr_util, only: assert_eq
19      USE suphec_m, ONLY: rcpd, rd, retv, rlstt, rlvtt, rtt      USE suphec_m, ONLY: rcpd, rd, retv, rlstt, rlvtt, rtt
20      USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2      USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2
# Line 90  contains Line 90  contains
90      sl = merge(RLSTT, RLVTT, tsurf < RTT)      sl = merge(RLSTT, RLVTT, tsurf < RTT)
91    
92      ! Q      ! Q
93      oq = 1. - (beta * coef * peqBcoef * dtime)      oq = 1. - beta * coef * peqBcoef * dtime
94      mq = beta * coef * (peqAcoef - qsat + dq_s_dt * tsurf) / oq      mq = beta * coef * (peqAcoef - qsat + dq_s_dt * tsurf) / oq
95      nq = beta * coef * (- 1. * dq_s_dt) / oq      nq = beta * coef * (- 1. * dq_s_dt) / oq
96    
Legend:
Removed from v.207  
changed lines
  Added in v.266
  ViewVC Help
Powered by ViewVC 1.1.21