Changeset 53 for LMDZ.3.3/trunk/libf/phylmd/physiq.F

Timestamp:

Feb 18, 2000, 12:37:13 PM (24 years ago)

Author:

lmdzadmin

Message:

Correction sur le calendrier. LF

File:

• LMDZ.3.3/trunk/libf/phylmd/physiq.F (modified) (9 diffs)

LMDZ.3.3/trunk/libf/phylmd/physiq.F

@@ -56 +56 @@
 #include "clesphys.h"
 #include "control.h"
+#include "temps.h"
 c======================================================================
       LOGICAL check ! Verifier la conservation du modele en eau
@@ -282 +283 @@
       REAL yv1(klon)            ! vents dans la premiere couche V
       LOGICAL offline           ! Controle du stockage ds "physique"
-      PARAMETER (offline=.FALSE.)
+      PARAMETER (offline=.true.)
+      INTEGER physid
       REAL pfrac_impa(klon,klev)! Produits des coefs lessivage impaction
       save pfrac_impa
@@ -671 +673 @@
       IF (ok_journe) THEN
 c
-C         CALL ymds2ju(1900, 1, 1, 0.0, zjulian)
-         CALL ymds2ju(anneeref, 1, 1, 0.0, zjulian)
-         zjulian = zjulian + dayref
+         CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
+         zjulian = zjulian + day_ini
 c
          CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon)
@@ -824 +825 @@
          ndex3d = 0
 c
-Cess         i = NINT(zout/zsto)
-Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
-Cess        CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
-c
-Cess         i = NINT(zout/zsto)
-Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
-Cess        CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
-c
       ENDIF ! fin de test sur ok_journe
 c
       IF (ok_mensuel) THEN
 c
-c         CALL ymds2ju(1900, 1, 1, 0.0, zjulian)
-         CALL ymds2ju(anneeref, 1, 1, 0.0, zjulian)
-         zjulian = zjulian + dayref
+         CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
+         zjulian = zjulian + day_ini
 c
          CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon)
@@ -1172 +1164 @@
          ndex3d = 0
 c
-Cess         i = NINT(zout/zsto)
-Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
-Cess        CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
-C
-Cess         i = NINT(zout/zsto)
-Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
-Cess        CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
-c
       ENDIF ! fin de test sur ok_mensuel
 c
@@ -1185 +1169 @@
       IF (ok_instan) THEN
 c
-c         CALL ymds2ju(1900, 1, 1, 0.0, zjulian)
-         CALL ymds2ju(anneeref, 1, 1, 0.0, zjulian)
-         zjulian = zjulian + dayref
+         CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
+         zjulian = zjulian + day_ini
 c
          CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon)
@@ -1253 +1236 @@
          ndex2d = 0
          ndex3d = 0
-c
-Cess         i = NINT(zout/zsto)
-Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
-Cess        CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
-c
-Cess         i = NINT(zout/zsto)
-Cess         CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
-Cess        CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d)
 c
       ENDIF
@@ -1943 +1918 @@
 C la physique s'applique
 C
-       write(*,*) 'Phytrac= '
       call phytrac (rnpb,
      I                   debut,
@@ -1955 +1929 @@
      O                   tr_seri)
 
-       write(*,*) 'OFFLINE= ', offline
       IF (offline) THEN
-       write(*,*) 'OFFLINE= ', offline
-      call phystoke (
-     I                   nlon,nlev,pdtphys,
+
+         call phystokenc (
+     I                   nlon,nlev,pdtphys,rlon,rlat,
      I                   pmfu, pmfd, pen_u, pde_u, pen_d, pde_d,
      I                   ycoefh,yu1,yv1,ftsol,pctsrf,
-     I                   frac_impa, frac_nucl)
+     I                   frac_impa, frac_nucl,
+     I                   pphis,paire,dtime,itap,
+     O                   physid)
 
       ENDIF