Changeset 1576


Ignore:
Timestamp:
Jul 27, 2016, 11:14:41 AM (8 years ago)
Author:
emillour
Message:

All models:
Clean up the dynamics/physics interface to converge with LMDZ5;
get rid of tracerdyn parameter (which was supposed to send information
about wether tracers should be advected or not in the Mars and Generic
models).
EM

Location:
trunk
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • trunk/DOC/chantiers/commit_importants.log

    r1575 r1576  
    18591859- grid_atob_m.F90: fix for some zoomed grid interpolation cases
    18601860
     1861**********************
     1862**** commit_v1576 ****
     1863**********************
     1864Ehouarn: Clean up the dynamics/physics interface to converge with LMDZ5;
     1865get rid of tracerdyn parameter (which was supposed to send information
     1866about wether tracers should be advected or not in the Mars and Generic
     1867models).
     1868

  • trunk/LMDZ.COMMON/libf/dynphy_lonlat/calfis.F

    r1549 r1576  
    183183      LOGICAL,SAVE :: firstcal=.true., debut=.true.
    184184!      REAL rdayvrai
    185 
    186       LOGICAL tracerdyn ! for generic/mars physics call ; possibly to get rid of
    187185
    188186! For Titan only right now:
    … …  
    707705     &                   zufi,zvfi,zrfi,ztfi,zqfi,
    708706     &                   flxwfi,pducov,
    709      &                   zdufi,zdvfi,zdtfi,zdqfi,zdpsrf,
    710      &                   tracerdyn)
     707     &                   zdufi,zdvfi,zdtfi,zdqfi,zdpsrf)
    711708
    712709!      if (planet_type.eq."earth") then

  • trunk/LMDZ.COMMON/libf/dynphy_lonlat/calfis_p.F

    r1572 r1576  
    253253      INTEGER :: offset 
    254254      INTEGER :: jjb,jje
    255 
    256       LOGICAL tracerdyn ! for generic/mars physics call ; possibly to get rid of
    257255
    258256! For Titan only right now:
    … …  
    945943     &                   flxwfi_omp,pducov,
    946944     &                   zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,
    947      &                   zdpsrf_omp,tracerdyn)
     945     &                   zdpsrf_omp)
    948946
    949947!      if (planet_type=="earth") then

  • trunk/LMDZ.COMMON/libf/dynphy_lonlat/phytitan/callphysiq_mod.F90

    r1549 r1576  
    1616                       zufi_omp,zvfi_omp,zrfi_omp,ztfi_omp,zqfi_omp,      &
    1717                       flxwfi_omp,pducov,                                 &
    18                        zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp,&
    19                        tracerdyn)
     18                       zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp)
    2019
    2120  USE control_mod, ONLY: planet_type
    … …  
    5251  REAL,INTENT(OUT) :: zdqfi_omp(klon,llm,nqtot) ! tendency on tracers
    5352  REAL,INTENT(OUT) :: zdpsrf_omp(klon) ! tendency on surface pressure
    54   LOGICAL,INTENT(OUT) :: tracerdyn
    5553
    5654  ! Local variables

  • trunk/LMDZ.COMMON/libf/dynphy_lonlat/phyvenus/callphysiq_mod.F90

    r1549 r1576  
    1616                       zufi_omp,zvfi_omp,zrfi_omp,ztfi_omp,zqfi_omp,      &
    1717                       flxwfi_omp,pducov,                                 &
    18                        zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp,&
    19                        tracerdyn)
     18                       zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp)
    2019
    2120  USE control_mod, ONLY: planet_type
    … …  
    5251  REAL,INTENT(OUT) :: zdqfi_omp(klon,llm,nqtot) ! tendency on tracers
    5352  REAL,INTENT(OUT) :: zdpsrf_omp(klon) ! tendency on surface pressure
    54   LOGICAL,INTENT(OUT) :: tracerdyn
    5553
    5654  ! Local variables

  • trunk/LMDZ.GENERIC/README

    r1573 r1576  
    12381238package (directory phystd) must be done in the related phystd/iniphysiq
    12391239routine.
     1240
     1241== 27/07/2016 == EM
     1242- further cleanup in the dynamics/physics interface: stop sending information
     1243  about tracer advection (tracerdyn) back from physics to dynamics as this
     1244  input parameter can be read and set in the dynamics.

  • trunk/LMDZ.GENERIC/libf/dyn3d/gcm.F

    r1543 r1576  
    11      PROGRAM gcm
    22
     3      use ioipsl_getincom, only: getin
    34      use infotrac, only: iniadvtrac, nqtot, iadv
    45      use sponge_mod,only: callsponge,mode_sponge,sponge
    … …  
    129130      character*80 abort_message
    130131
    131       LOGICAL tracer
    132           data tracer/.true./
    133 !      INTEGER nq
    134 
     132! flag to run with or without tracer transport (read from run.def)
     133      logical tracer
     134     
    135135C Calendrier
    136136      LOGICAL true_calendar
    … …  
    177177c-----------------------------------------------------------------------
    178178      CALL defrun_new( 99, .TRUE. )
     179      WRITE(*,*) ""
     180      WRITE(*,*) "Run with or without tracer transport ?"
     181      tracer=.true. ! default value
     182      call getin("tracer",tracer)
     183      WRITE(*,*)" tracer = ",tracer
    179184
    180185! Initialize tracers
    … …  
    479484        CALL calfis( nqtot, lafin ,rdayvrai,rday_ecri,time  ,
    480485     $                 ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
    481      $     du,dv,dteta,dq,w, dufi,dvfi,dhfi,dqfi,dpfi,tracer)
     486     $     du,dv,dteta,dq,w, dufi,dvfi,dhfi,dqfi,dpfi)
    482487
    483488

  • trunk/LMDZ.GENERIC/libf/dynphy_lonlat/calfis.F

    r1549 r1576  
    22     $            pucov,pvcov,pteta,pq,pmasse,pps,pp,ppk,pphis,pphi,
    33     $            pducov,pdvcov,pdteta,pdq,pw,
    4      $            pdufi,pdvfi,pdhfi,pdqfi,pdpsfi,tracer )
     4     $            pdufi,pdvfi,pdhfi,pdqfi,pdpsfi )
    55c
    66c    Auteur :  P. Le Van, F. Hourdin
    … …  
    5656c        pdtsfi         tendency for the surface temperature
    5757c
    58 c        tracer         Call tracer in  gcm.F ? (decided in callphys.def)
    59 c
    6058c=======================================================================
    6159c
    … …  
    106104      REAL pdqfi(iip1,jjp1,llm,nqtot)
    107105      REAL pdpsfi(iip1,jjp1)
    108       logical tracer
    109106
    110107c    Local variables :
    … …  
    441438     ,     flxwfi,
    442439C - sorties
    443      s     zdufi, zdvfi, zdtfi, zdqfi,zdpsrf,tracer)
     440     s     zdufi, zdvfi, zdtfi, zdqfi,zdpsrf)
    444441
    445442

  • trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/callphysiq_mod.F90

    r1549 r1576  
    1616                       zufi_omp,zvfi_omp,zrfi_omp,ztfi_omp,zqfi_omp,      &
    1717                       flxwfi_omp,pducov,                                 &
    18                        zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp,&
    19                        tracerdyn)
     18                       zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp)
    2019
    2120  USE control_mod, ONLY: planet_type
    … …  
    5251  REAL,INTENT(OUT) :: zdqfi_omp(klon,llm,nqtot) ! tendency on tracers
    5352  REAL,INTENT(OUT) :: zdpsrf_omp(klon) ! tendency on surface pressure
    54   LOGICAL,INTENT(OUT) :: tracerdyn
    5553
    5654  ! Local variables
    … …  
    9088              zdtfi_omp,      & ! pdt
    9189              zdqfi_omp,      & ! pdq
    92               zdpsrf_omp,     & ! pdpsrf
    93               tracerdyn)        ! tracerdyn (somewhat obsolete)
     90              zdpsrf_omp)       ! pdpsrf
    9491
    9592

  • trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

    r1549 r1576  
    108108      INTEGER ierr, aslun
    109109      REAL tmp1(0:llm),tmp2(0:llm)
    110       Logical  tracerdyn
    111110      integer :: nq !=1 ! number of tracers
    112111 
    … …  
    906905     ,     w,
    907906C - sorties
    908      s     du, dv, dtemp, dq,dpsurf,tracerdyn)
     907     s     du, dv, dtemp, dq,dpsurf)
    909908
    910909

  • trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

    r1549 r1576  
    1212                  pu,pv,pt,pq,             &
    1313                  flxw,                    &
    14                   pdu,pdv,pdt,pdq,pdpsrf,tracerdyn)
     14                  pdu,pdv,pdt,pdq,pdpsrf)
    1515 
    1616      use radinc_h, only : L_NSPECTI,L_NSPECTV,naerkind
    … …  
    132132!    pdq(ngrid,nlayer)        /
    133133!    pdpsrf(ngrid)             /
    134 !    tracerdyn                 call tracer in dynamical part of GCM ?
    135134!
    136135!
    … …  
    195194      real,intent(out) :: pdq(ngrid,nlayer,nq) ! Tracer tendencies (kg/kg_of_air/s).
    196195      real,intent(out) :: pdpsrf(ngrid)        ! Surface pressure tendency (Pa/s).
    197       logical,intent(out) :: tracerdyn         ! Signal to the dynamics to advect tracers or not.     
    198196
    199197! Local saved variables:
    … …  
    517515!        Initialize tracer names, indexes and properties.
    518516!        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    519          tracerdyn=tracer
    520517         IF (.NOT.ALLOCATED(noms)) ALLOCATE(noms(nq)) ! (because noms is an argument of physdem1 whether or not tracer is on)
    521518         if (tracer) then

  • trunk/LMDZ.MARS/README

    r1573 r1576  
    23122312package (directory phymars) must be done in the related phymars/iniphysiq
    23132313routine.
     2314
     2315== 27/07/2016 == EM
     2316- further cleanup in the dynamics/physics interface: stop sending information
     2317  about tracer advection (tracerdyn) back from physics to dynamics as this
     2318  input parameter can be read and set in the dynamics.

  • trunk/LMDZ.MARS/libf/dyn3d/gcm.F

    r1543 r1576  
    11      PROGRAM gcm
    22
     3      use ioipsl_getincom, only: getin
    34      use infotrac, only: iniadvtrac, nqtot, iadv
    45      use control_mod, only: day_step, iperiod, iphysiq, ndynstep,
    … …  
    129130      character*80 abort_message
    130131
    131       LOGICAL tracer
    132           data tracer/.true./
    133 !      INTEGER nq
    134 
     132! flag to run with or without tracer transport (read from run.def)
     133      logical tracer
     134     
    135135C Calendrier
    136136      LOGICAL true_calendar
    … …  
    161161c-----------------------------------------------------------------------
    162162      CALL defrun_new( 99, .TRUE. )
     163      WRITE(*,*) ""
     164      WRITE(*,*) "Run with or without tracer transport ?"
     165      tracer=.true. ! default value
     166      call getin("tracer",tracer)
     167      WRITE(*,*)" tracer = ",tracer
    163168
    164169! Initialize tracers
    … …  
    491496        CALL calfis( nqtot, lafin ,rdayvrai,rday_ecri,time  ,
    492497     $                 ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
    493      $     du,dv,dteta,dq,w, dufi,dvfi,dhfi,dqfi,dpfi,tracer)
     498     $     du,dv,dteta,dq,w, dufi,dvfi,dhfi,dqfi,dpfi)
    494499
    495500

  • trunk/LMDZ.MARS/libf/dynphy_lonlat/calfis.F

    r1549 r1576  
    22     $            pucov,pvcov,pteta,pq,pmasse,pps,pp,ppk,pphis,pphi,
    33     $            pducov,pdvcov,pdteta,pdq,pw,
    4      $            pdufi,pdvfi,pdhfi,pdqfi,pdpsfi,tracer )
     4     $            pdufi,pdvfi,pdhfi,pdqfi,pdpsfi )
    55c
    66c    Auteur :  P. Le Van, F. Hourdin
    … …  
    5656c        pdtsfi         tendency for the surface temperature
    5757c
    58 c        tracer         Call tracer in  gcm.F ? (decided in callphys.def)
    59 c
    6058c=======================================================================
    6159c
    … …  
    103101      REAL pdqfi(iip1,jjp1,llm,nq)
    104102      REAL pdpsfi(iip1,jjp1)
    105       logical tracer
    106103
    107104c    Local variables :
    … …  
    437434     ,     flxwfi,
    438435C - sorties
    439      s     zdufi, zdvfi, zdtfi, zdqfi,zdpsrf,tracer)
     436     s     zdufi, zdvfi, zdtfi, zdqfi,zdpsrf)
    440437
    441438

  • trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/callphysiq_mod.F90

    r1549 r1576  
    1616                       zufi_omp,zvfi_omp,zrfi_omp,ztfi_omp,zqfi_omp,      &
    1717                       flxwfi_omp,pducov,                                 &
    18                        zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp,&
    19                        tracerdyn)
     18                       zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp)
    2019
    2120  USE control_mod, ONLY: planet_type
    … …  
    5251  REAL,INTENT(OUT) :: zdqfi_omp(klon,llm,nqtot) ! tendency on tracers
    5352  REAL,INTENT(OUT) :: zdpsrf_omp(klon) ! tendency on surface pressure
    54   LOGICAL,INTENT(OUT) :: tracerdyn
    5553
    5654  ! Local variables
    … …  
    8987              zdtfi_omp,      & ! pdt
    9088              zdqfi_omp,      & ! pdq
    91               zdpsrf_omp,     & ! pdpsrf
    92               tracerdyn)        ! tracerdyn (somewhat obsolete)
     89              zdpsrf_omp)       ! pdpsrf
    9390
    9491

  • trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

    r1550 r1576  
    109109      INTEGER ierr, aslun
    110110      REAL tmp1(0:nlayer),tmp2(0:nlayer)
    111       Logical  tracerdyn
    112111      integer :: nq=1 ! number of tracers
    113112      real :: latitude(1), longitude(1), cell_area(1)
    … …  
    753752     ,     w,
    754753C - outputs
    755      s     du, dv, dtemp, dq,dpsurf,tracerdyn)
     754     s     du, dv, dtemp, dq,dpsurf)
    756755!      write(*,*) "testphys1d apres q", q(1,:)
    757756

  • trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

    r1549 r1576  
    1212     $            ,pu,pv,pt,pq
    1313     $            ,flxw
    14      $            ,pdu,pdv,pdt,pdq,pdpsrf,tracerdyn)
     14     $            ,pdu,pdv,pdt,pdq,pdpsrf)
    1515
    1616      use tracer_mod, only: noms, mmol, igcm_co2, igcm_n2,
    … …  
    146146c    pdq(ngrid,nlayer,nq)    |
    147147c    pdpsrf(ngrid)           |
    148 c    tracerdyn                 call tracer in dynamical part of GCM ?
    149148
    150149c
    … …  
    191190      REAL,INTENT(out) :: pdq(ngrid,nlayer,nq) ! tracer tendencies (../kg/s)
    192191      REAL,INTENT(out) :: pdpsrf(ngrid) ! surface pressure tendency (Pa/s)
    193       LOGICAL,INTENT(out) :: tracerdyn ! signal to the dynamics to advect tracers or not
    194192
    195193
    … …  
    434432c        initialize tracers
    435433c        ~~~~~~~~~~~~~~~~~~
    436          tracerdyn=tracer
    437434         IF (tracer) THEN
    438435            CALL initracer(ngrid,nq,qsurf)
Note: See TracChangeset for help on using the changeset viewer.