Changeset 879 for LMDZ4/trunk/libf/phylmd/write_histmth.h

Timestamp:

Jan 17, 2008, 10:20:32 PM (16 years ago)

Author:

Laurent Fairhead

Message:

Suite de la bascule vers une physique avec thermiques, nouvelle convection, poche froide ...
LF

File:

• LMDZ4/trunk/libf/phylmd/write_histmth.h (modified) (7 diffs)

LMDZ4/trunk/libf/phylmd/write_histmth.h

@@ -9 +9 @@
          itau_w = itau_phy + itap
 
+c  temperature tendency due to moist convective processes
+       DO l=1, klev
+       DO i=1, klon
+        zx_tmp_fi3d(i,l)=d_t_con(i,l)/pdtphys
+       ENDDO !i
+       ENDDO !l
+c
+
+cIM: 101003 : K/30min ==> K/s
+      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
+cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
+      CALL histwrite_phy(nid_mth,"dtthe",itau_w,zx_tmp_fi3d)
+c
+cIM: 101003 : K/30min ==> K/s
+      zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
+cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
+      CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
+c
 c
       IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
@@ -44 +62 @@
 c ENSEMBLES BEG
 cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
-      CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
+c     CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
 c
 cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
-      CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
+c     CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
 c
 c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
@@ -103 +121 @@
 c
 cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
-      CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
+c     CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
 c
       DO i = 1, klon
@@ -968 +986 @@
       CALL histwrite_phy(nid_mth,"dtlw",itau_w,zx_tmp_fi3d)
 c
-c  temperature tendency due to moist convective processes
-       DO l=1, klev
-       DO i=1, klon
-        zx_tmp_fi3d(i,l)=d_t_con(i,l)/pdtphys
-       ENDDO !i
-       ENDDO !l
-c
-cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
-      CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
+
+cIM: 101003 : K/30min ==> K/s
+      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
+cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
+      CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
+
+c
+      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
+cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
+      CALL histwrite_phy(nid_mth,"dqthe",itau_w,zx_tmp_fi3d)
 c
 
@@ -1078 +1097 @@
 c ENSEMBLES BEG
 cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
-      CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
+c     CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
 c
 cym     CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
-      CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
+c     CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
 c
 c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
@@ -1137 +1156 @@
 c
 cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
-      CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
+c     CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
 c
       DO i = 1, klon
@@ -1806 +1825 @@
 c
 cIM: 101003 : K/30min ==> K/s
-      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
-cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
-      CALL histwrite_phy(nid_mth,"dtthe",itau_w,zx_tmp_fi3d)
 c
       zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys