/[lmdze]/trunk/Sources/phylmd/clmain.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/clmain.f

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

revision 149 by guez, Thu Jun 18 12:23:44 2015 UTC revision 175 by guez, Fri Feb 5 16:02:34 2016 UTC
# Line 6  contains Line 6  contains
6    
7    SUBROUTINE clmain(dtime, itap, pctsrf, pctsrf_new, t, q, u, v, jour, rmu0, &    SUBROUTINE clmain(dtime, itap, pctsrf, pctsrf_new, t, q, u, v, jour, rmu0, &
8         co2_ppm, ts, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ftsoil, qsol, &         co2_ppm, ts, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ftsoil, qsol, &
9         paprs, pplay, snow, qsurf, evap, albe, alblw, fluxlat, rain_fall, &         paprs, pplay, snow, qsurf, evap, falbe, fluxlat, rain_fall, snow_f, &
10         snow_f, solsw, sollw, fder, rlat, rugos, debut, agesno, rugoro, d_t, &         solsw, sollw, fder, rlat, rugos, debut, agesno, rugoro, d_t, d_q, d_u, &
11         d_q, d_u, d_v, d_ts, flux_t, flux_q, flux_u, flux_v, cdragh, cdragm, &         d_v, d_ts, flux_t, flux_q, flux_u, flux_v, cdragh, cdragm, q2, &
12         q2, dflux_t, dflux_q, ycoefh, zu1, zv1, t2m, q2m, u10m, v10m, pblh, &         dflux_t, dflux_q, ycoefh, zu1, zv1, t2m, q2m, u10m, v10m, pblh, capcl, &
13         capcl, oliqcl, cteicl, pblt, therm, trmb1, trmb2, trmb3, plcl, &         oliqcl, cteicl, pblt, therm, trmb1, trmb2, trmb3, plcl, fqcalving, &
14         fqcalving, ffonte, run_off_lic_0, flux_o, flux_g, tslab)         ffonte, run_off_lic_0)
15    
16      ! From phylmd/clmain.F, version 1.6, 2005/11/16 14:47:19      ! From phylmd/clmain.F, version 1.6, 2005/11/16 14:47:19
17      ! Author: Z. X. Li (LMD/CNRS), date: 1993/08/18      ! Author: Z. X. Li (LMD/CNRS), date: 1993/08/18
# Line 73  contains Line 73  contains
73      REAL snow(klon, nbsrf)      REAL snow(klon, nbsrf)
74      REAL qsurf(klon, nbsrf)      REAL qsurf(klon, nbsrf)
75      REAL evap(klon, nbsrf)      REAL evap(klon, nbsrf)
76      REAL albe(klon, nbsrf)      REAL, intent(inout):: falbe(klon, nbsrf)
     REAL alblw(klon, nbsrf)  
77    
78      REAL fluxlat(klon, nbsrf)      REAL fluxlat(klon, nbsrf)
79    
# Line 85  contains Line 84  contains
84      ! solid water mass flux (kg/m2/s), positive down      ! solid water mass flux (kg/m2/s), positive down
85    
86      REAL, INTENT(IN):: solsw(klon, nbsrf), sollw(klon, nbsrf)      REAL, INTENT(IN):: solsw(klon, nbsrf), sollw(klon, nbsrf)
87      REAL fder(klon)      REAL, intent(in):: fder(klon)
88      REAL, INTENT(IN):: rlat(klon) ! latitude en degr\'es      REAL, INTENT(IN):: rlat(klon) ! latitude en degr\'es
89    
90      REAL rugos(klon, nbsrf)      REAL rugos(klon, nbsrf)
# Line 150  contains Line 149  contains
149      !           hauteur de neige, en kg/m2/s      !           hauteur de neige, en kg/m2/s
150      REAL run_off_lic_0(klon)      REAL run_off_lic_0(klon)
151    
     REAL flux_o(klon), flux_g(klon)  
     !IM "slab" ocean  
     ! flux_g---output-R-  flux glace (pour OCEAN='slab  ')  
     ! flux_o---output-R-  flux ocean (pour OCEAN='slab  ')  
   
     REAL tslab(klon)  
     ! tslab-in/output-R temperature du slab ocean (en Kelvin)  
     ! uniqmnt pour slab  
   
152      ! Local:      ! Local:
153    
     REAL y_flux_o(klon), y_flux_g(klon)  
     real ytslab(klon)  
154      REAL y_fqcalving(klon), y_ffonte(klon)      REAL y_fqcalving(klon), y_ffonte(klon)
155      real y_run_off_lic_0(klon)      real y_run_off_lic_0(klon)
156    
# Line 172  contains Line 160  contains
160    
161      REAL yts(klon), yrugos(klon), ypct(klon), yz0_new(klon)      REAL yts(klon), yrugos(klon), ypct(klon), yz0_new(klon)
162      REAL yalb(klon)      REAL yalb(klon)
     REAL yalblw(klon)  
163      REAL yu1(klon), yv1(klon)      REAL yu1(klon), yv1(klon)
164      ! on rajoute en output yu1 et yv1 qui sont les vents dans      ! on rajoute en output yu1 et yv1 qui sont les vents dans
165      ! la premiere couche      ! la premiere couche
# Line 187  contains Line 174  contains
174      REAL ysnow_f(klon)      REAL ysnow_f(klon)
175      ! solid water mass flux (kg/m2/s), positive down      ! solid water mass flux (kg/m2/s), positive down
176    
     REAL ysollw(klon), ysolsw(klon)  
177      REAL yfder(klon)      REAL yfder(klon)
178      REAL yrugm(klon), yrads(klon), yrugoro(klon)      REAL yrugm(klon), yrads(klon), yrugoro(klon)
179    
# Line 282  contains Line 268  contains
268      yts = 0.      yts = 0.
269      ysnow = 0.      ysnow = 0.
270      yqsurf = 0.      yqsurf = 0.
     yalb = 0.  
     yalblw = 0.  
271      yrain_f = 0.      yrain_f = 0.
272      ysnow_f = 0.      ysnow_f = 0.
273      yfder = 0.      yfder = 0.
     ysolsw = 0.  
     ysollw = 0.  
274      yrugos = 0.      yrugos = 0.
275      yu1 = 0.      yu1 = 0.
276      yv1 = 0.      yv1 = 0.
# Line 350  contains Line 332  contains
332               i = ni(j)               i = ni(j)
333               ypct(j) = pctsrf(i, nsrf)               ypct(j) = pctsrf(i, nsrf)
334               yts(j) = ts(i, nsrf)               yts(j) = ts(i, nsrf)
              ytslab(i) = tslab(i)  
335               ysnow(j) = snow(i, nsrf)               ysnow(j) = snow(i, nsrf)
336               yqsurf(j) = qsurf(i, nsrf)               yqsurf(j) = qsurf(i, nsrf)
337               yalb(j) = albe(i, nsrf)               yalb(j) = falbe(i, nsrf)
              yalblw(j) = alblw(i, nsrf)  
338               yrain_f(j) = rain_fall(i)               yrain_f(j) = rain_fall(i)
339               ysnow_f(j) = snow_f(i)               ysnow_f(j) = snow_f(i)
340               yagesno(j) = agesno(i, nsrf)               yagesno(j) = agesno(i, nsrf)
341               yfder(j) = fder(i)               yfder(j) = fder(i)
              ysolsw(j) = solsw(i, nsrf)  
              ysollw(j) = sollw(i, nsrf)  
342               yrugos(j) = rugos(i, nsrf)               yrugos(j) = rugos(i, nsrf)
343               yrugoro(j) = rugoro(i)               yrugoro(j) = rugoro(i)
344               yu1(j) = u1lay(i)               yu1(j) = u1lay(i)
345               yv1(j) = v1lay(i)               yv1(j) = v1lay(i)
346               yrads(j) = ysolsw(j) + ysollw(j)               yrads(j) = solsw(i, nsrf) + sollw(i, nsrf)
347               ypaprs(j, klev+1) = paprs(i, klev+1)               ypaprs(j, klev+1) = paprs(i, klev+1)
348               y_run_off_lic_0(j) = run_off_lic_0(i)               y_run_off_lic_0(j) = run_off_lic_0(i)
349               yu10mx(j) = u10m(i, nsrf)               yu10mx(j) = u10m(i, nsrf)
# Line 480  contains Line 458  contains
458            CALL clqh(dtime, itap, jour, debut, rlat, knon, nsrf, ni(:knon), &            CALL clqh(dtime, itap, jour, debut, rlat, knon, nsrf, ni(:knon), &
459                 pctsrf, ytsoil, yqsol, rmu0, co2_ppm, yrugos, yrugoro, yu1, &                 pctsrf, ytsoil, yqsol, rmu0, co2_ppm, yrugos, yrugoro, yu1, &
460                 yv1, coefh(:knon, :), yt, yq, yts, ypaprs, ypplay, ydelp, &                 yv1, coefh(:knon, :), yt, yq, yts, ypaprs, ypplay, ydelp, &
461                 yrads, yalb, yalblw, ysnow, yqsurf, yrain_f, ysnow_f, yfder, &                 yrads, yalb(:knon), ysnow, yqsurf, yrain_f, ysnow_f, yfder, &
462                 ysolsw, yfluxlat, pctsrf_new, yagesno, y_d_t, y_d_q, &                 yfluxlat, pctsrf_new, yagesno(:knon), y_d_t, y_d_q, &
463                 y_d_ts(:knon), yz0_new, y_flux_t, y_flux_q, y_dflux_t, &                 y_d_ts(:knon), yz0_new, y_flux_t, y_flux_q, y_dflux_t, &
464                 y_dflux_q, y_fqcalving, y_ffonte, y_run_off_lic_0, y_flux_o, &                 y_dflux_q, y_fqcalving, y_ffonte, y_run_off_lic_0)
                y_flux_g)  
465    
466            ! calculer la longueur de rugosite sur ocean            ! calculer la longueur de rugosite sur ocean
467            yrugm = 0.            yrugm = 0.
# Line 520  contains Line 497  contains
497    
498            evap(:, nsrf) = -flux_q(:, 1, nsrf)            evap(:, nsrf) = -flux_q(:, 1, nsrf)
499    
500            albe(:, nsrf) = 0.            falbe(:, nsrf) = 0.
           alblw(:, nsrf) = 0.  
501            snow(:, nsrf) = 0.            snow(:, nsrf) = 0.
502            qsurf(:, nsrf) = 0.            qsurf(:, nsrf) = 0.
503            rugos(:, nsrf) = 0.            rugos(:, nsrf) = 0.
# Line 529  contains Line 505  contains
505            DO j = 1, knon            DO j = 1, knon
506               i = ni(j)               i = ni(j)
507               d_ts(i, nsrf) = y_d_ts(j)               d_ts(i, nsrf) = y_d_ts(j)
508               albe(i, nsrf) = yalb(j)               falbe(i, nsrf) = yalb(j)
              alblw(i, nsrf) = yalblw(j)  
509               snow(i, nsrf) = ysnow(j)               snow(i, nsrf) = ysnow(j)
510               qsurf(i, nsrf) = yqsurf(j)               qsurf(i, nsrf) = yqsurf(j)
511               rugos(i, nsrf) = yz0_new(j)               rugos(i, nsrf) = yz0_new(j)
# Line 637  contains Line 612  contains
612                  q2(i, k, nsrf) = yq2(j, k)                  q2(i, k, nsrf) = yq2(j, k)
613               END DO               END DO
614            END DO            END DO
           !IM "slab" ocean  
           IF (nsrf == is_oce) THEN  
              DO j = 1, knon  
                 ! on projette sur la grille globale  
                 i = ni(j)  
                 IF (pctsrf_new(i, is_oce)>epsfra) THEN  
                    flux_o(i) = y_flux_o(j)  
                 ELSE  
                    flux_o(i) = 0.  
                 END IF  
              END DO  
           END IF  
   
           IF (nsrf == is_sic) THEN  
              DO j = 1, knon  
                 i = ni(j)  
                 ! On pond\`ere lorsque l'on fait le bilan au sol :  
                 IF (pctsrf_new(i, is_sic)>epsfra) THEN  
                    flux_g(i) = y_flux_g(j)  
                 ELSE  
                    flux_g(i) = 0.  
                 END IF  
              END DO  
   
           END IF  
615         end IF if_knon         end IF if_knon
616      END DO loop_surface      END DO loop_surface
617    
Legend:
Removed from v.149  
changed lines
  Added in v.175
  ViewVC Help
Powered by ViewVC 1.1.21