Index: trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F	(revision 3910)
+++ trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F	(revision 3916)
@@ -10,5 +10,5 @@
      &                     zmea, zstd, zsig, zgam, zthe,
      &                     emissiv, emisice, iceradius,
-     &                     dtemisice
+     &                     dtemisice, n2frac
       use comdiurn_h, only: sinlat, coslat, sinlon, coslon
       use comsoil_h, only: nsoilmx, layer, mlayer, inertiedat, volcapa
@@ -851,4 +851,8 @@
       write(*,*) " inertia = ",inertiedat(1,1)
 
+c  Initialisation n2frac
+c  --------------------------------------
+      n2frac(1)=1. ! default value for n2frac
+
 ! Initialize soil properties and temperature
 ! ------------------------------------------
@@ -920,5 +924,6 @@
       call physdem1("startfi.nc",nsoilmx,1,llm,nq,
      &                dtphys,time,
-     &                tsurf,tsoil,inertiedat,emis,albedodat,q2,qsurf)
+     &                tsurf,tsoil,inertiedat,emis,albedodat,q2,qsurf,
+     &                n2frac)
 
 c=======================================================================
Index: trunk/LMDZ.PLUTO/libf/phypluto/phyetat0_mod.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/phyetat0_mod.F90	(revision 3910)
+++ trunk/LMDZ.PLUTO/libf/phypluto/phyetat0_mod.F90	(revision 3916)
@@ -148,9 +148,10 @@
   call get_field(nid_start,"n2frac",n2frac,found)
   if (.not.found) then
-    call abort_physic(modname,"Failed loading <n2frac>",1)
+    n2frac(:)=1.
+    !call abort_physic(modname,"Failed loading <n2frac>",1)
   endif
 else
   ! If no startfi file, use fixed value to test
-  n2frac(:)=1
+  n2frac(:)=1.
 endif ! of if (startphy_file)
 if (is_master) write(*,*) "phyetat0: n2 fraction of ice  <n2frac> range:", &