Changeset 97 for trunk/libf/phyvenus/physiq.F

Timestamp:

Mar 22, 2011, 5:25:44 PM (14 years ago)

Author:

slebonnois

Message:

Serie de modifs SL pour homogeneisation des phytitan et phyvenus
Ca touche aussi aux liens phy/dyn (surtout a propos de clesphy0),
a verifier avec les autres, donc...

File:

• trunk/libf/phyvenus/physiq.F (modified) (10 diffs)

trunk/libf/phyvenus/physiq.F

@@ -5 +5 @@
       SUBROUTINE physiq (nlon,nlev,nqmax,
      .            debut,lafin,rjourvrai,gmtime,pdtphys,
-     .            paprs,pplay,ppk,pphi,pphis,presnivs,clesphy0,
+     .            paprs,pplay,ppk,pphi,pphis,presnivs,
      .            u,v,t,qx,
      .            omega,
@@ -12 +12 @@
       USE ioipsl
       USE histcom
+      USE infotrac
+      USE control_mod
       IMPLICIT none
 c======================================================================
@@ -71 +73 @@
 #include "dimsoil.h"
 #include "clesphys.h"
-#include "control.h"
 #include "temps.h"
 #include "comgeomphy.h"
-#include "advtrac.h"
 #include "iniprint.h"
 #include "raddim.h"
@@ -157 +157 @@
 #endif
 
-      INTEGER        longcles
-      PARAMETER    ( longcles = 20 )
-      REAL clesphy0( longcles      )
-c
 c Variables propres a la physique
 c
@@ -385 +381 @@
       REAL zx_tmp_2d(iim,jjmp1),zx_tmp_3d(iim,jjmp1,klev)
       REAL zx_lon(iim,jjmp1), zx_lat(iim,jjmp1)
-      REAL*4 z4_tmp_2d(iim,jjmp1),z4_tmp_3d(iim,jjmp1,klev)
-      REAL*4 z4_lon(iim,jjmp1), z4_lat(iim,jjmp1)
-      REAL*4 zpresnivs(klev)
-c     REAL z4_tmp_2d(iim,jjmp1),z4_tmp_3d(iim,jjmp1,klev)
-c     REAL z4_lon(iim,jjmp1), z4_lat(iim,jjmp1)
-c     REAL zpresnivs(klev)
 
       INTEGER nid_day, nid_mth, nid_ins
@@ -396 +386 @@
 c
       INTEGER nhori, nvert, idayref
-      REAL*4 zsto, zout, zsto1, zsto2, z4julian, zero
-c     REAL zsto, zout, zsto1, zsto2, z4julian, zero
+      REAL zsto, zout, zsto1, zsto2, zero
       parameter (zero=0.0e0)
       real zjulian
@@ -439 +428 @@
 c
 #include "YOMCST.h"
-#include "YOETHF.h"
-#include "FCTTRE.h"
 
 c======================================================================
@@ -489 +476 @@
      .       rlatd,rlond,ftsol,ftsoil,
      .       falbe, solsw, sollw,
-     .       dlw,radsol,clesphy0,
+     .       dlw,radsol,
      .      zmea,zstd,zsig,zgam,zthe,zpic,zval,
      .       tabcntr0,
@@ -577 +564 @@
      .                   iflag_ajs
 c
-         ecrit_mth = NINT(RDAY/dtime *ecritphy)  ! tous les ecritphy jours
+         ecrit_mth = NINT(RDAY/dtime*ecritphy)  ! tous les ecritphy jours
          IF (ok_mensuel) THEN
          WRITE(lunout,*)'La frequence de sortie mensuelle est de ', 
      .                   ecrit_mth
          ENDIF
+
          ecrit_day = NINT(RDAY/dtime *1.0)  ! tous les jours
          IF (ok_journe) THEN
@@ -587 +575 @@
      .                   ecrit_day
          ENDIF
-         ecrit_ins = NINT(RDAY/dtime/48.)  ! a chaque pas de temps ==> PB. dans time_counter pour 1mois
-         ecrit_ins = NINT(RDAY/dtime/12.)  ! toutes les deux heures
+
+         ecrit_ins = NINT(RDAY/dtime/*ecritphy)  ! Fraction de jour reglable
          IF (ok_instan) THEN
          WRITE(lunout,*)'La frequence de sortie instant. est de ',