Changeset 457 for LMDZ.3.3

Timestamp:

Mar 7, 2003, 5:57:53 PM (21 years ago)

Author:

lmdzadmin

Message:

Correction bug tsurf_new=RTT - 1.8 apres fonte et passage routine soil
Correction bug tsurf_new=RTT sic & lic
IM/JLD/LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90 (modified) (16 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

@@ -67 +67 @@
       & evap, fluxsens, fluxlat, dflux_l, dflux_s, &              
       & tsol_rad, tsurf_new, alb_new, alblw, emis_new, &
-      & z0_new, pctsrf_new, agesno)
+!IM cf. JLD      & z0_new, pctsrf_new, agesno)
+      & z0_new, pctsrf_new, agesno,fqcalving,ffonte)
 
 
@@ -165 +166 @@
   real, dimension(klon), intent(IN) :: sollw, sollwdown, swnet, swdown
   real, dimension(klon), intent(IN) :: ps, albedo
+!IM cf LF
+! real, dimension(klon), intent(INOUT) :: tsurf
+! real, dimension(klon), intent(IN) :: p1lay
   real, dimension(klon), intent(IN) :: tsurf, p1lay
   REAL, DIMENSION(klon), INTENT(INOUT) :: radsol,fder
@@ -188 +192 @@
   real, dimension(klon), intent(INOUT):: agesno
 
+! Flux thermique utiliser pour fondre la neige
+!jld a rajouter   real, dimension(klon), intent(INOUT):: ffonte
+  real, dimension(klon), intent(INOUT):: ffonte
+! Flux d'eau "perdue" par la surface et nécessaire pour que limiter la
+! hauteur de neige, en kg/m2/s
+!jld a rajouter   real, dimension(klon), intent(INOUT):: fqcalving
+  real, dimension(klon), intent(INOUT):: fqcalving
+
 ! Local
   character (len = 20),save :: modname = 'interfsurf_hq'
@@ -206 +218 @@
   real, DIMENSION(klon):: zfra
   logical              :: cumul = .false.
+!IM cf LF
+  integer :: i
 
   if (check) write(*,*) 'Entree ', modname
@@ -334 +348 @@
      &   radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, &
      &   petAcoef, peqAcoef, petBcoef, peqBcoef, &
-     &   tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
+     &   tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, &
+!IM cf JLD
+     &   fqcalving,ffonte)
 
 
@@ -494 +510 @@
              &  tsurf_new, pctsrf_new)
 
+!IM cf LF
+      DO i = 1, knon
+       IF (pctsrf_new(i,nisurf) < EPSFRA) then
+          snow(i) = 0.0
+!IM cf LF/JLD         tsurf(i) = RTT - 1.8
+          tsurf_new(i) = RTT - 1.8
+          IF (soil_model) tsoil(i,:) = RTT -1.8
+        endif
+      enddo
+
       CALL calbeta(dtime, nisurf, knon, snow, qsol, beta, capsol, dif_grnd)
       
       IF (soil_model) THEN 
-         CALL soil(dtime, nisurf, knon,snow, tsurf, tsoil,soilcap, soilflux)
+!IM cf LF/JLD        CALL soil(dtime, nisurf, knon,snow, tsurf, tsoil,soilcap, soilflux)
+         CALL soil(dtime, nisurf, knon,snow, tsurf_new, tsoil,soilcap, soilflux)
          cal(1:knon) = RCPD / soilcap(1:knon)
          radsol(1:knon) = radsol(1:knon)  + soilflux(1:knon)
@@ -517 +544 @@
          &   tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
 
+!IM cf JP 12.02.03
+!     DO i = 1, knon
+!      IF (pctsrf_new(i,nisurf) < EPSFRA) then
+!         snow(i) = 0.0      
+!         tsurf_new(i) = RTT - 1.8
+!         IF (soil_model) tsoil(i,:) = RTT -1.8
+!       endif
+!     enddo
+
     IF (ocean /= 'couple') THEN
       CALL fonte_neige( klon, knon, nisurf, dtime, &
@@ -523 +559 @@
              &   radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, &
              &   petAcoef, peqAcoef, petBcoef, peqBcoef, &
-             &   tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
+             &   tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, &
+!IM cf JLD
+             &   fqcalving,ffonte)
 
 !     calcul albedo
@@ -592 +630 @@
      &   radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, &
      &   petAcoef, peqAcoef, petBcoef, peqBcoef, &
-     &   tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
+     &   tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, &
+!IM cf JLD
+     &   fqcalving,ffonte)
 
 !
@@ -2439 +2479 @@
      & radsol, dif_grnd, t1lay, q1lay, u1lay, v1lay, &
      & petAcoef, peqAcoef, petBcoef, peqBcoef, &
-     & tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
+     & tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, &
+!IM cf JLD
+     & fqcalving,ffonte)
 
 ! Routine de traitement de la fonte de la neige dans le cas du traitement
@@ -2477 +2519 @@
 #include "FCTTRE.inc"
 #include "indicesol.inc"
+!IM cf JLD
+!#include "YOMCST.inc"
 
 ! Parametres d'entree
@@ -2493 +2537 @@
   real, dimension(klon), intent(INOUT):: tsurf_new, evap, fluxsens, fluxlat
   real, dimension(klon), intent(INOUT):: dflux_s, dflux_l
+! Flux thermique utiliser pour fondre la neige
+  real, dimension(klon), intent(INOUT):: ffonte
+! Flux d'eau "perdue" par la surface et necessaire pour que limiter la
+! hauteur de neige, en kg/m2/s
+  real, dimension(klon), intent(INOUT):: fqcalving
 
 ! Variables locales
+! Masse maximum de neige (kg/m2). Au dessus de ce seuil, la neige
+! en exces "s'ecoule" (calving)
+  real, parameter :: snow_max=1.
   integer :: i
   real, dimension(klon) :: zx_mh, zx_nh, zx_oh
@@ -2507 +2559 @@
   real, parameter :: t_grnd = 271.35, t_coup = 273.15
 !! PB temporaire en attendant mieux pour le modele de neige
+! REAL, parameter :: chasno = RLMLT/(2.3867E+06*0.15)
   REAL, parameter :: chasno = 3.334E+05/(2.3867E+06*0.15)
 !
@@ -2585 +2638 @@
 ! Y'a-t-il fonte de neige?
 !
+  ffonte=0.
   do i = 1, knon
     neige_fond = ((snow(i) > epsfra .OR. nisurf == is_sic .OR. nisurf == is_lic) &
@@ -2590 +2644 @@
     if (neige_fond) then
       fq_fonte = MIN( MAX((tsurf_new(i)-RTT )/chasno,0.0),snow(i))
+      ffonte(i) = fq_fonte * RLMLT/dtime
       snow(i) = max(0., snow(i) - fq_fonte)
       bil_eau_s(i) = bil_eau_s(i) + fq_fonte 
       tsurf_new(i) = tsurf_new(i) - fq_fonte * chasno  
-      IF (nisurf == is_sic .OR. nisurf == is_lic ) tsurf_new(i) = RTT -1.8
+!IM cf JLD OK     
+      IF (nisurf == is_sic .OR. nisurf == is_lic ) tsurf_new(i) = RTT
       d_ts(i) = tsurf_new(i) - tsurf(i)
 !      zx_h_ts(i) = tsurf_new(i) * RCPD * zx_pkh(i)
@@ -2611 +2667 @@
 !!$      fq_fonte = bilan_f / zx_sl(i)
     endif
-
+!
+!   s'il y a une hauteur trop importante de neige, elle s'coule
+    fqcalving(i) = max(0., snow(i) - snow_max)/dtime
+    snow(i)=min(snow(i),snow_max)
+!
     IF (nisurf == is_ter) then
       qsol(i) = qsol(i) + bil_eau_s(i)