Changeset 655 for LMDZ4/trunk/libf/phylmd/physiq.F

Timestamp:

Jun 22, 2005, 2:51:34 PM (19 years ago)

Author:

lmdzadmin

Message:

Correction dans l'appel de chemhook_end (INCA) AC
MAF

File:

• LMDZ4/trunk/libf/phylmd/physiq.F (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -2894 +2894 @@
 #include "calcul_divers.h"
 c
-c=============================================================
-c
-c Convertir les incrementations en tendances
-c
-      DO k = 1, klev
-      DO i = 1, klon
-         d_u(i,k) = ( u_seri(i,k) - u(i,k) ) / dtime
-         d_v(i,k) = ( v_seri(i,k) - v(i,k) ) / dtime
-         d_t(i,k) = ( t_seri(i,k)-t(i,k) ) / dtime
-         d_qx(i,k,ivap) = ( q_seri(i,k) - qx(i,k,ivap) ) / dtime
-         d_qx(i,k,iliq) = ( ql_seri(i,k) - qx(i,k,iliq) ) / dtime
-      ENDDO
-      ENDDO
-c
-      IF (nqmax.GE.3) THEN
-      DO iq = 3, nqmax
-      DO  k = 1, klev
-      DO  i = 1, klon
-         d_qx(i,k,iq) = ( tr_seri(i,k,iq-2) - qx(i,k,iq) ) / dtime
-      ENDDO
-      ENDDO
-      ENDDO
-      ENDIF
-c
-cIM rajout diagnostiques bilan KP pour analyse MJO par Jun-Ichi Yano
-c#include "write_bilKP_ins.h"
-c#include "write_bilKP_ave.h"
-c
-c Sauvegarder les valeurs de t et q a la fin de la physique:
-c
-      DO k = 1, klev
-      DO i = 1, klon
-         t_ancien(i,k) = t_seri(i,k)
-         q_ancien(i,k) = q_seri(i,k)
-      ENDDO
-      ENDDO
-c
-c 22.03.04 BEG
-c=============================================================
-c   Ecriture des sorties
-c=============================================================
-#ifdef CPP_IOIPSL
-
-#ifdef histhf
-#include "write_histhf.h"
-#endif
-
-#ifdef histday
-#include "write_histday.h"
-#include "write_histday_seri.h"
-#endif
-
-#ifdef histmth
-#include "write_histmth.h"
-#endif
-
-#ifdef histins
-#include "write_histins.h"
-#endif
-
-#ifdef histREGDYN
-#include "write_histREGDYN.h"
-#endif
-
-#ifdef histISCCP
-#include "write_histISCCP.h"
-#endif
-
-
-#ifdef histmthNMC
-#include "write_histmthNMC.h"
-#endif
-
-#endif
-
 #ifdef INCA
 #ifdef INCAINFO
@@ -2997 +2922 @@
            WRITE(lunout,*)'OK.'
 #endif
+#endif
+
+c=============================================================
+c
+c Convertir les incrementations en tendances
+c
+      DO k = 1, klev
+      DO i = 1, klon
+         d_u(i,k) = ( u_seri(i,k) - u(i,k) ) / dtime
+         d_v(i,k) = ( v_seri(i,k) - v(i,k) ) / dtime
+         d_t(i,k) = ( t_seri(i,k)-t(i,k) ) / dtime
+         d_qx(i,k,ivap) = ( q_seri(i,k) - qx(i,k,ivap) ) / dtime
+         d_qx(i,k,iliq) = ( ql_seri(i,k) - qx(i,k,iliq) ) / dtime
+      ENDDO
+      ENDDO
+c
+      IF (nqmax.GE.3) THEN
+      DO iq = 3, nqmax
+      DO  k = 1, klev
+      DO  i = 1, klon
+         d_qx(i,k,iq) = ( tr_seri(i,k,iq-2) - qx(i,k,iq) ) / dtime
+      ENDDO
+      ENDDO
+      ENDDO
+      ENDIF
+c
+cIM rajout diagnostiques bilan KP pour analyse MJO par Jun-Ichi Yano
+c#include "write_bilKP_ins.h"
+c#include "write_bilKP_ave.h"
+c
+c Sauvegarder les valeurs de t et q a la fin de la physique:
+c
+      DO k = 1, klev
+      DO i = 1, klon
+         t_ancien(i,k) = t_seri(i,k)
+         q_ancien(i,k) = q_seri(i,k)
+      ENDDO
+      ENDDO
+c
+c 22.03.04 BEG
+c=============================================================
+c   Ecriture des sorties
+c=============================================================
+#ifdef CPP_IOIPSL
+
+#ifdef histhf
+#include "write_histhf.h"
+#endif
+
+#ifdef histday
+#include "write_histday.h"
+#include "write_histday_seri.h"
+#endif
+
+#ifdef histmth
+#include "write_histmth.h"
+#endif
+
+#ifdef histins
+#include "write_histins.h"
+#endif
+
+#ifdef histREGDYN
+#include "write_histREGDYN.h"
+#endif
+
+#ifdef histISCCP
+#include "write_histISCCP.h"
+#endif
+
+
+#ifdef histmthNMC
+#include "write_histmthNMC.h"
+#endif
+
 #endif