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

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

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

revision 170 by guez, Wed Apr 29 15:47:56 2015 UTC revision 171 by guez, Tue Sep 29 19:48:59 2015 UTC
# Line 4  module calcul_fluxs_m Line 4  module calcul_fluxs_m
4    
5  contains  contains
6    
7    SUBROUTINE calcul_fluxs(nisurf, dtime, tsurf, p1lay, cal, beta, coef1lay, &    SUBROUTINE calcul_fluxs(dtime, tsurf, p1lay, cal, beta, coef1lay, ps, &
8         ps, qsurf, radsol, dif_grnd, t1lay, q1lay, u1lay, v1lay, petAcoef, &         qsurf, radsol, dif_grnd, t1lay, q1lay, u1lay, v1lay, petAcoef, &
9         peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, fluxlat, fluxsens, &         peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, fluxlat, fluxsens, &
10         dflux_s, dflux_l)         dflux_s, dflux_l)
11    
# Line 15  contains Line 15  contains
15      ! L. Fairhead April 2000      ! L. Fairhead April 2000
16    
17      USE abort_gcm_m, ONLY: abort_gcm      USE abort_gcm_m, ONLY: abort_gcm
     USE indicesol, ONLY: is_ter  
18      USE fcttre, ONLY: dqsatl, dqsats, foede, foeew, qsatl, qsats, thermcep      USE fcttre, ONLY: dqsatl, dqsats, foede, foeew, qsatl, qsats, thermcep
19      USE interface_surf, ONLY: run_off      USE indicesol, ONLY: is_ter
20      use nr_util, only: assert_eq      use nr_util, only: assert_eq
21      USE suphec_m, ONLY: rcpd, rd, retv, rkappa, rlstt, rlvtt, rtt      USE suphec_m, ONLY: rcpd, rd, retv, rkappa, rlstt, rlvtt, rtt
22      USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2      USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2
23    
     integer, intent(IN):: nisurf ! surface a traiter  
24      real, intent(IN):: dtime      real, intent(IN):: dtime
25      real, intent(IN):: tsurf(:) ! (knon) temperature de surface      real, intent(IN):: tsurf(:) ! (knon) temperature de surface
26      real, intent(IN):: p1lay(:) ! (knon) pression 1er niveau (milieu de couche)      real, intent(IN):: p1lay(:) ! (knon) pression 1er niveau (milieu de couche)
# Line 42  contains Line 40  contains
40      ! coefficients A de la résolution de la couche limite pour t et q      ! coefficients A de la résolution de la couche limite pour t et q
41    
42      real, intent(IN):: petBcoef(:), peqBcoef(:) ! (knon)      real, intent(IN):: petBcoef(:), peqBcoef(:) ! (knon)
43      ! petBcoef coeff. B de la resolution de la CL pour t      ! coeff. B de la resolution de la CL pour t et q
     ! peqBcoef coeff. B de la resolution de la CL pour q  
44    
45      real, intent(OUT):: tsurf_new(:) ! (knon) température au sol      real, intent(OUT):: tsurf_new(:) ! (knon) température au sol
46      real, intent(OUT):: evap(:), fluxlat(:), fluxsens(:) ! (knon)      real, intent(OUT):: evap(:) ! (knon)
47      ! fluxlat flux de chaleur latente  
48      ! fluxsens flux de chaleur sensible      real, intent(OUT):: fluxlat(:), fluxsens(:) ! (knon)
49        ! flux de chaleur latente et sensible
50    
51      real, intent(OUT):: dflux_s(:), dflux_l(:) ! (knon)      real, intent(OUT):: dflux_s(:), dflux_l(:) ! (knon)
52      ! Dérivées des flux dF/dTs (W m-2 K-1)      ! dérivées des flux de chaleurs sensible et latente par rapport à
53      ! dflux_s derivee du flux de chaleur sensible / Ts      ! Ts (W m-2 K-1)
     ! dflux_l derivee du flux de chaleur latente / Ts  
54    
55      ! Local:      ! Local:
56      integer i      integer i
# Line 74  contains Line 72  contains
72           size(evap), size(fluxlat), size(fluxsens), size(dflux_s), &           size(evap), size(fluxlat), size(fluxsens), size(dflux_s), &
73           size(dflux_l)/), "calcul_fluxs knon")           size(dflux_l)/), "calcul_fluxs knon")
74    
     if (size(run_off) /= knon .AND. nisurf == is_ter) then  
        print *, 'Bizarre, le nombre de points continentaux'  
        print *, 'a change entre deux appels. J''arrete.'  
        call abort_gcm('calcul_fluxs', 'Pb run_off', 1)  
     endif  
   
75      ! Traitement humidite du sol      ! Traitement humidite du sol
76    
77      IF (thermcep) THEN      IF (thermcep) THEN
Legend:
Removed from v.170  
changed lines
  Added in v.171
  ViewVC Help
Powered by ViewVC 1.1.21