Index: LMDZ5/trunk/libf/phy_common/iophys.F90
===================================================================
--- LMDZ5/trunk/libf/phy_common/iophys.F90	(revision 2632)
+++ LMDZ5/trunk/libf/phy_common/iophys.F90	(revision 2632)
@@ -0,0 +1,132 @@
+      subroutine iophys_ecrit(nom,lllm,titre,unite,px)
+      USE dimphy
+      USE mod_phys_lmdz_para
+      USE mod_grid_phy_lmdz
+      IMPLICIT NONE
+
+
+
+!  Ecriture de variables diagnostiques au choix dans la physique 
+!  dans un fichier NetCDF nomme  'diagfi'. Ces variables peuvent etre
+!  3d (ex : temperature), 2d (ex : temperature de surface), ou
+!  0d (pour un scalaire qui ne depend que du temps : ex : la longitude
+!  solaire)
+!  (ou encore 1d, dans le cas de testphys1d, pour sortir une colonne)
+!  La periode d'ecriture est donnee par 
+!  "ecritphy " regle dans le fichier de controle de run :  run.def
+!
+!    writediagfi peut etre appele de n'importe quelle subroutine
+!    de la physique, plusieurs fois. L'initialisation et la creation du
+!    fichier se fait au tout premier appel.
+!
+! WARNING : les variables dynamique (u,v,t,q,ps)
+!  sauvees par writediagfi avec une
+! date donnee sont legerement differentes que dans le fichier histoire car 
+! on ne leur a pas encore ajoute de la dissipation et de la physique !!!
+! IL est  RECOMMANDE d'ajouter les tendance physique a ces variables
+! avant l'ecriture dans diagfi (cf. physiq.F)
+!  
+! Modifs: Aug.2010 Ehouarn: enforce outputs to be real*4
+!
+!  parametres (input) :
+!  ----------
+!      unit : unite logique du fichier de sortie (toujours la meme)
+!      nom  : nom de la variable a sortir (chaine de caracteres)
+!      titre: titre de la variable (chaine de caracteres)
+!      unite : unite de la variable (chaine de caracteres)
+!      px : variable a sortir (real 0, 1, 2, ou 3d)
+!
+!=================================================================
+
+#include "dimensions.h"
+#include "paramet.h"
+#include "netcdf.inc"
+#include "iotd.h"
+
+
+! Arguments on input:
+      integer lllm
+      character (len=*) :: nom,titre,unite
+      integer imjmax
+      parameter (imjmax=100000)
+      real px(klon_omp,lllm)
+      real xglo(klon_glo,lllm)
+      real zx(iim,jjp1,lllm)
+
+
+      CALL Gather(px,xglo)
+!$OMP MASTER
+      IF (is_mpi_root) THEN       
+        CALL Grid1Dto2D_glo(xglo,zx)
+        call iotd_ecrit(nom,lllm,titre,unite,zx)
+      ENDIF
+!$OMP END MASTER
+
+      return
+      end
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+      SUBROUTINE iophys_ini
+      USE mod_phys_lmdz_para
+      USE vertical_layers_mod, ONLY: presnivs
+      USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
+      USE dimphy, ONLY: klev
+
+      IMPLICIT NONE
+
+!=======================================================================
+!
+!   Auteur:  L. Fairhead  ,  P. Le Van, Y. Wanherdrick, F. Forget
+!   -------
+!
+!   Objet:
+!   ------
+!
+!   'Initialize' the diagfi.nc file: write down dimensions as well
+!   as time-independent fields (e.g: geopotential, mesh area, ...)
+!
+!=======================================================================
+!-----------------------------------------------------------------------
+!   Declarations:
+!   -------------
+
+real pi
+
+!   Arguments:
+!   ----------
+
+!$OMP MASTER
+    IF (is_mpi_root) THEN       
+pi=2.*asin(1.)
+call iotd_ini('phys.nc   ', &
+!iim,jjp1,llm,rlonv(1:iim)*180./pi,rlatu*180./pi,presnivs)
+size(lon_reg),size(lat_reg),klev,lon_reg(:)*180./pi,lat_reg*180./pi,presnivs)
+    ENDIF
+!$OMP END MASTER
+
+      END
+
+#ifdef und
+      SUBROUTINE gr_fi_ecrit(nfield,nlon,iim,jjmp1,fi,ecrit)
+      IMPLICIT none
+
+!=======================================================================
+      INTEGER nfield,nlon,iim,jjmp1, jjm
+      REAL fi(nlon,nfield), ecrit(iim*jjmp1,nfield)
+
+      INTEGER i, n, ig
+
+      jjm = jjmp1 - 1
+      DO n = 1, nfield
+         DO i=1,iim
+            ecrit(i,n) = fi(1,n)
+            ecrit(i+jjm*iim,n) = fi(nlon,n)
+         ENDDO
+         DO ig = 1, nlon - 2
+           ecrit(iim+ig,n) = fi(1+ig,n)
+         ENDDO
+      ENDDO
+      RETURN
+      END
+
+#endif
Index: LMDZ5/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90	(revision 2630)
+++ LMDZ5/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90	(revision 2632)
@@ -1007,12 +1007,9 @@
 
 
-      do it=1,nbtr
-         write(str2,'(i2.2)') it
-         call iophys_ecrit('TRA'//str2,klev,'SOURCE','',tr_seri(:,:,it))
-      enddo
-
 
       if (debutphy) then
+#ifdef IOPHYS_DUST
          CALL iophys_ini
+#endif
          nbreg_ind=1
          nbreg_bb=1
@@ -1022,4 +1019,11 @@
       endif
 
+
+#ifdef IOPHYS_DUST
+      do it=1,nbtr
+         write(str2,'(i2.2)') it
+         call iophys_ecrit('TRA'//str2,klev,'SOURCE','',tr_seri(:,:,it))
+      enddo
+#endif
 
   
@@ -2611,4 +2615,5 @@
 !      print *,'tr_seri = ',SUM(tr_seri(:,:,3)),MINVAL(tr_seri(:,:,3)),
 !     .                                         MAXVAL(tr_seri(:,:,3))
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
@@ -2616,4 +2621,5 @@
          call iophys_ecrit('fav'//str2,1,'SOURCE','',source_tr(:,it))
       enddo
+#endif
 
       do it=1,nbtr
@@ -2643,4 +2649,5 @@
                                                      source_tr,flux_tr)    
 
+#ifdef IOPHYS_DUST
 
       do it=1,nbtr
@@ -2649,4 +2656,5 @@
          call iophys_ecrit('fap'//str2,1,'SOURCE','',source_tr(:,it))
       enddo
+#endif
 
       IF (lminmax) THEN
@@ -2667,4 +2675,5 @@
 !======================================================================
 !
+#ifdef IOPHYS_DUST
 !
       print *,'INPUT TO PRECUREMISSION'
@@ -2705,4 +2714,5 @@
          call iophys_ecrit('flux_sparam_ind',1,'flux_sparam_ind','',flux_sparam_ind)
          call iophys_ecrit('flux_sparam_bb',1,'flux_sparam_bb','',flux_sparam_bb)
+#endif
 
 
@@ -2738,4 +2748,5 @@
 !                      EMISSIONS OF FINE AEROSOLS
 !=======================================================================
+#ifdef IOPHYS_DUST
 !
       do it=1,nbtr
@@ -2744,4 +2755,5 @@
          call iophys_ecrit('fpr'//str2,1,'SOURCE','',flux_tr(:,it))
       enddo
+#endif
 
       CALL finemission(zdz,pdtphys,zalt,kminbc,kmaxbc,                     & 
@@ -2784,4 +2796,5 @@
 
 
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
@@ -2789,4 +2802,5 @@
          call iophys_ecrit('f'//str2,1,'SOURCE','',flux_tr(:,it))
       enddo
+#endif
 !
 !
@@ -2842,8 +2856,10 @@
       ENDIF
 
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
          call iophys_ecrit('TRC'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
       CALL deposition(vdep_oce,vdep_sic,vdep_ter,vdep_lic,pctsrf,      &
@@ -2882,8 +2898,10 @@
 !======================================================================
 
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
          call iophys_ecrit('TRD'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 
@@ -2951,4 +2969,5 @@
 !======================================================================
 
+#ifdef IOPHYS_DUST
       print*,'iflag_conv=',iflag_conv
       call iophys_ecrit('coefh',klev,'coefh','',coefh)
@@ -2960,4 +2979,5 @@
          call iophys_ecrit('TRE'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 
@@ -3153,8 +3173,10 @@
 !
 !=======================================================================
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
          call iophys_ecrit('TRF'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 
@@ -3245,8 +3267,10 @@
 !
 
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
          call iophys_ecrit('TRG'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 
@@ -3418,8 +3442,10 @@
 !======================================================================
 !
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
          call iophys_ecrit('TRH'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 
@@ -3599,4 +3625,5 @@
 !=======================================================================
 !     
+#ifdef IOPHYS_DUST
       call iophys_ecrit('da',klev,'da','',da)
       call iophys_ecrit('phi',klev,'phi','',phi)
@@ -3623,4 +3650,5 @@
          call iophys_ecrit('TRI'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 
@@ -3736,8 +3764,10 @@
 !                         COMPUTING THE BURDEN
 !=======================================================================
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
          call iophys_ecrit('TRJ'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 !   
@@ -4078,8 +4108,10 @@
 !  Stockage sur bande histoire
 !======================================================================
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
          call iophys_ecrit('TRK'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 
Index: LMDZ5/trunk/libf/phymar/iophys.F90
===================================================================
--- LMDZ5/trunk/libf/phymar/iophys.F90	(revision 2630)
+++ 	(revision )
@@ -1,132 +1,0 @@
-      subroutine iophys_ecrit(nom,lllm,titre,unite,px)
-      USE dimphy
-      USE mod_phys_lmdz_para
-      USE mod_grid_phy_lmdz
-      IMPLICIT NONE
-
-
-
-!  Ecriture de variables diagnostiques au choix dans la physique 
-!  dans un fichier NetCDF nomme  'diagfi'. Ces variables peuvent etre
-!  3d (ex : temperature), 2d (ex : temperature de surface), ou
-!  0d (pour un scalaire qui ne depend que du temps : ex : la longitude
-!  solaire)
-!  (ou encore 1d, dans le cas de testphys1d, pour sortir une colonne)
-!  La periode d'ecriture est donnee par 
-!  "ecritphy " regle dans le fichier de controle de run :  run.def
-!
-!    writediagfi peut etre appele de n'importe quelle subroutine
-!    de la physique, plusieurs fois. L'initialisation et la creation du
-!    fichier se fait au tout premier appel.
-!
-! WARNING : les variables dynamique (u,v,t,q,ps)
-!  sauvees par writediagfi avec une
-! date donnee sont legerement differentes que dans le fichier histoire car 
-! on ne leur a pas encore ajoute de la dissipation et de la physique !!!
-! IL est  RECOMMANDE d'ajouter les tendance physique a ces variables
-! avant l'ecriture dans diagfi (cf. physiq.F)
-!  
-! Modifs: Aug.2010 Ehouarn: enforce outputs to be real*4
-!
-!  parametres (input) :
-!  ----------
-!      unit : unite logique du fichier de sortie (toujours la meme)
-!      nom  : nom de la variable a sortir (chaine de caracteres)
-!      titre: titre de la variable (chaine de caracteres)
-!      unite : unite de la variable (chaine de caracteres)
-!      px : variable a sortir (real 0, 1, 2, ou 3d)
-!
-!=================================================================
-
-#include "dimensions.h"
-#include "paramet.h"
-#include "netcdf.inc"
-#include "iotd.h"
-
-
-! Arguments on input:
-      integer lllm
-      character (len=*) :: nom,titre,unite
-      integer imjmax
-      parameter (imjmax=100000)
-      real px(klon_omp,lllm)
-      real xglo(klon_glo,lllm)
-      real zx(iim,jjp1,lllm)
-
-
-      CALL Gather(px,xglo)
-!$OMP MASTER
-      IF (is_mpi_root) THEN       
-        CALL Grid1Dto2D_glo(xglo,zx)
-        call iotd_ecrit(nom,lllm,titre,unite,zx)
-      ENDIF
-!$OMP END MASTER
-
-      return
-      end
-
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-      SUBROUTINE iophys_ini
-      USE mod_phys_lmdz_para
-      USE vertical_layers_mod, ONLY: presnivs
-      USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
-      USE dimphy, ONLY: klev
-
-      IMPLICIT NONE
-
-!=======================================================================
-!
-!   Auteur:  L. Fairhead  ,  P. Le Van, Y. Wanherdrick, F. Forget
-!   -------
-!
-!   Objet:
-!   ------
-!
-!   'Initialize' the diagfi.nc file: write down dimensions as well
-!   as time-independent fields (e.g: geopotential, mesh area, ...)
-!
-!=======================================================================
-!-----------------------------------------------------------------------
-!   Declarations:
-!   -------------
-
-real pi
-
-!   Arguments:
-!   ----------
-
-!$OMP MASTER
-    IF (is_mpi_root) THEN       
-pi=2.*asin(1.)
-call iotd_ini('phys.nc   ', &
-!iim,jjp1,llm,rlonv(1:iim)*180./pi,rlatu*180./pi,presnivs)
-size(lon_reg),size(lat_reg),klev,lon_reg(:)*180./pi,lat_reg*180./pi,presnivs)
-    ENDIF
-!$OMP END MASTER
-
-      END
-
-#ifdef und
-      SUBROUTINE gr_fi_ecrit(nfield,nlon,iim,jjmp1,fi,ecrit)
-      IMPLICIT none
-
-!=======================================================================
-      INTEGER nfield,nlon,iim,jjmp1, jjm
-      REAL fi(nlon,nfield), ecrit(iim*jjmp1,nfield)
-
-      INTEGER i, n, ig
-
-      jjm = jjmp1 - 1
-      DO n = 1, nfield
-         DO i=1,iim
-            ecrit(i,n) = fi(1,n)
-            ecrit(i+jjm*iim,n) = fi(nlon,n)
-         ENDDO
-         DO ig = 1, nlon - 2
-           ecrit(iim+ig,n) = fi(1+ig,n)
-         ENDDO
-      ENDDO
-      RETURN
-      END
-
-#endif