Changeset 521 for LMDZ.3.3

Timestamp:

Apr 30, 2004, 10:56:51 AM (20 years ago)

Author:

lmdzadmin

Message:

Suppression d'une sortie POLDER (qui venait de Johannes)
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (modified) (12 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

@@ -145 +145 @@
 c
 c
-      LOGICAL ok_polder ! sortir échantillonné de manière POLDER
-      save ok_polder
 c======================================================================
 c
@@ -956 +954 @@
       !        ok_ade=F -AIE=topswai-topsw
 
-      ! For POLDER swath
-      INTEGER pyr, pmo, pday    ! Year, month and day
-      INTEGER poldermask(klon)  ! POLDER swath mask (0 or 1)
-
       REAL aerindex(klon)       ! POLDER aerosol index
      
@@ -1184 +1178 @@
 #include "ini_histmth.h"
 
+
+#define histmthNMC
 #undef histmthNMC
-cccccccc#define histmthNMC
 #ifdef histmthNMC
 #include "ini_histmthNMC.h"
@@ -1191 +1186 @@
 
 #include "ini_histins.h"
-#include "ini_histpol.h"
-
+
+#undef histREGDYN
 #ifdef histREGDYN
 #include "ini_histREGDYN.h"
 #endif
 
+#undef histISCCP
 #ifdef histISCCP
 #include "ini_histISCCP.h"
@@ -1427 +1423 @@
 C Repartition du longwave par sous-surface linearisee
 Cn
+
        DO nsrf = 1, nbsrf
        DO i = 1, klon
@@ -1440 +1437 @@
 
 
-cIM   CALL clmain(dtime,itap,date0,pctsrf,
       CALL clmain(dtime,itap,date0,pctsrf,pctsrf_new,
      e            t_seri,q_seri,u_seri,v_seri,
@@ -1462 +1458 @@
 CXXX Incrementation des flux
 CXXX
+
       zxfluxt=0.
       zxfluxq=0.
@@ -1656 +1653 @@
           ELSE ! ok_cvl
 
-c          print*,'Avant conema OUI'
           CALL conema3 (dtime,
      .        paprs,pplay,t_seri,q_seri,
@@ -1668 +1664 @@
      .        ,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr
      .        ,clwcon0)
-          print*,'Apres conema3 '
 
           ENDIF ! ok_cvl
@@ -1724 +1719 @@
 c     CALL homogene(paprs, q_seri, d_q_con, u_seri,v_seri,
 c    .              d_u_con, d_v_con)
+
       DO k = 1, klev
         DO i = 1, klon
@@ -1836 +1832 @@
 c   -------------
       do k=1,klev
-         ratqss(:,k)=ratqsbas+(ratqshaut-ratqsbas)*
-     s   min((paprs(:,1)-pplay(:,k))/(paprs(:,1)-30000.),1.)        
+cIM RAJOUT boucle do=i
+       do i=1, klon
+cIM         ratqss(:,k)=ratqsbas+(ratqshaut-ratqsbas)*
+cIM     s   min((paprs(:,1)-pplay(:,k))/(paprs(:,1)-30000.),1.) 
+         ratqss(i,k)=ratqsbas+(ratqshaut-ratqsbas)*
+     s   min((paprs(i,1)-pplay(i,k))/(paprs(i,1)-30000.),1.) 
+cIM      print*,' IMratqs STABLE i, k',i,k,ratqss(i,k)
+       enddo 
       enddo
 
@@ -2760 +2762 @@
      .      t_ancien, q_ancien, rnebcon, ratqs, clwcon)
       ENDIF
+      
 
       RETURN