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

Timestamp:

Apr 13, 2001, 4:58:01 PM (23 years ago)

Author:

lmdzadmin

Message:

Rajout date0 dans physiq pour interface avec ORCHIDEE
Fermeture des fichiers restart par restclo dans abort_gcm
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

@@ -94 +94 @@
 ccc      PARAMETER (soil_model=.FALSE.)
       logical ok_veget
-      parameter (ok_veget = .false.)
-c      parameter (ok_veget = .true.)
+      parameter (ok_veget = .true.)
+c      parameter (ok_veget = .false.)
 c======================================================================
 c Dans les versions precedentes, l'eau liquide nuageuse utilisee dans
@@ -118 +118 @@
 c
       LOGICAL ok_instan ! sortir le fichier instantane
-      PARAMETER (ok_instan=.true.)
+      PARAMETER (ok_instan=.false.)
 c
       LOGICAL ok_region ! sortir le fichier regional
@@ -499 +499 @@
       character*80 abort_message
       logical ok_sync
-
+      real date0
 
 c
@@ -1525 +1525 @@
 c Verifier les temperatures
 c
-
       CALL hgardfou(t_seri,ftsol,'debutphy')
 c
@@ -1593 +1592 @@
 
       fder = 0.
-
-      CALL clmain(dtime,itap,pctsrf,
+      date0 = day_ini
+
+      CALL clmain(dtime,itap,date0,pctsrf,
      e            t_seri,q_seri,u_seri,v_seri,
      e            julien, rmu0, 
@@ -1607 +1607 @@
      s            dsens, devap,
      s            ycoefh,yu1,yv1) 
+
 c
 C§§§ PB
@@ -1633 +1634 @@
 c         evap(i) = - fluxq(i,1) ! flux d'evaporation au sol
          evap(i) = - zxfluxq(i,1) ! flux d'evaporation au sol
+C LF test signe flux
+         sens(i) = zxfluxt(i,1)
+         evap(i) = zxfluxq(i,1)
          fder(i) = dsens(i) + devap(i)
       ENDDO
@@ -1738 +1742 @@
       CALL abort
       ENDIF
+
       CALL homogene(paprs, q_seri, d_q_con, u_seri,v_seri,
      .              d_u_con, d_v_con)
       DO k = 1, klev
-      DO i = 1, klon
+        DO i = 1, klon
          t_seri(i,k) = t_seri(i,k) + d_t_con(i,k)
          q_seri(i,k) = q_seri(i,k) + d_q_con(i,k)
          u_seri(i,k) = u_seri(i,k) + d_u_con(i,k)
          v_seri(i,k) = v_seri(i,k) + d_v_con(i,k)
-      ENDDO
+        ENDDO
       ENDDO
       IF (check) THEN