Changeset 97 for trunk/libf/phyvenus


Ignore:
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...

Location:
trunk/libf/phyvenus
Files:
26 edited

Legend:

Unmodified
Added
Removed

  • trunk/libf/phyvenus/YOMCST.h

    r3 r97  
    11!
    2 ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/YOMCST.h,v 1.1.1.1 2004/05/19 12:53:08 lmdzadmin Exp $
     2!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
     3!                 veillez  n'utiliser que des ! pour les commentaires
     4!                 et  bien positionner les & des lignes de continuation
     5!                 (les placer en colonne 6 et en colonne 73)
    36!
    47! A1.0 Fundamental constants
    … …  
    710      REAL RDAY,REA,REPSM,RSIYEA,RSIDAY,ROMEGA
    811! A1.1.bis Constantes concernant l'orbite de la Terre:
    9       REAL R_ecc, R_peri, R_incl
     12      REAL :: R_ecc, R_peri, R_incl
    1013! A1.2 Geoide
    1114      REAL RA,RG,R1SA
    1215! A1.3 Radiation
    13 !     REAL RSIGMA,RI0
     16!     REAL :: RSIGMA,RI0
    1417      REAL RSIGMA
    1518! A1.4 Thermodynamic gas phase
    … …  
    2730      REAL RALPD,RBETD,RGAMD
    2831!
    29 !    S      ,RSIGMA,RI0
    30       COMMON/YOMCST/RPI   ,RCLUM ,RHPLA ,RKBOL ,RNAVO
    31      S      ,RDAY  ,REA   ,REPSM ,RSIYEA,RSIDAY,ROMEGA
    32      s      ,R_ecc, R_peri, R_incl
    33      S      ,RA    ,RG    ,R1SA
    34      S      ,RSIGMA
    35      S      ,R     ,RMD   ,RMV   ,RD    ,RV    ,RCPD
    36      S      ,RCPV  ,RCVD  ,RCVV  ,RKAPPA,RETV
    37      S      ,RCW   ,RCS
    38      S      ,RLVTT ,RLSTT ,RLMLT ,RTT   ,RATM
    39      S      ,RESTT ,RALPW ,RBETW ,RGAMW ,RALPS ,RBETS ,RGAMS
    40      S      ,RALPD ,RBETD ,RGAMD
    41 !    ------------------------------------------------------------------
     32      COMMON/YOMCST/RPI ,RCLUM, RHPLA, RKBOL, RNAVO ,RDAY  ,REA         &
     33     & ,REPSM ,RSIYEA,RSIDAY,ROMEGA , R_ecc, R_peri, R_incl,            &
     34     & ,RA    ,RG ,R1SA                                                 &
     35     & ,RSIGMA,R ,RMD   ,RMV   ,RD    ,RV    ,RCPD ,RCPV,RCVD           &
     36     & ,RCVV  ,RKAPPA,RETV ,RCW   ,RCS ,RLVTT ,RLSTT ,RLMLT ,RTT ,RATM  &
     37     & ,RESTT ,RALPW ,RBETW ,RGAMW ,RALPS ,RBETS ,RGAMS ,RALPD ,RBETD   &
     38     & ,RGAMD

  • trunk/libf/phyvenus/clcdrag.F90

    r3 r97  
    3030! ================================================================= c
    3131!
    32 #include "YOMCST.inc"
    33 #include "YOETHF.inc"
    34 #include "clesphys.inc"
     32#include "YOMCST.h"
     33#include "clesphys.h"
    3534!
    3635! Quelques constantes et options:

  • trunk/libf/phyvenus/clesphys.h

    r3 r97  
    11!
    2 ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/clesphys.h,v 1.1.1.1 2004/05/19 12:53:08 lmdzadmin Exp $
     2!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
     3!                 veillez  n'utiliser que des ! pour les commentaires
     4!                 et  bien positionner les & des lignes de continuation
     5!                 (les placer en colonne 6 et en colonne 73)
    36!
    4 c..include cles_phys.h
    5 c
    6        LOGICAL cycle_diurne,soil_model,new_oliq
     7!
     8!..include cles_phys.h
     9!
     10       LOGICAL cycle_diurne,soil_model
    711       LOGICAL ok_orodr,ok_orolf,ok_gw_nonoro
    8        INTEGER nbapp_rad, iflag_con, iflag_ajs
     12       INTEGER nbapp_rad, nbapp_chim, iflag_con, iflag_ajs
     13       REAL    ecritphy
    914       REAL    solaire
    10 
    11 c Parametres pour PBL:
    12        real z0,lmixmin
    13        REAL ksta
     15       REAL    z0, lmixmin
     16       REAL    ksta, inertie
    1417       LOGICAL ok_kzmin
    15 
    16 c surface
    17        real inertie
    18 
    19 cIM lev_histhf  : niveau sorties 6h
    20 cIM lev_histday : niveau sorties journalieres
    21 cIM lev_histmth : niveau sorties mensuelles
    2218       INTEGER lev_histhf, lev_histday, lev_histmth
    2319
    24        COMMON/clesphys/cycle_diurne, soil_model, new_oliq,
    25      ,     ok_orodr, ok_orolf, ok_gw_nonoro,
    26      ,     nbapp_rad, iflag_con
    27      ,     , iflag_ajs, solaire, z0, lmixmin, ksta
    28      ,     , ok_kzmin, lev_histhf, lev_histday, lev_histmth
    29      ,     , inertie
     20       COMMON/clesphys/cycle_diurne, soil_model,                        &
     21     &     ok_orodr, ok_orolf, ok_gw_nonoro, nbapp_rad, nbapp_chim      &
     22     &     , ecritphy                                                   &
     23     &     , iflag_con, iflag_ajs, solaire, z0, lmixmin, ksta           &
     24     &     , ok_kzmin, lev_histhf, lev_histday, lev_histmth             &
     25     &     , inertie
     26

  • trunk/libf/phyvenus/clmain.F

    r3 r97  
    134134c
    135135#include "YOMCST.h"
    136 #include "YOETHF.h"
    137 #include "FCTTRE.h"
    138136      REAL u1lay(klon), v1lay(klon)
    139137      REAL delp(klon,klev)
    … …  
    499497#include "dimphy.h"
    500498#include "YOMCST.h"
    501 #include "YOETHF.h"
    502 #include "FCTTRE.h"
    503499#include "dimsoil.h"
    504500#include "iniprint.h"
    … …  
    994990      PARAMETER (zxli=.FALSE.)
    995991c
    996 #include "YOETHF.h"
    997 #include "FCTTRE.h"
    998992      fsta(x) = 1.0 / (1.0+10.0*x*(1+8.0*x))
    999993      fins(x) = SQRT(1.0-18.0*x)

  • trunk/libf/phyvenus/clmain.classic

    r3 r97  
    134134c
    135135#include "YOMCST.h"
    136 #include "YOETHF.h"
    137 #include "FCTTRE.h"
    138136      REAL u1lay(klon), v1lay(klon)
    139137      REAL delp(klon,klev)
    … …  
    499497#include "dimphy.h"
    500498#include "YOMCST.h"
    501 #include "YOETHF.h"
    502 #include "FCTTRE.h"
    503499#include "dimsoil.h"
    504500#include "iniprint.h"
    … …  
    994990      PARAMETER (zxli=.FALSE.)
    995991c
    996 #include "YOETHF.h"
    997 #include "FCTTRE.h"
    998992      fsta(x) = 1.0 / (1.0+10.0*x*(1+8.0*x))
    999993      fins(x) = SQRT(1.0-18.0*x)

  • trunk/libf/phyvenus/clmain.simple

    r3 r97  
    128128c
    129129#include "YOMCST.h"
    130 #include "YOETHF.h"
    131 #include "FCTTRE.h"
    132130      REAL u1lay(klon), v1lay(klon)
    133131      REAL delp(klon,klev)
    … …  
    370368#include "dimphy.h"
    371369#include "YOMCST.h"
    372 #include "YOETHF.h"
    373 #include "FCTTRE.h"
    374370#include "dimsoil.h"
    375371#include "iniprint.h"
    … …  
    865861      PARAMETER (zxli=.FALSE.)
    866862c
    867 #include "YOETHF.h"
    868 #include "FCTTRE.h"
    869863      fsta(x) = 1.0 / (1.0+10.0*x*(1+8.0*x))
    870864      fins(x) = SQRT(1.0-18.0*x)

  • trunk/libf/phyvenus/conf_phys.F90

    r3 r97  
    1111   implicit none
    1212
    13 #include "YOMCST.inc"
    14 #include "clesphys.inc"
     13#include "YOMCST.h"
     14#include "clesphys.h"
    1515#include "compbl.h"
    16 !
    17 ! Configuration de la "physique" de LMDZ a l'aide de la fonction
    18 ! GETIN de IOIPSL
    19 !
    20 ! LF 05/2001
    21 !
    22 !--- Ca lit le physiq.def ---
    23 
    24 !
     16
    2517! ok_journe:  sorties journalieres
    2618! ok_mensuel: sorties mensuelles
    … …  
    3426! Local
    3527  integer              :: numout = 6
    36   real*4               :: zzz
    37 ! real                 :: zzz
     28
     29!
     30! Configuration de la "physique" de LMDZ a l'aide de la fonction
     31! GETIN de IOIPSL
     32!
     33! LF 05/2001
     34!
     35!--- Ca lit le physiq.def ---
     36
     37!******************* parametres anciennement lus dans gcm.def
     38
     39!Config  Key  = cycle_diurne
     40!Config  Desc = Cycle ddiurne
     41!Config  Def  = y
     42!Config  Help = Cette option permet d'eteidre le cycle diurne.
     43!Config         Peut etre util pour accelerer le code !
     44       cycle_diurne = .TRUE.
     45       CALL getin('cycle_diurne',cycle_diurne)
     46
     47!Config  Key  = soil_model
     48!Config  Desc = Modele de sol
     49!Config  Def  = y
     50!Config  Help = Choix du modele de sol (Thermique ?)
     51!Config         Option qui pourait un string afin de pouvoir
     52!Config         plus de choix ! Ou meme une liste d'options !
     53       soil_model = .TRUE.
     54       CALL getin('soil_model',soil_model)
     55
     56!Config  Key  = ok_orodr
     57!Config  Desc = Oro drag
     58!Config  Def  = y
     59!Config  Help = GW drag orographie
     60!Config         
     61       ok_orodr = .TRUE.
     62       CALL getin('ok_orodr',ok_orodr)
     63
     64!Config  Key  =  ok_orolf
     65!Config  Desc = Oro lift
     66!Config  Def  = n
     67!Config  Help = GW lift orographie (pas utilise)
     68       ok_orolf = .TRUE.
     69       CALL getin('ok_orolf', ok_orolf)
     70
     71!Config  Key  = ok_gw_nonoro
     72!Config  Desc = Gravity waves parameterization
     73!Config  Def  = n
     74!Config  Help = GW drag non-orographique
     75       ok_gw_nonoro = .FALSE.
     76       CALL getin('ok_gw_nonoro',ok_gw_nonoro)
     77
     78!Config  Key  = nbapp_rad
     79!Config  Desc = Frequence d'appel au rayonnement
     80!Config  Def  = 12
     81!Config  Help = Nombre  d'appels des routines de rayonnements
     82!Config         par jour.
     83       nbapp_rad = 12
     84       CALL getin('nbapp_rad',nbapp_rad)
     85
     86!Config  Key  = nbapp_chim
     87!Config  Desc = Frequence d'appel a la chimie
     88!Config  Def  = 1
     89!Config  Help = Nombre  d'appels des routines de chimie
     90!Config         par jour.
     91       nbapp_chim = 1
     92       CALL getin('nbapp_chim',nbapp_chim)
     93
     94!Config  Key  = iflag_con
     95!Config  Desc = Flag de convection
     96!Config  Def  = 0
     97!Config  Help = Flag  pour la convection les options suivantes existent :
     98!Config         0 : ajsec simple (VENUS, TITAN)
     99!Config         1 pour LMD,
     100!Config         2 pour Tiedtke,
     101!Config         3 pour CCM(NCAR) 
     102       iflag_con = 0
     103       CALL getin('iflag_con',iflag_con)
     104
     105!******************* fin parametres anciennement lus dans gcm.def
    38106
    39107!Config Key  = OK_journe
    … …  
    63131  ok_instan = .false.
    64132  call getin('OK_instan', ok_instan)
     133!
     134!Config  Key  = ecritphy
     135!Config  Desc = Frequence d'ecriture dans histmth et histins
     136!Config  Def  = 1
     137!Config  Help = frequence de l'ecriture du fichier histmth et histins
     138!Config         en jours.
     139!
     140       ecritphy = 1.
     141       CALL getin('ecritphy', ecritphy)
    65142!
    66143!
    … …  
    83160! VENUS
    84161! R_ecc = 0.006787
    85   zzz   = 0.0
    86   call getin('R_ecc', zzz)
    87   R_ecc = zzz
     162  R_ecc   = 0.0
     163  call getin('R_ecc', R_ecc)
    88164!!
    89165!Config Key  = R_peri
    … …  
    93169!               
    94170! VENUS
    95   zzz = 0.
    96   call getin('R_peri', zzz)
    97   R_peri = zzz
     171  R_peri = 0.
     172  call getin('R_peri', R_peri)
    98173!!
    99174!Config Key  = R_incl
    … …  
    103178!               
    104179! VENUS
    105   zzz = 0.0
    106   call getin('R_incl', zzz)
    107   R_incl = zzz
     180  R_incl = 0.0
     181  call getin('R_incl', R_incl)
    108182!
    109183!Config Key  = solaire
    … …  
    113187!Config Help =
    114188!
    115   zzz = 2620.
    116     call getin('solaire', zzz)
    117   solaire = zzz
     189  solaire = 2620.
     190    call getin('solaire', solaire)
    118191!
    119192!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    140213!Config Help =
    141214!
    142   zzz = 1.0e-7
    143   call getin('ksta',zzz)
    144   ksta = zzz
     215  ksta = 1.0e-7
     216  call getin('ksta',ksta)
    145217
    146218!
    … …  
    150222!Config Help =
    151223!
    152   zzz = 1.0e-2
    153   call getin('z0',zzz)
    154   z0 = zzz
     224  z0 = 1.0e-2
     225  call getin('z0',z0)
    155226
    156227!
    … …  
    160231!Config Help =
    161232!
    162   zzz = 35.
    163   call getin('lmixmin',zzz)
    164   lmixmin = zzz
    165 
    166 !
     233  lmixmin = 35.
     234  call getin('lmixmin',lmixmin)
     235
    167236!
    168237!Config Key  = ok_kzmin
    … …  
    186255!Config Key  = inertie
    187256!Config Desc =
    188 !Config Def  = 30000.
    189 !Config Help =
    190 !
    191   zzz = 30000.
    192   call getin('inertie',zzz)
    193   inertie = zzz
     257!Config Def  = 2000.
     258!Config Help =
     259!
     260  inertie = 2000.
     261  call getin('inertie',inertie)
    194262!
    195263!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    239307  write(numout,*)' ##############################################'
    240308  write(numout,*)' Configuration des parametres de la physique: '
     309  write(numout,*)' cycle_diurne = ', cycle_diurne
     310  write(numout,*)' soil_model = ', soil_model
     311  write(numout,*)' ok_orodr = ', ok_orodr
     312  write(numout,*)' ok_orolf = ', ok_orolf
     313  write(numout,*)' ok_gw_nonoro = ', ok_gw_nonoro
     314  write(numout,*)' nbapp_rad = ', nbapp_rad
     315  write(numout,*)' nbapp_chim = ', nbapp_chim
     316  write(numout,*)' iflag_con = ', iflag_con
    241317  write(numout,*)' Sortie journaliere = ', ok_journe
    242318  write(numout,*)' Sortie mensuelle = ', ok_mensuel
    243319  write(numout,*)' Sortie instantanee = ', ok_instan
     320  write(numout,*)' frequence sorties = ', ecritphy 
    244321  write(numout,*)' Sortie bilan d''energie, if_ebil =', if_ebil
    245322  write(numout,*)' Excentricite = ',R_ecc

  • trunk/libf/phyvenus/diagphy.F

    r3 r97  
    77     $    , d_etp_tot, d_qt_tot, d_ec_tot
    88     $    , fs_bound, fq_bound)
     9
     10! ATTENTION !! PAS DU TOUT A JOUR POUR VENUS OU TITAN...
     11
    912C======================================================================
    1013C
    … …  
    5053#include "dimphy.h"
    5154#include "YOMCST.h"
    52 #include "YOETHF.h"
    5355C
    5456C     Input variables
    … …  
    211213#include "dimphy.h"
    212214#include "YOMCST.h"
    213 #include "YOETHF.h"
    214215C
    215216c     Input variables

  • trunk/libf/phyvenus/ini_histday.h

    r3 r97  
    99c
    1010         idayref = day_ref
    11          z4julian = zjulian
    12          CALL ymds2ju(annee_ref, 1, idayref, zero, z4julian)
     11         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
    1312c
    1413         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
    … …  
    1817         ENDDO
    1918         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)
    20          z4_lon = zx_lon
    21          z4_lat = zx_lat
    22          zpresnivs=presnivs
    23          CALL histbeg("histday", iim,z4_lon(:,1), jjmp1,z4_lat(1,:),
    24      .                 1,iim,1,jjmp1, itau_phy, z4julian, zsto,
     19
     20c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
     21c    d'ou le -1.* devant zx_lon et zx_lat       
     22         CALL histbeg("histday.nc", iim,-1.*zx_lon(:,1), jjmp1,
     23     .                 -1.*zx_lat(1,:),
     24     .                 1,iim,1,jjmp1, itau_phy, zjulian, zsto,
    2525     .                 nhori, nid_day)
    2626         write(*,*)'Journee ', itau_phy, zjulian
    2727         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
    28      .                 klev, zpresnivs, nvert)
     28     .                 klev, presnivs, nvert)
    2929
    3030c-------------------------------------------------------

  • trunk/libf/phyvenus/ini_histins.h

    r3 r97  
    99c
    1010         idayref = day_ref
    11          z4julian = zjulian
    12          CALL ymds2ju(annee_ref, 1, idayref, zero, z4julian)
     11         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
    1312c
    1413         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
    … …  
    1817         ENDDO
    1918         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)
    20          z4_lon = zx_lon
    21          z4_lat = zx_lat
    22          zpresnivs=presnivs
    23          CALL histbeg("histins", iim,z4_lon(:,1), jjmp1,z4_lat(1,:),
    24      .                 1,iim,1,jjmp1, itau_phy, z4julian, zsto1,
     19         
     20c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
     21c    d'ou le -1.* devant zx_lon et zx_lat       
     22         CALL histbeg("histins.nc", iim,-1.*zx_lon(:,1), jjmp1,
     23     .                 -1.*zx_lat(1,:),
     24     .                 1,iim,1,jjmp1, itau_phy, zjulian, zsto1,
    2525     .                 nhori, nid_ins)
    2626         write(*,*)'Inst ', itau_phy, zjulian
    2727         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb",
    28      .                 klev, zpresnivs, nvert)
     28     .                 klev, presnivs, nvert)
    2929
    3030c-------------------------------------------------------

  • trunk/libf/phyvenus/ini_histmth.h

    r3 r97  
    1414c
    1515         idayref = day_ref
    16          z4julian = zjulian
    17          CALL ymds2ju(annee_ref, 1, idayref, zero, z4julian)
     16         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
    1817c
    1918         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
    … …  
    2524         
    2625c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
    27          z4_lon = -1.*zx_lon
    28          z4_lat = -1.*zx_lat
    29          
    30          zpresnivs=presnivs
    31          CALL histbeg("histmth.nc", iim,z4_lon(:,1), jjmp1,z4_lat(1,:),
    32      .                 1,iim,1,jjmp1, itau_phy, z4julian, zsto,
     26c    d'ou le -1.* devant zx_lon et zx_lat       
     27         CALL histbeg("histmth.nc", iim,-1.*zx_lon(:,1), jjmp1,
     28     .                 -1.*zx_lat(1,:),
     29     .                 1,iim,1,jjmp1, itau_phy, zjulian, zsto,
    3330     .                 nhori, nid_mth)
    3431         write(*,*)'Mensuel ', itau_phy, zjulian
    3532         CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
    36      .                 klev, zpresnivs, nvert)
     33     .                 klev, presnivs, nvert)
    3734
    3835c-------------------------------------------------------

  • trunk/libf/phyvenus/interface_surf.F90

    r3 r97  
    2525  END INTERFACE
    2626
    27 #include "YOMCST.inc"
     27#include "YOMCST.h"
    2828
    2929  CONTAINS
    … …  
    210210!
    211211
    212 #include "YOETHF.inc"
    213 #include "FCTTRE.inc"
    214 
    215212! Parametres d'entree
    216213  integer, intent(IN) :: knon, klon

  • trunk/libf/phyvenus/phyetat0.F

    r3 r97  
    77     .            rlat,rlon, tsol,tsoil,
    88     .           albe, solsw, sollw,
    9      .           fder,radsol,clesphy0,
     9     .           fder,radsol,
    1010     .    zmea, zstd, zsig, zgam, zthe, zpic, zval, 
    1111     .           tabcntr0,
    … …  
    4444      LOGICAL ancien_ok
    4545
    46       INTEGER        longcles
    47       PARAMETER    ( longcles = 20 )
    48       REAL clesphy0( longcles )
    49 c
    5046      REAL xmin, xmax
    5147c
    … …  
    8884         ENDDO
    8985c
    90          cycle_diurne   = .FALSE.
    91          soil_model     = .FALSE.
    92          new_oliq       = .FALSE.
    93          ok_orodr       = .FALSE.
    94          ok_orolf       = .FALSE.
    95          ok_gw_nonoro   = .FALSE.
    96 
    97 
    98          IF( clesphy0(1).NE.tab_cntrl( 5 ) )  THEN
    99              tab_cntrl( 5 ) = clesphy0(1)
    100          ENDIF
    101 
    102          IF( clesphy0(2).NE.tab_cntrl( 6 ) )  THEN
    103              tab_cntrl( 6 ) = clesphy0(2)
    104          ENDIF
    105 
    106          IF( clesphy0(3).NE.tab_cntrl( 7 ) )  THEN
    107              tab_cntrl( 7 ) = clesphy0(3)
    108          ENDIF
    109 
    110          IF( clesphy0(4).NE.tab_cntrl( 8 ) )  THEN
    111              tab_cntrl( 8 ) = clesphy0(4)
    112          ENDIF
    113 
    114          IF( clesphy0(5).NE.tab_cntrl( 9 ) )  THEN
    115              tab_cntrl( 9 ) = clesphy0( 5 )
    116          ENDIF
    117 
    118          IF( clesphy0(6).NE.tab_cntrl( 10 ) )  THEN
    119              tab_cntrl( 10 ) = clesphy0( 6 )
    120          ENDIF
    121 
    122          IF( clesphy0(7).NE.tab_cntrl( 11 ) )  THEN
    123              tab_cntrl( 11 ) = clesphy0( 7 )
    124          ENDIF
    125 
    126          IF( clesphy0(8).NE.tab_cntrl( 12 ) )  THEN
    127              tab_cntrl( 12 ) = clesphy0( 8 )
    128          ENDIF
    129 
    13086
    13187         dtime        = tab_cntrl(1)
    13288         radpas       = tab_cntrl(2)
    133          iflag_con    = tab_cntrl(5)
    134          nbapp_rad    = tab_cntrl(6)
    135 
    136 
    137          cycle_diurne    = .FALSE.
    138          soil_model      = .FALSE.
    139          new_oliq        = .FALSE.
    140          ok_orodr        = .FALSE.
    141          ok_orolf        = .FALSE.
    142          ok_gw_nonoro    = .FALSE.
    143 
    144          IF( tab_cntrl( 7) .EQ. 1. )    cycle_diurne  = .TRUE.
    145          IF( tab_cntrl( 8) .EQ. 1. )       soil_model = .TRUE.
    146          IF( tab_cntrl( 9) .EQ. 1. )         new_oliq = .TRUE.
    147          IF( tab_cntrl(10) .EQ. 1. )         ok_orodr = .TRUE.
    148          IF( tab_cntrl(11) .EQ. 1. )         ok_orolf = .TRUE.
    149          IF( tab_cntrl(12) .EQ. 1. )     ok_gw_nonoro = .TRUE.
    15089
    15190      ENDIF

  • trunk/libf/phyvenus/phyredem.F

    r3 r97  
    1010     .    zmea, zstd, zsig, zgam, zthe, zpic, zval,
    1111     .           t_ancien)
     12
     13      USE control_mod
    1214      IMPLICIT none
    1315c======================================================================
    … …  
    2022#include "dimsoil.h"
    2123#include "clesphys.h"
    22 #include "control.h"
    2324#include "temps.h"
    2425c======================================================================

  • trunk/libf/phyvenus/physiq.F

    r3 r97  
    55      SUBROUTINE physiq (nlon,nlev,nqmax,
    66     .            debut,lafin,rjourvrai,gmtime,pdtphys,
    7      .            paprs,pplay,ppk,pphi,pphis,presnivs,clesphy0,
     7     .            paprs,pplay,ppk,pphi,pphis,presnivs,
    88     .            u,v,t,qx,
    99     .            omega,
    … …  
    1212      USE ioipsl
    1313      USE histcom
     14      USE infotrac
     15      USE control_mod
    1416      IMPLICIT none
    1517c======================================================================
    … …  
    7173#include "dimsoil.h"
    7274#include "clesphys.h"
    73 #include "control.h"
    7475#include "temps.h"
    7576#include "comgeomphy.h"
    76 #include "advtrac.h"
    7777#include "iniprint.h"
    7878#include "raddim.h"
    … …  
    157157#endif
    158158
    159       INTEGER        longcles
    160       PARAMETER    ( longcles = 20 )
    161       REAL clesphy0( longcles      )
    162 c
    163159c Variables propres a la physique
    164160c
    … …  
    385381      REAL zx_tmp_2d(iim,jjmp1),zx_tmp_3d(iim,jjmp1,klev)
    386382      REAL zx_lon(iim,jjmp1), zx_lat(iim,jjmp1)
    387       REAL*4 z4_tmp_2d(iim,jjmp1),z4_tmp_3d(iim,jjmp1,klev)
    388       REAL*4 z4_lon(iim,jjmp1), z4_lat(iim,jjmp1)
    389       REAL*4 zpresnivs(klev)
    390 c     REAL z4_tmp_2d(iim,jjmp1),z4_tmp_3d(iim,jjmp1,klev)
    391 c     REAL z4_lon(iim,jjmp1), z4_lat(iim,jjmp1)
    392 c     REAL zpresnivs(klev)
    393383
    394384      INTEGER nid_day, nid_mth, nid_ins
    … …  
    396386c
    397387      INTEGER nhori, nvert, idayref
    398       REAL*4 zsto, zout, zsto1, zsto2, z4julian, zero
    399 c     REAL zsto, zout, zsto1, zsto2, z4julian, zero
     388      REAL zsto, zout, zsto1, zsto2, zero
    400389      parameter (zero=0.0e0)
    401390      real zjulian
    … …  
    439428c
    440429#include "YOMCST.h"
    441 #include "YOETHF.h"
    442 #include "FCTTRE.h"
    443430
    444431c======================================================================
    … …  
    489476     .       rlatd,rlond,ftsol,ftsoil,
    490477     .       falbe, solsw, sollw,
    491      .       dlw,radsol,clesphy0,
     478     .       dlw,radsol,
    492479     .      zmea,zstd,zsig,zgam,zthe,zpic,zval,
    493480     .       tabcntr0,
    … …  
    577564     .                   iflag_ajs
    578565c
    579          ecrit_mth = NINT(RDAY/dtime *ecritphy)  ! tous les ecritphy jours
     566         ecrit_mth = NINT(RDAY/dtime*ecritphy)  ! tous les ecritphy jours
    580567         IF (ok_mensuel) THEN
    581568         WRITE(lunout,*)'La frequence de sortie mensuelle est de ',
    582569     .                   ecrit_mth
    583570         ENDIF
     571
    584572         ecrit_day = NINT(RDAY/dtime *1.0)  ! tous les jours
    585573         IF (ok_journe) THEN
    … …  
    587575     .                   ecrit_day
    588576         ENDIF
    589          ecrit_ins = NINT(RDAY/dtime/48.)  ! a chaque pas de temps ==> PB. dans time_counter pour 1mois
    590          ecrit_ins = NINT(RDAY/dtime/12.)  ! toutes les deux heures
     577
     578         ecrit_ins = NINT(RDAY/dtime/*ecritphy)  ! Fraction de jour reglable
    591579         IF (ok_instan) THEN
    592580         WRITE(lunout,*)'La frequence de sortie instant. est de ',

  • trunk/libf/phyvenus/phytrac.F

    r3 r97  
    2626
    2727      USE ioipsl
     28      USE infotrac
     29      USE control_mod
    2830
    2931      IMPLICIT none
    … …  
    4244#include "temps.h"
    4345#include "paramet.h"
    44 #include "control.h"
    4546#include "comgeomphy.h"
    46 #include "advtrac.h"
    4747c======================================================================
    4848
    … …  
    124124c   k,i,it,tr_seri(klon,klev,nqmax),pplay(nlon,nlev),
    125125      integer nqCO_OCS
    126       real zprof(klev,nqmx)
     126      real zprof(klev,nqtot)
    127127c      real pzero,gamma
    128128c      parameter (pzero=85000.)
    129129c      parameter (gamma=5000.)
    130130      REAL alpha
    131       real deltatr(klon,klev,nqmx) ! ecart au profil de ref zprof
    132       real tau(klev,nqmx)          ! temps de relaxation vers le profil (s)
     131      real deltatr(klon,klev,nqtot) ! ecart au profil de ref zprof
     132      real tau(klev,nqtot)          ! temps de relaxation vers le profil (s)
    133133      save zprof,tau
    134134c======================================================================

  • trunk/libf/phyvenus/radlwsw.1mat

    r3 r97  
    3838#include "dimphy.h"
    3939#include "raddim.h"
    40 #include "YOETHF.h"
    4140c
    4241      real rmu0(klon), fract(klon), dist

  • trunk/libf/phyvenus/radlwsw.F

    r3 r97  
    4242#include "dimphy.h"
    4343#include "raddim.h"
    44 #include "YOETHF.h"
    4544#include "YOMCST.h"
    4645#include "clesphys.h"

  • trunk/libf/phyvenus/radlwsw.NewtonCool

    r3 r97  
    4141#include "dimphy.h"
    4242#include "raddim.h"
    43 #include "YOETHF.h"
    4443#include "YOMCST.h"
    4544#include "clesphys.h"

  • trunk/libf/phyvenus/radlwsw.multimat

    r3 r97  
    4242#include "dimphy.h"
    4343#include "raddim.h"
    44 #include "YOETHF.h"
    4544#include "YOMCST.h"
    4645#include "clesphys.h"

  • trunk/libf/phyvenus/suphec.F

    r3 r97  
    55C
    66#include "YOMCST.h"
    7 #include "YOETHF.h"
    87cIM cf. JLD
    98       LOGICAL firstcall
    … …  
    153152      WRITE(UNIT=6,FMT='(''     Rv/Rd-1  = '',e13.7)') RETV
    154153C
    155 C     ----------------------------------------------------------------
    156 C
    157 C*       6.    DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.
    158 C              ---------------------------------------------
    159 C
    160       RCW=RCPV
    161       WRITE(UNIT=6,FMT='('' *** Thermodynamic, liquid  ***'')')
    162       WRITE(UNIT=6,FMT='(''         Cw   = '',E13.7)') RCW
    163 C
    164 C     ----------------------------------------------------------------
    165 C
    166 C*       7.    DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.
    167 C              --------------------------------------------
    168 C
    169       RCS=RCPV
    170       WRITE(UNIT=6,FMT='('' *** thermodynamic, solid   ***'')')
    171       WRITE(UNIT=6,FMT='(''         Cs   = '',E13.7)') RCS
    172 C
    173 C     ----------------------------------------------------------------
    174 C
    175 C*       8.    DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.
    176 C              ----------------------------------------------------
    177 C
    178       RTT=273.16
    179       RLVTT=2.5008E+6
    180       RLSTT=2.8345E+6
    181       RLMLT=RLSTT-RLVTT
    182       RATM=100000.
    183       WRITE(UNIT=6,FMT='('' *** Thermodynamic, trans.  ***'')')
    184       WRITE(UNIT=6,FMT='('' Fusion point  = '',E13.7)') RTT
    185       WRITE(UNIT=6,FMT='(''        RLvTt  = '',E13.7)') RLVTT
    186       WRITE(UNIT=6,FMT='(''        RLsTt  = '',E13.7)') RLSTT
    187       WRITE(UNIT=6,FMT='(''        RLMlt  = '',E13.7)') RLMLT
    188       WRITE(UNIT=6,FMT='('' Normal press. = '',E13.7)') RATM
    189       WRITE(UNIT=6,FMT='('' Latent heat :  '')')
    190 C
    191 C     ----------------------------------------------------------------
    192 C
    193 C*       9.    SATURATED VAPOUR PRESSURE.
    194 C              --------------------------
    195 C
    196       RESTT=611.14
    197       RGAMW=(RCW-RCPV)/RV
    198       RBETW=RLVTT/RV+RGAMW*RTT
    199       RALPW=LOG(RESTT)+RBETW/RTT+RGAMW*LOG(RTT)
    200       RGAMS=(RCS-RCPV)/RV
    201       RBETS=RLSTT/RV+RGAMS*RTT
    202       RALPS=LOG(RESTT)+RBETS/RTT+RGAMS*LOG(RTT)
    203       RGAMD=RGAMS-RGAMW
    204       RBETD=RBETS-RBETW
    205       RALPD=RALPS-RALPW
    206 C
    207 C     ------------------------------------------------------------------
    208 c
    209 c calculer les constantes pour les fonctions thermodynamiques
    210 c
    211       RVTMP2=RCPV/RCPD-1.
    212       RHOH2O=RATM/100.
    213       R2ES=RESTT*RD/RV
    214       R3LES=17.269
    215       R3IES=21.875
    216       R4LES=35.86
    217       R4IES=7.66
    218       R5LES=R3LES*(RTT-R4LES)
    219       R5IES=R3IES*(RTT-R4IES)
    220 C
    221154      RETURN
    222155      END

  • trunk/libf/phyvenus/testphys1d.F

    r3 r97  
    5454      REAL gru,grv   ! prescribed "geostrophic" background wind
    5555      REAL temp(llm)   ! temperature at the middle of the layers
    56       REAL q(llm,nqmx) ! tracer mixing ratio (e.g. kg/kg)
     56      REAL q(llm,nqtot) ! tracer mixing ratio (e.g. kg/kg)
    5757      REAL tsoil(nsoilmx)   ! subsurface soik temperature (K)
    5858      REAL zlay(llm)   ! altitude estimee dans les couches (km)
    … …  
    6666      REAL dudyn(llm),dvdyn(llm),dtempdyn(llm)
    6767      REAL dpsurf   
    68       REAL dq(llm,nqmx)
    69       REAL dqdyn(llm,nqmx)
     68      REAL dq(llm,nqtot)
     69      REAL dqdyn(llm,nqtot)
    7070
    7171c   Various intermediate variables
    … …  
    287287c    ---------------------------
    288288
    289       DO iq=1,nqmx
     289      DO iq=1,nqtot
    290290        DO ilayer=1,nlayer
    291291           q(ilayer,iq) = 0.
    … …  
    375375c       --------------------
    376376
    377       CALL physiq (1,llm,nqmx,
     377      CALL physiq (1,llm,nqtot,
    378378     ,     firstcall,lastcall,
    379379     ,     day,time,dtphys,

  • trunk/libf/phyvenus/ustarhb.F

    r3 r97  
    3636      LOGICAL check(klon)   ! True=>chk if Richardson no.>critcal
    3737c
    38 #include "YOETHF.h"
    39 #include "FCTTRE.h"
    4038      DO i = 1, knon
    4139        zx_alf1 = 1.0

  • trunk/libf/phyvenus/write_histday.h

    r3 r97  
    77      ndex3d = 0
    88      zx_tmp_2d = 0.
    9       z4_tmp_2d = 0.
    109      zx_tmp_3d = 0.
    11       z4_tmp_3d = 0.
    1210      zx_tmp_fi2d=0.
    1311      zx_tmp_fi3d=0.
    … …  
    2624      i = NINT(zout/zsto)
    2725      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
    28       z4_tmp_2d=zx_tmp_2d
    29       CALL histwrite(nid_day,"phis",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     26      CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    3027C
    3128      i = NINT(zout/zsto)
    3229      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
    33       z4_tmp_2d=zx_tmp_2d
    34       CALL histwrite(nid_day,"aire",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     30      CALL histwrite(nid_day,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    3531c
    3632      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ftsol,zx_tmp_2d)
    37       z4_tmp_2d=zx_tmp_2d
    38       CALL histwrite(nid_day,"tsol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     33      CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    3934c
    4035      DO i = 1, klon
    … …  
    4237      ENDDO
    4338      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    44       z4_tmp_2d=zx_tmp_2d
    45       CALL histwrite(nid_day,"psol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     39      CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    4640c
    4741      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
    48       z4_tmp_2d=zx_tmp_2d
    49       CALL histwrite(nid_day,"ue",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     42      CALL histwrite(nid_day,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    5043c
    5144      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
    52       z4_tmp_2d=zx_tmp_2d
    53       CALL histwrite(nid_day,"ve",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     45      CALL histwrite(nid_day,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    5446c
    5547      ENDIF !lev_histday.GE.1
    … …  
    6153c
    6254      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    63       z4_tmp_3d=zx_tmp_3d
    64       CALL histwrite(nid_day,"temp",itau_w,z4_tmp_3d,
     55      CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
    6556     .                                   iim*jjmp1*klev,ndex3d)
    6657c
    6758      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    68       z4_tmp_3d=zx_tmp_3d
    69       CALL histwrite(nid_day,"pres",itau_w,z4_tmp_3d,
     59      CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
    7060     .                                   iim*jjmp1*klev,ndex3d)
    7161c
    7262      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    73       z4_tmp_3d=zx_tmp_3d
    74       CALL histwrite(nid_day,"geop",itau_w,z4_tmp_3d,
     63      CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
    7564     .                                   iim*jjmp1*klev,ndex3d)
    7665c
    7766      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    78       z4_tmp_3d=zx_tmp_3d
    79       CALL histwrite(nid_day,"vitu",itau_w,z4_tmp_3d,
     67      CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
    8068     .                                   iim*jjmp1*klev,ndex3d)
    8169c
    8270      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    83       z4_tmp_3d=zx_tmp_3d
    84       CALL histwrite(nid_day,"vitv",itau_w,z4_tmp_3d,
     71      CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
    8572     .                                   iim*jjmp1*klev,ndex3d)
    8673c
    8774      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    88       z4_tmp_3d=zx_tmp_3d
    89       CALL histwrite(nid_day,"vitw",itau_w,z4_tmp_3d,
     75      CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
    9076     .                                   iim*jjmp1*klev,ndex3d)
    9177c
    … …  
    10086c
    10187      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
    102       z4_tmp_2d=zx_tmp_2d
    103       CALL histwrite(nid_day,"tops",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     88      CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    10489c
    10590      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    106       z4_tmp_2d=zx_tmp_2d
    107       CALL histwrite(nid_day,"topl",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     91      CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    10892c
    10993      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
    110       z4_tmp_2d=zx_tmp_2d
    111       CALL histwrite(nid_day,"sols",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     94      CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    11295c
    11396      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
    114       z4_tmp_2d=zx_tmp_2d
    115       CALL histwrite(nid_day,"soll",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     97      CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    11698c
    11799c 3D
    … …  
    119101      zx_tmp_fi3d(1:klon,1:klev)=swnet(1:klon,1:klev)
    120102      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    121       z4_tmp_3d=zx_tmp_3d
    122       CALL histwrite(nid_day,"SWnet",itau_w,z4_tmp_3d,
     103      CALL histwrite(nid_day,"SWnet",itau_w,zx_tmp_3d,
    123104     .                                   iim*jjmp1*klev,ndex3d)
    124105c
    125106      zx_tmp_fi3d(1:klon,1:klev)=lwnet(1:klon,1:klev)
    126107      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    127       z4_tmp_3d=zx_tmp_3d
    128       CALL histwrite(nid_day,"LWnet",itau_w,z4_tmp_3d,
     108      CALL histwrite(nid_day,"LWnet",itau_w,zx_tmp_3d,
    129109     .                                   iim*jjmp1*klev,ndex3d)
    130110c
    … …  
    136116c en K/s     
    137117      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
    138       z4_tmp_3d=zx_tmp_3d
    139       CALL histwrite(nid_day,"dtdyn",itau_w,z4_tmp_3d,
     118      CALL histwrite(nid_day,"dtdyn",itau_w,zx_tmp_3d,
    140119     .                                   iim*jjmp1*klev,ndex3d)
    141120c
    142121c en K/s     
    143122c      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
    144 c      z4_tmp_3d=zx_tmp_3d
    145 c      CALL histwrite(nid_day,"dtphy",itau_w,z4_tmp_3d,
     123c      CALL histwrite(nid_day,"dtphy",itau_w,zx_tmp_3d,
    146124c     .                                   iim*jjmp1*klev,ndex3d)
    147125c
    148126c en K/s     
    149127      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,d_t_vdf,zx_tmp_3d)
    150       z4_tmp_3d=zx_tmp_3d
    151       CALL histwrite(nid_day,"dtvdf",itau_w,z4_tmp_3d,
     128      CALL histwrite(nid_day,"dtvdf",itau_w,zx_tmp_3d,
    152129     .                                   iim*jjmp1*klev,ndex3d)
    153130c
    154131c en K/s     
    155132      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,d_t_ajs,zx_tmp_3d)
    156       z4_tmp_3d=zx_tmp_3d
    157       CALL histwrite(nid_day,"dtajs",itau_w,z4_tmp_3d,
     133      CALL histwrite(nid_day,"dtajs",itau_w,zx_tmp_3d,
    158134     .                                   iim*jjmp1*klev,ndex3d)
    159135c
    … …  
    161137      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
    162138      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    163       z4_tmp_3d=zx_tmp_3d
    164       CALL histwrite(nid_day,"dtswr",itau_w,z4_tmp_3d,
     139      CALL histwrite(nid_day,"dtswr",itau_w,zx_tmp_3d,
    165140     .                                   iim*jjmp1*klev,ndex3d)
    166141c
    … …  
    168143      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
    169144      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    170       z4_tmp_3d=zx_tmp_3d
    171       CALL histwrite(nid_day,"dtlwr",itau_w,z4_tmp_3d,
     145      CALL histwrite(nid_day,"dtlwr",itau_w,zx_tmp_3d,
    172146     .                                   iim*jjmp1*klev,ndex3d)
    173147c en K/s     
    174148c     zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
    175149c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    176 c     z4_tmp_3d=zx_tmp_3d
    177 c     CALL histwrite(nid_day,"dtec",itau_w,z4_tmp_3d,
     150c     CALL histwrite(nid_day,"dtec",itau_w,zx_tmp_3d,
    178151c    .                                   iim*jjmp1*klev,ndex3d)
    179152c
    180153c en (m/s)/s     
    181154      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_dyn, zx_tmp_3d)
    182       z4_tmp_3d=zx_tmp_3d
    183       CALL histwrite(nid_day,"dudyn",itau_w,z4_tmp_3d,
     155      CALL histwrite(nid_day,"dudyn",itau_w,zx_tmp_3d,
    184156     .                                   iim*jjmp1*klev,ndex3d)
    185157c
    186158c en (m/s)/s     
    187159      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
    188       z4_tmp_3d=zx_tmp_3d
    189       CALL histwrite(nid_day,"duvdf",itau_w,z4_tmp_3d,
     160      CALL histwrite(nid_day,"duvdf",itau_w,zx_tmp_3d,
    190161     .                                   iim*jjmp1*klev,ndex3d)
    191162c
    192163c en (m/s)/s     
    193164      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_ajs, zx_tmp_3d)
    194       z4_tmp_3d=zx_tmp_3d
    195       CALL histwrite(nid_day,"duajs",itau_w,z4_tmp_3d,
     165      CALL histwrite(nid_day,"duajs",itau_w,zx_tmp_3d,
    196166     .                                   iim*jjmp1*klev,ndex3d)
    197167c
    198168c en (m/s)/s     
    199169      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d)
    200       z4_tmp_3d=zx_tmp_3d
    201       CALL histwrite(nid_day,"dugwo",itau_w,z4_tmp_3d,
     170      CALL histwrite(nid_day,"dugwo",itau_w,zx_tmp_3d,
    202171     .                                   iim*jjmp1*klev,ndex3d)
    203172c
    204173c en (m/s)/s     
    205174      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_hin, zx_tmp_3d)
    206       z4_tmp_3d=zx_tmp_3d
    207       CALL histwrite(nid_day,"dugwno",itau_w,z4_tmp_3d,
     175      CALL histwrite(nid_day,"dugwno",itau_w,zx_tmp_3d,
    208176     .                                   iim*jjmp1*klev,ndex3d)
    209177c
    … …  
    211179c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
    212180c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
    213 c     z4_tmp_3d=-1.*zx_tmp_3d
    214 c     CALL histwrite(nid_day,"dvvdf",itau_w,z4_tmp_3d,
     181c     zx_tmp_3d=-1.*zx_tmp_3d
     182c     CALL histwrite(nid_day,"dvvdf",itau_w,zx_tmp_3d,
    215183c    .                                   iim*jjmp1*klev,ndex3d)
    216184c
    … …  
    224192         IF (iq.LE.99) THEN
    225193      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
    226       z4_tmp_3d=zx_tmp_3d
    227194      WRITE(str2,'(i2.2)') iq
    228       CALL histwrite(nid_day,"trac"//str2,itau_w,z4_tmp_3d,
     195      CALL histwrite(nid_day,"trac"//str2,itau_w,zx_tmp_3d,
    229196     .                                   iim*jjmp1*klev,ndex3d)
    230197         ELSE
    … …  
    236203c
    237204c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxu , zx_tmp_2d)
    238 c      z4_tmp_2d=zx_tmp_2d
    239205c      CALL histwrite(nid_day,"taux_",itau_w,
    240 c    $      z4_tmp_2d,iim*jjmp1,ndex2d)
     206c    $      zx_tmp_2d,iim*jjmp1,ndex2d)
    241207c     
    242208c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d)
    243 c      z4_tmp_2d=zx_tmp_2d
    244209c      CALL histwrite(nid_day,"tauy_",itau_w,
    245 c    $      z4_tmp_2d,iim*jjmp1,ndex2d)
     210c    $      zx_tmp_2d,iim*jjmp1,ndex2d)
    246211c
    247212c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    248 c      z4_tmp_2d=zx_tmp_2d
    249 c     CALL histwrite(nid_day,"cdrm",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     213c     CALL histwrite(nid_day,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    250214c
    251215c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    252 c      z4_tmp_2d=zx_tmp_2d
    253 c     CALL histwrite(nid_day,"cdrh",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     216c     CALL histwrite(nid_day,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    254217c
    255218      ENDIF !lev_histday.GE.5

  • trunk/libf/phyvenus/write_histins.h

    r3 r97  
    77      ndex3d = 0
    88      zx_tmp_2d = 0.
    9       z4_tmp_2d = 0.
    109      zx_tmp_3d = 0.
    11       z4_tmp_3d = 0.
    1210      zx_tmp_fi2d=0.
    1311      zx_tmp_fi3d=0.

  • trunk/libf/phyvenus/write_histmth.h

    r3 r97  
    77      ndex3d = 0
    88      zx_tmp_2d = 0.
    9       z4_tmp_2d = 0.
    109      zx_tmp_3d = 0.
    11       z4_tmp_3d = 0.
    1210      zx_tmp_fi2d=0.
    1311      zx_tmp_fi3d=0.
    … …  
    2523      i = NINT(zout/zsto)
    2624      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
    27       z4_tmp_2d=zx_tmp_2d
    28       CALL histwrite(nid_mth,"phis",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     25      CALL histwrite(nid_mth,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    2926C
    3027      i = NINT(zout/zsto)
    3128      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
    32       z4_tmp_2d=zx_tmp_2d
    33       CALL histwrite(nid_mth,"aire",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     29      CALL histwrite(nid_mth,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    3430c
    3531      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ftsol,zx_tmp_2d)
    36       z4_tmp_2d=zx_tmp_2d
    37       CALL histwrite(nid_mth,"tsol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     32      CALL histwrite(nid_mth,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    3833c
    3934      DO i = 1, klon
    … …  
    4136      ENDDO
    4237      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    43       z4_tmp_2d=zx_tmp_2d
    44       CALL histwrite(nid_mth,"psol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     38      CALL histwrite(nid_mth,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    4539c
    4640c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
    47 c     z4_tmp_2d=zx_tmp_2d
    48 c     CALL histwrite(nid_mth,"ue",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     41c     CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    4942c
    5043c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
    5144c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
    52 c     z4_tmp_2d=-1.*zx_tmp_2d
    53 c     CALL histwrite(nid_mth,"ve",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     45c     zx_tmp_2d=-1.*zx_tmp_2d
     46c     CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    5447c
    5548      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    56       z4_tmp_2d=zx_tmp_2d
    57       CALL histwrite(nid_mth,"cdragh",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     49      CALL histwrite(nid_mth,"cdragh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    5850c
    5951      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    60       z4_tmp_2d=zx_tmp_2d
    61       CALL histwrite(nid_mth,"cdragm",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     52      CALL histwrite(nid_mth,"cdragm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    6253c
    6354      ENDIF !lev_histmth.GE.1
    … …  
    6960c
    7061      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    71       z4_tmp_3d=zx_tmp_3d
    72       CALL histwrite(nid_mth,"temp",itau_w,z4_tmp_3d,
     62      CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
    7363     .                                   iim*jjmp1*klev,ndex3d)
    7464c
    7565      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    76       z4_tmp_3d=zx_tmp_3d
    77       CALL histwrite(nid_mth,"pres",itau_w,z4_tmp_3d,
     66      CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
    7867     .                                   iim*jjmp1*klev,ndex3d)
    7968c
    8069      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    81       z4_tmp_3d=zx_tmp_3d
    82       CALL histwrite(nid_mth,"geop",itau_w,z4_tmp_3d,
     70      CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
    8371     .                                   iim*jjmp1*klev,ndex3d)
    8472c
    8573      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    86       z4_tmp_3d=zx_tmp_3d
    87       CALL histwrite(nid_mth,"vitu",itau_w,z4_tmp_3d,
     74      CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
    8875     .                                   iim*jjmp1*klev,ndex3d)
    8976c
    9077      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    9178c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
    92       z4_tmp_3d=-1.*zx_tmp_3d
    93       CALL histwrite(nid_mth,"vitv",itau_w,z4_tmp_3d,
     79      zx_tmp_3d=-1.*zx_tmp_3d
     80      CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
    9481     .                                   iim*jjmp1*klev,ndex3d)
    9582c
    9683      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    97       z4_tmp_3d=zx_tmp_3d
    98       CALL histwrite(nid_mth,"vitw",itau_w,z4_tmp_3d,
     84      CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
    9985     .                                   iim*jjmp1*klev,ndex3d)
    10086c
    10187c en (m/s)/s     
    10288      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_dyn, zx_tmp_3d)
    103       z4_tmp_3d=zx_tmp_3d
    104       CALL histwrite(nid_mth,"dudyn",itau_w,z4_tmp_3d,
     89      CALL histwrite(nid_mth,"dudyn",itau_w,zx_tmp_3d,
    10590     .                                   iim*jjmp1*klev,ndex3d)
    10691c
    10792c en (m/s)/s     
    10893      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
    109       z4_tmp_3d=zx_tmp_3d
    110       CALL histwrite(nid_mth,"duvdf",itau_w,z4_tmp_3d,
     94      CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
    11195     .                                   iim*jjmp1*klev,ndex3d)
    11296c
    11397c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, mang, zx_tmp_3d)
    114 c     z4_tmp_3d=zx_tmp_3d
    115 c     CALL histwrite(nid_mth,"mang",itau_w,z4_tmp_3d,
     98c     CALL histwrite(nid_mth,"mang",itau_w,zx_tmp_3d,
    11699c    .                                   iim*jjmp1*klev,ndex3d)
    117100c
    118101      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ycoefh, zx_tmp_3d)
    119       z4_tmp_3d=zx_tmp_3d
    120       CALL histwrite(nid_mth,"Kz",itau_w,z4_tmp_3d,
     102      CALL histwrite(nid_mth,"Kz",itau_w,zx_tmp_3d,
    121103     .                                   iim*jjmp1*klev,ndex3d)
    122104c
    … …  
    126108          IF (iq.LE.99) THEN
    127109       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
    128        z4_tmp_3d=zx_tmp_3d
    129110       WRITE(str2,'(i2.2)') iq
    130        CALL histwrite(nid_mth,tnom(iq),itau_w,z4_tmp_3d,
     111       CALL histwrite(nid_mth,tnom(iq),itau_w,zx_tmp_3d,
    131112     .                                   iim*jjmp1*klev,ndex3d)
    132113          ELSE
    … …  
    138119c
    139120      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
    140       z4_tmp_2d=zx_tmp_2d
    141       CALL histwrite(nid_mth,"tops",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     121      CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    142122c
    143123      ENDIF !lev_histmth.GE.2
    … …  
    151131c
    152132      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
    153       z4_tmp_2d=zx_tmp_2d
    154       CALL histwrite(nid_mth,"topl",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     133      CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    155134c
    156135      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
    157       z4_tmp_2d=zx_tmp_2d
    158       CALL histwrite(nid_mth,"sols",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     136      CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    159137c
    160138      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
    161       z4_tmp_2d=zx_tmp_2d
    162       CALL histwrite(nid_mth,"soll",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     139      CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    163140c
    164141c 3D
    … …  
    166143      zx_tmp_fi3d(1:klon,1:klev)=swnet(1:klon,1:klev)
    167144      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    168       z4_tmp_3d=zx_tmp_3d
    169       CALL histwrite(nid_mth,"SWnet",itau_w,z4_tmp_3d,
     145      CALL histwrite(nid_mth,"SWnet",itau_w,zx_tmp_3d,
    170146     .                                   iim*jjmp1*klev,ndex3d)
    171147c
    172148      zx_tmp_fi3d(1:klon,1:klev)=lwnet(1:klon,1:klev)
    173149      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    174       z4_tmp_3d=zx_tmp_3d
    175       CALL histwrite(nid_mth,"LWnet",itau_w,z4_tmp_3d,
     150      CALL histwrite(nid_mth,"LWnet",itau_w,zx_tmp_3d,
    176151     .                                   iim*jjmp1*klev,ndex3d)
    177152c
    178153      zx_tmp_fi3d(1:klon,1:klev)=fluxt(1:klon,1:klev)
    179154      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    180       z4_tmp_3d=zx_tmp_3d
    181       CALL histwrite(nid_mth,"fluxvdf",itau_w,z4_tmp_3d,
     155      CALL histwrite(nid_mth,"fluxvdf",itau_w,zx_tmp_3d,
    182156     .                                   iim*jjmp1*klev,ndex3d)
    183157c
    184158      zx_tmp_fi3d(1:klon,1:klev)=flux_dyn(1:klon,1:klev)
    185159      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    186       z4_tmp_3d=zx_tmp_3d
    187       CALL histwrite(nid_mth,"fluxdyn",itau_w,z4_tmp_3d,
     160      CALL histwrite(nid_mth,"fluxdyn",itau_w,zx_tmp_3d,
    188161     .                                   iim*jjmp1*klev,ndex3d)
    189162c
    190163      zx_tmp_fi3d(1:klon,1:klev)=flux_ajs(1:klon,1:klev)
    191164      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    192       z4_tmp_3d=zx_tmp_3d
    193       CALL histwrite(nid_mth,"fluxajs",itau_w,z4_tmp_3d,
     165      CALL histwrite(nid_mth,"fluxajs",itau_w,zx_tmp_3d,
    194166     .                                   iim*jjmp1*klev,ndex3d)
    195167c
    196168c     zx_tmp_fi3d(1:klon,1:klev)=flux_ec(1:klon,1:klev)
    197169c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    198 c     z4_tmp_3d=zx_tmp_3d
    199 c     CALL histwrite(nid_mth,"fluxec",itau_w,z4_tmp_3d,
     170c     CALL histwrite(nid_mth,"fluxec",itau_w,zx_tmp_3d,
    200171c    .                                   iim*jjmp1*klev,ndex3d)
    201172c
    … …  
    207178c en K/s     
    208179      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
    209       z4_tmp_3d=zx_tmp_3d
    210       CALL histwrite(nid_mth,"dtdyn",itau_w,z4_tmp_3d,
     180      CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
    211181     .                                   iim*jjmp1*klev,ndex3d)
    212182c
    213183c en K/s     
    214184c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
    215 c     z4_tmp_3d=zx_tmp_3d
    216 c     CALL histwrite(nid_mth,"dtphy",itau_w,z4_tmp_3d,
     185c     CALL histwrite(nid_mth,"dtphy",itau_w,zx_tmp_3d,
    217186c    .                                   iim*jjmp1*klev,ndex3d)
    218187c
    219188c en K/s     
    220189      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,d_t_vdf,zx_tmp_3d)
    221       z4_tmp_3d=zx_tmp_3d
    222       CALL histwrite(nid_mth,"dtvdf",itau_w,z4_tmp_3d,
     190      CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
    223191     .                                   iim*jjmp1*klev,ndex3d)
    224192c
    225193c en K/s     
    226194      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,d_t_ajs,zx_tmp_3d)
    227       z4_tmp_3d=zx_tmp_3d
    228       CALL histwrite(nid_mth,"dtajs",itau_w,z4_tmp_3d,
     195      CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
    229196     .                                   iim*jjmp1*klev,ndex3d)
    230197c
    … …  
    232199      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
    233200      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    234       z4_tmp_3d=zx_tmp_3d
    235       CALL histwrite(nid_mth,"dtswr",itau_w,z4_tmp_3d,
     201      CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
    236202     .                                   iim*jjmp1*klev,ndex3d)
    237203c
    … …  
    239205      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
    240206      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    241       z4_tmp_3d=zx_tmp_3d
    242       CALL histwrite(nid_mth,"dtlwr",itau_w,z4_tmp_3d,
     207      CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
    243208     .                                   iim*jjmp1*klev,ndex3d)
    244209c en K/s     
    245210c     zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
    246211c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    247 c     z4_tmp_3d=zx_tmp_3d
    248 c     CALL histwrite(nid_mth,"dtec",itau_w,z4_tmp_3d,
     212c     CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
    249213c    .                                   iim*jjmp1*klev,ndex3d)
    250214c
    251215c en (m/s)/s     
    252216      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_ajs, zx_tmp_3d)
    253       z4_tmp_3d=zx_tmp_3d
    254       CALL histwrite(nid_mth,"duajs",itau_w,z4_tmp_3d,
     217      CALL histwrite(nid_mth,"duajs",itau_w,zx_tmp_3d,
    255218     .                                   iim*jjmp1*klev,ndex3d)
    256219c
    257220c en (m/s)/s     
    258221      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d)
    259       z4_tmp_3d=zx_tmp_3d
    260       CALL histwrite(nid_mth,"dugwo",itau_w,z4_tmp_3d,
     222      CALL histwrite(nid_mth,"dugwo",itau_w,zx_tmp_3d,
    261223     .                                   iim*jjmp1*klev,ndex3d)
    262224c
    263225c en (m/s)/s     
    264226      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_hin, zx_tmp_3d)
    265       z4_tmp_3d=zx_tmp_3d
    266       CALL histwrite(nid_mth,"dugwno",itau_w,z4_tmp_3d,
     227      CALL histwrite(nid_mth,"dugwno",itau_w,zx_tmp_3d,
    267228     .                                   iim*jjmp1*klev,ndex3d)
    268229c
    … …  
    270231c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
    271232c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
    272 c     z4_tmp_3d=-1.*zx_tmp_3d
    273 c     CALL histwrite(nid_mth,"dvvdf",itau_w,z4_tmp_3d,
     233c     zx_tmp_3d=-1.*zx_tmp_3d
     234c     CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
    274235c    .                                   iim*jjmp1*klev,ndex3d)
    275236c
    … …  
    281242c
    282243c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxu , zx_tmp_2d)
    283 c      z4_tmp_2d=zx_tmp_2d
    284244c      CALL histwrite(nid_mth,"taux_",itau_w,
    285 c    $      z4_tmp_2d,iim*jjmp1,ndex2d)
     245c    $      zx_tmp_2d,iim*jjmp1,ndex2d)
    286246c     
    287247c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d)
    288 c      z4_tmp_2d=zx_tmp_2d
    289248c      CALL histwrite(nid_mth,"tauy_",itau_w,
    290 c    $      z4_tmp_2d,iim*jjmp1,ndex2d)
     249c    $      zx_tmp_2d,iim*jjmp1,ndex2d)
    291250c
    292251c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    293 c      z4_tmp_2d=zx_tmp_2d
    294 c     CALL histwrite(nid_mth,"cdrm",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     252c     CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    295253c
    296254c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    297 c      z4_tmp_2d=zx_tmp_2d
    298 c     CALL histwrite(nid_mth,"cdrh",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     255c     CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    299256c
    300257      ENDIF !lev_histmth.GE.5
Note: See TracChangeset for help on using the changeset viewer.