Changeset 204 for LMDZ.3.3/trunk/libf/phylmd/phystokenc.F

Timestamp:

Apr 13, 2001, 12:44:53 PM (23 years ago)

Author:

lmdz

Message:

Debogage du guidage et de la version debranchee et abandon de la version
debranchee non-netcdf FH/MAF
LF

File:

• LMDZ.3.3/trunk/libf/phylmd/phystokenc.F (modified) (11 diffs)

LMDZ.3.3/trunk/libf/phylmd/phystokenc.F

@@ -6 +6 @@
      I                   pt,pmfu, pmfd, pen_u, pde_u, pen_d, pde_d,
      I                   pcoefh,yu1,yv1,ftsol,pctsrf,
-     I                   frac_impa,frac_nucl,
+     I                   pfrac_impa,pfrac_nucl,
      I                   pphis,paire,dtime,itap)
       USE ioipsl
@@ -63 +63 @@
 c   ----------
 c
-      REAL frac_impa(klon,klev)
-      REAL frac_nucl(klon,klev)
+      REAL pfrac_impa(klon,klev)
+      REAL pfrac_nucl(klon,klev)
 c
 c Arguments necessaires pour les sources et puits de traceur
@@ -82 +82 @@
       REAL coefh(klon,klev) ! flux detraine dans le panache descendant
         REAL t(klon,klev)
+      REAL frac_impa(klon,klev)
+      REAL frac_nucl(klon,klev)
+      REAL rain(klon)
 
       REAL pyu1(klon),pyv1(klon)
@@ -95 +98 @@
       save t,mfu,mfd,en_u,de_u,en_d,de_d,coefh,dtcum
       save iadvtr,irec
+      save frac_impa,frac_nucl,rain
       save pyu1,pyv1,pftsol,ppsrf
 
@@ -133 +137 @@
          do k=1,klev
             do i=1,klon
+               frac_impa(i,k)=1.
+               frac_nucl(i,k)=1.
                mfu(i,k)=0.
                mfd(i,k)=0.
@@ -144 +150 @@
          enddo
          do i=1,klon
+            rain(i)=0.
             pyv1(i)=0.
             pyu1(i)=0.
@@ -159 +166 @@
       do k=1,klev
          do i=1,klon
+            frac_impa(i,k)=frac_impa(i,k)*pfrac_impa(i,k)
+            frac_nucl(i,k)=frac_nucl(i,k)*pfrac_nucl(i,k)
             mfu(i,k)=mfu(i,k)+pmfu(i,k)*pdtphys
             mfd(i,k)=mfd(i,k)+pmfd(i,k)*pdtphys
@@ -187 +196 @@
          do k=1,klev
             do i=1,klon
+c              frac_impa=frac_impa : c'est la fraction cumulee qu'on stoke
+c              frac_nucl=frac_nucl : c'est la fraction cumulee qu'on stoke
                mfu(i,k)=mfu(i,k)/dtcum
                mfd(i,k)=mfd(i,k)/dtcum
@@ -198 +209 @@
          enddo
          do i=1,klon
+            rain(i)=rain(i)/dtcum
             pyv1(i)=pyv1(i)/dtcum
             pyu1(i)=pyu1(i)/dtcum
@@ -226 +238 @@
 
 ccccc
+      print*,'AVANT ECRITURE'
          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t, zx_tmp_3d)
          CALL histwrite(physid,"t",itap,zx_tmp_3d,
      .                                   iim*(jjm+1)*klev,ndex)
+      print*,'APRES ECRITURE'
 
          CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfu, zx_tmp_3d)
@@ -279 +293 @@
          CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol4, zx_tmp_2d)
       CALL histwrite(physid,"ftsol4",itap,zx_tmp_2d,
+     .                                   iim*(jjm+1),ndex)
+
+        CALL gr_fi_ecrit(1,klon,iim,jjm+1, rain, zx_tmp_2d)
+      CALL histwrite(physid,"rain",itap,zx_tmp_2d,
      .                                   iim*(jjm+1),ndex)