Changeset 4795 for LMDZ6/branches


Ignore:
Timestamp:
Feb 1, 2024, 4:53:24 PM (4 months ago)
Author:
nfevrier
Message:

First save of N. Février's modifications

Location:
LMDZ6/branches/LMDZ_DECOUPLE
Files:
13 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/LMDZ_DECOUPLE/DefLists/field_def_lmdz.xml

    r3854 r4795  
    173173        <field id="flat_x"    long_name="flat off-wake"    unit="W/m2" />
    174174        <field id="flat_w"    long_name="flat within-wake"    unit="W/m2" />
    175         <field id="delta_tsurf"    long_name="Temperature difference (w-x)"    unit="K" enabled="false" /> <!-- ym --field disabled-- delta_tsurf is not conform to the associated grid (the is an arry containing all sub_surface) -->
     175        <field id="delta_tsurf"    long_name="Temperature difference (w-x)"    unit="K" enabled="false" /> <!-- ym field disabled delta_tsurf is not conform to the associated grid (the is an arry containing all sub_surface) -->
    176176        <field id="cdragh_x"    long_name="cdragh off-wake"    unit="-" />
    177177        <field id="cdragh_w"    long_name="cdragh within-wake"    unit="-" />
    … …  
    713713        <field id="cool_volc"    long_name="LW cooling rate from volcano"    unit="K/s" />
    714714        <field id="pvap"    long_name="pvap intermediary variable" unit="-">pres*ovap*461.5 / (287.04*(1.+ (10.9491/18.0153)*ovap)) </field>
     715        <!-- Nicolas Février -->
     716        <field id="dAIBCM_bcscav"    long_name="Below-cloud scavenging of AIBCM"    unit="kg/kg/s" />
     717        <field id="dAIBCM_insc"    long_name="In-cloud scavenging of AIBCM"    unit="kg/kg/s" />
     718        <field id="dAIBCM_evapls"    long_name="Reevaporation of rain for AIBCM"    unit="kg/kg/s" />
     719        <field id="dAIBCM_con"    long_name="Convective scavenging of AIBCM"    unit="kg/kg/s" />
     720        <field id="dAIBCM_vdf"    long_name="Boundary layer of AIBCM"    unit="kg/kg/s" />
     721
     722        <field id="dASBCM_bcscav"    long_name="Below-cloud scavenging of ASBCM"    unit="kg/kg/s" />
     723        <field id="dASBCM_insc"    long_name="In-cloud scavenging of ASBCM"    unit="kg/kg/s" />
     724        <field id="dASBCM_evapls"    long_name="Reevaporation of rain for ASBCM"    unit="kg/kg/s" />
     725        <field id="dASBCM_con"    long_name="Convective scavenging of ASBCM"    unit="kg/kg/s" />
     726        <field id="dASBCM_vdf"    long_name="Boundary layer of ASBCM"    unit="kg/kg/s" />
     727
     728        <field id="dASAPp1a_bcscav"    long_name="Below-cloud scavenging of ASAPp1a"    unit="kg/kg/s" />
     729        <field id="dASAPp1a_insc"    long_name="In-cloud scavenging of ASAPp1a"    unit="kg/kg/s" />
     730        <field id="dASAPp1a_evapls"    long_name="Reevaporation of rain for ASAPp1a"    unit="kg/kg/s" />
     731        <field id="dASAPp1a_con"    long_name="Convective scavenging of ASAPp1a"    unit="kg/kg/s" />
     732        <field id="dASAPp1a_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     733
     734        <field id="dASAPp2a_bcscav"    long_name="Below-cloud scavenging of ASAPp2a"    unit="kg/kg/s" />
     735        <field id="dASAPp2a_insc"    long_name="In-cloud scavenging of ASAPp2a"    unit="kg/kg/s" />
     736        <field id="dASAPp2a_evapls"    long_name="Reevaporation of rain for ASAPp2a"    unit="kg/kg/s" />
     737        <field id="dASAPp2a_con"    long_name="Convective scavenging of ASAPp2a"    unit="kg/kg/s" />
     738        <field id="dASAPp2a_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     739
     740        <field id="dASARp1a_bcscav"    long_name="Below-cloud scavenging of ASARp1a"    unit="kg/kg/s" />
     741        <field id="dASARp1a_insc"    long_name="In-cloud scavenging of ASARp1a"    unit="kg/kg/s" />
     742        <field id="dASARp1a_evapls"    long_name="Reevaporation of rain for ASARp1a"    unit="kg/kg/s" />
     743        <field id="dASARp1a_con"    long_name="Convective scavenging of ASARp1a"    unit="kg/kg/s" />
     744        <field id="dASARp1a_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     745
     746        <field id="dASARp2a_bcscav"    long_name="Below-cloud scavenging of ASARp2a"    unit="kg/kg/s" />
     747        <field id="dASARp2a_insc"    long_name="In-cloud scavenging of ASARp2a"    unit="kg/kg/s" />
     748        <field id="dASARp2a_evapls"    long_name="Reevaporation of rain for ASARp2a"    unit="kg/kg/s" />
     749        <field id="dASARp2a_con"    long_name="Convective scavenging of ASARp2a"    unit="kg/kg/s" />
     750        <field id="dASARp2a_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     751
     752        <field id="dAIPOMM_bcscav"    long_name="Below-cloud scavenging of AIPOMM"    unit="kg/kg/s" />
     753        <field id="dAIPOMM_insc"    long_name="In-cloud scavenging of AIPOMM"    unit="kg/kg/s" />
     754        <field id="dAIPOMM_evapls"    long_name="Reevaporation of rain for AIPOMM"    unit="kg/kg/s" />
     755        <field id="dAIPOMM_con"    long_name="Convective scavenging of AIPOMM"    unit="kg/kg/s" />
     756        <field id="dAIPOMM_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     757
     758        <field id="dASPOMM_bcscav"    long_name="Below-cloud scavenging of ASPOMM"    unit="kg/kg/s" />
     759        <field id="dASPOMM_insc"    long_name="In-cloud scavenging of ASPOMM"    unit="kg/kg/s" />
     760        <field id="dASPOMM_evapls"    long_name="Reevaporation of rain for ASPOMM"    unit="kg/kg/s" />
     761        <field id="dASPOMM_con"    long_name="Convective scavenging of ASPOMM"    unit="kg/kg/s" />
     762        <field id="dASPOMM_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     763
     764        <field id="dCIDUSTM_bcscav"    long_name="Below-cloud scavenging of CIDUSTM"    unit="kg/kg/s" />
     765        <field id="dCIDUSTM_insc"    long_name="In-cloud scavenging of CIDUSTM"    unit="kg/kg/s" />
     766        <field id="dCIDUSTM_evapls"    long_name="Reevaporation of rain for CIDUSTM"    unit="kg/kg/s" />
     767        <field id="dCIDUSTM_con"    long_name="Convective scavenging of CIDUSTM"    unit="kg/kg/s" />
     768        <field id="dCIDUSTM_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     769
     770        <field id="dASSSM_bcscav"    long_name="Below-cloud scavenging of ASSSM"    unit="kg/kg/s" />
     771        <field id="dASSSM_insc"    long_name="In-cloud scavenging of ASSSM"    unit="kg/kg/s" />
     772        <field id="dASSSM_evapls"    long_name="Reevaporation of rain for ASSSM"    unit="kg/kg/s" />
     773        <field id="dASSSM_con"    long_name="Convective scavenging of ASSSM"    unit="kg/kg/s" />
     774        <field id="dASSSM_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     775
     776        <field id="dCSSSM_bcscav"    long_name="Below-cloud scavenging of CSSSM"    unit="kg/kg/s" />
     777        <field id="dCSSSM_insc"    long_name="In-cloud scavenging of CSSSM"    unit="kg/kg/s" />
     778        <field id="dCSSSM_evapls"    long_name="Reevaporation of rain for CSSSM"    unit="kg/kg/s" />
     779        <field id="dCSSSM_con"    long_name="Convective scavenging of CSSSM"    unit="kg/kg/s" />
     780        <field id="dCSSSM_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     781
     782        <field id="dSSSSM_bcscav"    long_name="Below-cloud scavenging of SSSSM"    unit="kg/kg/s" />
     783        <field id="dSSSSM_insc"    long_name="In-cloud scavenging of SSSSM"    unit="kg/kg/s" />
     784        <field id="dSSSSM_evapls"    long_name="Reevaporation of rain for SSSSM"    unit="kg/kg/s" />
     785        <field id="dSSSSM_con"    long_name="Convective scavenging of SSSSM"    unit="kg/kg/s" />
     786        <field id="dSSSSM_vdf"    long_name="Boundary layer"    unit="kg/kg/s" />
     787
     788
    715789    </field_group>
    716790

  • LMDZ6/branches/LMDZ_DECOUPLE/DefLists/file_def_histday_lmdz.xml

    r3854 r4795  
    615615                <field field_ref="l_mixmin_sic" level="10" />
    616616                <field field_ref="ozone_daylight" level="10" />
     617                <!-- Nicolas Février -->
     618                <field field_ref="dAIBCM_bcscav" level="5" />
     619                <field field_ref="dAIBCM_insc" level="5" />
     620                <field field_ref="dAIBCM_evapls" level="5" />
     621                <field field_ref="dAIBCM_con" level="5" />
     622                <field field_ref="dAIBCM_vdf" level="5" />
     623
     624                <field field_ref="dASBCM_bcscav" level="5" />
     625                <field field_ref="dASBCM_insc" level="5" />
     626                <field field_ref="dASBCM_evapls" level="5" />
     627                <field field_ref="dASBCM_con" level="5" />
     628                <field field_ref="dASBCM_vdf" level="5" />
     629
     630                <field field_ref="dASAPp1a_bcscav" level="5" />
     631                <field field_ref="dASAPp1a_insc" level="5" />
     632                <field field_ref="dASAPp1a_evapls" level="5" />
     633                <field field_ref="dASAPp1a_con" level="5" />
     634                <field field_ref="dASAPp1a_vdf" level="5" />
     635
     636                <field field_ref="dASAPp2a_bcscav" level="5" />
     637                <field field_ref="dASAPp2a_insc" level="5" />
     638                <field field_ref="dASAPp2a_evapls" level="5" />
     639                <field field_ref="dASAPp2a_con" level="5" />
     640                <field field_ref="dASAPp2a_vdf" level="5" />
     641
     642                <field field_ref="dASARp1a_bcscav" level="5" />
     643                <field field_ref="dASARp1a_insc" level="5" />
     644                <field field_ref="dASARp1a_evapls" level="5" />
     645                <field field_ref="dASARp1a_con" level="5" />
     646                <field field_ref="dASARp1a_vdf" level="5" />
     647
     648                <field field_ref="dASARp2a_bcscav" level="5" />
     649                <field field_ref="dASARp2a_insc" level="5" />
     650                <field field_ref="dASARp2a_evapls" level="5" />
     651                <field field_ref="dASARp2a_con" level="5" />
     652                <field field_ref="dASARp2a_vdf" level="5" />
     653
     654                <field field_ref="dAIPOMM_bcscav" level="5" />
     655                <field field_ref="dAIPOMM_insc" level="5" />
     656                <field field_ref="dAIPOMM_evapls" level="5" />
     657                <field field_ref="dAIPOMM_con" level="5" />
     658                <field field_ref="dAIPOMM_vdf" level="5" />
     659
     660                <field field_ref="dASPOMM_bcscav" level="5" />
     661                <field field_ref="dASPOMM_insc" level="5" />
     662                <field field_ref="dASPOMM_evapls" level="5" />
     663                <field field_ref="dASPOMM_con" level="5" />
     664                <field field_ref="dASPOMM_vdf" level="5" />
     665
     666                <field field_ref="dCIDUSTM_bcscav" level="5" />
     667                <field field_ref="dCIDUSTM_insc" level="5" />
     668                <field field_ref="dCIDUSTM_evapls" level="5" />
     669                <field field_ref="dCIDUSTM_con" level="5" />
     670                <field field_ref="dCIDUSTM_vdf" level="5" />
     671
     672                <field field_ref="dASSSM_bcscav" level="5" />
     673                <field field_ref="dASSSM_insc" level="5" />
     674                <field field_ref="dASSSM_evapls" level="5" />
     675                <field field_ref="dASSSM_con" level="5" />
     676                <field field_ref="dASSSM_vdf" level="5" />
     677
     678                <field field_ref="dCSSSM_bcscav" level="5" />
     679                <field field_ref="dCSSSM_insc" level="5" />
     680                <field field_ref="dCSSSM_evapls" level="5" />
     681                <field field_ref="dCSSSM_con" level="5" />
     682                <field field_ref="dCSSSM_vdf" level="5" />
     683
     684                <field field_ref="dSSSSM_bcscav" level="5" />
     685                <field field_ref="dSSSSM_insc" level="5" />
     686                <field field_ref="dSSSSM_evapls" level="5" />
     687                <field field_ref="dSSSSM_con" level="5" />
     688                <field field_ref="dSSSSM_vdf" level="5" />
     689
    617690            </field_group>
    618691            <field_group operation="average" grid_ref="grid_out_presinter">
    … …  
    640713                <field field_ref="kz_max" level="4" operation="maximum" />
    641714                <field field_ref="evu" level="10" />
     715
    642716            </field_group>
    643717           

  • LMDZ6/branches/LMDZ_DECOUPLE/DefLists/file_def_histmth_lmdz.xml

    r3854 r4795  
    550550                <field field_ref="ozone" level="3" />
    551551                <field field_ref="ozone_daylight" level="2" />
     552
     553                <field field_ref="dAIBCM_bcscav" level="5" />
     554                <field field_ref="dAIBCM_insc" level="5" />
     555                <field field_ref="dAIBCM_evapls" level="5" />
     556                <field field_ref="dAIBCM_con" level="5" />
     557                <field field_ref="dAIBCM_vdf" level="5" />
     558
     559                <field field_ref="dASBCM_bcscav" level="5" />
     560                <field field_ref="dASBCM_insc" level="5" />
     561                <field field_ref="dASBCM_evapls" level="5" />
     562                <field field_ref="dASBCM_con" level="5" />
     563                <field field_ref="dASBCM_vdf" level="5" />
     564
     565                <field field_ref="dASAPp1a_bcscav" level="5" />
     566                <field field_ref="dASAPp1a_insc" level="5" />
     567                <field field_ref="dASAPp1a_evapls" level="5" />
     568                <field field_ref="dASAPp1a_con" level="5" />
     569                <field field_ref="dASAPp1a_vdf" level="5" />
     570
     571                <field field_ref="dASAPp2a_bcscav" level="5" />
     572                <field field_ref="dASAPp2a_insc" level="5" />
     573                <field field_ref="dASAPp2a_evapls" level="5" />
     574                <field field_ref="dASAPp2a_con" level="5" />
     575                <field field_ref="dASAPp2a_vdf" level="5" />
     576
     577                <field field_ref="dASARp1a_bcscav" level="5" />
     578                <field field_ref="dASARp1a_insc" level="5" />
     579                <field field_ref="dASARp1a_evapls" level="5" />
     580                <field field_ref="dASARp1a_con" level="5" />
     581                <field field_ref="dASARp1a_vdf" level="5" />
     582
     583                <field field_ref="dASARp2a_bcscav" level="5" />
     584                <field field_ref="dASARp2a_insc" level="5" />
     585                <field field_ref="dASARp2a_evapls" level="5" />
     586                <field field_ref="dASARp2a_con" level="5" />
     587                <field field_ref="dASARp2a_vdf" level="5" />
     588
     589                <field field_ref="dAIPOMM_bcscav" level="5" />
     590                <field field_ref="dAIPOMM_insc" level="5" />
     591                <field field_ref="dAIPOMM_evapls" level="5" />
     592                <field field_ref="dAIPOMM_con" level="5" />
     593                <field field_ref="dAIPOMM_vdf" level="5" />
     594
     595                <field field_ref="dASPOMM_bcscav" level="5" />
     596                <field field_ref="dASPOMM_insc" level="5" />
     597                <field field_ref="dASPOMM_evapls" level="5" />
     598                <field field_ref="dASPOMM_con" level="5" />
     599
     600                <field field_ref="dCIDUSTM_bcscav" level="5" />
     601                <field field_ref="dCIDUSTM_insc" level="5" />
     602                <field field_ref="dCIDUSTM_evapls" level="5" />
     603                <field field_ref="dCIDUSTM_con" level="5" />
     604                <field field_ref="dCIDUSTM_vdf" level="5" />
     605
     606                <field field_ref="dASSSM_bcscav" level="5" />
     607                <field field_ref="dASSSM_insc" level="5" />
     608                <field field_ref="dASSSM_evapls" level="5" />
     609                <field field_ref="dASSSM_con" level="5" />
     610                <field field_ref="dASSSM_vdf" level="5" />
     611
     612                <field field_ref="dCSSSM_bcscav" level="5" />
     613                <field field_ref="dCSSSM_insc" level="5" />
     614                <field field_ref="dCSSSM_evapls" level="5" />
     615                <field field_ref="dCSSSM_con" level="5" />
     616                <field field_ref="dCSSSM_vdf" level="5" />
     617
     618                <field field_ref="dSSSSM_bcscav" level="5" />
     619                <field field_ref="dSSSSM_insc" level="5" />
     620                <field field_ref="dSSSSM_evapls" level="5" />
     621                <field field_ref="dSSSSM_con" level="5" />
     622                <field field_ref="dSSSSM_vdf" level="5" />
     623
    552624                <field field_ref="upwd" level="7" />
    553625                <field field_ref="ep" level="7" />
    … …  
    664736                <field field_ref="cool_volc" level="1" />
    665737                <field field_ref="heat_volc" level="1" />
     738
    666739            </field_group>
    667740            <field_group operation="average" grid_ref="grid_out_presinter">

  • LMDZ6/branches/LMDZ_DECOUPLE/arch/arch-X64_IRENE.fcm

    r3435 r4795  
    1414%MPI_FFLAGS
    1515%OMP_FFLAGS          -qopenmp
    16 %BASE_LD             -i4 -r8 -auto -L/ccc/products/mkl-17.0.4.196/intel--17.0.4.196__openmpi--2.0.2/default/17.0.4.196/mkl/lib/intel64 -lmkl_intel_lp64 -lmkl_core -lmkl_sequential -lpthread
     16%BASE_LD -i4 -r8 -auto -L/ccc/products/mkl-20.0.0/system/default/20.0.0/mkl/lib/intel64 -lmkl_intel_lp64 -lmkl_core -lmkl_sequential -lpthread
    1717%MPI_LD
    1818%OMP_LD              -qopenmp

  • LMDZ6/branches/LMDZ_DECOUPLE/libf/dyn3dmem/dynredem_loc.F90

    r2622 r4795  
    1111  USE infotrac
    1212  USE netcdf, ONLY: NF90_CREATE, NF90_DEF_DIM, NF90_INQ_VARID, NF90_GLOBAL,    &
    13                     NF90_CLOSE,  NF90_PUT_ATT, NF90_UNLIMITED, NF90_CLOBBER
     13                    NF90_CLOSE,  NF90_PUT_ATT, NF90_UNLIMITED, NF90_CLOBBER, &
     14                    NF90_64BIT_OFFSET
    1415  USE dynredem_mod, ONLY: cre_var, put_var, err, modname, fil
    1516  USE comvert_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt,pa,preff, &
    … …  
    110111
    111112!--- File creation
    112   CALL err(NF90_CREATE(fichnom,NF90_CLOBBER,nid))
     113  CALL err(NF90_CREATE(fichnom,IOR(NF90_CLOBBER,NF90_64BIT_OFFSET),nid))
    113114
    114115!--- Some global attributes

  • LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/Dust/inscav_spl.F

    r2630 r4795  
    11      SUBROUTINE inscav_spl(pdtime,it,masse,henry,kk,qliq,
    22     .                   flxr,flxs,zrho,zdz,t,x,
    3      .                   his_dh)
     3     .                   dx)
    44      USE dimphy
    55      IMPLICIT NONE

  • LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/lsc_scav.F90

    r2351 r4795  
    77                    oliq,flxr,flxs,rneb,beta_fisrt,  &
    88                    beta_v1,pplay,paprs,t,tr_seri,d_tr_insc,          &
    9                     d_tr_bcscav,d_tr_evap,qPrls)
     9                    d_tr_bcscav,d_tr_evap,qPrls,beta_st,j)
    1010  USE ioipsl
    1111  USE dimphy
    … …  
    1313  USE mod_phys_lmdz_para
    1414  USE traclmdz_mod
    15   USE infotrac_phy,ONLY : nbtr
     15  USE infotrac_phy,ONLY : nbtr,solsym
    1616  USE iophy
    1717  IMPLICIT NONE
    … …  
    4747  REAL,DIMENSION(klon,klev,nbtr),INTENT(OUT)  :: d_tr_bcscav  ! tendance de traceur
    4848  REAL,DIMENSION(klon,klev,nbtr),INTENT(OUT)  :: d_tr_evap
     49  ! Nicolas Février
     50  REAL,DIMENSION(klon,klev,12),INTENT(OUT)       :: beta_st
     51  INTEGER,INTENT(IN)                             :: j           ! indice parmi les 12 traceurs de WETDEP
     52
    4953  REAL,DIMENSION(klon,nbtr),INTENT(OUT)       :: qPrls      !jyg: concentration tra dans pluie LS a la surf.
    5054  REAL :: dxin,dxev                              ! tendance temporaire de traceur incloud
    … …  
    111115      R_r=0.001            !  mean raindrop radius (m)
    112116      R_s=0.001            !  mean snow crystal radius (m)
    113       frac_fine_scav=0.7
    114       frac_coar_scav=0.7
     117      !frac_fine_scav=0.7
     118      !frac_coar_scav=0.7
    115119!     frac_aer=0.5 ~ droplet size shrinks by evap
    116120      frac_aer=0.5
    … …  
    133137   print*,'R_r',R_r
    134138   print*,'R_s',R_s
    135    print*,'frac_fine_scav',frac_fine_scav
    136    print*,'frac_coar_scav',frac_coar_scav
     139   !print*,'frac_fine_scav',frac_fine_scav
     140   !print*,'frac_coar_scav',frac_coar_scav
    137141   print*,'frac_aer ev',frac_aer
    138142
    … …  
    142146!
    143147! initialization
     148
     149  ! Nicolas Février
     150  ! Différentiation BC soluble/insoluble
     151
     152  IF ((TRIM(solsym(it))=='AIBCM') .OR. (TRIM(solsym(it))=='AIPOMM') .OR. (TRIM(solsym(it))=='CIDUSTM')) THEN
     153     frac_fine_scav=0.2
     154     frac_coar_scav=0.2
     155     !PRINT *,'insoluble okidoki'
     156  ELSE
     157     frac_fine_scav=0.7
     158     frac_coar_scav=0.7
     159  ENDIF
     160
    144161  dxin=0.
    145162  dxev=0.
    … …  
    223240      his_dh(i)=his_dh(i)-dxin*zmass(i,k)/pdtime !  kg/m2/s
    224241      d_tr_insc(i,k,it)=dxin
     242      ! Nicolas Février
     243      beta_st(i,k,j)=beta_st(i,k,j)-(exp(-scav(i,k)*beta*pdtime)-1.)/pdtime
    225244
    226245!  below-cloud impaction
    … …  
    239258     his_dh(i)=his_dh(i)-dxbc(i,k)*zmass(i,k)/pdtime !  kg/m2/s
    240259     d_tr_bcscav(i,k,it)=dxbc(i,k)
     260     ! Nicolas Février
     261     beta_st(i,k,j)=beta_st(i,k,j)+3./4.*(water+ice)
    241262    ENDIF
    242263
    … …  
    279300   ENDDO  !! do k
    280301
     302! Debug beta_st
     303    print *,"BETA LSC_SCAV",beta_st(5000,35,7)
     304
     305
    281306!jyg (20130114)
    282307   DO i = 1,klon
    … …  
    298323!     ENDDO
    299324
     325!print*, "TAGADATSOINTSOIN"
     326
    300327END SUBROUTINE lsc_scav

  • LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/minmaxqfi.F90

    r2311 r4795  
    22! $Id$
    33!
    4 SUBROUTINE minmaxqfi(zq,qmin,qmax,comment)
     4SUBROUTINE minmaxqfi(zq,qmin,qmax,comment,it)
    55  USE dimphy
    66  USE print_control_mod, ONLY: prt_level
     7  USE infotrac_phy, ONLY: nbtr, type_trac, conv_flg, solsym
    78  IMPLICIT NONE
    89
    … …  
    1112  REAL,INTENT(IN)                         :: qmin,qmax
    1213  CHARACTER(LEN=*),INTENT(IN)             :: comment
     14  INTEGER,INTENT(IN)                      :: it
    1315
    1416! Local 
    … …  
    2426        ENDIF
    2527     ENDDO
    26      IF (jbad.GT.0) THEN
     28     IF ((jbad.GT.0) .AND. ((TRIM(solsym(it))=='AIBCM') .OR. (TRIM(solsym(it))=='ASBCM'))) THEN ! Nicolas Février : alerte seulement si le BC est concerné
    2729        WRITE(*,*)comment, "k = ", k, "jbad > 0"
    28         if (prt_level >= 1) then
     30        if (prt_level >= 1) then ! Nicolas Février on print tout
    2931           DO i = 1, jbad
    3032              WRITE(*,*) "i,q=", jadrs(i), zq(jadrs(i),k)
    … …  
    3335     ENDIF
    3436  ENDDO
    35  
     37
    3638END SUBROUTINE minmaxqfi

  • LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/phys_output_mod.F90

    r3853 r4795  
    143143    REAL, DIMENSION(NSW,2)          :: spbnds_sun !bounds of spectband
    144144
    145     WRITE(lunout,*) 'Debut phys_output_mod.F90'
     145    WRITE(lunout,*) 'DEBUG - ANNE Debut phys_output_mod.F90'
    146146    ! Initialisations (Valeurs par defaut
    147147
    … …  
    312312  ecrit_files(9) = ecrit_files(3)
    313313
     314  WRITE(lunout,*) 'DEBUG - ANNE nfiles=',nfiles
     315
    314316  DO iff=1,nfiles
    315317
    … …  
    498500            DO iq=nqo+1,nqtot
    499501            iiq=niadv(iq)
     502           
     503            !WRITE(lunout,*) 'UWUWU nom à utiliser : ','d'//trim(tname(iq))//'_bcscav'
     504            WRITE(lunout,*) 'DEBUG - ANNE / iq=',iq,'/ d//trim(tname(iq))//_vdf=','d'//trim(tname(iq))//'_vdf',' / Tendance tracer //ttext(iiq)=','Tendance tracer '//ttext(iiq)
     505
    500506            o_trac(iq-nqo) = ctrl_out((/ 1, 5, 5, 5, 10, 10, 11, 11, 11, 11 /), &
    501507                           tname(iiq),'Tracer '//ttext(iiq), "-",  &
    502508                           (/ '', '', '', '', '', '', '', '', '', '' /))
    503             o_dtr_vdf(iq-nqo) = ctrl_out((/ 4, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
     509            o_dtr_vdf(iq-nqo) = ctrl_out((/ 5, 5, 5, 5, 10, 10, 11, 11, 11, 11 /), &
    504510                              'd'//trim(tname(iq))//'_vdf',  &
    505511                              'Tendance tracer '//ttext(iiq), "-" , &
    … …  
    526532                                     (/ '', '', '', '', '', '', '', '', '', '' /))
    527533
    528             o_dtr_insc(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
     534            o_dtr_insc(iq-nqo) = ctrl_out((/ 5, 5, 5, 5, 10, 10, 11, 11, 11, 11 /), &
    529535                               'd'//trim(tname(iq))//'_insc', &
    530536                               'Tendance tracer '//ttext(iiq), "-", &
    531537                               (/ '', '', '', '', '', '', '', '', '', '' /))
    532538
    533             o_dtr_bcscav(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
     539            o_dtr_bcscav(iq-nqo) = ctrl_out((/ 5, 5, 5, 5, 10, 10, 11, 11, 11, 11 /), &
    534540                                 'd'//trim(tname(iq))//'_bcscav', &
    535541                                 'Tendance tracer '//ttext(iiq), "-", &
    536542                                 (/ '', '', '', '', '', '', '', '', '', '' /))
    537543
    538             o_dtr_evapls(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
     544            o_dtr_evapls(iq-nqo) = ctrl_out((/ 5, 5, 5, 5, 10, 10, 11, 11, 11, 11 /), &
    539545                                 'd'//trim(tname(iq))//'_evapls', &
    540546                                 'Tendance tracer '//ttext(iiq), "-", &
    … …  
    566572                                 (/ '', '', '', '', '', '', '', '', '', '' /))
    567573
    568             o_dtr_dry(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
    569                               'cum'//'d'//trim(tname(iq))//'_dry', &
     574            o_dtr_dry(iq-nqo) = ctrl_out((/ 5, 5, 5, 5, 10, 10, 11, 11, 11, 11 /), &
     575                              'd'//trim(tname(iq))//'_dry', &
    570576                              'tracer tendency dry deposition'//ttext(iiq), "-", &
    571577                              (/ '', '', '', '', '', '', '', '', '', '' /))
    … …  
    579585
    580586   ENDDO !  iff
     587
     588   WRITE(lunout,*)'DEBUG - ANNE / on est sorti de la boucle sur iff'
    581589
    582590    ! Updated write frequencies due to phys_out_filetimesteps.

  • LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/phys_output_write_mod.F90

    r3842 r4795  
    66  USE phytrac_mod, ONLY : d_tr_cl, d_tr_th, d_tr_cv, d_tr_lessi_impa, &
    77       d_tr_lessi_nucl, d_tr_insc, d_tr_bcscav, d_tr_evapls, d_tr_ls,  &
    8        d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav, flux_tr_dry
     8       d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav, flux_tr_dry,list_aero_wetdep
    99
    1010  ! Author: Abderrahmane IDELKADI (original include file)
    … …  
    446446    REAL,DIMENSION(klon,klev) :: z, dz
    447447    REAL,DIMENSION(klon)      :: zrho, zt
     448    ! Nicolas Février
     449    INTEGER :: j
    448450
    449451    ! On calcul le nouveau tau:
    … …  
    23492351#endif
    23502352!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     2353
     2354       ! Nicolas Février
     2355       ! On rajoute type_trac == 'inca'
    23512356       IF (iflag_phytrac == 1 ) then
    2352          IF (type_trac == 'lmdz' .OR. type_trac == 'coag') THEN
    2353            DO iq=nqo+1, nqtot
     2357         IF (type_trac == 'lmdz' .OR. type_trac == 'coag' .OR. type_trac == 'inca') THEN
     2358           ! Nicolas Février
     2359           ! On ne fait la boucle que sur les aérosols voulus
     2360           !DO iq=nqo+1, nqtot
     2361           DO j = 1,size(list_aero_wetdep)
     2362           !DO iq=6
     2363           !iq=list_aero_wetdep(6)+2
     2364             iq=list_aero_wetdep(j)+3
    23542365             !--3D fields
    2355              CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
     2366             !CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
    23562367             CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
    2357              CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
     2368             !CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
    23582369             CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo))
    2359              CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))
    2360              CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))
     2370             !CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))
     2371             !CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))
    23612372             CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo))
    23622373             CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo))
    23632374             CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo))
    2364              CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))
    2365              CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))
    2366              CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))
    2367              CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))
    2368              CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo))
     2375             !CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))
     2376             !CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))
     2377             !CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))
     2378             !CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))
     2379             !/CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo))
    23692380            !--2D fields
    2370              CALL histwrite_phy(o_dtr_dry(iq-nqo), flux_tr_dry(:,iq-nqo))
     2381             !CALL histwrite_phy(o_dtr_dry(iq-nqo), flux_tr_dry(:,iq-nqo))
    23712382             zx_tmp_fi2d=0.
    23722383             IF (vars_defined) THEN
    … …  
    23752386                ENDDO
    23762387             ENDIF
    2377              CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
     2388             !CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
    23782389           ENDDO !--iq
    23792390         ENDIF   !--type_trac

  • LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/phytrac_mod.F90

    r3649 r4795  
    4545  REAL,DIMENSION(:,:),ALLOCATABLE,SAVE      :: d_tr_dry ! Td depot sec/traceur (1st layer),ALLOCATABLE,SAVE  jyg
    4646  REAL,DIMENSION(:,:),ALLOCATABLE,SAVE      :: flux_tr_dry ! depot sec/traceur (surface),ALLOCATABLE,SAVE    jyg
     47  !Nicolas Février
     48  INTEGER,DIMENSION(12),SAVE :: list_aero_wetdep
     49  REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: beta_st,beta_cv
     50!$OMP THREADPRIVATE(list_aero_wetdep,beta_st,beta_cv)
     51
    4752
    4853!$OMP THREADPRIVATE(qPa,qMel,qTrdi,dtrcvMA,d_tr_th,d_tr_lessi_impa,d_tr_lessi_nucl)
    … …  
    7176       ALLOCATE(d_tr_th(klon,klev,nbtr))
    7277       ALLOCATE(d_tr_lessi_impa(klon,klev,nbtr),d_tr_lessi_nucl(klon,klev,nbtr))
     78       ! Nicolas Février
     79       ALLOCATE(beta_st(klon,klev,12),beta_cv(klon,klev,12))
    7380
    7481
    … …  
    295302    !=======================================================================================
    296303
    297     INTEGER :: i, k, it
     304    INTEGER :: i, k, it, j ! j rajouté par Nicolas Février
    298305    INTEGER :: nsplit
    299306
    … …  
    361368!$OMP THREADPRIVATE(ccntrAA_in,ccntrENV_in,coefcoli_in)
    362369!$OMP THREADPRIVATE(convscav)
     370
     371
     372    ! Nicolas Février
     373
     374    INTEGER,SAVE :: id_ASAPp1a,id_ASAPp2a,id_ASARp1a,id_ASARp2a,id_ASBCM,id_AIBCM,id_CIDUSTM,id_AIPOMM,id_ASPOMM,id_ASSSM,id_CSSSM,id_SSSSM
     375!$OMP THREADPRIVATE(id_ASAPp1a,id_ASAPp2a,id_ASARp1a,id_ASARp2a,id_ASBCM,id_AIBCM,id_CIDUSTM,id_AIPOMM,id_ASPOMM,id_ASSSM,id_CSSSM,id_SSSSM)
     376
    363377    !RomP <<<
    364378    !######################################################################
    365379    !                    -- INITIALIZATION --
    366380    !######################################################################
     381
     382    ! Nicolas Février
     383    ! On force iflag_lscav à 4
     384    !iflag_lscav=4
    367385
    368386    DO k=1,klev
    … …  
    374392       END DO
    375393    END DO
     394
     395   ! Nicolas Février
     396   DO j=1, 12
     397      DO k=1,klev
     398         DO i=1,klon
     399              beta_st(i,k,j)=0.
     400              beta_cv(i,k,j)=0.
     401          END DO
     402       END DO
     403    END DO
     404
    376405
    377406    DO it=1, nbtr
    … …  
    452481       !
    453482!$OMP MASTER
    454        iflag_lscav_omp=1
     483       !iflag_lscav_omp=4 !Nicolas Février
    455484       call getin('iflag_lscav', iflag_lscav_omp)
    456485       ccntrAA_omp=1
    … …  
    495524          source(:,:)=init_source(:,:)
    496525          CALL tracinca_init(aerosol,lessivage)
     526          ! Nicolas Février
     527          ! On met le lessivage à TRUE tout en désactivant le lessivage de INCA
     528          ! correspondant pour ne pas faire le lessivage deux fois.
     529          ! Exemple : commenter CALL WETDEP dans aerosolmain.F90
     530          ! lessivage=.TRUE.
    497531       CASE('repr')
    498532          source(:,:)=0.
    … …  
    614648       write(lunout,*)  'flag_cvltr    = ', flag_cvltr
    615649
    616       IF (lessivage .AND. type_trac .EQ. 'inca') THEN
    617           CALL abort_physic('phytrac', 'lessivage=T config_inca=inca impossible',1)
    618 !          STOP
    619       ENDIF
     650      ! IF (lessivage .AND. type_trac .EQ. 'inca') THEN
     651      !    CALL abort_physic('phytrac', 'lessivage=T config_inca=inca impossible',1)
     652!     !     STOP
     653      ! ENDIF
    620654       !
    621655    ENDIF ! of IF (debutphy)
    … …  
    680714
    681715    END SELECT
     716
     717! Nicolas Février
     718
     719    !IF (aerosol_couple) THEN   !--we get aerosols from tr_seri array from INCA
     720      !
     721      !--copy fields from INCA tr_seri
     722      !--convert to ug m-3 unit for consistency with offline fields
     723      !
     724      DO it=1,nbtr
     725         SELECT CASE(trim(solsym(it)))
     726           CASE ("ASAPp1a")
     727             id_ASAPp1a = it
     728             aerosol(it)=.TRUE.
     729           CASE ("ASAPp2a")
     730             id_ASAPp2a = it
     731             aerosol(it)=.TRUE.
     732           CASE ("ASARp1a")
     733             id_ASARp1a = it
     734             aerosol(it)=.TRUE.
     735           CASE ("ASARp2a")
     736             id_ASARp2a = it
     737             aerosol(it)=.TRUE.
     738           CASE ("ASBCM")
     739             id_ASBCM = it
     740             aerosol(it)=.TRUE.
     741           CASE ("AIBCM")
     742             id_AIBCM = it
     743             aerosol(it)=.TRUE.
     744           CASE ("CIDUSTM")
     745             id_CIDUSTM = it
     746             aerosol(it)=.TRUE.
     747           CASE ("AIPOMM")
     748             id_AIPOMM = it
     749             aerosol(it)=.TRUE.
     750           CASE ("ASPOMM")
     751             id_ASPOMM = it
     752             aerosol(it)=.TRUE.
     753           CASE ("ASSSM")
     754             id_ASSSM = it
     755             aerosol(it)=.TRUE.
     756           CASE ("CSSSM")
     757             id_CSSSM = it
     758             aerosol(it)=.TRUE.
     759           CASE ("SSSSM")
     760             id_SSSSM = it
     761             aerosol(it)=.TRUE.
     762         END SELECT
     763      ENDDO
     764    !END IF
     765
     766    list_aero_wetdep=(/ id_ASAPp1a,id_ASAPp2a,id_ASARp1a,id_ASARp2a,id_ASBCM,id_AIBCM,id_CIDUSTM,id_AIPOMM,id_ASPOMM,id_ASSSM,id_CSSSM,id_SSSSM /)
     767
     768    !WRITE(lunout,*) 'id_AIBCM=',id_AIBCM
     769    !WRITE(lunout,*) 'id_ASBCM=',id_ASBCM
     770    !WRITE(lunout,*) 'shape_lmdz tr_seri=',shape(tr_seri)
     771
    682772    !======================================================================
    683773    !       -- Calcul de l'effet de la convection --
    684774    !======================================================================
     775
     776    ! Nicolas Février
     777
     778    !WRITE(lunout,*) 'CHECK1 TR_SERI AIBCM',SUM(tr_seri(:,:,id_AIBCM))
     779    !WRITE(lunout,*) 'CHECK1 TR_SERI ASBCM',SUM(tr_seri(:,:,id_ASBCM))
    685780
    686781    IF (iflag_con_trac==1) THEN
    … …  
    699794          ELSE   
    700795             !--nouveau transport convectif de Emanuel
     796             
     797             ! Nicolas Février
     798
     799             IF ( ANY( list_aero_wetdep==it ) ) THEN
     800                flag_cvltr(it)=.TRUE.
     801                ccntrAA(it) =ccntrAA_in
     802                ccntrENV(it)=ccntrENV_in
     803                coefcoli(it)=coefcoli_in
     804             END IF
    701805
    702806             IF (flag_cvltr(it)) THEN
    … …  
    744848          END DO
    745849
    746           CALL minmaxqfi(tr_seri(:,:,it),0.,1.e33,'convection it = '//solsym(it))
     850          CALL minmaxqfi(tr_seri(:,:,it),0.,1.e33,'convection it = '//solsym(it),it)  ! Nicolas Février it rajouté en argument
     851
     852          !WRITE(lunout,*) 'CHECK2 TR_SERI AIBCM',SUM(tr_seri(:,:,id_AIBCM))
     853          !WRITE(lunout,*) 'CHECK2 TR_SERI ASBCM',SUM(tr_seri(:,:,id_ASBCM))
    747854
    748855       END DO ! nbtr
    … …  
    856963
    857964       DO it=1, nbtr
    858           !
     965          !WRITE(lunout,*) 'BOUM'
    859966          IF( pbl_flg(it) /= 0 ) THEN
     967             !WRITE(lunout,*) 'BADABOUM'
    860968             !
    861969             CALL cltrac(pdtphys, coefh,t_seri,       &
    … …  
    887995       !
    888996    ELSE IF (iflag_vdf_trac==0) THEN
     997
     998    !WRITE(lunout,*) 'TRALALA'
    889999       !
    8901000       !   Injection of source in the first model layer
    … …  
    8971007       !
    8981008    ELSE IF (iflag_vdf_trac==-1) THEN
     1009
     1010    !WRITE(lunout,*) 'BOFBOFBOF'
    8991011       !
    9001012       ! Nothing happens
    … …  
    9111023    !   POUR INCA le lessivage est fait directement dans INCA
    9121024    !======================================================================
     1025   
     1026    ! Nicolas Février
     1027
     1028    ! id_ASAPp1a,id_ASAPp2a,id_ASARp1a,id_ASARp2a,id_ASBCM,id_AIBCM,id_CIDUSTM,id_AIPOMM,id_ASPOMM,id_ASSSM,id_CSSSM,id_SSSSM
     1029
     1030!    IF (aerosol_couple) THEN   !--we get aerosols from tr_seri array from INCA
     1031!      !
     1032!      !--copy fields from INCA tr_seri
     1033!      !--convert to ug m-3 unit for consistency with offline fields
     1034!      !
     1035!      DO it=1,nbtr
     1036!         SELECT CASE(trim(solsym(it)))
     1037!           CASE ("ASAPp1a")
     1038!             id_ASAPp1a = it
     1039!             aerosol(it)=.TRUE.
     1040!           CASE ("ASAPp2a")
     1041!             id_ASAPp2a = it
     1042!             aerosol(it)=.TRUE.
     1043!           CASE ("ASARp1a")
     1044!             id_ASARp1a = it
     1045!             aerosol(it)=.TRUE.
     1046!           CASE ("ASARp2a")
     1047!             id_ASARp2a = it
     1048!             aerosol(it)=.TRUE.
     1049!           CASE ("ASBCM")
     1050!             id_ASBCM = it
     1051!             aerosol(it)=.TRUE.
     1052!           CASE ("AIBCM")
     1053!             id_AIBCM = it
     1054!             aerosol(it)=.TRUE.
     1055!           CASE ("CIDUSTM")
     1056!             id_CIDUSTM = it
     1057!             aerosol(it)=.TRUE.
     1058!           CASE ("AIPOMM")
     1059!             id_AIPOMM = it
     1060!             aerosol(it)=.TRUE.
     1061!           CASE ("ASPOMM")
     1062!             id_ASPOMM = it
     1063!             aerosol(it)=.TRUE.
     1064!           CASE ("ASSSM")
     1065!             id_ASSSM = it
     1066!             aerosol(it)=.TRUE.
     1067!           CASE ("CSSSM")
     1068!             id_CSSSM = it
     1069!             aerosol(it)=.TRUE.
     1070!           CASE ("SSSSM")
     1071!             id_SSSSM = it
     1072!             aerosol(it)=.TRUE.
     1073!         END SELECT
     1074!      ENDDO
     1075!    END IF
     1076
     1077!    list_aero_wetdep=(/ id_ASAPp1a,id_ASAPp2a,id_ASARp1a,id_ASARp2a,id_ASBCM,id_AIBCM,id_CIDUSTM,id_AIPOMM,id_ASPOMM,id_ASSSM,id_CSSSM,id_SSSSM /)
     1078   
     1079    lessivage = .TRUE.
     1080
     1081    !WRITE(lunout,*) 'CHECK3 TR_SERI AIBCM',SUM(tr_seri(:,:,id_AIBCM))
     1082    !WRITE(lunout,*) 'CHECK3 TR_SERI ASBCM',SUM(tr_seri(:,:,id_ASBCM))
     1083
    9131084
    9141085    IF (lessivage) THEN
    … …  
    9221093       ! choix du lessivage
    9231094       !
     1095
     1096       !WRITE(lunout,*) 'Nico iflag_lscav=',iflag_lscav
     1097
    9241098       IF (iflag_lscav .EQ. 3 .OR. iflag_lscav .EQ. 4) THEN
    9251099          ! ********  Olivier Boucher version (3) possibly with modified ql_incl (4)
    9261100          !
    927           DO it = 1, nbtr
     1101          DO j = 1,size(list_aero_wetdep)
     1102             it=list_aero_wetdep(j)
    9281103
    9291104             IF (aerosol(it)) THEN
    … …  
    9321107             ! Liu (2001) proposed to use 1.5e-3 kg/kg
    9331108
    934 !jyg<
    935 !!             CALL lsc_scav(pdtphys,it,iflag_lscav,ql_incl,prfl,psfl,rneb,beta_fisrt,  &
     1109             !WRITE(lunout,*) 'On_lessive_id = ',it
     1110
    9361111             CALL lsc_scav(pdtphys,it,iflag_lscav,aerosol,ql_incl,prfl,psfl,rneb,beta_fisrt,  &
    937 !>jyg
    9381112                  beta_v1,pplay,paprs,t_seri,tr_seri,d_tr_insc,   &
    939                   d_tr_bcscav,d_tr_evapls,qPrls)
     1113                  d_tr_bcscav,d_tr_evapls,qPrls,beta_st,j)
    9401114
    9411115             !large scale scavenging tendency
    942              DO k = 1, klev
    943                 DO i = 1, klon
    944                    d_tr_ls(i,k,it)=d_tr_insc(i,k,it)+d_tr_bcscav(i,k,it)+d_tr_evapls(i,k,it)
    945                    tr_seri(i,k,it)=tr_seri(i,k,it)+d_tr_ls(i,k,it)
    946                 ENDDO
    947              ENDDO
    948              CALL minmaxqfi(tr_seri(:,:,it),0.,1.e33,'lsc scav it = '//solsym(it))
     1116             ! Nicolas Février : commenté pour ne pas faire le lessivage ici. Il est à la place fait dans SETHET.
     1117             !DO k = 1, klev
     1118             !   DO i = 1, klon
     1119             !      d_tr_ls(i,k,it)=d_tr_insc(i,k,it)+d_tr_bcscav(i,k,it)+d_tr_evapls(i,k,it)
     1120             !      tr_seri(i,k,it)=tr_seri(i,k,it)+d_tr_ls(i,k,it)
     1121             !   ENDDO
     1122             !ENDDO
     1123             !CALL minmaxqfi(tr_seri(:,:,it),0.,1.e33,'lsc scav it = '//solsym(it),it) ! Nicolas Février it rajouté en argument
    9491124             ENDIF
    9501125
    9511126          END DO  !tr
     1127
     1128! Debug beta_st
     1129    print *,"BETA PHYTRAC",beta_st(5000,35,7)
     1130
    9521131
    9531132#ifdef CPP_StratAer
    … …  
    10881267       ENDIF  !  iflag_lscav . EQ. 1, 2, 3 or 4
    10891268       !
     1269    !write(lunout,*) 'THATSALLFOLKS'
    10901270    ENDIF !  lessivage
     1271
     1272    !WRITE(lunout,*) 'CHECK4 TR_SERI AIBCM',SUM(tr_seri(:,:,id_AIBCM))
     1273    !WRITE(lunout,*) 'CHECK4 TR_SERI ASBCM',SUM(tr_seri(:,:,id_ASBCM))
     1274
    10911275
    10921276
    … …  
    11041288            tau_aero, piz_aero, cg_aero,        ccm,       &
    11051289            rfname,                                        &
    1106             tr_seri,  source)     
     1290            tr_seri,  source,beta_st)     
    11071291       
    11081292       
    11091293    ENDIF
     1294
     1295    !WRITE(lunout,*) 'CHECK5 TR_SERI AIBCM',SUM(tr_seri(:,:,id_AIBCM))
     1296    !WRITE(lunout,*) 'CHECK5 TR_SERI ASBCM',SUM(tr_seri(:,:,id_ASBCM))
     1297
     1298
     1299    !write(lunout,*) 'MERCI INCA POUR LES SERVICES'
    11101300    !=============================================================
    11111301    !   Ecriture des sorties

  • LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/radlwsw_m.F90

    r3666 r4795  
    873873         ZTOPLWAIAERO,ZSOLLWAIAERO, &
    874874         ZLWADAERO, & !--NL
    875          ok_ade, ok_aie, ok_volcan, flag_aerosol,flag_aerosol_strat, flag_aer_feedback) ! flags aerosols
     875         ok_ade, ok_aie, ok_volcan, flag_aerosol,flag_aerosol_strat, .FALSE.) !Nicolas Février
    876876           
    877877!        print *,'RADLWSW: apres RECMWF'

  • LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/tracinca_mod.F90

    r2784 r4795  
    4242       tau_aero, piz_aero, cg_aero,        ccm,       &
    4343       rfname,                                        &
    44        tr_seri,  source)     
     44       tr_seri,  source,beta_st)     
    4545
    4646!========================================================
    … …  
    111111    REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf ! Pourcentage de sol f(nature du sol)
    112112
     113! Nicolas Février
     114! Lessivage de LMDZ à envoyer vers INCA
     115    REAL,DIMENSION(klon,klev,12),INTENT(IN) :: beta_st
    113116
    114117  ! InOutput argument
    … …  
    145148    ENDIF
    146149
     150    !PRINT*, 'IL EST PASSE PAR ICI'
     151
    147152    CALL aerosolmain(                    &
    148153         aerosol_couple,tr_seri,pdtphys, &
    … …  
    153158         rfname,ccm,lafin)
    154159#endif
     160
     161    !PRINT*, 'IL REPASSERA PAR LA'
     162
     163! Debug beta_st
     164    print *,"BETA TRACINCA",beta_st(5000,35,7)
    155165
    156166
    … …  
    190200         nbp_lon,    & !nx
    191201         nbp_lat,    & !ny
    192          source )
     202         source,     &
     203         beta_st)
     204
    193205#endif
    194206   
     207    !PRINT*, 'AVEC DE LA CHANCE PEUT-ETRE ICI AUSSI'
     208
    195209    CALL VTe(VTinca)
    196210    CALL VTb(VTphysiq)
Note: See TracChangeset for help on using the changeset viewer.