Changeset 918 for LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/physiq.F

Timestamp:

Feb 27, 2008, 4:24:26 PM (17 years ago)

Author:

Laurent Fairhead

Message:

• correction du bug ISCCP (n'ecrire ptop que quand ISCCP est appelé)
• petite inversion de boucle dans isccp_cloud_types.F pour aller + vite
• "CFisation" d'un certain nombre d'unités pour les hist*
• les suggestions de JL pour rugoro

SD
LF

File:

• LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/physiq.F (modified) (4 diffs)

LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/physiq.F

@@ -47 +47 @@
 c   CLEFS CPP POUR LES IO
 c   =====================
-c#define histhf
+#define histhf
 #define histday
 #define histmth
 c#define histins
-c#define histmthNMC
-c#define histISCCP
+#define histmthNMC
+#define histISCCP
 c======================================================================
 c    modif   ( P. Le Van ,  12/10/98 )
@@ -1487 +1487 @@
 cym Attention pbase pas initialise dans concvl !!!!
         pbase(:)=0
-        
         first=.false.
       endif 
@@ -1675 +1674 @@
              rugoro(i) = MAX(1.0e-05, zstd(i)*zsig(i)/2.0)
            ENDDO
+
+cFRED           WRITE(lunout,*)'FRED rugoro ...'
+cFRED           WRITE(lunout,*) rugoro (416), rugoro (417), rugoro(418)
+cFRED     .                    , rugoro(419)
+cFRED           WRITE(lunout,*) rugoro (421), rugoro (422), rugoro(423)
+cFRED           WRITE(lunout,*) rugoro (424), rugoro (425), rugoro(426)
+cFRED           WRITE(lunout,*) rugoro (427), rugoro (428), rugoro(429)
+cFRED           WRITE(lunout,*) rugoro (430), rugoro (431), rugoro(567)
+cFRED     .                    , rugoro(568)
+cFRED           WRITE(lunout,*) rugoro (574), rugoro (575), rugoro(576)
+cFRED           WRITE(lunout,*) rugoro (577), rugoro (698), rugoro(699)
+cFRED           WRITE(lunout,*) rugoro (711), rugoro (712), rugoro(713)
+
+          WRITE(lunout,*)'FRED rugoro =0.'
+          rugoro=0.
+
            CALL SUGWD(klon,klev,paprs,pplay)
            DO i=1,klon
@@ -2594 +2609 @@
 #ifdef histISCCP
       IF (ok_isccp) THEN
+c SD initialisation
+         fq_isccp(:,:,:,:)= 0.
+         totalcldarea(:,:)= 0.
+         meanptop(:,:)    = 0.
+         meantaucld(:,:)  = 0.
+         nbsunlit(:,:,:)  = 0.
+c
 cIM appel simulateur toutes les  NINT(freq_ISCCP/dtime) heures
        IF (MOD(itap,NINT(freq_ISCCP/dtime)).EQ.0) THEN
 #include "calcul_simulISCCP.h"
-       ENDIF !(MOD(itap,NINT(freq_ISCCP/dtime))
+       ENDIF !(MOD(itap,NINT(freq_ISCCP/dtime)))
       ENDIF !ok_isccp
 #endif