Ignore:
Timestamp:
Jul 21, 2024, 1:07:18 PM (7 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
Files:
53 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/bilan_dyn.F

    r5082 r5093  
    2626c====================================================================
    2727c
    28 c   Sous-programme consacre � des diagnostics dynamiques de base
     28c   Sous-programme consacre à des diagnostics dynamiques de base
    2929c
    3030c
    … …  
    8989      real ww
    9090
    91 c   variables dynamiques interm�diaires
     91c   variables dynamiques intermédiaires
    9292      REAL vcont(iip1,jjm,llm),ucont(iip1,jjp1,llm)
    9393      REAL ang(iip1,jjp1,llm),unat(iip1,jjp1,llm)
    … …  
    9797      REAL bern(iip1,jjp1,llm)
    9898
    99 c   champ contenant les scalaires advect�s.
     99c   champ contenant les scalaires advectés.
    100100      real Q(iip1,jjp1,llm,nQ)
    101101   
    102 c   champs cumul�s
     102c   champs cumulés
    103103      real ps_cum(iip1,jjp1)
    104104      real masse_cum(iip1,jjp1,llm)
    … …  
    308308c   ----------------------------
    309309
    310 c   �nergie cin�tique
     310c   énergie cinétique
    311311      ucont(:,:,:)=0
    312312      CALL covcont(llm,ucov,vcov,ucont,vcont)
    313313      CALL enercin(vcov,ucov,vcont,ucont,ecin)
    314314
    315 c   moment cin�tique
     315c   moment cinétique
    316316      do l=1,llm
    317317         ang(:,:,l)=ucov(:,:,l)+constang(:,:)
    … …  
    373373      enddo
    374374
    375 c    flux m�ridien
     375c    flux méridien
    376376c    -------------
    377377      do iQ=1,nQ
    … …  
    440440 
    441441c=====================================================================
    442 c   Transport m�ridien
     442c   Transport méridien
    443443c=====================================================================
    444444

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/guide_mod.F90

    r5088 r5093  
    9999    CALL getpar('guide_teta',.false.,guide_teta,'guidage de T par Teta')
    100100
    101     CALL getpar('guide_add',.false.,guide_add,'for�age constant?')
     101    CALL getpar('guide_add',.false.,guide_add,'forçage constant?')
    102102    CALL getpar('guide_zon',.false.,guide_zon,'guidage moy zonale')
    103103    if (guide_zon .and. abs(grossismx - 1.) > 0.01) &
    … …  
    121121
    122122
    123 ! Sauvegarde du for�age
     123! Sauvegarde du forçage
    124124    CALL getpar('guide_sav',.false.,guide_sav,'sauvegarde guidage')
    125125    CALL getpar('iguide_sav',4,iguide_sav,'freq. sauvegarde guidage')
    … …  
    17191719        call nf95_put_var(nid, varid_alpha_q, alpha_q)
    17201720! --------------------------------------------------------------------
    1721 ! Cr�ation des variables sauvegard�es
     1721! Création des variables sauvegardées
    17221722! --------------------------------------------------------------------
    17231723        ierr = nf90_redef(nid)

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j.F

    r5086 r5093  
    5757c=====================================================================
    5858      if (lnew) then
    59 c   on r�initialise les r�indicages et les poids
     59c   on réinitialise les réindicages et les poids
    6060c=====================================================================
    6161
    … …  
    9696c
    9797c ...   Modif . P. Le Van    ( 20/01/98) ....
    98 c       Modif Fr�d�ric Hourdin (3/01/02)
     98c       Modif Frédéric Hourdin (3/01/02)
    9999
    100100        IF(pgcm(i,lb(i))==0.OR.

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j1.F

    r5086 r5093  
    5757c=====================================================================
    5858      if (lnew) then
    59 c   on r�initialise les r�indicages et les poids
     59c   on réinitialise les réindicages et les poids
    6060c=====================================================================
    6161
    … …  
    9696c
    9797c ...   Modif . P. Le Van    ( 20/01/98) ....
    98 c       Modif Fr�d�ric Hourdin (3/01/02)
     98c       Modif Frédéric Hourdin (3/01/02)
    9999
    100100        IF(pgcm(i,lb(i))==0.OR.

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/top_bound.F

    r5082 r5093  
    2727c   ------
    2828c
    29 c   Dissipation lin�aire (ex top_bound de la physique)
     29c   Dissipation linéaire (ex top_bound de la physique)
    3030c
    3131c=======================================================================

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/comdissnew.h

    r1952 r5093  
    55!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
    66!                 veillez à n'utiliser que des ! pour les commentaires
    7 !                 et à bien positionner les & des lignes de continuation 
     7!                 et à bien positionner les & des lignes de continuation
    88!                 (les placer en colonne 6 et en colonne 73)
    99!

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.F

    r5086 r5093  
    2424c
    2525c=======================================================================
    26 c    Discretisation verticale en coordonn�e hybride (ou sigma)
     26c    Discretisation verticale en coordonnée hybride (ou sigma)
    2727c
    2828c=======================================================================
    … …  
    208208c     Calcul au milieu des couches :
    209209c     WARNING : le choix de placer le milieu des couches au niveau de
    210 c     pression interm�diaire est arbitraire et pourrait etre modifi�.
     210c     pression intermédiaire est arbitraire et pourrait etre modifié.
    211211c     Le calcul du niveau pour la derniere couche
    212212c     (on met la meme distance (en log pression)  entre P(llm)
    213213c     et P(llm -1) qu'entre P(llm-1) et P(llm-2) ) est
    214 c     Specifique.  Ce choix est sp�cifi� ici ET dans exner_milieu.F
     214c     Specifique.  Ce choix est spécifié ici ET dans exner_milieu.F
    215215
    216216      DO l = 1, llm-1
    … …  
    284284c     L'objectif est de calculer newsig telle que
    285285c       (1 -pa/preff)*exp(1-1./newsig**2)+(pa/preff)*newsig = sig
    286 c     Cela ne se r�soud pas analytiquement:
    287 c     => on r�soud par iterration bourrine
     286c     Cela ne se résoud pas analytiquement:
     287c     => on résoud par iterration bourrine
    288288c     ----------------------------------------------
    289289c     Information  : where exp(1-1./x**2) become << x
    … …  
    320320              newsig=(X2+newsig)*0.5
    321321          end if
    322 c         Test : on arete lorsque on approxime sig � moins de 0.01 m pr�s
     322c         Test : on arete lorsque on approxime sig à moins de 0.01 m près
    323323c         (en pseudo altitude) :
    324324          IF(abs(10.*log(F))<1.E-5) goto 999

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpre.F

    r5086 r5093  
    4141
    4242c       CALCUL DE LA PRESSION DE SURFACE
    43 c       Les coefficients ap et bp sont pass�s en common
    44 c       Calcul de la pression au sol en mb optimis�e pour
     43c       Les coefficients ap et bp sont passés en common
     44c       Calcul de la pression au sol en mb optimisée pour
    4545c       la vectorialisation
    4646                   
    … …  
    6767c RECONSTRUCTION DES CHAMPS CONTRAVARIANTS
    6868c Le programme ppm3d travaille avec les composantes
    69 c de vitesse et pas les flux, on doit donc passer de l'un � l'autre
    70 c Dans le m�me temps, on fait le changement d'orientation du vent en v
     69c de vitesse et pas les flux, on doit donc passer de l'un à l'autre
     70c Dans le même temps, on fait le changement d'orientation du vent en v
    7171      do l=1,llm
    7272          do j=1,jjm
    … …  
    9999     
    100100c INVERSION DES NIVEAUX
    101 c le programme ppm3d travaille avec une 3�me coordonn�e invers�e par rapport
     101c le programme ppm3d travaille avec une 3ème coordonnée inversée par rapport
    102102c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
    103 c On passe donc des niveaux du LMDZ � ceux de Lin
     103c On passe donc des niveaux du LMDZ à ceux de Lin
    104104     
    105105      do l=1,llm+1

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ppm3d.F

    r5079 r5093  
    6868      implicit none
    6969
    70 c     rajout de d�clarations
     70c     rajout de déclarations
    7171c      integer Jmax,kmax,ndt0,nstep,k,j,i,ic,l,js,jn,imh,iad,jad,krd
    7272c      integer iu,iiu,j2,jmr,js0,jt
    … …  
    769769      DO k=1,NLAY
    770770      DO I=1,IMR
    771 c     j=1 c'est le p�le Sud, j=JNP c'est le p�le Nord
     771c     j=1 c'est le pôle Sud, j=JNP c'est le pôle Nord
    772772      Q(I,  2,k,IC) = Q(I,  1,k,IC)
    773773      Q(I,JMR,k,IC) = Q(I,JNP,k,IC)

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/traceurpole.F

    r2622 r5093  
    2828
    2929   
    30 c On impose une seule valeur au pôle Sud j=jjm+1=jjp1       
     30c On impose une seule valeur au pôle Sud j=jjm+1=jjp1
    3131      sommemasses=0
    3232      sommemqs=0

  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/bilan_dyn_loc.F

    r5082 r5093  
    3030c====================================================================
    3131c
    32 c   Sous-programme consacre � des diagnostics dynamiques de base
     32c   Sous-programme consacre à des diagnostics dynamiques de base
    3333c
    3434c
    … …  
    8888      real ww
    8989
    90 c   variables dynamiques interm�diaires
     90c   variables dynamiques intermédiaires
    9191      REAL,SAVE,ALLOCATABLE :: vcont(:,:,:),ucont(:,:,:)
    9292      REAL,SAVE,ALLOCATABLE :: ang(:,:,:),unat(:,:,:)
    … …  
    9696      REAL,SAVE,ALLOCATABLE :: bern(:,:,:)
    9797
    98 c   champ contenant les scalaires advect�s.
     98c   champ contenant les scalaires advectés.
    9999      real,SAVE,ALLOCATABLE :: Q(:,:,:,:)
    100100   
    101 c   champs cumul�s
     101c   champs cumulés
    102102      real,SAVE,ALLOCATABLE ::  ps_cum(:,:)
    103103      real,SAVE,ALLOCATABLE ::  masse_cum(:,:,:)
    … …  
    370370      jje=jj_end
    371371   
    372 c   �nergie cin�tique
     372c   énergie cinétique
    373373!      ucont(:,jjb:jje,:)=0
    374374
    … …  
    382382      CALL enercin_loc(vcov,ucov,vcont,ucont,ecin)
    383383
    384 c   moment cin�tique
     384c   moment cinétique
    385385!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    386386      do l=1,llm
    … …  
    490490      enddo
    491491
    492 c    flux m�ridien
     492c    flux méridien
    493493c    -------------
    494494      do iQ=1,nQ
    … …  
    642642
    643643c=====================================================================
    644 c   Transport m�ridien
     644c   Transport méridien
    645645c=====================================================================
    646646

  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/conf_gcm.F90

    r5082 r5093  
    121121
    122122  !Config  Key  = prt_level
    123   !Config  Desc = niveau d'impressions de d�bogage
     123  !Config  Desc = niveau d'impressions de débogage
    124124  !Config  Def  = 0
    125   !Config  Help = Niveau d'impression pour le d�bogage
     125  !Config  Help = Niveau d'impression pour le débogage
    126126  !Config         (0 = minimum d'impression)
    127127  prt_level = 0
    … …  
    921921
    922922     !Config  Key  = use_mpi_alloc
    923      !Config  Desc = Utilise un buffer MPI en m�moire globale
     923     !Config  Desc = Utilise un buffer MPI en mémoire globale
    924924     !Config  Def  = false
    925925     !Config  Help = permet d'activer l'utilisation d'un buffer MPI
    926      !Config         en m�moire globale a l'aide de la fonction MPI_ALLOC.
    927      !Config         Cela peut am�liorer la bande passante des transferts MPI
     926     !Config         en mémoire globale a l'aide de la fonction MPI_ALLOC.
     927     !Config         Cela peut améliorer la bande passante des transferts MPI
    928928     !Config         d'un facteur 2 
    929929     use_mpi_alloc=.FALSE.
    … …  
    933933     !Config  Desc = activation de la version strato
    934934     !Config  Def  = .FALSE.
    935      !Config  Help = active la version stratosph�rique de LMDZ de F. Lott
     935     !Config  Help = active la version stratosphérique de LMDZ de F. Lott
    936936
    937937     ok_strato=.FALSE.

  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/guide_loc_mod.F90

    r5088 r5093  
    100100    CALL getpar('guide_teta',.false.,guide_teta,'guidage de T par Teta')
    101101
    102     CALL getpar('guide_add',.false.,guide_add,'for�age constant?')
     102    CALL getpar('guide_add',.false.,guide_add,'foréage constant?')
    103103    CALL getpar('guide_zon',.false.,guide_zon,'guidage moy zonale')
    104104    if (guide_zon .and. abs(grossismx - 1.) > 0.01) &
    … …  
    121121    CALL getpar('plim_guide_BL',85000.,plim_guide_BL,'BL top presnivs value')
    122122
    123 ! Sauvegarde du for�age
     123! Sauvegarde du forçage
    124124    CALL getpar('guide_sav',.false.,guide_sav,'sauvegarde guidage')
    125125    CALL getpar('iguide_sav',4,iguide_sav,'freq. sauvegarde guidage')
    … …  
    22692269        call nf95_put_var(nid, varid_alpha_q, zq)
    22702270! --------------------------------------------------------------------
    2271 ! Cr�ation des variables sauvegard�es
     2271! Création des variables sauvegardées
    22722272! --------------------------------------------------------------------
    22732273        ierr = nf90_redef(nid)

  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/parallel_lmdz.F90

    r5089 r5093  
    222222!Config  Desc = taille des blocs openmp
    223223!Config  Def  = 1
    224 !Config  Help = defini la taille des packets d'it�ration openmp
     224!Config  Help = defini la taille des packets d'itération openmp
    225225!Config         distribue a chaque tache lors de l'entree dans une
    226226!Config         boucle parallelisee

  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/times.F90

    r4600 r5093  
    194194      endif
    195195     
    196     ENDIF  ! using_mp�
     196    ENDIF  ! using_mpi
    197197  end subroutine allgather_timer_average
    198198 

  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/top_bound_loc.F

    r5082 r5093  
    2727c   ------
    2828c
    29 c   Dissipation lin�aire (ex top_bound de la physique)
     29c   Dissipation linéaire (ex top_bound de la physique)
    3030c
    3131c=======================================================================

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_calfis_loc.F90

    r5090 r5093  
    187187  ! Introduction du splitting (FH)
    188188  ! Question pour Yann :
    189   ! J'ai �t� surpris au d�but que les tableaux zufi_omp, zdufi_omp n'co soitent
     189  ! J'ai été surpris au début que les tableaux zufi_omp, zdufi_omp n'co soitent
    190190  ! en SAVE. Je crois comprendre que c'est parce que tu voulais qu'il
    191191  ! soit allocatable (plutot par exemple que de passer une dimension
    192   ! d�pendant du process en argument des routines) et que, du coup,
    193   ! le SAVE �vite d'avoir � refaire l'allocation � chaque appel.
     192  ! dépendant du process en argument des routines) et que, du coup,
     193  ! le SAVE évite d'avoir à refaire l'allocation à chaque appel.
    194194  ! Tu confirmes ?
    195   ! J'ai suivi le m�me principe pour les zdufic_omp
     195  ! J'ai suivi le même principe pour les zdufic_omp
    196196  ! Mais c'est surement bien que tu controles.
    197197  !

  • LMDZ6/branches/Amaury_dev/libf/misc/iniprint.h

    r2352 r5093  
    44!
    55! gestion des impressions de sorties et de débogage
    6 ! lunout:    unité du fichier dans lequel se font les sorties 
     6! lunout:    unité du fichier dans lequel se font les sorties
    77!                           (par defaut 6, la sortie standard)
    88! prt_level: niveau d'impression souhaité (0 = minimum)

  • LMDZ6/branches/Amaury_dev/libf/misc/write_field.F90

    r5090 r5093  
    282282                      //trim(int2str(pos+offset))      &   
    283283                      //'," ---> ",g22.16," | ")'
    284 ! d�pent de l'impl�mention, sur compaq, c'est necessaire
     284! dépend de l'implémention, sur compaq, c'est necessaire
    285285!            Pos=Pos+ColumnSize
    286286          endif

  • LMDZ6/branches/Amaury_dev/libf/phy_common/mod_phys_lmdz_mpi_transfert.F90

    r5081 r5093  
    6565!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    6666
    67 !! -- Les chaine de charact�re -- !!
     67!! -- Les chaine de charactère -- !!
    6868
    6969  SUBROUTINE bcast_mpi_c(var1)

  • LMDZ6/branches/Amaury_dev/libf/phy_common/mod_phys_lmdz_omp_transfert.F90

    r5058 r5093  
    109109!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    110110
    111 !! -- Les chaine de charact�re -- !!
     111!! -- Les chaine de charactère -- !!
    112112
    113113  SUBROUTINE bcast_omp_c(var)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/add_wake_tend.F90

    r5081 r5093  
    11SUBROUTINE add_wake_tend(zddeltat, zddeltaq, zds, zdas, zddensw, zddensaw, zoccur, text, abortphy)
    22!===================================================================
    3 ! Ajoute les tendances li�es aux diverses parametrisations physiques aux
     3! Ajoute les tendances liées aux diverses parametrisations physiques aux
    44! variables d'etat des poches froides.
    55!===================================================================

  • LMDZ6/branches/Amaury_dev/libf/phylmd/calcul_fluxs_mod.F90

    r5082 r5093  
    5555!   dflux_s      derivee du flux de chaleur sensible / Ts
    5656!   dflux_l      derivee du flux de chaleur latente  / Ts
    57 !   sens_prec_liq flux sensible li� aux echanges de precipitations liquides
     57!   sens_prec_liq flux sensible lié aux echanges de precipitations liquides
    5858!   sens_prec_sol                                    precipitations solides
    59 !   lat_prec_liq  flux latent li� aux echanges de precipitations liquides
     59!   lat_prec_liq  flux latent lié aux echanges de precipitations liquides
    6060!   lat_prec_sol                                  precipitations solides
    6161

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/cosp_output_write_mod.F90

    r5082 r5093  
    457457
    458458    ! ug On récupère le type écrit de la structure:
    459     !       Assez moche, �| refaire si meilleure méthode...
     459    !       Assez moche, à refaire si meilleure méthode...
    460460    IF (INDEX(var%cosp_typeecrit(iff), "once") > 0) THEN
    461461       typeecrit = 'once'
    … …  
    559559
    560560    ! ug On récupère le type écrit de la structure:
    561     !       Assez moche, �| refaire si meilleure méthode...
     561    !       Assez moche, à refaire si meilleure méthode...
    562562    IF (INDEX(var%cosp_typeecrit(iff), "once") > 0) THEN
    563563       typeecrit = 'once'

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/lidar_simulator.F90

    r5087 r5093  
    690690! Lum_norm=f(tetaS,tau_cloud) derived from adding-doubling calculations
    691691!        valid ONLY ABOVE OCEAN (albedo_sfce=5%)
    692 !        valid only in one viewing direction (theta_v=30�, phi_s-phi_v=320�)
     692!        valid only in one viewing direction (theta_v=30°, phi_s-phi_v=320°)
    693693!        based on adding-doubling radiative transfer computation
    694694!        for tau values (0 to 100) and for tetas values (0 to 80)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/cosp_output_write_mod.F90

    r5082 r5093  
    474474
    475475    ! ug On récupère le type écrit de la structure:
    476     !       Assez moche, �| refaire si meilleure méthode...
     476    !       Assez moche, à refaire si meilleure méthode...
    477477    IF (INDEX(var%cosp_typeecrit(iff), "once") > 0) THEN
    478478       typeecrit = 'once'
    … …  
    576576
    577577    ! ug On récupère le type écrit de la structure:
    578     !       Assez moche, �| refaire si meilleure méthode...
     578    !       Assez moche, à refaire si meilleure méthode...
    579579    IF (INDEX(var%cosp_typeecrit(iff), "once") > 0) THEN
    580580       typeecrit = 'once'

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/parasol.F90

    r5082 r5093  
    8181    ! Lum_norm=f(PARASOL_SZA,tau_cloud) derived from adding-doubling calculations
    8282    !        valid ONLY ABOVE OCEAN (albedo_sfce=5%)
    83     !        valid only in one viewing direction (theta_v=30�, phi_s-phi_v=320�)
     83    !        valid only in one viewing direction (theta_v=30°, phi_s-phi_v=320°)
    8484    !        based on adding-doubling radiative transfer computation
    8585    !        for PARASOL_TAU values (0 to 100) and for PARASOL_SZA values (0 to 80)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/lmdz_cosp_output_write_mod.F90

    r5082 r5093  
    670670
    671671    ! ug On récupère le type écrit de la structure:
    672     !       Assez moche, �| refaire si meilleure méthode...
     672    !       Assez moche, à refaire si meilleure méthode...
    673673    IF (INDEX(var%cosp_typeecrit(iff), "once") > 0) THEN
    674674       typeecrit = 'once'
    … …  
    766766
    767767    ! ug On récupère le type écrit de la structure:
    768     !       Assez moche, �| refaire si meilleure méthode...
     768    !       Assez moche, à refaire si meilleure méthode...
    769769    IF (INDEX(var%cosp_typeecrit(iff), "once") > 0) THEN
    770770       typeecrit = 'once'

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/parasol.F90

    r5082 r5093  
    8181    ! Lum_norm=f(PARASOL_SZA,tau_cloud) derived from adding-doubling calculations
    8282    !        valid ONLY ABOVE OCEAN (albedo_sfce=5%)
    83     !        valid only in one viewing direction (theta_v=30�, phi_s-phi_v=320�)
     83    !        valid only in one viewing direction (theta_v=30°, phi_s-phi_v=320°)
    8484    !        based on adding-doubling radiative transfer computation
    8585    !        for PARASOL_TAU values (0 to 100) and for PARASOL_SZA values (0 to 80)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cv30_routines.F90

    r5086 r5093  
    30823082
    30833083  ! fraction deau condensee dans les melanges convertie en precip : epm
    3084   ! et eau condens�e pr�cipit�e dans masse d'air satur� : l_m*dM_m/dzdz.dzdz
     3084  ! et eau condensée précipitée dans masse d'air saturé : l_m*dM_m/dzdz.dzdz
    30853085  DO j = 1, nam1
    30863086    DO k = 1, j - 1
    … …  
    32773277
    32783278        ! On fait varier epmax en fn de la cape
    3279         ! Il faut donc recalculer ep, et hp qui a d�j� �t� calcul� et
    3280         ! qui en d�pend
    3281         ! Toutes les autres variables fn de ep sont calcul�es plus bas.
     3279        ! Il faut donc recalculer ep, et hp qui a déjà été calculé et
     3280        ! qui en dépend
     3281        ! Toutes les autres variables fn de ep sont calculées plus bas.
    32823282
    32833283        INCLUDE "cvthermo.h"

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cv3_routines.F90

    r5086 r5093  
    11011101
    11021102!JAM--------------------------------------------------------------------
    1103 ! Calcul de la quantit� d'eau sous forme de glace
     1103! Calcul de la quantité d'eau sous forme de glace
    11041104! --------------------------------------------------------------------
    11051105  INTEGER nl, len
    … …  
    30723072!!---end jyg---
    30733073
    3074 ! --------retour � la formulation originale d'Emanuel.
     3074! --------retour à la formulation originale d'Emanuel.
    30753075        IF (cvflag_ice) THEN
    30763076
    … …  
    30823082
    30833083!JAM  Attention: evap=sigt*E
    3084 !    Modification: evap devient l'�vaporation en milieu de couche
    3085 !    car n�cessaire dans cv3_yield
    3086 !    Du coup, il faut modifier pas mal d'�quations...
     3084!    Modification: evap devient l'évaporation en milieu de couche
     3085!    car nécessaire dans cv3_yield
     3086!    Du coup, il faut modifier pas mal d'équations...
    30873087!    et l'expression de afac qui devient afac1
    30883088!    revap=sqrt((prec(i+1)+prec(i))/2)
    … …  
    31033103!JYG    Dans sa formulation originale, Emanuel calcule l'evaporation par:
    31043104! c             evap(il,i)=sigt*afac*revap
    3105 ! ce qui n'est pas correct. Dans cv_routines, la formulation a �t� modifiee.
     3105! ce qui n'est pas correct. Dans cv_routines, la formulation a été modifiee.
    31063106! Ici,l'evaporation evap est simplement calculee par l'equation de
    31073107! conservation.
    … …  
    49024902    ENDDO       ! k
    49034903
    4904 ! 14/01/15 AJ delta n'a rien � faire l�...                                                 
     4904! 14/01/15 AJ delta n'a rien à faire là...
    49054905    DO il = 1, ncum                                                  ! cld
    49064906!!      IF (wa(il,i)>0.0 .AND. iflag(il)<=1) &                         ! cld
    … …  
    49184918
    49194919! IM cf. FH
    4920 ! 14/01/15 AJ ne correspond pas � ce qui a �t� cod� par JYG et SB           
     4920! 14/01/15 AJ ne correspond pas à ce qui a été codé par JYG et SB
    49214921                                                         
    49224922      IF (iflag_clw==0) THEN                                         ! cld
    … …  
    49914991
    49924992! fraction deau condensee dans les melanges convertie en precip : epm
    4993 ! et eau condens�e pr�cipit�e dans masse d'air satur� : l_m*dM_m/dzdz.dzdz
     4993! et eau condensée précipitée dans masse d'air saturé : l_m*dM_m/dzdz.dzdz
    49944994  DO j = 1, nl
    49954995    DO k = 1, nl
    … …  
    51415141
    51425142        ! On fait varier epmax en fn de la cape
    5143         ! Il faut donc recalculer ep, et hp qui a d�j� �t� calcul� et
    5144         ! qui en d�pend
    5145         ! Toutes les autres variables fn de ep sont calcul�es plus bas.
     5143        ! Il faut donc recalculer ep, et hp qui a déjà été calculé et
     5144        ! qui en dépend
     5145        ! Toutes les autres variables fn de ep sont calculées plus bas.
    51465146
    51475147  include "cvthermo.h"
    … …  
    51785178
    51795179        ! il faut calculer la cape: on fait un calcule simple car tant qu'on ne
    5180         ! connait pas ep, on ne connait pas les m�langes, ddfts etc... qui sont
     5180        ! connait pas ep, on ne connait pas les mélanges, ddfts etc... qui sont
    51815181        ! necessaires au calcul de la cape dans la nouvelle physique
    51825182       

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cv_driver.F90

    r3492 r5093  
    568568                ,cape,ep,hp,icb,inb,clw,nk,t,h,lv &
    569569                ,epmax_diag)
    570         ! on écrase ep et recalcule hp
     570        ! on écrase ep et recalcule hp
    571571    END IF
    572572
    … …  
    699699! si icvflag_Tpa=0, alors la fraction de glace dans l'ascendance adiabatique est
    700700  ! fonction de la temperature de l'environnement et la temperature de l'ascendance est
    701   ! calculee en deux itérations, une en supposant qu'il n'y a pas de glace et l'autre 
     701  ! calculee en deux itérations, une en supposant qu'il n'y a pas de glace et l'autre
    702702  ! en ajoutant la glace (ancien schéma d'Arnaud Jam).
    703703! si icvflag_Tpa=1, alors la fraction de glace dans l'ascendance adiabatique est

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cvltr.F90

    r5082 r5093  
    4949  REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: sij        ! fraction dair de lenv
    5050  REAL,DIMENSION(klon,klev),INTENT(IN)      :: wght_cvfd  ! weights of the layers feeding convection
    51   REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: elij       ! contenu en eau condens�e sp�cifique/conc deau condens�e massique
     51  REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: elij       ! contenu en eau condensée spécifique/conc deau condensée massique
    5252  REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: epmlmMm    ! eau condensee precipitee dans mel masse dair sat
    5353  REAL,DIMENSION(klon,klev),INTENT(IN)      :: eplaMm    ! eau condensee precipitee dans aa masse dair sat
    5454
    55   REAL,DIMENSION(klon,klev),INTENT(IN)      :: clw        ! contenu en eau condens�e dans lasc adiab
     55  REAL,DIMENSION(klon,klev),INTENT(IN)      :: clw        ! contenu en eau condensée dans lasc adiab
    5656  REAL,DIMENSION(klon),INTENT(IN)           :: sigd
    5757  INTEGER,DIMENSION(klon),INTENT(IN)        :: icb,inb
    … …  
    8080
    8181! RomP ! les variables sont nettoyees des valeurs aberrantes
    82   REAL,DIMENSION(klon,klev)         :: Pa, Pm  ! pluie AA et m�langes, var temporaire
     82  REAL,DIMENSION(klon,klev)         :: Pa, Pm  ! pluie AA et mélanges, var temporaire
    8383  REAL,DIMENSION(klon,klev)         :: pmflxs,pmflxr ! pmflxrIN,pmflxsIN sans valeur aberante
    8484  REAL,DIMENSION(klon,klev)         :: mp            ! flux de masse
    … …  
    155155! On prend la moyenne des precip entre le niveau i+1 et i
    156156! I=3/4* (P(1+1)+P(i))/2 / (sigd*r*rho_l)
    157 !  1000kg/m3= densit� de l'eau
     157!  1000kg/m3= densité de l'eau
    158158! 0.75e-3 = 3/4 /1000
    159 ! Par la suite, I est tout le temps multipli� par sig_d pour avoir l'impaction sur la surface de la maille
    160 ! on le n�glige ici pour simplifier le code
     159! Par la suite, I est tout le temps multiplié par sig_d pour avoir l'impaction sur la surface de la maille
     160! on le néglige ici pour simplifier le code
    161161  do j=1,klev-1
    162162   do i=1,klon
    … …  
    201201  END DO
    202202
    203 ! suppression des valeurs tr�s faibles (~1e-320)
     203! suppression des valeurs très faibles (~1e-320)
    204204! multiplication de levaporation pour lavoir par unite de temps
    205205! et par unite de surface de la maille
    … …  
    468468     kappa(i,j)=1.
    469469     if(j==1) then
    470        qDi(i,j,it)=qDi(i,j+1,it) !orig tr(i,j,it)   ! mp(1)=0 donc tout vient de la couche sup�rieure
     470       qDi(i,j,it)=qDi(i,j+1,it) !orig tr(i,j,it)   ! mp(1)=0 donc tout vient de la couche supérieure
    471471     elseif(mp(i,j+1)>mp(i,j).and.mp(i,j+1)>1.e-10) then
    472472       qDi(i,j,it)=qDi(i,j+1,it)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/mod_1D_amma_read.F90

    r5088 r5093  
    389389       endif
    390390       if (annee_ref==2006 .and. day1<day_ini_amma) then
    391         print*,'AMMA a d�but� le 10 juillet 2006',day1,day_ini_amma
     391        print*,'AMMA a débuté le 10 juillet 2006',day1,day_ini_amma
    392392        print*,'Changer dayref dans run.def'
    393393        stop

  • LMDZ6/branches/Amaury_dev/libf/phylmd/hbtm_mod.F90

    r3774 r5093  
    181181    ! ++ r4 = 35.86             7.66           Kelvin
    182182    ! ++  q_sat = eps*e_sat/(p-(1-eps)*e_sat)
    183     ! ++ deriv� :
     183    ! ++ derivé :
    184184    ! ++ =========
    185185    ! ++                   r3*(Tf-r4)*q_sat(T,p)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/init_be.F90

    r5082 r5093  
    6262  WRITE(*,*)'PASSAGE init_be ...'
    6363
    64 ! la source est maintenant d�finie independemment de la valeur de klev.
     64! la source est maintenant définie independemment de la valeur de klev.
    6565!!! Source actuellement definie pour klev = 19 et klev >= 39
    6666!!  IF (klev /= 19 .AND. klev<39) CALL abort_physic("init_be","Source du be7 necessite klev=19 ou klev>=39",1)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_tso.F

    r3900 r5093  
    623623c #EU      TsisSV(ikl,isl)  = max(TaT_SV(ikl)-15.,TsisSV(ikl,isl))
    624624
    625           !XF 18/11/2018 to avoid ST reaching 70�C!!
     625          !XF 18/11/2018 to avoid ST reaching 70°C!!
    626626          !It is an error compensation but does not work over tundra
    627627 

  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_alp.F90

    r5087 r5093  
    2929!   Nicolas Rochetin et Jean-Yves Grandpeix
    3030!         pour la fermeture stochastique. 2012
    31 !   Fr�d�ric Hourdin :
     31!   Frédéric Hourdin :
    3232!         netoyage informatique. 2022
    3333!   
    … …  
    8181      INTEGER ig,k,l
    8282      integer nsrf
    83       real rhobarz0(ngrid)                    ! Densit� au LCL
     83      real rhobarz0(ngrid)                    ! Densité au LCL
    8484      logical ok_lcl(ngrid)                   ! Existence du LCL des thermiques
    8585      integer klcl(ngrid)                     ! Niveau du LCL
    8686      real interp(ngrid)                      ! Coef d'interpolation pour le LCL
    8787!--Triggering
    88       real, parameter :: su_cst=4e4              ! Surface unite: celle d'un updraft �l�mentaire
     88      real, parameter :: su_cst=4e4              ! Surface unite: celle d'un updraft élémentaire
    8989      real, parameter :: hcoef=1             ! Coefficient directeur pour le calcul de s2
    90       real, parameter :: hmincoef=0.3        ! Coefficient directeur pour l'ordonn�e � l'origine pour le calcul de s2
    91       real, parameter :: eps1=0.3            ! Fraction de surface occup�e par la population 1 : eps1=n1*s1/(fraca0*Sd)
    92       real, dimension(ngrid) :: hmin         ! Ordonn�e � l'origine pour le calcul de s2
     90      real, parameter :: hmincoef=0.3        ! Coefficient directeur pour l'ordonnée à l'origine pour le calcul de s2
     91      real, parameter :: eps1=0.3            ! Fraction de surface occupée par la population 1 : eps1=n1*s1/(fraca0*Sd)
     92      real, dimension(ngrid) :: hmin         ! Ordonnée à l'origine pour le calcul de s2
    9393      real, dimension(ngrid) :: zmax_moy     ! Hauteur moyenne des thermiques : zmax_moy = zlcl + 0.33 (zmax-zlcl)
    9494      real, parameter :: zmax_moy_coef=0.33
    … …  
    100100      real, dimension(ngrid,nlay) :: pbl_tke_max       ! Profil de TKE moyenne
    101101      real, dimension(ngrid) :: pbl_tke_max0           ! TKE moyenne au LCL
    102       real, dimension(ngrid,nlay) :: w_ls              ! Vitesse verticale grande �chelle (m/s)
    103       real, parameter :: coef_m=1.            ! On consid�re un rendement pour alp_bl_fluct_m
    104       real, parameter :: coef_tke=1.          ! On consid�re un rendement pour alp_bl_fluct_tke
     102      real, dimension(ngrid,nlay) :: w_ls              ! Vitesse verticale grande échelle (m/s)
     103      real, parameter :: coef_m=1.            ! On considère un rendement pour alp_bl_fluct_m
     104      real, parameter :: coef_tke=1.          ! On considère un rendement pour alp_bl_fluct_tke
    105105      real :: zdp
    106106      real, dimension(ngrid) :: alp_int,dp_int
    … …  
    171171!!jyg fin
    172172
    173 !------------Calcul des propri�t�s du thermique au LCL
     173!------------Calcul des propriétés du thermique au LCL
    174174  IF ( (iflag_trig_bl>=1) .or. (iflag_clos_bl>=1) ) THEN
    175175
    … …  
    198198   enddo
    199199
    200 !-----Calcul de la TKE transport�e par les thermiques : therm_tke_max
     200!-----Calcul de la TKE transportée par les thermiques : therm_tke_max
    201201   call thermcell_tke_transport(ngrid,nlay,ptimestep,fm0,entr0,  &  ! in
    202202             rg,pplev,therm_tke_max)                               ! out
    203203!   print *,' thermcell_tke_transport -> '   !!jyg
    204204
    205 !-----Calcul des profils verticaux de TKE hors thermiques : env_tke_max, et de la vitesse verticale grande �chelle : W_ls
     205!-----Calcul des profils verticaux de TKE hors thermiques : env_tke_max, et de la vitesse verticale grande échelle : W_ls
    206206   do l=1,nlay
    207207    do ig=1,ngrid
    208      pbl_tke_max(ig,l)=fraca(ig,l)*therm_tke_max(ig,l)+(1.-fraca(ig,l))*env_tke_max(ig,l)         !  Recalcul de TKE moyenne apr�s transport de TKE_TH
    209      env_tke_max(ig,l)=(pbl_tke_max(ig,l)-fraca(ig,l)*therm_tke_max(ig,l))/(1.-fraca(ig,l))       !  Recalcul de TKE dans  l'environnement apr�s transport de TKE_TH
    210      w_ls(ig,l)=-1.*omega(ig,l)/(RG*rhobarz(ig,l))                                                !  Vitesse verticale de grande �chelle
     208     pbl_tke_max(ig,l)=fraca(ig,l)*therm_tke_max(ig,l)+(1.-fraca(ig,l))*env_tke_max(ig,l)         !  Recalcul de TKE moyenne après transport de TKE_TH
     209     env_tke_max(ig,l)=(pbl_tke_max(ig,l)-fraca(ig,l)*therm_tke_max(ig,l))/(1.-fraca(ig,l))       !  Recalcul de TKE dans  l'environnement après transport de TKE_TH
     210     w_ls(ig,l)=-1.*omega(ig,l)/(RG*rhobarz(ig,l))                                                !  Vitesse verticale de grande échelle
    211211    enddo
    212212   enddo
    … …  
    253253
    254254
    255 !-----Epaisseur du nuage (depth) et d�termination de la queue du spectre de panaches (n2,s2) et du panache le plus gros (s_max)
     255!-----Epaisseur du nuage (depth) et détermination de la queue du spectre de panaches (n2,s2) et du panache le plus gros (s_max)
    256256   do ig=1,ngrid
    257257     zmax_moy(ig)=zlcl(ig)+zmax_moy_coef*(zmax(ig)-zlcl(ig))
    … …  
    328328  enddo
    329329
    330 !-----S�curit� ALP infinie
     330!-----Sécurité ALP infinie
    331331  do ig=1,ngrid
    332332   if (fraca0(ig)>0.98) alp_bl_stat(ig)=2.

  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_dry.F90

    r5087 r5093  
    9090! 1. le flux de masse en haut  f_star(ig,l+1)
    9191! 2. la temperature potentielle virtuelle dans la couche ztva(ig,l)
    92 ! 3. la vitesse au carr� en haut zw2(ig,l+1)
     92! 3. la vitesse au carré en haut zw2(ig,l+1)
    9393!------------------------------------------------------------------------
    9494

  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_dv2.F90

    r5087 r5093  
    113113               zf2=1.
    114114
    115 !  la premi�re fois on multiplie le coefficient de freinage
     115!  la première fois on multiplie le coefficient de freinage
    116116!  par le module du vent dans la couche en dessous.
    117117!  Mais pourquoi donc ???

  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_flux2.F90

    r5087 r5093  
    174174! FH Version en cours de test;
    175175! par rapport a thermcell_flux, on fait une grande boucle sur "l"
    176 ! et on modifie le flux avec tous les contr�les appliques d'affilee
     176! et on modifie le flux avec tous les contrôles appliques d'affilee
    177177! pour la meme couche
    178178! Momentanement, on duplique le calcule du flux pour pouvoir comparer

  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_height.F90

    r5087 r5093  
    3939      enddo
    4040
    41 ! On traite le cas particulier qu'il faudrait �viter ou le thermique
     41! On traite le cas particulier qu'il faudrait éviter ou le thermique
    4242! atteind le haut du modele ...
    4343      do ig=1,ngrid

  • LMDZ6/branches/Amaury_dev/libf/phylmd/o3_chem_m.F90

    r5086 r5093  
    2424
    2525    integer, intent(in):: julien ! jour julien, 1 <= julien <= 360
    26     real, intent(in):: gmtime ! heure de la journ�e en fraction de jour
     26    real, intent(in):: gmtime ! heure de la journée en fraction de jour
    2727    real, intent(in):: t_seri(:, :) ! (klon, nbp_lev) temperature, in K
    2828
    … …  
    5959    real earth_long
    6060    ! (longitude vraie de la Terre dans son orbite solaire, par
    61     ! rapport au point vernal (21 mars), en degr�s)
     61    ! rapport au point vernal (21 mars), en degrés)
    6262
    6363    real pmu0(klon) ! mean of cosine of solar zenith angle during "pdtphys"

  • LMDZ6/branches/Amaury_dev/libf/phylmd/phyaqua_mod.F90

    r5090 r5093  
    942942      if (type_profil==20) then
    943943      print*,'Profile SST 20'
    944 !     Méthode 13 "Qmax2K" plateau réel �|  l'Equateur augmenté +2K
     944!     Méthode 13 "Qmax2K" plateau réel é|  l'Equateur augmenté +2K
    945945
    946946      do j=1,klon
    … …  
    951951      if (type_profil==21) then
    952952      print*,'Profile SST 21'
    953 !     Méthode 13 "Qmax2K" plateau réel �|  l'Equateur augmenté +2K
     953!     Méthode 13 "Qmax2K" plateau réel é|  l'Equateur augmenté +2K
    954954      do j=1,klon
    955955        phy_sst(j,i)=252.+55.*(1-sin(rlatd(j))**2)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/phys_output_var_mod.F90

    r5087 r5093  
    164164 
    165165  REAL, SAVE, ALLOCATABLE:: tkt(:) ! (klon)
    166   ! �paisseur (m) de la couche de diffusion thermique (microlayer)
     166  ! épaisseur (m) de la couche de diffusion thermique (microlayer)
    167167  ! cool skin thickness
    168168
    169169  REAL, SAVE, ALLOCATABLE:: tks(:) ! (klon)
    170   ! �paisseur (m) de la couche de diffusion de masse (microlayer)
     170  ! épaisseur (m) de la couche de diffusion de masse (microlayer)
    171171 
    172172  REAL, SAVE, ALLOCATABLE:: taur(:) ! (klon) momentum flux due to rain, in Pa

  • LMDZ6/branches/Amaury_dev/libf/phylmd/plevel.F90

    r5088 r5093  
    7474  ! =====================================================================
    7575  IF (lnew) THEN
    76     ! on r�nitialise les r�ndicages et les poids
     76    ! on rénitialise les réndicages et les poids
    7777    ! =====================================================================
    7878
    … …  
    108108
    109109      ! ...   Modif . P. Le Van    ( 20/01/98) ....
    110       ! Modif Fr��ic Hourdin (3/01/02)
     110      ! Modif Frédéric Hourdin (3/01/02)
    111111
    112112      aist(i) = log(pgcm(i,lb(i))/pres)/log(pgcm(i,lb(i))/pgcm(i,lt(i)))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/qsatseaw_1D.F90

    r4722 r5093  
    4848!!      Zeng, X., Zhao, M., and Dickinson, R. E., 1998 : Intercomparaison of bulk
    4949!!      aerodynamic algorithm for the computation of sea surface fluxes using
    50 !!      TOGA COARE and TAO data. Journal of Climate, vol 11, n�10, pp 2628--2644
     50!!      TOGA COARE and TAO data. Journal of Climate, vol 11, n°10, pp 2628--2644
    5151!!
    5252!!

  • LMDZ6/branches/Amaury_dev/libf/phylmd/slab_heat_transp_mod.F90

    r5082 r5093  
    409409      ELSE
    410410          alpha(:)=alpham
    411           ! Tsurf-Tdeep ~ 10� in the Tropics
     411          ! Tsurf-Tdeep ~ 10° in the Tropics
    412412      ENDIF
    413413

  • LMDZ6/branches/Amaury_dev/libf/phylmd/surf_land_orchidee_mod.F90

    r5082 r5093  
    254254       IF ( .NOT. ALLOCATED(albedo_keep)) THEN
    255255!ym          ALLOCATE(albedo_keep(klon))
    256 !ym bizarre que non allou� en knon precedement
     256!ym bizarre que non alloué en knon precedement
    257257          ALLOCATE(albedo_keep(knon))
    258258          ALLOCATE(zlev(knon))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/surf_ocean_mod.F90

    r5087 r5093  
    142142
    143143    REAL, intent(out):: tkt(:) ! (knon)
    144     ! �paisseur (m) de la couche de diffusion thermique (microlayer)
     144    ! épaisseur (m) de la couche de diffusion thermique (microlayer)
    145145    ! cool skin thickness
    146146
    147147    REAL, intent(out):: tks(:) ! (knon)
    148     ! �paisseur (m) de la couche de diffusion de masse (microlayer)
     148    ! épaisseur (m) de la couche de diffusion de masse (microlayer)
    149149
    150150    REAL, intent(out):: taur(:) ! (knon)

  • LMDZ6/branches/Amaury_dev/libf/phylmdiso/add_wake_tend.F90

    r5082 r5093  
    55        )
    66!===================================================================
    7 ! Ajoute les tendances li�es 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   ! D�claration 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 impos�e.
     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 r�partition
     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               D�finitions 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.