Changeset 759 for LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histrac.h

Timestamp:

May 22, 2007, 3:55:22 PM (17 years ago)

Author:

lsce

Message:

correction de bug pour le couplage avec le modele inca - ACo

File:

• LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histrac.h (modified) (8 diffs)

LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histrac.h

@@ -62 +62 @@
 #ifdef INCA
 #ifdef INCAINFO
+#ifdef INCA_NMHC
       DO it=1, phtcnt
       WRITE(str2,'(i2.2)') it
@@ -85 +86 @@
      .                                   iim*(jjm+1)*klev,ndex3d)
       ENDDO
+#endif
+#ifdef INCA_CH4
+       DO it=1, phtcnt
+      WRITE(str2,'(i2.2)') it
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,it),
+     .     zx_tmp_3d)
+      CALL histwrite(nid_tra,"j"//str2,itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+      ENDDO
+
+      DO it=1, hetcnt
+      WRITE(str2,'(i2.2)') it
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,it),
+     .     zx_tmp_3d)
+      CALL histwrite(nid_tra,"w"//str2,itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+      ENDDO
+
+      DO it=1, extcnt
+      WRITE(str2,'(i2.2)') it
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,it),
+     .     zx_tmp_3d)
+      CALL histwrite(nid_tra,"ext"//str2,itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+      ENDDO
+#endif
 
       DO it=1, nfs
@@ -94 +121 @@
       ENDDO
 #else
+#ifdef INCA_NMHC
       CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2),
      .     zx_tmp_3d)
@@ -124 +152 @@
      .                                   iim*(jjm+1)*klev,ndex3d)
 #endif
-      DO it=1, grpcnt
-
-      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
-      zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it)
-      CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d,
-     .                                   iim*(jjm+1)*klev,ndex3d)
-      ENDDO
+#ifdef INCA_CH4
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2),
+     .     zx_tmp_3d)
+      CALL histwrite(nid_tra,"jO3",itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,4),
+     .     zx_tmp_3d)
+      CALL histwrite(nid_tra,"jNO2",itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,13),
+     .     zx_tmp_3d)
+      CALL histwrite(nid_tra,"jH2O2",itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,1),
+     .     zx_tmp_3d)
+      CALL histwrite(nid_tra,"wHNO3",itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,krates(1,1,1),
+     .     zx_tmp_3d)
+      CALL histwrite(nid_tra,"kN2O5",itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,1),
+     .     zx_tmp_3d)
+      CALL histwrite(nid_tra,"LghtNO",itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+#endif
+#endif
+!      DO it=1, grpcnt
+!
+!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
+!      zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it)
+!      CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d,
+!     .                                   iim*(jjm+1)*klev,ndex3d)
+!      ENDDO
 #endif
 
@@ -228 +287 @@
 #endif
        else
-#ifdef INCA
+#ifdef INCA_CH4
+       if ( solsym(it) == 'OX' ) then
+       zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3)
+       end if
+#endif
+#ifdef INCA_NMHC
        if ( solsym(it) == 'OX' ) then
        zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3)
@@ -264 +328 @@
 
 #ifdef INCA
-!#ifdef INCA_CH4
       CALL gr_fi_ecrit(1, klon,iim,jjm+1, o3_col(1), zx_tmp_2d)
       CALL histwrite(nid_tra,"O3_column",itau_w,zx_tmp_2d,
@@ -363 +426 @@
 #endif
 
+#ifdef INCA_CH4
+      DO it=1, grpcnt
+
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
+      zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it)
+      CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+      ENDDO
+
+#endif
+
+
 #ifdef INCA_NMHC
       CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_basprod(1,1),
@@ -398 +473 @@
       CALL histwrite(nid_tra,"CO_loss",itau_w,zx_tmp_3d,
      .                                   iim*(jjm+1)*klev,ndex)
+
+      DO it=1, grpcnt
+
+      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
+      zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it)
+      CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d,
+     .                                   iim*(jjm+1)*klev,ndex3d)
+      ENDDO
+
 
 #endif