Ignore:
Timestamp:
Jul 18, 2013, 3:32:27 PM (11 years ago)
Author:
Ehouarn Millour
Message:

Déplacement de nombreuses variables de physiq.F vers phys_local_var_mod.
UG
................................
Moving of numerous vars from physiq.F to phys_local_var_mod.
UG

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r1791 r1797  
    173173      PARAMETER (ok_region=.FALSE.)
    174174c======================================================================
    175       real weak_inversion(klon),dthmin(klon)
    176175      real seuil_inversion
    177176      save seuil_inversion
    … …  
    184183      REAL zz,znum,zden
    185184      REAL wmax_th(klon)
    186       REAL zmax_th(klon)
    187185      REAL tau_overturning_th(klon)
    188186
    … …  
    220218      REAL u(klon,klev)
    221219      REAL v(klon,klev)
    222       REAL t(klon,klev),theta(klon,klev),thetal(klon,klev)
     220      REAL t(klon,klev),thetal(klon,klev)
    223221c thetal: ligne suivante a decommenter si vous avez les fichiers     MPL 20130625
    224222c fth_fonctions.F90 et parkind1.F90
    … …  
    227225      REAL qx(klon,klev,nqtot)
    228226      REAL flxmass_w(klon,klev)
    229       REAL omega(klon,klev) ! vitesse verticale en Pa/s
    230227      REAL d_u(klon,klev)
    231228      REAL d_v(klon,klev)
    … …  
    242239      real clw(klon,klev),elij(klon,klev,klev)
    243240      real epmlmMm(klon,klev,klev),eplaMm(klon,klev)
    244       real wdtrainA(klon,klev),wdtrainM(klon,klev)
    245241! RomP <<<
    246242!IM definition dynamique o_trac dans phys_output_open
    … …  
    281277c
    282278c
    283 c prw: precipitable water
    284       real prw(klon)
    285279
    286280      REAL convliq(klon,klev)  ! eau liquide nuageuse convective
    … …  
    295289c flwp, fiwp = Liquid Water Path & Ice Water Path (kg/m2)
    296290c flwc, fiwc = Liquid Water Content & Ice Water Content (kg/kg)
    297       REAL flwp(klon), fiwp(klon)
    298       REAL flwc(klon,klev), fiwc(klon,klev)
    299291      REAL flwp_c(klon), fiwp_c(klon)
    300292      REAL flwc_c(klon,klev), fiwc_c(klon,klev)
    … …  
    302294      REAL flwc_s(klon,klev), fiwc_s(klon,klev)
    303295
    304       REAL evap_pot(klon,nbsrf)
    305296
    306297cIM ISCCP simulator v3.4
    … …  
    528519c$OMP THREADPRIVATE(itap)
    529520c
    530       real slp(klon) ! sea level pressure
    531 c
    532       REAL fevap(klon,nbsrf)
    533       REAL fluxlat(klon,nbsrf)
    534 c
    535       REAL qsol(klon)
    536521      REAL,save ::  solarlong0
    537522c$OMP THREADPRIVATE(solarlong0)
    … …  
    544529c
    545530      INTEGER igwd,idx(klon),itest(klon)
    546 c
    547       REAL agesno(klon,nbsrf)
    548531c
    549532c      REAL,allocatable,save :: run_off_lic_0(:)
    … …  
    565548
    566549      REAL mip(klon,klev)  ! mass flux shed by the adiab ascent at each level
    567       REAL Vprecip(klon,klev+1)   ! precipitation vertical profile
    568550c
    569551      REAL wape_prescr, fip_prescr
    … …  
    589571c$OMP THREADPRIVATE(wake_s_min_lsp)
    590572
    591       real ale_wake(klon)
    592       real alp_wake(klon)
    593573
    594574      real ok_wk_lsp(klon)
    … …  
    609589      REAL wake_dtPBL(klon,klev)      ! Wake : differential heating (wake - unpertubed) PBL
    610590      REAL wake_dqPBL(klon,klev)      ! Wake : differential moistening (wake - unpertubed) PBL
    611       REAL wake_omg(klon,klev)        ! Wake : velocity difference (wake - unpertubed)
    612591      REAL wake_ddeltat(klon,klev),wake_ddeltaq(klon,klev)
    613592      REAL wake_dp_deltomg(klon,klev) ! Wake : gradient vertical de wake_omg
    … …  
    615594c
    616595cpourquoi y'a pas de save??
    617       REAL wake_h(klon)               ! Wake : hauteur de la poche froide
    618596c
    619597      INTEGER wake_k(klon)            ! Wake sommet
    … …  
    645623cRR:fin declarations poches froides
    646624c=======================================================================================================
    647 
    648       REAL zw2(klon,klev+1)
    649       REAL fraca(klon,klev+1)       
     625       
    650626      REAL ztv(klon,klev),ztva(klon,klev)
    651627      REAL zpspsk(klon,klev)
    … …  
    665641      real env_tke_max0(klon)                                ! TKE dans l'environnement au LCL
    666642
    667 c---Spectre de thermiques de type 2 au LCL
    668       real n2(klon),s2(klon)
    669       real ale_bl_stat(klon)
    670 
    671643c---D\'eclenchement stochastique
    672644      integer :: tau_trig(klon)
    673       real proba_notrig(klon)
    674       real random_notrig(klon)
    675645
    676646c--------Statistical Boundary Layer Closure: ALP_BL--------
    … …  
    680650      real env_tke_max(klon,klev)                            ! Profil de TKE dans l'environnement
    681651
    682 c---Fermeture statistique
    683       real alp_bl_det(klon)                                     ! ALP d\'terministe du thermique unique
    684       real alp_bl_fluct_m(klon)                                 ! ALP li\'ee aux fluctuations de flux de masse sous-nuageux
    685       real alp_bl_fluct_tke(klon)                               ! ALP li\'ee aux fluctuations d'\'energie cin\'etique sous-nuageuse
    686       real alp_bl_conv(klon)                                    ! ALP li\'ee \`a grande \'echelle
    687       real alp_bl_stat(klon)                                    ! ALP totale
    688652
    689653ccc fin nrlmd le 10/04/2012
    … …  
    697661c Variables locales:
    698662c
    699       REAL cdragh(klon) ! drag coefficient pour T and Q
    700       REAL cdragm(klon) ! drag coefficient pour vent
    701663cAA
    702664cAA  Pour phytrac
    703665      REAL u1(klon)             ! vents dans la premiere couche U
    704666      REAL v1(klon)             ! vents dans la premiere couche V
    705 
    706       REAL zxffonte(klon), zxfqcalving(klon),zxfqfonte(klon)
    707667
    708668c@$$      LOGICAL offline           ! Controle du stockage ds "physique"
    … …  
    713673! RomP >>>
    714674      REAL beta_prec_fisrt(klon,klev) ! taux de conv de l'eau cond (fisrt)
    715       REAL beta_prec(klon,klev)       ! taux de conv de l'eau cond (utilise)
    716675! RomP <<<
    717676      INTEGER       :: iii
    … …  
    722681c
    723682cIM 050204 END
    724       REAL evap(klon), devap(klon) ! evaporation et sa derivee
    725       REAL sens(klon), dsens(klon) ! chaleur sensible et sa derivee
    726 
    727       REAL bils(klon) ! bilan de chaleur au sol
    728 
    729       REAL wfbilo(klon,nbsrf) ! bilan d'eau, pour chaque
    730 C                             ! type de sous-surface et pondere par la fraction
    731       REAL wfbils(klon,nbsrf) ! bilan de chaleur au sol, pour chaque
    732 C                             ! type de sous-surface et pondere par la fraction
    733       REAL slab_wfbils(klon)  ! bilan de chaleur au sol pour le cas de slab, sur les points d'ocean
    734 
    735       REAL fder(klon)         
    736       REAL ve(klon) ! integr. verticale du transport meri. de l'energie
    737       REAL vq(klon) ! integr. verticale du transport meri. de l'eau
    738       REAL ue(klon) ! integr. verticale du transport zonal de l'energie
    739       REAL uq(klon) ! integr. verticale du transport zonal de l'eau
    740 c
    741       REAL frugs(klon,nbsrf)
    742       REAL zxrugs(klon) ! longueur de rugosite
     683      REAL devap(klon) ! evaporation et sa derivee
     684      REAL dsens(klon) ! chaleur sensible et sa derivee
     685
    743686c
    744687c Conditions aux limites
    … …  
    804747      REAL diafra(klon,klev)  ! fraction nuageuse
    805748      REAL cldliq(klon,klev)  ! eau liquide nuageuse
    806       REAL cldfra(klon,klev)  ! fraction nuageuse
    807       REAL cldtau(klon,klev)  ! epaisseur optique
    808       REAL cldemi(klon,klev)  ! emissivite infrarouge
    809749c
    810750CXXX PB
    811751      REAL fluxq(klon,klev, nbsrf)   ! flux turbulent d'humidite
    812       REAL fluxt(klon,klev, nbsrf)   ! flux turbulent de chaleur
    813       REAL fluxu(klon,klev, nbsrf)   ! flux turbulent de vitesse u
    814       REAL fluxv(klon,klev, nbsrf)   ! flux turbulent de vitesse v
    815752c
    816753      REAL zxfluxt(klon, klev)
    … …  
    818755      REAL zxfluxu(klon, klev)
    819756      REAL zxfluxv(klon, klev)
    820 CXXX
    821 c
    822       REAL fsollw(klon, nbsrf)   ! bilan flux IR pour chaque sous surface
    823       REAL fsolsw(klon, nbsrf)   ! flux solaire absorb. pour chaque sous surface
     757
    824758c Le rayonnement n'est pas calcule tous les pas, il faut donc
    825759c                      sauvegarder les sorties du rayonnement
    … …  
    834768      REAL conv_q(klon,klev) ! convergence de l'humidite (kg/kg/s)
    835769      REAL conv_t(klon,klev) ! convergence de la temperature(K/s)
    836 c
    837       REAL cldl(klon),cldm(klon),cldh(klon) !nuages bas, moyen et haut
    838       REAL cldt(klon),cldq(klon) !nuage total, eau liquide integree
    839 c
    840       REAL zxtsol(klon), zxqsurf(klon), zxsnow(klon), zxfluxlat(klon)
     770
     771c
     772      REAL zxsnow(klon)
    841773      REAL zxsnow_dummy(klon)
    842774c
    … …  
    857789      REAL t_coup
    858790      PARAMETER (t_coup=234.0)
    859 c
    860       REAL zphi(klon,klev)
     791
    861792cym A voir plus tard !!
    862793cym      REAL zx_relief(iim,jjmp1)
    … …  
    864795c
    865796c Grandeurs de sorties
    866       REAL s_pblh(klon), s_lcl(klon), s_capCL(klon)
    867       REAL s_oliqCL(klon), s_cteiCL(klon), s_pblt(klon)
    868       REAL s_therm(klon), s_trmb1(klon), s_trmb2(klon)
     797      REAL s_capCL(klon)
     798      REAL s_oliqCL(klon), s_cteiCL(klon)
     799      REAL s_trmb1(klon), s_trmb2(klon)
    869800      REAL s_trmb3(klon)
    870801cKE43
    871802c Variables locales pour la convection de K. Emanuel (sb):
    872 c
    873       REAL upwd(klon,klev)      ! saturated updraft mass flux
    874       REAL dnwd(klon,klev)      ! saturated downdraft mass flux
    875       REAL dnwd0(klon,klev)     ! unsaturated downdraft mass flux
     803
    876804      REAL tvp(klon,klev)       ! virtual temp of lifted parcel
    877       REAL plcl(klon)           ! Lifting Condensation Level
    878       REAL plfc(klon)           ! Level of Free Convection
    879       REAL wbeff(klon)          ! saturated updraft velocity at LFC
    880805      CHARACTER*40 capemaxcels  !max(CAPE)
    881806
    … …  
    902827c eva: evaporation de l'eau liquide nuageuse
    903828c vdf: couche limite (Vertical DiFfusion)
    904       REAL rneb(klon,klev)
    905829
    906830! tendance nulles
    … …  
    914838cIM 081204 END
    915839c
    916       REAL pmfu(klon,klev), pmfd(klon,klev)
    917840      REAL pen_u(klon,klev), pen_d(klon,klev)
    918841      REAL pde_u(klon,klev), pde_d(klon,klev)
    919842      INTEGER kcbot(klon), kctop(klon), kdtop(klon)
    920       REAL pmflxr(klon,klev+1), pmflxs(klon,klev+1)
    921       REAL prfl(klon,klev+1), psfl(klon,klev+1)
    922 c
    923       REAL rain_lsc(klon)
    924       REAL snow_lsc(klon)
    925843c
    926844      REAL ratqsc(klon,klev)
    … …  
    935853      logical ok_newmicro
    936854      save ok_newmicro
    937       real ref_liq(klon,klev), ref_ice(klon,klev)
    938855c$OMP THREADPRIVATE(ok_newmicro)
    939856      save fact_cldcon,facttemps
    … …  
    973890c Variables locales pour effectuer les appels en serie
    974891c
    975       REAL zx_rh(klon,klev)
    976892cIM RH a 2m (la surface)
    977       REAL rh2m(klon), qsat2m(klon)
    978       REAL tpot(klon), tpote(klon)
    979893      REAL Lheat
    980894
    … …  
    10911005      REAL ZRCPD
    10921006c-jld ec_conser
    1093       REAL t2m(klon,nbsrf)  ! temperature a 2m
    10941007      REAL q2m(klon,nbsrf)  ! humidite a 2m
    10951008
    10961009cIM: t2m, q2m, ustar, u10m, v10m et t2mincels, t2maxcels
    1097       REAL zt2m(klon), zq2m(klon)             !temp., hum. 2m moyenne s/ 1 maille
    1098       REAL zustar(klon),zu10m(klon), zv10m(klon)  ! u* et vents a 10m moyennes s/1 maille
    10991010      CHARACTER*40 t2mincels, t2maxcels       !t2m min., t2m max
    11001011      CHARACTER*40 tinst, tave, typeval
    11011012      REAL cldtaupi(klon,klev)  ! Cloud optical thickness for pre-industrial (pi) aerosols
    11021013
    1103       REAL re(klon, klev)       ! Cloud droplet effective radius
    1104       REAL fl(klon, klev)  ! denominator of re
    11051014
    11061015      REAL re_top(klon), fl_top(klon) ! CDR at top of liquid water clouds
    … …  
    12281137! Utilise notamment en 1D mais peut etre active egalement en 3D
    12291138! en imposant la valeur de igout.
    1230 c======================================================================
    1231 
     1139c======================================================================d
    12321140      if (prt_level.ge.1) then
    12331141          igout=klon/2+1/klon
    … …  
    12501158
    12511159c======================================================================
    1252 
    1253 cym => necessaire pour iflag_con != 2   
    1254       pmfd(:,:) = 0.
    1255       pen_u(:,:) = 0.
    1256       pen_d(:,:) = 0.
    1257       pde_d(:,:) = 0.
    1258       pde_u(:,:) = 0.
    1259       aam=0.
    1260 
    1261       torsfc=0.
    1262       forall (k=1: llm) zmasse(:, k) = (paprs(:, k)-paprs(:, k+1)) / rg
    12631160
    12641161      if (first) then
    … …  
    13011198
    13021199      endif  ! first
     1200
     1201cym => necessaire pour iflag_con != 2   
     1202      pmfd(:,:) = 0.
     1203      pen_u(:,:) = 0.
     1204      pen_d(:,:) = 0.
     1205      pde_d(:,:) = 0.
     1206      pde_u(:,:) = 0.
     1207      aam=0.
     1208
     1209      torsfc=0.
     1210      forall (k=1: llm) zmasse(:, k) = (paprs(:, k)-paprs(:, k+1)) / rg
     1211
     1212
    13031213
    13041214       modname = 'physiq'
Note: See TracChangeset for help on using the changeset viewer.