Changeset 1035 for LMDZ4/trunk/libf

Timestamp:

Oct 28, 2008, 3:45:24 PM (16 years ago)

Author:

lmdzadmin

Message:

Bug physiq: RDAY a la place de 86400. pour calcul tendances dues au rayonnement (FC)
Ajout test coherence flags Tiedtke et simulateur ISCCP
FC/IM

File:

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

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -1254 +1254 @@
      .   iflag_coupl,iflag_clos,iflag_wake
       print*,'CYCLE_DIURNE', cycle_diurne
-
-c
+c
+      IF (iflag_con.EQ.2.AND.iflag_cldcon.NE.-1) THEN
+         abort_message = 'Tiedtke convection needs iflag_cldcon=-1 '
+         CALL abort_gcm (modname,abort_message,1)
+      ENDIF
+c
+      IF(ok_isccp.AND.iflag_con.LE.2) THEN
+         abort_message = 'iflag_con >= 3 for KE and ISCCP simulator'
+         CALL abort_gcm (modname,abort_message,1)
+      ENDIF
 c
 c Initialiser les compteurs:
@@ -2683 +2691 @@
       ENDIF
 
-c
 cIM calcul nuages par le simulateur ISCCP
 c
 #ifdef histISCCP
       IF (ok_isccp) THEN
+c
+c
 cIM appel simulateur toutes les  NINT(freq_ISCCP/dtime) heures
        IF (MOD(itap,NINT(freq_ISCCP/dtime)).EQ.0) THEN
@@ -2700 +2709 @@
 
       ENDIF
-
 c
 c 2. NUAGES STARTIFORMES
@@ -2941 +2949 @@
       DO i = 1, klon
          t_seri(i,k) = t_seri(i,k)
-     .               + (heat(i,k)-cool(i,k)) * dtime/86400.
+     .               + (heat(i,k)-cool(i,k)) * dtime/RDAY
       ENDDO
       ENDDO