Changeset 1779 for LMDZ5/trunk


Ignore:
Timestamp:
Jul 5, 2013, 10:19:53 AM (11 years ago)
Author:
Laurent Fairhead
Message:

Suite de modifs pour le 1d
MPL

Follow-up on modifications for 1d cases
MPL

Location:
LMDZ5/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/calltherm.F90

    r1776 r1779  
    1717     &      ,alp_bl_conv,alp_bl_stat &
    1818!!! fin nrlmd le 10/04/2012
    19      &                    )
     19     &      ,zqla,ztva )
    2020
    2121      USE dimphy
    … …  
    6363      real zqla(klon,klev)
    6464      real zqta(klon,klev)
    65       real ztv(klon,klev)
     65      real ztv(klon,klev),ztva(klon,klev)
    6666      real zpspsk(klon,klev)
    6767      real ztla(klon,klev)
    … …  
    257257     &      ,alp_bl_conv,alp_bl_stat &
    258258!!! fin nrlmd le 10/04/2012
    259      &                        )
     259     &      ,ztva )
    260260           if (prt_level.gt.10) write(lunout,*)'Apres thermcell_main OK'
    261261         else

  • LMDZ5/trunk/libf/phylmd/physiq.F

    r1764 r1779  
    3535      use netcdf95, only: nf95_close
    3636cIM for NMC files
    37       use netcdf, only: nf90_fill_real
     37c     use netcdf, only: nf90_fill_real
     38      use netcdf
    3839      use mod_phys_lmdz_mpi_data, only: is_mpi_root
    3940      USE aero_mod
    … …  
    5051      USE CFMIP_point_locations
    5152      IMPLICIT none
    52 c======================================================================
    53 c
    54 c Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
    55 c
    56 c Objet: Moniteur general de la physique du modele
    57 cAA      Modifications quant aux traceurs :
    58 cAA                  -  uniformisation des parametrisations ds phytrac
    59 cAA                  -  stockage des moyennes des champs necessaires
    60 cAA                     en mode traceur off-line
    61 c======================================================================
    62 c   CLEFS CPP POUR LES IO
    63 c   =====================
     53!>======================================================================
     54!!
     55!! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
     56!!
     57!! Objet: Moniteur general de la physique du modele
     58!!AA      Modifications quant aux traceurs :
     59!!AA                  -  uniformisation des parametrisations ds phytrac
     60!!AA                  -  stockage des moyennes des champs necessaires
     61!!AA                     en mode traceur off-line
     62!!======================================================================
     63!!   CLEFS CPP POUR LES IO
     64!!   =====================
    6465#define histNMC
    6566c#define histISCCP
    66 c======================================================================
    67 c    modif   ( P. Le Van ,  12/10/98 )
    68 c
    69 c  Arguments:
    70 c
    71 c nlon----input-I-nombre de points horizontaux
    72 c nlev----input-I-nombre de couches verticales, doit etre egale a klev
    73 c debut---input-L-variable logique indiquant le premier passage
    74 c lafin---input-L-variable logique indiquant le dernier passage
    75 c jD_cur       -R-jour courant a l'appel de la physique (jour julien)
    76 c jH_cur       -R-heure courante a l'appel de la physique (jour julien)
    77 c pdtphys-input-R-pas d'integration pour la physique (seconde)
    78 c paprs---input-R-pression pour chaque inter-couche (en Pa)
    79 c pplay---input-R-pression pour le mileu de chaque couche (en Pa)
    80 c pphi----input-R-geopotentiel de chaque couche (g z) (reference sol)
    81 c pphis---input-R-geopotentiel du sol
    82 c presnivs-input_R_pressions approximat. des milieux couches ( en PA)
    83 c u-------input-R-vitesse dans la direction X (de O a E) en m/s
    84 c v-------input-R-vitesse Y (de S a N) en m/s
    85 c t-------input-R-temperature (K)
    86 c qx------input-R-humidite specifique (kg/kg) et d'autres traceurs
    87 c d_t_dyn-input-R-tendance dynamique pour "t" (K/s)
    88 c d_q_dyn-input-R-tendance dynamique pour "q" (kg/kg/s)
    89 c flxmass_w -input-R- flux de masse verticale
    90 c d_u-----output-R-tendance physique de "u" (m/s/s)
    91 c d_v-----output-R-tendance physique de "v" (m/s/s)
    92 c d_t-----output-R-tendance physique de "t" (K/s)
    93 c d_qx----output-R-tendance physique de "qx" (kg/kg/s)
    94 c d_ps----output-R-tendance physique de la pression au sol
    95 cIM
    96 c PVteta--output-R-vorticite potentielle a des thetas constantes
    97 c======================================================================
     67!!======================================================================
     68!!    modif   ( P. Le Van ,  12/10/98 )
     69!!
     70!!  Arguments:
     71!!
     72!! nlon----input-I-nombre de points horizontaux
     73!! nlev----input-I-nombre de couches verticales, doit etre egale a klev
     74!! debut---input-L-variable logique indiquant le premier passage
     75!! lafin---input-L-variable logique indiquant le dernier passage
     76!! jD_cur       -R-jour courant a l'appel de la physique (jour julien)
     77!! jH_cur       -R-heure courante a l'appel de la physique (jour julien)
     78!! pdtphys-input-R-pas d'integration pour la physique (seconde)
     79!! paprs---input-R-pression pour chaque inter-couche (en Pa)
     80!! pplay---input-R-pression pour le mileu de chaque couche (en Pa)
     81!! pphi----input-R-geopotentiel de chaque couche (g z) (reference sol)
     82!! pphis---input-R-geopotentiel du sol
     83!! presnivs-input_R_pressions approximat. des milieux couches ( en PA)
     84!! u-------input-R-vitesse dans la direction X (de O a E) en m/s
     85!! v-------input-R-vitesse Y (de S a N) en m/s
     86!! t-------input-R-temperature (K)
     87!! qx------input-R-humidite specifique (kg/kg) et d'autres traceurs
     88!! d_t_dyn-input-R-tendance dynamique pour "t" (K/s)
     89!! d_q_dyn-input-R-tendance dynamique pour "q" (kg/kg/s)
     90!! flxmass_w -input-R- flux de masse verticale
     91!! d_u-----output-R-tendance physique de "u" (m/s/s)
     92!! d_v-----output-R-tendance physique de "v" (m/s/s)
     93!! d_t-----output-R-tendance physique de "t" (K/s)
     94!! d_qx----output-R-tendance physique de "qx" (kg/kg/s)
     95!! d_ps----output-R-tendance physique de la pression au sol
     96!!IM
     97!! PVteta--output-R-vorticite potentielle a des thetas constantes
     98!!======================================================================
    9899#include "dimensions.h"
    99100      integer jjmp1
    … …  
    218219      REAL u(klon,klev)
    219220      REAL v(klon,klev)
    220       REAL t(klon,klev),theta(klon,klev)
     221      REAL t(klon,klev),theta(klon,klev),thetal(klon,klev)
     222c thetal: ligne suivante a decommenter si vous avez les fichiers     MPL 20130625
     223c fth_fonctions.F90 et parkind1.F90
     224c sinon thetal=theta
     225c     REAL fth_thetae,fth_thetav,fth_thetal
    221226      REAL qx(klon,klev,nqtot)
    222227      REAL flxmass_w(klon,klev)
    … …  
    642647      REAL zw2(klon,klev+1)
    643648      REAL fraca(klon,klev+1)       
    644       REAL ztv(klon,klev) 
     649      REAL ztv(klon,klev),ztva(klon,klev)
    645650      REAL zpspsk(klon,klev)
    646       REAL ztla(klon,klev)
     651      REAL ztla(klon,klev),zqla(klon,klev)
    647652      REAL zthl(klon,klev)
    648653
    … …  
    655660      real w0(klon)                                          ! Vitesse des thermiques au LCL
    656661      real w_conv(klon)                                      ! Vitesse verticale de grande \'echelle au LCL
    657       real tke0(klon,klev+1)                              ! TKE au début du pas de temps
     662      real tke0(klon,klev+1)                                 ! TKE au début du pas de temps
    658663      real therm_tke_max0(klon)                              ! TKE dans les thermiques au LCL
    659664      real env_tke_max0(klon)                                ! TKE dans l'environnement au LCL
    … …  
    10571062     .           prof2d_av = 3, prof3d_av = 4)
    10581063      character*30 nom_fichier
    1059       character*10 varname
     1064      character*40 varname
    10601065      character*40 vartitle
    10611066      character*20 varunits
    … …  
    26682673     s      ,alp_bl_conv,alp_bl_stat
    26692674ccc fin nrlmd le 10/04/2012
    2670      s                )
     2675     s      ,zqla,ztva )
    26712676
    26722677ccc nrlmd le 10/04/2012
    … …  
    34513456         solsw=0.
    34523457         radsol=0.
     3458         swup=0.    ! MPL 27102011 pour les fichiers AMMA_profiles et AMMA_scalars
     3459         swup0=0.
     3460         swdn=0.
     3461         swdn0=0.
     3462         lwup=0.
     3463         lwup0=0.
     3464         lwdn=0.
     3465         lwdn0=0.
    34533466      END IF
    34543467
    … …  
    39633976cJYG/IM theta en fin de pas de temps de physique
    39643977        theta(i,k)=t_seri(i,k)*(100000./pplay(i,k))**(RD/RCPD)
     3978c thetal: 2 lignes suivantes a decommenter si vous avez les fichiers     MPL 20130625
     3979c fth_fonctions.F90 et parkind1.F90
     3980c sinon thetal=theta
     3981c       thetal(i,k)=fth_thetal(pplay(i,k),t_seri(i,k),q_seri(i,k),
     3982c    :         ql_seri(i,k))
     3983        thetal(i,k)=theta(i,k)
    39653984      ENDDO
    39663985      ENDDO

  • LMDZ5/trunk/libf/phylmd/thermcell_main.F90

    r1738 r1779  
    1919     &                  ,alp_bl_conv,alp_bl_stat &
    2020!!! fin nrlmd le 10/04/2012
    21      &                         )
     21     &                  ,ztva  )
    2222
    2323      USE dimphy
Note: See TracChangeset for help on using the changeset viewer.