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

Timestamp:

Nov 26, 2001, 2:47:57 PM (23 years ago)

Author:

lmdzadmin

Message:

Tag version 0 qui marche en couple/force
LF

File:

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

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

@@ -371 +371 @@
       REAL albsollw(klon)
       SAVE albsollw                 ! albedo du sol total
+      REAL albsol1(klon)
+      SAVE albsol1                 ! albedo du sol total
+      REAL albsollw1(klon)
+      SAVE albsollw1                 ! albedo du sol total
 
       REAL wo(klon,klev)
@@ -576 +580 @@
 c
       INTEGER nhori, nvert
-      REAL zsto, zout, zjulian
+      REAL zsto, zout
+      real zjulian
+      save zjulian
 
       character*20 modname
@@ -728 +734 @@
 ccc         ecrit_ins = NINT(86400./dtime *0.25)  ! 4 fois par jour
          ecrit_ins = NINT(86400./dtime/48.)  ! a chaque pas de temps
+         ecrit_ins = NINT(86400./dtime/12.)  ! toutes les deux heures
          IF (ok_instan) THEN
          PRINT*, 'La frequence de sortie instant. est de ', ecrit_ins
@@ -750 +757 @@
       endif        
 c
+c
+c Gestion calendrier
+
+         CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
+         zjulian = zjulian + day_ini
+
 c
       IF (ok_journe) THEN
@@ -1807 +1820 @@
       fder = dlw
 
-      CALL clmain(dtime,itap,date0,pctsrf,
+      CALL clmain(dtime,itap,zjulian,pctsrf,
      e            t_seri,q_seri,u_seri,v_seri,
      e            julien, rmu0, 
@@ -1971 +1984 @@
 c$$$     .        ,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,
 c$$$     .        pori_con,plcl_con,dtma_con,dtlcl_con)
-          if (1.eq.1) then ! vectorise
+          if (0.eq.1) then ! vectorise
           CALL conemav (dtime,paprs,pplay,t_seri,q_seri,
      .        u_seri,v_seri,tr_seri,nbtr,
@@ -2197 +2210 @@
      .               + falblw(i,is_sic) * pctsrf(i,is_sic)
       ENDDO
+!      if (debut) then
+!        albsol1 = albsol
+!        albsollw1 = albsollw
+!      endif
+!      albsol = albsol1
+!      albsollw = albsollw1
       CALL radlwsw ! nouveau rayonnement (compatible Arpege-IFS)
      e            (dist, rmu0, fract, co2_ppm, solaire,
@@ -2331 +2350 @@
       IF (iflag_con.EQ.4) THEN
 c           on ajoute les tendances calculees par KE43
+c$$$ OM on onhibe la convection sur les traceurs
         DO iq=1, nqmax-2 ! Sandrine a -3 ???
-        DO k = 1, nlev
-        DO i = 1, klon
-          tr_seri(i,k,iq) = tr_seri(i,k,iq) + d_tr(i,k,iq)
-        ENDDO
-        ENDDO
+c$$$ OM on inhibe la convection sur les traceur
+c$$$        DO k = 1, nlev
+c$$$        DO i = 1, klon
+c$$$          tr_seri(i,k,iq) = tr_seri(i,k,iq) + d_tr(i,k,iq)
+c$$$        ENDDO
+c$$$        ENDDO
         WRITE(iqn,'(i2.2)') iq
         CALL minmaxqfi(tr_seri(1,1,iq),0.,1.e33,'couche lim iq='//iqn)