Changeset 2632

Timestamp:

Sep 23, 2016, 10:15:55 AM (8 years ago)

Author:

fhourdin

Message:

Ajustement sur la version dusts.

Location:

LMDZ5/trunk/libf

Files:

• phy_common/iophys.F90 (added)
• phylmd/Dust/phytracr_spl_mod.F90 (modified) (23 diffs)
• phymar/iophys.F90 (deleted)

LMDZ5/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90

@@ -1007 +1007 @@
 
 
-      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 +1019 @@
       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 +2615 @@
 !      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 +2621 @@
          call iophys_ecrit('fav'//str2,1,'SOURCE','',source_tr(:,it))
       enddo
+#endif
 
       do it=1,nbtr
@@ -2643 +2649 @@
                                                      source_tr,flux_tr)    
 
+#ifdef IOPHYS_DUST
 
       do it=1,nbtr
@@ -2649 +2656 @@
          call iophys_ecrit('fap'//str2,1,'SOURCE','',source_tr(:,it))
       enddo
+#endif
 
       IF (lminmax) THEN
@@ -2667 +2675 @@
 !======================================================================
 !
+#ifdef IOPHYS_DUST
 !
       print *,'INPUT TO PRECUREMISSION'
@@ -2705 +2714 @@
          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 +2748 @@
 !                      EMISSIONS OF FINE AEROSOLS
 !=======================================================================
+#ifdef IOPHYS_DUST
 !
       do it=1,nbtr
@@ -2744 +2755 @@
          call iophys_ecrit('fpr'//str2,1,'SOURCE','',flux_tr(:,it))
       enddo
+#endif
 
       CALL finemission(zdz,pdtphys,zalt,kminbc,kmaxbc,                     & 
@@ -2784 +2796 @@
 
 
+#ifdef IOPHYS_DUST
       do it=1,nbtr
          write(str2,'(i2.2)') it
@@ -2789 +2802 @@
          call iophys_ecrit('f'//str2,1,'SOURCE','',flux_tr(:,it))
       enddo
+#endif
 !
 !
@@ -2842 +2856 @@
       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 +2898 @@
 !======================================================================
 
+#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 +2969 @@
 !======================================================================
 
+#ifdef IOPHYS_DUST
       print*,'iflag_conv=',iflag_conv
       call iophys_ecrit('coefh',klev,'coefh','',coefh)
@@ -2960 +2979 @@
          call iophys_ecrit('TRE'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 
@@ -3153 +3173 @@
 !
 !=======================================================================
+#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 +3267 @@
 !
 
+#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 +3442 @@
 !======================================================================
 !
+#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 +3625 @@
 !=======================================================================
 !     
+#ifdef IOPHYS_DUST
       call iophys_ecrit('da',klev,'da','',da)
       call iophys_ecrit('phi',klev,'phi','',phi)
@@ -3623 +3650 @@
          call iophys_ecrit('TRI'//str2,klev,'SOURCE','',tr_seri(:,:,it))
       enddo
+#endif
 
 
@@ -3736 +3764 @@
 !                         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 +4108 @@
 !  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