Initialisation de ammonium a partir de obs.nml au lieu de bio_fasham.nml

src/extras/bio/bio_fasham.F90

@@ -299,8 +299,8 @@
 !      cc(n,i)=n_initial
 !      cc(n,i)=nit(i)					!CHG3
       cc(n,i)=nprof(i)                                  !CHG3
-!      cc(n,i)=aprof(i)                                  !CHG5
-      cc(a,i)=a_initial                                !CHG5
+      cc(n,i)=aprof(i)                                  !CHG5
+!      cc(a,i)=a_initial                                !CHG5
       cc(l,i)=l_initial
    end do
 

src/observations/observations.F90

@@ -44,10 +44,10 @@
    REALTYPE, public, dimension(:), allocatable   :: tprof
 
 !  'observed' nitrate profile
-   REALTYPE, public, dimension(:), allocatable   :: nprof,aprof         !CHG3
+   REALTYPE, public, dimension(:), allocatable   :: nprof         !CHG3
 
 !  'observed' ammonium profile
-!   REALTYPE, public, dimension(:), allocatable   :: aprof         !CHG5
+   REALTYPE, public, dimension(:), allocatable   :: aprof         !CHG5
 
 !  'observed' oxygen profile
    REALTYPE, public, dimension(:), allocatable   :: o2_prof
@@ -504,7 +504,7 @@
    NRelaxTau = 0.
 
    allocate(ARelaxTau(0:nlev),stat=rc)                                  !CHG5
-   if (rc /= 0) STOP 'init_observations: Error allocating (NRelaxTau)'
+   if (rc /= 0) STOP 'init_observations: Error allocating (ARelaxTau)'
    ARelaxTau = 0.
 
    allocate(uprof(0:nlev),stat=rc)