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

Timestamp:

Dec 15, 2003, 6:50:41 PM (20 years ago)

Author:

lmdzadmin

Message:

Phasage avec la version de Ionela
IM/LF

File:

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

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

@@ -694 +694 @@
      alb_new(1 : knon)  = alb_neig(1 : knon)*zfra(1:knon) + &
     &                     0.6 * (1.0-zfra(1:knon))
-!!     alb_new(1 : knon)  = 0.6
+!IM cf FH/GK     alb_new(1 : knon)  = 0.6
+!       alb_new(1 : knon)  = 0.82
+!IM cf JLD/ GK
+!IM: 211003 Ksta0.77      alb_new(1 : knon)  = 0.77
+!IM: KstaTER0.8 & LMD_ARMIP5    alb_new(1 : knon)  = 0.8
+!IM: KstaTER0.77 & LMD_ARMIP6    
+        alb_new(1 : knon)  = 0.77
+
 !
 ! Rugosite
@@ -1144 +1151 @@
 
 !IM cf. JP +++
-    albedo_keep(:) = (albedo_out(:,1)+albedo_out(:,2))/2.
+!IM BUG BUG BUG albedo_keep(:) = (albedo_out(:,1)+albedo_out(:,2))/2.
+    albedo_keep(1:knon) = (albedo_out(1:knon,1)+albedo_out(1:knon,2))/2.
 !IM cf. JP ---
 
@@ -2599 +2607 @@
 ! Masse maximum de neige (kg/m2). Au dessus de ce seuil, la neige
 ! en exces "s'ecoule" (calving)
-  real, parameter :: snow_max=1.
+!  real, parameter :: snow_max=1.
+!IM cf JLD/GK
+  real, parameter :: snow_max=3000.
   integer :: i
   real, dimension(klon) :: zx_mh, zx_nh, zx_oh
@@ -2614 +2624 @@
 ! REAL, parameter :: chasno = RLMLT/(2.3867E+06*0.15)
   REAL, parameter :: chasno = 3.334E+05/(2.3867E+06*0.15)
+!IM cf JLD/ GKtest
+  REAL, parameter :: chaice = 3.334E+05/(2.3867E+06*0.15)
+! fin GKtest
 !
   logical, save         :: check = .FALSE.
@@ -2702 +2715 @@
       tsurf_new(i) = tsurf_new(i) - fq_fonte * chasno  
 !IM cf JLD OK     
-      IF (nisurf == is_sic .OR. nisurf == is_lic ) tsurf_new(i) = RTT
+!IM cf JLD/ GKtest fonte aussi pour la glace
+!      IF (nisurf == is_sic .OR. nisurf == is_lic ) tsurf_new(i) = RTT
+      IF (nisurf == is_sic .OR. nisurf == is_lic ) THEN
+        fq_fonte = MAX((tsurf_new(i)-RTT )/chaice,0.0)
+        ffonte(i) = ffonte(i) + fq_fonte * RLMLT/dtime
+        bil_eau_s(i) = bil_eau_s(i) + fq_fonte
+        tsurf_new(i) = RTT
+      ENDIF
+! fin GKtest
       d_ts(i) = tsurf_new(i) - tsurf(i)
 !      zx_h_ts(i) = tsurf_new(i) * RCPD * zx_pkh(i)
@@ -2729 +2750 @@
       run_off(i) = run_off(i) + MAX(qsol(i) - max_eau_sol, 0.0)
       qsol(i) = MIN(qsol(i), max_eau_sol) 
-    else
-      run_off(i) = run_off(i) + MAX(bil_eau_s(i), 0.0)
+!IM : 0601003 else
+!IM: run_off(i)  
+!IM : 061003   run_off(i) = run_off(i) + MAX(bil_eau_s(i), 0.0)
     endif
   enddo