Ignore:
Timestamp:
Apr 13, 2015, 10:21:09 AM (10 years ago)
Author:
Laurent Fairhead
Message:

Merged trunk changes 2216:2237 into testing branch

Location:
LMDZ5/branches/testing
Files:
1 deleted
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/branches/testing

  • LMDZ5/branches/testing/libf/dyn3d/calfis.F

    r2056 r2258  
    2121     $                  pdq,
    2222     $                  flxw,
    23      $                  clesphy0,
    2423     $                  pdufi,
    2524     $                  pdvfi,
    … …  
    131130      REAL,INTENT(OUT) :: pdpsfi(iip1,jjp1) ! tendency on surface pressure (Pa/s)
    132131
    133       INTEGER,PARAMETER :: longcles = 20
    134       REAL,INTENT(IN) :: clesphy0( longcles ) ! unused
    135 
    136132
    137133c    Local variables :
    … …  
    465461     .             zphis,
    466462     .             presnivs,
    467      .             clesphy0,
    468463     .             zufi,
    469464     .             zvfi,

  • LMDZ5/branches/testing/libf/dyn3d/ce0l.F90

    r1999 r2258  
    4141#include "temps.h"
    4242#include "logic.h"
    43   INTEGER, PARAMETER            :: longcles=20
    44   REAL,    DIMENSION(longcles)  :: clesphy0
    4543  REAL,    DIMENSION(iip1,jjp1) :: masque
    4644  CHARACTER(LEN=15)             :: calnd
    4745  REAL,    DIMENSION(iip1,jjp1) :: phis ! geopotentiel au sol
    4846!-------------------------------------------------------------------------------
    49   CALL conf_gcm( 99, .TRUE. , clesphy0 )
     47  CALL conf_gcm( 99, .TRUE. )
    5048
    5149  CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))

  • LMDZ5/branches/testing/libf/dyn3d/conf_gcm.F90

    r2160 r2258  
    22! $Id$
    33
    4 SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
     4SUBROUTINE conf_gcm( tapedef, etatinit )
    55
    66  USE control_mod
    … …  
    2323  !     etatinit  :     = TRUE   , on ne  compare pas les valeurs des para-
    2424  !     -metres  du zoom  avec  celles lues sur le fichier start .
    25   !      clesphy0 :  sortie  .
    26 
    27   LOGICAL etatinit
    28   INTEGER tapedef
    29 
    30   INTEGER        longcles
    31   PARAMETER(     longcles = 20 )
    32   REAL clesphy0( longcles )
     25
     26  LOGICAL,INTENT(IN) :: etatinit
     27  INTEGER,INTENT(IN) :: tapedef
    3328
    3429  !   Declarations :
    … …  
    4136  include "temps.h"
    4237  include "comconst.h"
    43 
    44   ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
    45   ! include "clesphys.h"
    4638  include "iniprint.h"
    4739
    … …  
    381373  ip_ebil_dyn = 0
    382374  CALL getin('ip_ebil_dyn',ip_ebil_dyn)
    383 
    384   DO i = 1, longcles
    385      clesphy0(i) = 0.
    386   ENDDO
    387375
    388376  !cc  ....   P. Le Van , ajout  le 7/03/95 .pour le zoom ...
    … …  
    734722     !Config  Help = extension en longitude  de la zone du zoom 
    735723     !Config         ( fraction de la zone totale)
    736      dzoomx = 0.0
     724     dzoomx = 0.2
    737725     CALL getin('dzoomx',dzoomx)
     726     call assert(dzoomx < 1, "conf_gcm: dzoomx must be < 1")
    738727
    739728     !Config  Key  = dzoomy
    … …  
    742731     !Config  Help = extension en latitude de la zone  du zoom 
    743732     !Config         ( fraction de la zone totale)
    744      dzoomy = 0.0
     733     dzoomy = 0.2
    745734     CALL getin('dzoomy',dzoomy)
     735     call assert(dzoomy< 1, "conf_gcm: dzoomy must be < 1")
    746736
    747737     !Config  Key  = taux

  • LMDZ5/branches/testing/libf/dyn3d/gcm.F

    r2160 r2258  
    2626! Only INCA needs these informations (from the Earth's physics)
    2727      USE indice_sol_mod
     28      USE mod_phys_lmdz_para, ONLY : klon_mpi_para_nb
    2829#endif
    2930
    … …  
    3334! dynamique -> physique pour l'initialisation
    3435#ifdef CPP_PHYS
    35       USE dimphy
    36       USE comgeomphy
    37       USE mod_phys_lmdz_para, ONLY : klon_mpi_para_nb
     36!      USE dimphy
     37!      USE comgeomphy
    3838#endif
    3939!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    9090!#include "indicesol.h"
    9191#endif
    92       INTEGER         longcles
    93       PARAMETER     ( longcles = 20 )
    94       REAL  clesphy0( longcles )
    95       SAVE  clesphy0
    96 
    97 
    9892
    9993      REAL zdtvr
    … …  
    125119      LOGICAL first
    126120
    127       LOGICAL call_iniphys
    128       data call_iniphys/.true./
     121!      LOGICAL call_iniphys
     122!      data call_iniphys/.true./
    129123
    130124c+jld variables test conservation energie
    … …  
    147141      REAL :: heure
    148142
    149 
    150 c-----------------------------------------------------------------------
    151 c    variables pour l'initialisation de la physique :
    152 c    ------------------------------------------------
    153       INTEGER ngridmx
    154       PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm   )
    155       REAL zcufi(ngridmx),zcvfi(ngridmx)
    156       REAL latfi(ngridmx),lonfi(ngridmx)
    157       REAL airefi(ngridmx)
    158       SAVE latfi, lonfi, airefi
    159 
    160143c-----------------------------------------------------------------------
    161144c   Initialisations:
    … …  
    175158c  ---------------------------------------
    176159c
    177 ! Ehouarn: dump possibility of using defrun
    178 !#ifdef CPP_IOIPSL
    179       CALL conf_gcm( 99, .TRUE. , clesphy0 )
     160      CALL conf_gcm( 99, .TRUE.)
    180161      if (mod(iphysiq, iperiod) /= 0) call abort_gcm("conf_gcm",
    181162     s "iphysiq must be a multiple of iperiod", 1)
    182 !#else
    183 !      CALL defrun( 99, .TRUE. , clesphy0 )
    184 !#endif
    185163
    186164!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    199177#ifdef CPP_PHYS
    200178      CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
    201       call InitComgeomphy
     179!      call InitComgeomphy ! now done in iniphysiq
    202180#endif
    203181!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    215193        call ioconf_calendar('noleap')
    216194        write(lunout,*)'CALENDRIER CHOISI: Terrestre a 365 jours/an'
    217       else if (calend == 'earth_366d') then
     195      else if (calend == 'gregorian') then
    218196        call ioconf_calendar('gregorian')
    219197        write(lunout,*)'CALENDRIER CHOISI: Terrestre bissextile'
    … …  
    316294
    317295C
    318 C on remet le calendrier à zero si demande
     296C on remet le calendrier \`a zero si demande
    319297c
    320298      IF (start_time /= starttime) then
    321299        WRITE(lunout,*)' GCM: Attention l''heure de depart lue dans le'
    322      &,' fichier restart ne correspond pas à celle lue dans le run.def'
     300     &,' fichier restart ne correspond pas a celle lue dans le run.def'
    323301        IF (raz_date == 1) then
    324302          WRITE(lunout,*)'Je prends l''heure lue dans run.def'
    … …  
    428406c   -------------------------------
    429407
    430       IF (call_iniphys.and.(iflag_phys==1.or.iflag_phys>=100)) THEN
    431          latfi(1)=rlatu(1)
    432          lonfi(1)=0.
    433          zcufi(1) = cu(1)
    434          zcvfi(1) = cv(1)
    435          DO j=2,jjm
    436             DO i=1,iim
    437                latfi((j-2)*iim+1+i)= rlatu(j)
    438                lonfi((j-2)*iim+1+i)= rlonv(i)
    439                zcufi((j-2)*iim+1+i) = cu((j-1)*iip1+i)
    440                zcvfi((j-2)*iim+1+i) = cv((j-1)*iip1+i)
    441             ENDDO
    442          ENDDO
    443          latfi(ngridmx)= rlatu(jjp1)
    444          lonfi(ngridmx)= 0.
    445          zcufi(ngridmx) = cu(ip1jm+1)
    446          zcvfi(ngridmx) = cv(ip1jm-iim)
    447          CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)
    448          WRITE(lunout,*)
    449      .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
     408      IF ((iflag_phys==1).or.(iflag_phys>=100)) THEN
    450409! Physics:
    451410#ifdef CPP_PHYS
    452          CALL iniphysiq(ngridmx,llm,daysec,day_ini,dtphys/nsplit_phys,
    453      &                latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp,
     411         CALL iniphysiq(iim,jjm,llm,daysec,day_ini,dtphys/nsplit_phys,
     412     &                rlatu,rlonv,aire,cu,cv,rad,g,r,cpp,
    454413     &                iflag_phys)
    455414#endif
    456          call_iniphys=.false.
    457       ENDIF ! of IF (call_iniphys.and.(iflag_phys.eq.1))
     415      ENDIF ! of IF ((iflag_phys==1).or.(iflag_phys>=100))
    458416
    459417c  numero de stockage pour les fichiers de redemarrage:
    … …  
    530488
    531489
    532       CALL leapfrog(ucov,vcov,teta,ps,masse,phis,q,clesphy0,
    533      .              time_0)
     490      CALL leapfrog(ucov,vcov,teta,ps,masse,phis,q,time_0)
    534491
    535492      END

  • LMDZ5/branches/testing/libf/dyn3d/leapfrog.F

    r2056 r2258  
    44c
    55c
    6       SUBROUTINE leapfrog(ucov,vcov,teta,ps,masse,phis,q,clesphy0,
    7      &                    time_0)
     6      SUBROUTINE leapfrog(ucov,vcov,teta,ps,masse,phis,q,time_0)
    87
    98
    … …  
    7069#include "academic.h"
    7170
    72 ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
    73 ! #include "clesphys.h"
    74 
    75       INTEGER,PARAMETER :: longcles = 20
    76       REAL,INTENT(IN) :: clesphy0( longcles ) ! not used
    7771      REAL,INTENT(IN) :: time_0 ! not used
    7872
    … …  
    446440     $               ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
    447441     $               du,dv,dteta,dq,
    448      $               flxw,
    449      $               clesphy0, dufi,dvfi,dtetafi,dqfi,dpfi  )
     442     $               flxw,dufi,dvfi,dtetafi,dqfi,dpfi  )
    450443
    451444c      ajout des tendances physiques:
Note: See TracChangeset for help on using the changeset viewer.