Changeset 1055 for LMDZ4/trunk

Timestamp:

Nov 14, 2008, 5:57:39 PM (16 years ago)

Author:

lmdzadmin

Message:

Correction sorties niveaux pression
AI/IM

Location:

LMDZ4/trunk/libf/phylmd

Files:

• phys_output_mod.F90 (modified) (2 diffs)
• phys_output_write.h (modified) (1 diff)
• physiq.F (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/phys_output_mod.F90

@@ -594 +594 @@
 !      on ecrit u v t q a 850 700 500 200 au niv 3
 
+   zstophym(iff) = ecrit_files(iff)
+   type_ecri(1) = 'inst(X)'
+   type_ecri(2) = 'inst(X)'
+   type_ecri(3) = 'inst(X)'
+   type_ecri(4) = 'inst(X)'
         DO k=1, nlevSTD
      IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
@@ -606 +611 @@
      ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
        ENDDO
+   zstophym(iff) = dtime
+   type_ecri(1) = 'ave(X)'
+   type_ecri(2) = 'ave(X)'
+   type_ecri(3) = 'ave(X)'
+   type_ecri(4) = 'inst(X)'
 
  CALL histdef2d(iff,flag_t_oce_sic,"t_oce_sic","Temp mixte oce-sic","K")

LMDZ4/trunk/libf/phylmd/phys_output_write.h

@@ -567 +567 @@
 ! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
        IF (flag_ulevsSTD(iff)<=lev_files(iff)) THEN
-       CALL histwrite_phy(nid_files(iff),"u"//bb2,  
-     & itau_w,ulevSTD(:,k))
+       CALL histwrite_phy(nid_files(iff),"u"//bb2,
+     & itau_w,uwriteSTD(:,k,iff))
        ENDIF
 
        IF (flag_vlevsSTD(iff)<=lev_files(iff)) THEN
       CALL histwrite_phy(nid_files(iff),"v"//bb2,  
-     & itau_w,vlevSTD(:,k))
+     & itau_w,vwriteSTD(:,k,iff))
        ENDIF
 
        IF (flag_wlevsSTD(iff)<=lev_files(iff)) THEN
-      CALL histwrite_phy(nid_files(iff),"w"//bb2,  
-     & itau_w,wlevSTD(:,k))
+      CALL histwrite_phy(nid_files(iff),"w"//bb2,
+     & itau_w,wwriteSTD(:,k,iff))
        ENDIF
 
        IF (flag_philevsSTD(iff)<=lev_files(iff)) THEN
       CALL histwrite_phy(nid_files(iff),
-     $                  "phi"//bb2,  
-     &               itau_w,philevSTD(:,k))
+     $                  "phi"//bb2,
+     &               itau_w,phiwriteSTD(:,k,iff))
        ENDIF
 
        IF (flag_qlevsSTD(iff)<=lev_files(iff)) THEN
-      CALL histwrite_phy(nid_files(iff),"q"//bb2,  
-     & itau_w, qlevSTD(:,k))
+      CALL histwrite_phy(nid_files(iff),"q"//bb2,
+     & itau_w, qwriteSTD(:,k,iff))
        ENDIF
 
        IF (flag_tlevsSTD(iff)<=lev_files(iff)) THEN
-      CALL histwrite_phy(nid_files(iff),"t"//bb2,  
-     & itau_w, tlevSTD(:,k))
+      CALL histwrite_phy(nid_files(iff),"t"//bb2,
+     & itau_w, twriteSTD(:,k,iff))
        ENDIF
 

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -263 +263 @@
       real ulevSTD(klon,nlevSTD), vlevSTD(klon,nlevSTD)
       real wlevSTD(klon,nlevSTD) 
+
+      real twriteSTD(klon,nlevSTD,nfiles)
+      real qwriteSTD(klon,nlevSTD,nfiles)
+      real rhwriteSTD(klon,nlevSTD,nfiles)
+      real phiwriteSTD(klon,nlevSTD,nfiles)
+      real uwriteSTD(klon,nlevSTD,nfiles)
+      real vwriteSTD(klon,nlevSTD,nfiles)
+      real wwriteSTD(klon,nlevSTD,nfiles)
 c
 c nout : niveau de output des variables a une pression donnee
@@ -3315 +3323 @@
 c
 #include "calcul_STDlev.h"
+      twriteSTD(:,:,1)=tsumSTD(:,:,2)
+      qwriteSTD(:,:,1)=qsumSTD(:,:,2)
+      rhwriteSTD(:,:,1)=rhsumSTD(:,:,2)
+      phiwriteSTD(:,:,1)=phisumSTD(:,:,2)
+      uwriteSTD(:,:,1)=usumSTD(:,:,2)
+      vwriteSTD(:,:,1)=vsumSTD(:,:,2)
+      wwriteSTD(:,:,1)=wsumSTD(:,:,2)
+
+      twriteSTD(:,:,2)=tsumSTD(:,:,1)
+      qwriteSTD(:,:,2)=qsumSTD(:,:,1)
+      rhwriteSTD(:,:,2)=rhsumSTD(:,:,1)
+      phiwriteSTD(:,:,2)=phisumSTD(:,:,1)
+      uwriteSTD(:,:,2)=usumSTD(:,:,1)
+      vwriteSTD(:,:,2)=vsumSTD(:,:,1)
+      wwriteSTD(:,:,2)=wsumSTD(:,:,1)
+
+      twriteSTD(:,:,3)=tlevSTD(:,:)
+      qwriteSTD(:,:,3)=qlevSTD(:,:)
+      rhwriteSTD(:,:,3)=rhlevSTD(:,:)
+      phiwriteSTD(:,:,3)=philevSTD(:,:)
+      uwriteSTD(:,:,3)=ulevSTD(:,:)
+      vwriteSTD(:,:,3)=vlevSTD(:,:)
+      wwriteSTD(:,:,3)=wlevSTD(:,:)
+
+      twriteSTD(:,:,4)=tlevSTD(:,:)
+      qwriteSTD(:,:,4)=qlevSTD(:,:)
+      rhwriteSTD(:,:,4)=rhlevSTD(:,:)
+      phiwriteSTD(:,:,4)=philevSTD(:,:)
+      uwriteSTD(:,:,4)=ulevSTD(:,:)
+      vwriteSTD(:,:,4)=vlevSTD(:,:)
+      wwriteSTD(:,:,4)=wlevSTD(:,:)
 c
 c slp sea level pressure