Ignore:
Timestamp:
Aug 17, 2006, 5:41:51 PM (18 years ago)
Author:
Laurent Fairhead
Message:

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/V3_test/libf/phylmd/fisrtilp.F

    r559 r704  
    1010
    1111c
     12      USE dimphy
    1213      IMPLICIT none
    1314c======================================================================
     
    1819c======================================================================
    1920c======================================================================
    20 #include "dimensions.h"
    21 #include "dimphy.h"
     21cym#include "dimensions.h"
     22cym#include "dimphy.h"
    2223#include "YOMCST.h"
    2324#include "tracstoke.h"
     
    9091      LOGICAL appel1er
    9192      SAVE appel1er
     93c$OMP THREADPRIVATE(appel1er)
    9294c
    9395c---------------------------------------------------------------
     
    99101      REAL a_tr_sca(4)
    100102      save a_tr_sca
     103c$OMP THREADPRIVATE(a_tr_sca)
    101104c
    102105c Variables intermediaires
     
    110113      REAL zlh_solid(klon), zm_solid
    111114cIM
    112       INTEGER klevm1
     115cym      INTEGER klevm1
    113116c---------------------------------------------------------------
    114117c
     
    210213c
    211214cIM : klevm1
    212       klevm1=klev-1
     215cym      klevm1=klev-1
    213216      DO 9999 k = klev, 1, -1
    214217c
     
    269272     .                            /RG/dtime
    270273
    271 c pour la glace, on réévapore toute la précip dans la couche du dessous
     274c pour la glace, on r�vapore toute la pr�ip dans la couche du dessous
    272275c la glace venant de la couche du dessus est simplement dans la couche
    273276c du dessous.
     
    322325c   zqn   : eau totale dans le nuage
    323326c   zcond : eau condensee moyenne dans la maille.
    324 c           on prend en compte le réchauffement qui diminue la partie condensee
     327c           on prend en compte le rhauffement qui diminue la partie condensee
    325328c
    326329c   Version avec les raqts
     
    427430            zct(i)=1./cld_tau_lsc
    428431         endif
    429 c  quantité d'eau à élminier.
     432c  quantit�d'eau ��minier.
    430433         zchau(i) = zct(i)*dtime/FLOAT(ninter) * zoliq(i)
    431434     .         *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl(i))**2)) *(1.-zfice(i))
Note: See TracChangeset for help on using the changeset viewer.