Ignore:
Timestamp:
Jul 21, 2024, 1:07:18 PM (4 months ago)
Author:
abarral
Message:

Use latest FCM source (2021.05.0) [Note: we still use the legacy FCM1 build system]
Correct UTF8 encoding of french chars


Compil OK (tested: oldrad/rrtm/ecrad, para/seq/1D)
Convergence (ref r5063) bench 33x OK oldrad orch2.0 (tested: para/seq)

Location:
LMDZ6/branches/Amaury_dev/libf/phylmdiso
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmdiso/add_wake_tend.F90

    r5082 r5093  
    55        )
    66!===================================================================
    7 ! Ajoute les tendances lies aux diverses parametrisations physiques aux
     7! Ajoute les tendances liées aux diverses parametrisations physiques aux
    88! variables d'etat des poches froides.
    99!===================================================================
  • LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90

    r5086 r5093  
    127127  ! wdens_ref: initial number of wakes per unit area (3D) or per
    128128  ! unit length (2D), at the beginning of each time step
    129   ! Tgw    : 1 sur la p�riode de onde de gravit�
    130   ! Cgw    : vitesse de propagation de onde de gravit
     129  ! Tgw    : 1 sur la période de onde de gravité
     130  ! Cgw    : vitesse de propagation de onde de gravité
    131131  ! LL     : distance entre 2 poches
    132132
    133133  ! -------------------------------------------------------------------------
    134   ! Dclaration de variables
     134  ! Déclaration de variables
    135135  ! -------------------------------------------------------------------------
    136136
     
    196196  ! -------------------
    197197
    198   ! Variables fixer
     198  ! Variables à fixer
    199199  INTEGER, SAVE                                         :: igout
    200200  !$OMP THREADPRIVATE(igout)
     
    383383  ! Configuration de coefgw,stark,wdens (22/02/06 by YU Jingmei)
    384384
    385   ! coefgw : Coefficient pour les ondes de gravit
     385  ! coefgw : Coefficient pour les ondes de gravité
    386386  ! stark : Coefficient k dans Cstar=k*sqrt(2*WAPE)
    387   ! wdens : Densit surfacique de poche froide
     387  ! wdens : Densité surfacique de poche froide
    388388  ! -------------------------------------------------------------------------
    389389
     
    10881088
    10891089    ! cc nrlmd   Ajout d'un recalcul de wdens dans le cas d'un entrainement
    1090     ! n�gatif de ktop � kupper --------
    1091     ! cc           On calcule pour cela une densit wdens0 pour laquelle on
     1090    ! négatif de ktop à kupper --------
     1091    ! cc           On calcule pour cela une densité wdens0 pour laquelle on
    10921092    ! aurait un entrainement nul ---
    10931093    !jyg<
     
    10961096    ! des descentes unsaturees. Nous faisons alors l'hypothese que la
    10971097    ! convection profonde cree directement de nouvelles poches, sans passer
    1098     ! par les thermiques. La nouvelle valeur de wdens est alors impose.
     1098    ! par les thermiques. La nouvelle valeur de wdens est alors imposée.
    10991099
    11001100    DO i = 1, klon
     
    11951195      DO i = 1, klon
    11961196        IF (wk_adv(i)) THEN
    1197           ! cc nrlmd          Introduction du taux de mortalit des poches et
     1197          ! cc nrlmd          Introduction du taux de mortalité des poches et
    11981198          ! test sur sigmaw_max=0.4
    11991199          ! cc         d_sigmaw(i) = gfl(i)*Cstar(i)*dtimesub
     
    13081308
    13091309    ! c      DO i=1,klon
    1310     ! c        print*,'Pente entre 0 et kupper (r�f�rence)'
     1310    ! c        print*,'Pente entre 0 et kupper (référence)'
    13111311    ! c     $           ,omg(i,kupper(i)+1)/(pupper(i)-ph(i,1))
    13121312    ! c        print*,'Pente entre ktop et kupper'
     
    16021602
    16031603
    1604           ! Coefficient de rpartition
     1604          ! Coefficient de répartition
    16051605
    16061606          crep(i, k) = crep_sol*(ph(i,kupper(i))-ph(i,k))/ &
     
    16461646!
    16471647
    1648           ! cc nrlmd          Prise en compte du taux de mortalit
    1649           ! cc               Dfinitions de entr, detr
     1648          ! cc nrlmd          Prise en compte du taux de mortalité
     1649          ! cc               Définitions de entr, detr
    16501650!jyg<
    16511651!!            detr(i, k) = 0.
     
    16641664
    16651665
    1666           ! ajout d'un effet onde de gravit -Tgw(k)*deltatw(k) 03/02/06 YU
     1666          ! ajout d'un effet onde de gravité -Tgw(k)*deltatw(k) 03/02/06 YU
    16671667          ! Jingmei
    16681668
     
    26012601        discrim = b*b - 4.*a*c
    26022602        ! print*, 'x, a, b, c, discrim', x, a, b, c, discrim
    2603         IF (a+b>=0.) THEN !! Condition suffisante pour la positivit de ovap
     2603        IF (a+b>=0.) THEN !! Condition suffisante pour la positivité de ovap
    26042604          alpha1(i) = 1.
    26052605        ELSE
  • LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyaqua_mod.F90

    r5090 r5093  
    963963      if (type_profil==20) then
    964964      print*,'Profile SST 20'
    965 !     Méthode 13 "Qmax2K" plateau réel |  l'Equateur augmenté +2K
     965!     Méthode 13 "Qmax2K" plateau réel é|  l'Equateur augmenté +2K
    966966
    967967      do j=1,klon
     
    972972      if (type_profil==21) then
    973973      print*,'Profile SST 21'
    974 !     Méthode 13 "Qmax2K" plateau réel |  l'Equateur augmenté +2K
     974!     Méthode 13 "Qmax2K" plateau réel é|  l'Equateur augmenté +2K
    975975      do j=1,klon
    976976        phy_sst(j,i)=252.+55.*(1-sin(rlatd(j))**2)
Note: See TracChangeset for help on using the changeset viewer.