Index: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/albedo.F
===================================================================
--- LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/albedo.F	(revision 898)
+++ LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/albedo.F	(revision 899)
@@ -75,5 +75,5 @@
       ENDDO
       IF (srmu .NE. 0.0) THEN
-         albedo(i) = salb / srmu * fmagic
+         albedo(i) = salb / srmu * fmagic+pmagic
       ELSE ! nuit polaire (on peut prendre une valeur quelconque)
          albedo(i) = fmagic
@@ -120,5 +120,5 @@
       ENDDO
       IF (srmu .NE. 0.0) THEN
-         albedo(i) = salb / srmu * fmagic
+         albedo(i) = salb / srmu * fmagic+pmagic
       ELSE ! nuit polaire (on peut prendre une valeur quelconque)
          albedo(i) = fmagic
@@ -170,5 +170,5 @@
 c
          fauxo = ( 1.47 - ACOS( rmu0(i) ) )/0.15
-         albedo(i) = fmagic*( .03 + .630/( 1. + fauxo*fauxo))
+         albedo(i) = fmagic*( .03 + .630/( 1. + fauxo*fauxo))+pmagic
          albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
       ENDDO
@@ -181,5 +181,5 @@
          rmu0(i) = MAX(rmu0(i),0.0)
 cIM:orig albedo(i) = 0.058/(rmu0(i) + 0.30)
-         albedo(i) = fmagic * 0.058/(rmu0(i) + 0.30)
+         albedo(i) = fmagic * 0.058/(rmu0(i) + 0.30)+pmagic
          albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
       ENDDO
Index: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/clesphys.h
===================================================================
--- LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/clesphys.h	(revision 898)
+++ LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/clesphys.h	(revision 899)
@@ -31,6 +31,7 @@
 !IM ok_kzmin : clef calcul Kzmin dans la CL de surface cf FH
        LOGICAL ok_kzmin
-!IM fmagic : parametre pour regler l albedo sur ocean
-       REAL fmagic
+!IM, MAFo fmagic, pmagic : parametres - additionnel et multiplicatif - 
+!                          pour regler l albedo sur ocean
+       REAL fmagic, pmagic
 !IM lev_histhf  : niveau sorties 6h
 !IM lev_histday : niveau sorties journalieres
@@ -53,5 +54,6 @@
      &     , CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt                     &
      &     , top_height, overlap, cdmmax, cdhmax, ksta, ksta_ter        &
-     &     , ok_kzmin, fmagic, lev_histhf, lev_histday, lev_histmth     &
+     &     , ok_kzmin, fmagic, pmagic                                   &
+     &     , lev_histhf, lev_histday, lev_histmth                       &
      &     , type_run, ok_isccp, ok_regdyn                              &
      &     , lonmin_ins, lonmax_ins, latmin_ins, latmax_ins             &
Index: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/conf_phys.F90
===================================================================
--- LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/conf_phys.F90	(revision 898)
+++ LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/conf_phys.F90	(revision 899)
@@ -83,5 +83,5 @@
   REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp
   LOGICAL,SAVE :: ok_kzmin_omp
-  REAL, SAVE ::  fmagic_omp
+  REAL, SAVE ::  fmagic_omp, pmagic_omp
   INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
   CHARACTER*4, SAVE :: type_run_omp
@@ -603,10 +603,20 @@
 !
 !Config Key  = fmagic
-!Config Desc = 
+!Config Desc = additionnal multiplicator factor used for albedo
 !Config Def  = 1.
-!Config Help = Used in albedo.F
+!Config Help = additionnal multiplicator factor used in albedo.F
 !
   fmagic_omp = 1.
   call getin('fmagic',fmagic_omp)
+
+!
+!Config Key  = pmagic
+!Config Desc = additional factor used for albedo
+!Config Def  = 0.
+!Config Help = additional factor used in albedo.F
+!
+  pmagic_omp = 0.
+  call getin('pmagic',pmagic_omp)
+
 
 !Config Key = ok_lic_melt
@@ -822,4 +832,5 @@
     ok_kzmin = ok_kzmin_omp
     fmagic = fmagic_omp
+    pmagic = pmagic_omp
     iflag_pbl = iflag_pbl_omp
     lev_histhf = lev_histhf_omp
@@ -915,4 +926,5 @@
   write(numout,*)' ok_kzmin = ',ok_kzmin 
   write(numout,*)' fmagic = ',fmagic
+  write(numout,*)' pmagic = ',pmagic
   write(numout,*)' ok_ade = ',ok_ade
   write(numout,*)' ok_aie = ',ok_aie