Changeset 101 for trunk


Ignore:
Timestamp:
Mar 23, 2011, 4:23:31 PM (14 years ago)
Author:
slebonnois
Message:

SL: modifications pour arriver a compiler le gcm VENUS !
Ca marche !
A noter: modifs de makelmdz

Location:
trunk
Files:
4 added
58 edited

Legend:

Unmodified
Added
Removed

  • trunk/libf/dyn3d/ce0l.F90

    r97 r101  
    2727#endif
    2828  IMPLICIT NONE
     29#include "iniprint.h"
    2930#ifndef CPP_EARTH
    3031  WRITE(lunout,*)'limit_netcdf: Earth-specific routine, needs Earth physics'
    … …  
    3637#include "paramet.h"
    3738#include "indicesol.h"
    38 #include "iniprint.h"
    3939#include "temps.h"
    4040#include "logic.h"

  • trunk/libf/dyn3d/gcm.F

    r97 r101  
    2222! dynamique -> physique pour l'initialisation
    2323! Ehouarn: for now these only apply to Earth:
    24 #ifdef CPP_EARTH
     24#ifdef CPP_PHYS
    2525      USE dimphy
    2626      USE comgeomphy
     27#endif
     28#ifdef CPP_EARTH
    2729      USE mod_phys_lmdz_para, ONLY : klon_mpi_para_nb
    2830#endif
    … …  
    180182!#ifdef CPP_EARTH
    181183#ifdef CPP_PHYS
    182       CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
    183       call InitComgeomphy
     184      CALL init_phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
     185      call initcomgeomphy
    184186#endif
    185187!      endif
    … …  
    482484      else
    483485! A voir pour Titan et Venus
    484       write(lunout,*)'A VOIR POUR VENUS ET TITAN: separation en annees...'
     486      write(lunout,*)'A VOIR POUR VENUS/TITAN: separation en annees...'
    485487      endif ! planet_type
    486488

  • trunk/libf/dyn3d/leapfrog.F

    r97 r101  
    196196      real,external :: cpdet
    197197      real :: ztetaec(ip1jmp1,llm)
     198
     199c dummy: sinon cette routine n'est jamais compilee...
     200      if(1.eq.0) then
     201        CALL init_phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
     202      endif
    198203
    199204      itaufin   = nday*day_step

  • trunk/libf/dyn3d/logic.h

    r1 r101  
    1010     &  statcl,conser,apdiss,apdelq,saison,ecripar,fxyhypb,ysinus       &
    1111     &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
    12      &  ,ok_limit,ok_etat0
     12     &  ,ok_limit,ok_etat0,iflag_trac
    1313
    1414      LOGICAL purmats,forward,leapf,apphys,statcl,conser,               &
    1515     & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus                      &
    1616     &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
    17      &  ,ok_limit,ok_etat0
     17     &  ,ok_limit,ok_etat0,iflag_trac
    1818
    1919      INTEGER iflag_phys

  • trunk/libf/dyn3d/sortvarc.F

    r97 r101  
    55     $(itau,ucov,teta,ps,masse,pk,phis,vorpot,phi,bern,dp,time ,
    66     $ vcov )
     7
     8      use control_mod,only:resetvarc
    79      IMPLICIT NONE
    810

  • trunk/libf/dyn3dpar/ce0l.F90

    r97 r101  
    3030#endif
    3131  IMPLICIT NONE
     32#include "iniprint.h"
    3233#ifndef CPP_EARTH
    3334  WRITE(lunout,*)'limit_netcdf: Earth-specific routine, needs Earth physics'
    … …  
    3940#include "paramet.h"
    4041#include "indicesol.h"
    41 #include "iniprint.h"
    4242#include "temps.h"
    4343#include "logic.h"

  • trunk/libf/dyn3dpar/gcm.F

    r97 r101  
    188188!#ifdef CPP_EARTH
    189189#ifdef CPP_PHYS
    190         CALL Init_Phys_lmdz(iim,jjp1,llm,mpi_size,distrib_phys)
     190        CALL init_phys_lmdz(iim,jjp1,llm,mpi_size,distrib_phys)
    191191#endif
    192192!      endif ! of if (planet_type.eq."earth")
    … …  
    210210#ifdef CPP_PHYS
    211211c$OMP PARALLEL
    212       call InitComgeomphy
     212      call initcomgeomphy
    213213c$OMP END PARALLEL
    214214#endif
    … …  
    513513      else
    514514! A voir pour Titan et Venus
    515       write(lunout,*)'A VOIR POUR VENUS ET TITAN: separation en annees...'
     515      write(lunout,*)'A VOIR POUR VENUS/TITAN: separation en annees...'
    516516      endif ! planet_type
    517517

  • trunk/libf/dyn3dpar/leapfrog_p.F

    r97 r101  
    203203      LOGICAL :: ok_start_timer=.FALSE.
    204204      LOGICAL, SAVE :: firstcall=.TRUE.
     205
     206c dummy: sinon cette routine n'est jamais compilee...
     207      if(1.eq.0) then
     208        CALL init_phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
     209      endif
    205210
    206211c$OMP MASTER

  • trunk/libf/dyn3dpar/logic.h

    r1 r101  
    1010     &  statcl,conser,apdiss,apdelq,saison,ecripar,fxyhypb,ysinus       &
    1111     &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
    12      &  ,ok_limit,ok_etat0
     12     &  ,ok_limit,ok_etat0,iflag_trac
    1313
    1414      LOGICAL purmats,forward,leapf,apphys,statcl,conser,               &
    1515     & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus                      &
    1616     &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
    17      &  ,ok_limit,ok_etat0
     17     &  ,ok_limit,ok_etat0,iflag_trac
    1818
    1919      INTEGER iflag_phys

  • trunk/libf/dyn3dpar/sortvarc.F

    r97 r101  
    55     $(itau,ucov,teta,ps,masse,pk,phis,vorpot,phi,bern,dp,time ,
    66     $ vcov )
     7
     8      use control_mod,only:resetvarc
    79      IMPLICIT NONE
    810

  • trunk/libf/phyvenus/YOMCST.h

    r97 r101  
    1010      REAL RDAY,REA,REPSM,RSIYEA,RSIDAY,ROMEGA
    1111! A1.1.bis Constantes concernant l'orbite de la Terre:
    12       REAL :: R_ecc, R_peri, R_incl
     12      REAL R_ecc, R_peri, R_incl
    1313! A1.2 Geoide
    1414      REAL RA,RG,R1SA
    1515! A1.3 Radiation
    16 !     REAL :: RSIGMA,RI0
     16!     REAL RSIGMA,RI0
    1717      REAL RSIGMA
    1818! A1.4 Thermodynamic gas phase
    … …  
    3131!
    3232      COMMON/YOMCST/RPI ,RCLUM, RHPLA, RKBOL, RNAVO ,RDAY  ,REA         &
    33      & ,REPSM ,RSIYEA,RSIDAY,ROMEGA , R_ecc, R_peri, R_incl,            &
     33     & ,REPSM ,RSIYEA,RSIDAY,ROMEGA , R_ecc, R_peri, R_incl             &
    3434     & ,RA    ,RG ,R1SA                                                 &
    3535     & ,RSIGMA,R ,RMD   ,RMV   ,RD    ,RV    ,RCPD ,RCPV,RCVD           &

  • trunk/libf/phyvenus/aaam_bud.F

    r3 r101  
    66     i                   p, u, v)
    77c
     8      use dimphy
    89      implicit none
    910c======================================================================
    … …  
    7879
    7980#include "dimensions.h"
    80 #include "dimphy.h"
    8181c
    8282c ARGUMENTS

  • trunk/libf/phyvenus/ajsec.F

    r3 r101  
    55      SUBROUTINE ajsec(paprs, pplay, ppk, tfi, ufi, vfi, nq, qfi,
    66     .                             d_tfi, d_ufi, d_vfi, d_qfi)
     7
     8      use dimphy
    79      IMPLICIT none
    810c======================================================================
    … …  
    2527c======================================================================
    2628#include "dimensions.h"
    27 #include "dimphy.h"
    2829#include "YOMCST.h"
    2930      REAL paprs(klon,klev+1), pplay(klon,klev)
    … …  
    3536      REAL qfi(klon,klev,nq), d_qfi(klon,klev,nq)
    3637c
    37       INTEGER limbas, limhau ! les couches a ajuster
    38 ccc      PARAMETER (limbas=klev-3, limhau=klev)
    39       PARAMETER (limbas=1, limhau=klev)
     38      INTEGER,save :: limbas, limhau ! les couches a ajuster
    4039c
    4140      REAL zh(klon,klev)
    … …  
    5049c Initialisation:
    5150c
     51      limbas=1
     52      limhau=klev
     53
    5254      DO k = 1, klev
    5355      DO i = 1, klon

  • trunk/libf/phyvenus/ballon.F

    r3 r101  
    22     i                   temp, p, u, v, geop)
    33
     4      use dimphy
    45      implicit none
    56
    … …  
    9495
    9596#include "dimensions.h"
    96 #include "dimphy.h"
    9797#include "YOMCST.h"
    9898c
    … …  
    256256     .                       phi,lam,ubal,vbal,abal)
    257257
     258      use dimphy
    258259      implicit none
    259260
    … …  
    284285
    285286#include "dimensions.h"
    286 #include "dimphy.h"
    287287#include "YOMCST.h"
    288288c
    … …  
    355355      subroutine pos_implem(phi,lam,ubal,vbal,dt)
    356356
     357      use dimphy
    357358      implicit none
    358359
    … …  
    371372
    372373#include "dimensions.h"
    373 #include "dimphy.h"
    374374#include "YOMCST.h"
    375375c

  • trunk/libf/phyvenus/clmain.F

    r97 r101  
    3434      USE ioipsl
    3535      USE interface_surf
     36      use dimphy
    3637      IMPLICIT none
    3738c======================================================================
    … …  
    6566c======================================================================
    6667#include "dimensions.h"
    67 #include "dimphy.h"
    6868c$$$ PB ajout pour soil
    6969#include "dimsoil.h"
    … …  
    488488
    489489      USE interface_surf
     490      use dimphy
    490491
    491492      IMPLICIT none
    … …  
    495496c======================================================================
    496497#include "dimensions.h"
    497 #include "dimphy.h"
    498498#include "YOMCST.h"
    499499#include "dimsoil.h"
    … …  
    766766     e                  paprs,pplay,delp,
    767767     s                  d_ven,flux_v)
     768
     769      use dimphy
    768770      IMPLICIT none
    769771c======================================================================
    … …  
    789791c======================================================================
    790792#include "dimensions.h"
    791 #include "dimphy.h"
    792793#include "iniprint.h"
    793794      INTEGER knon
    … …  
    894895     .                  ts,u,v,t,
    895896     .                  pcfm, pcfh)
     897
     898      use dimphy
    896899      IMPLICIT none
    897900c======================================================================
    … …  
    914917c======================================================================
    915918#include "dimensions.h"
    916 #include "dimphy.h"
    917919#include "YOMCST.h"
    918920#include "iniprint.h"
    … …  
    949951      PARAMETER (prandtl=0.4)
    950952      INTEGER isommet ! le sommet de la couche limite
    951 c TEST VENUS
    952       PARAMETER (isommet=klev)
    953 c     PARAMETER (isommet=5)
    954953
    955954      LOGICAL tvirtu ! calculer Ri d'une maniere plus performante
    … …  
    995994      DATA appel1er /.TRUE./
    996995c
     996      isommet=klev
     997
    997998      IF (appel1er) THEN
    998999        if (prt_level > 9) THEN
    … …  
    11771178      SUBROUTINE coefkz2(knon, paprs, pplay,t,
    11781179     .                  pcfm, pcfh)
     1180
     1181      use dimphy
    11791182      IMPLICIT none
    11801183c======================================================================
    … …  
    11931196c======================================================================
    11941197#include "dimensions.h"
    1195 #include "dimphy.h"
    11961198#include "YOMCST.h"
    11971199#include "iniprint.h"

  • trunk/libf/phyvenus/clmain.classic

    r97 r101  
    3434      USE ioipsl
    3535      USE interface_surf
     36      use dimphy
    3637      IMPLICIT none
    3738c======================================================================
    … …  
    6566c======================================================================
    6667#include "dimensions.h"
    67 #include "dimphy.h"
    6868c$$$ PB ajout pour soil
    6969#include "dimsoil.h"
    … …  
    488488
    489489      USE interface_surf
     490      use dimphy
    490491
    491492      IMPLICIT none
    … …  
    495496c======================================================================
    496497#include "dimensions.h"
    497 #include "dimphy.h"
    498498#include "YOMCST.h"
    499499#include "dimsoil.h"
    … …  
    766766     e                  paprs,pplay,delp,
    767767     s                  d_ven,flux_v)
     768
     769      use dimphy
    768770      IMPLICIT none
    769771c======================================================================
    … …  
    789791c======================================================================
    790792#include "dimensions.h"
    791 #include "dimphy.h"
    792793#include "iniprint.h"
    793794      INTEGER knon
    … …  
    894895     .                  ts,u,v,t,
    895896     .                  pcfm, pcfh)
     897
     898      use dimphy
    896899      IMPLICIT none
    897900c======================================================================
    … …  
    914917c======================================================================
    915918#include "dimensions.h"
    916 #include "dimphy.h"
    917919#include "YOMCST.h"
    918920#include "iniprint.h"
    … …  
    11771179      SUBROUTINE coefkz2(knon, paprs, pplay,t,
    11781180     .                  pcfm, pcfh)
     1181
     1182      use dimphy
    11791183      IMPLICIT none
    11801184c======================================================================
    … …  
    11931197c======================================================================
    11941198#include "dimensions.h"
    1195 #include "dimphy.h"
    11961199#include "YOMCST.h"
    11971200#include "iniprint.h"

  • trunk/libf/phyvenus/clmain.simple

    r97 r101  
    2828c---------------------------------------------------------------
    2929      USE ioipsl
     30      use dimphy
    3031      IMPLICIT none
    3132c======================================================================
    … …  
    5960c======================================================================
    6061#include "dimensions.h"
    61 #include "dimphy.h"
    6262c$$$ PB ajout pour soil
    6363#include "dimsoil.h"
    … …  
    359359
    360360      USE interface_surf
     361      use dimphy
    361362
    362363      IMPLICIT none
    … …  
    366367c======================================================================
    367368#include "dimensions.h"
    368 #include "dimphy.h"
    369369#include "YOMCST.h"
    370370#include "dimsoil.h"
    … …  
    637637     e                  paprs,pplay,delp,
    638638     s                  d_ven,flux_v)
     639
     640      use dimphy
    639641      IMPLICIT none
    640642c======================================================================
    … …  
    660662c======================================================================
    661663#include "dimensions.h"
    662 #include "dimphy.h"
    663664#include "iniprint.h"
    664665      INTEGER knon
    … …  
    765766     .                  ts,u,v,t,
    766767     .                  pcfm, pcfh)
     768
     769      use dimphy
    767770      IMPLICIT none
    768771c======================================================================
    … …  
    785788c======================================================================
    786789#include "dimensions.h"
    787 #include "dimphy.h"
    788790#include "YOMCST.h"
    789791#include "iniprint.h"
    … …  
    10461048      SUBROUTINE coefkz2(knon, paprs, pplay,t,
    10471049     .                  pcfm, pcfh)
     1050
     1051      use dimphy
    10481052      IMPLICIT none
    10491053c======================================================================
    … …  
    10621066c======================================================================
    10631067#include "dimensions.h"
    1064 #include "dimphy.h"
    10651068#include "YOMCST.h"
    10661069#include "iniprint.h"

  • trunk/libf/phyvenus/cltrac.F

    r3 r101  
    44      SUBROUTINE cltrac(dtime,coef,t,tr,flux,paprs,pplay,delp,
    55     s                  d_tr)
     6
     7      use dimphy
    68      IMPLICIT none
    79c======================================================================
    … …  
    2628c======================================================================
    2729#include "dimensions.h"
    28 #include "dimphy.h"
    2930      REAL dtime
    3031      REAL coef(klon,klev)

  • trunk/libf/phyvenus/coefkzmin.F

    r3 r101  
    55     .   ,km,kn)
    66c      SUBROUTINE coefkzmin(ngrid,zlev,teta,ustar,km,kn)
     7
     8      use dimphy
    79      IMPLICIT NONE
    810
    911#include "dimensions.h"
    10 #include "dimphy.h"
    1112#include "YOMCST.h"
    1213
    … …  
    5758      integer l_mix,ngrid
    5859
    59 
    60       integer nlay,nlev
    61       PARAMETER (nlay=klev)
    62       PARAMETER (nlev=klev+1)
    6360
    6461      integer ig,k

  • trunk/libf/phyvenus/comcstVE.h

    r3 r101  
    22c INCLUDE comcstVE.h
    33
    4         integer nlve,nnuve,nbztopve,nbpsve,nbmat
    5         parameter (nlve=kflev) ! fichiers Vincent: same grid than GCM
     4        integer nnuve,nbztopve,nbpsve,nbmat
    65        parameter (nnuve=68)   ! fichiers Vincent et Bullock
    76        parameter (nbpsve=16)   ! number of psurf in Vincent's matrixes

  • trunk/libf/phyvenus/diagphy.F

    r97 r101  
    4848C======================================================================
    4949C
     50      use dimphy
    5051      implicit none
    5152
    5253#include "dimensions.h"
    53 #include "dimphy.h"
    5454#include "YOMCST.h"
    5555C
    … …  
    208208c======================================================================
    209209 
     210      use dimphy
    210211      IMPLICIT NONE
    211212C
    212213#include "dimensions.h"
    213 #include "dimphy.h"
    214214#include "YOMCST.h"
    215215C

  • trunk/libf/phyvenus/drag_noro.F

    r3 r101  
    99     s                   d_t, d_u, d_v)
    1010c
     11      use dimphy
    1112      IMPLICIT none
    1213
    1314#include "dimensions.h"
    1415#include "paramet.h"
    15 #include "dimphy.h"
    1616
    1717c======================================================================

  • trunk/libf/phyvenus/flott_gwd_ran.F

    r3 r101  
    1212C
    1313C
     14      use dimphy
    1415      implicit none
    1516
    1617#include "dimensions.h"
    1718#include "paramet.h"
    18 #include "dimphy.h"
    1919
    2020#include "YOEGWD.h"
    … …  
    4040
    4141c ON CONSERVE LA MEMOIRE un certain temps AVEC UN SAVE ET UN COEF
    42       real coef,d_u_sav(KLON,KLEV),d_v_sav(KLON,KLEV)
     42      real coef
    4343      parameter (coef = 0.986)
    44       save d_u_sav,d_v_sav
     44      real,save,allocatable :: d_u_sav(:,:),d_v_sav(:,:)
     45      LOGICAL firstcall
     46      SAVE firstcall
     47      DATA firstcall/.true./
    4548
    4649C  GRAVITY-WAVES SPECIFICATIONS
    … …  
    8588
    8689C  INITIALISATION
     90
     91      IF (firstcall) THEN
     92        allocate(d_u_sav(klon,klev),d_v_sav(klon,klev))
     93        firstcall=.false.
     94      ENDIF
    8795
    8896      BVSEC=1.E-6

  • trunk/libf/phyvenus/gwprofil.F

    r3 r101  
    3636C     ntop.                                         
    3737C
     38      use dimphy
    3839      IMPLICIT NONE
    3940
    4041#include "dimensions.h"
    4142#include "paramet.h"
    42 #include "dimphy.h"
    4343
    4444#include "YOMCST.h"

  • trunk/libf/phyvenus/gwstress.F

    r3 r101  
    5252c
    5353c-----------------------------------------------------------------------
     54      use dimphy
    5455      implicit none
    5556
    5657#include "dimensions.h"
    5758#include "paramet.h"
    58 #include "dimphy.h"
    5959
    6060#include "YOMCST.h"

  • trunk/libf/phyvenus/hgardfou.F

    r3 r101  
    33!
    44      SUBROUTINE hgardfou (t,tsol,text)
     5
     6      use dimphy
    57      IMPLICIT none
    68c======================================================================
    … …  
    810c======================================================================
    911#include "dimensions.h"
    10 #include "dimphy.h"
    1112#include "YOMCST.h"
    1213      REAL t(klon,klev), tsol(klon)

  • trunk/libf/phyvenus/ini_histmth.h

    r97 r101  
    122122             IF (iq.LE.99) THEN
    123123          WRITE(str2,'(i2.2)') iq
    124           CALL histdef(nid_mth, tnom(iq), ttext(iq), "ppm",
     124          CALL histdef(nid_mth, tname(iq), ttext(iq), "ppm",
    125125     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    126126     .                "ave(X)", zsto,zout)

  • trunk/libf/phyvenus/inifis.F

    r3 r101  
    77     $           plat,plon,parea,
    88     $           prad,pg,pr,pcpp)
     9
     10      use dimphy
    911      IMPLICIT NONE
    1012c
    … …  
    4244 
    4345#include "dimensions.h"
    44 #include "dimphy.h"
    4546
    4647      REAL prad,pg,pr,pcpp,punjours

  • trunk/libf/phyvenus/iniphysiq.F

    r3 r101  
    99     $           plat,plon,parea,pcu,pcv,
    1010     $           prad,pg,pr,pcpp)
     11
     12      use dimphy
     13      USE comgeomphy
    1114      IMPLICIT NONE
    1215c
    … …  
    4447 
    4548#include "dimensions.h"
    46 #include "dimphy.h"
    47 #include "comgeomphy.h"
    4849
    4950      REAL prad,pg,pr,pcpp,punjours

  • trunk/libf/phyvenus/load_ksi.F

    r3 r101  
    11      SUBROUTINE load_ksi(ksive)
    22     
     3      use dimphy
    34      IMPLICIT none
    45
    56#include "dimensions.h"
    6 #include "dimphy.h"
    7 #include "raddim.h"
    87#include "YOMCST.h"
    98#include "comcstVE.h"
    … …  
    3635
    3736c local variables
    38       integer i,j,ig,band,pve
     37      integer i,j,ig,band,pve,nlve
    3938      integer mat,Nb,m,Nmat,nl_init,mat0
    4039      parameter(nl_init=8)
    … …  
    4544      real   dlambda                  ! cm-1 
    4645
     46      nlve = kflev
     47
    4748c ------------------------
    4849c Loading the ksi file
    4950c ------------------------
    50  
    51        file = "ksi_global.txt"
     51
     52      file = "ksi_global.txt"
    5253      open(10,file=file)
    5354     

  • trunk/libf/phyvenus/load_psi.F

    r3 r101  
    33     S              temp, psimap, deltapsimap)
    44     
     5      use dimphy
    56      IMPLICIT none
    67
    78#include "dimensions.h"
    8 #include "dimphy.h"
    9 #include "raddim.h"
    109#include "YOMCST.h"
    1110#include "comcstVE.h"
    … …  
    4645
    4746c local variables
    48       integer i,j,ig,band
     47      integer i,j,ig,band,nlve
    4948      integer mat,m,mat0
    5049      character*100 file
    51       real   bplck(0:nlve+1,nnuve)    ! Planck luminances in table layers
    52       real   y(0:nlve,nnuve)          ! intermediaire Planck
    53       real   zdblay(0:nlve+1,nnuve)   ! gradient en temperature de planck
     50      real   bplck(0:kflev+1,nnuve)    ! Planck luminances in table layers
     51      real   y(0:kflev,nnuve)          ! intermediaire Planck
     52      real   zdblay(0:kflev+1,nnuve)   ! gradient en temperature de planck
    5453      real   ksi
    5554      real   factp,factz
     55
     56      nlve = kflev  ! (doit correspondre, pour bplck, y et zdblay)
    5657
    5758c -----------------------

  • trunk/libf/phyvenus/lw_venus_ve.F

    r3 r101  
    55     S              ZFLNET)
    66     
     7      use dimphy
    78      IMPLICIT none
    89
    910#include "dimensions.h"
    10 #include "dimphy.h"
    11 #include "raddim.h"
    1211#include "YOMCST.h"
    1312C

  • trunk/libf/phyvenus/lwi.F

    r3 r101  
    55c          differente des versions *.1mat 
    66c !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     7      use dimphy
    78      implicit none
    89
    910
    1011#include "dimensions.h"
    11 #include "dimphy.h"
    1212#include "comg1d.h"  !pour grads 1D
    1313#include "YOMCST.h"

  • trunk/libf/phyvenus/minmaxqfi.F

    r3 r101  
    33!
    44      SUBROUTINE minmaxqfi(zq,qmin,qmax,comment)
     5
     6      use dimphy
    57      IMPLICIT none
    68
    79#include "dimensions.h"
    8 #include "dimphy.h"
    910
    1011      CHARACTER*(*) comment

  • trunk/libf/phyvenus/orodrag.F

    r3 r101  
    88     r                 , pvom,pvol,pte )
    99     
     10      use dimphy
    1011      IMPLICIT NONE
    1112
    1213#include "dimensions.h"
    1314#include "paramet.h"
    14 #include "dimphy.h"
    1515
    1616c

  • trunk/libf/phyvenus/orosetup.F

    r3 r101  
    9595c
    9696c-----------------------------------------------------------------------
     97      use dimphy
    9798      implicit none
    9899
    99100#include "dimensions.h"
    100101#include "paramet.h"
    101 #include "dimphy.h"
    102102
    103103#include "YOMCST.h"

  • trunk/libf/phyvenus/phyetat0.F

    r97 r101  
    1111     .           tabcntr0,
    1212     .           t_ancien,ancien_ok)
     13
     14      use dimphy
    1315      IMPLICIT none
    1416c======================================================================
    … …  
    1719c======================================================================
    1820#include "dimensions.h"
    19 #include "dimphy.h"
    2021#include "netcdf.inc"
    2122#include "dimsoil.h"

  • trunk/libf/phyvenus/phyredem.F

    r97 r101  
    1212
    1313      USE control_mod
     14      use dimphy
    1415      IMPLICIT none
    1516c======================================================================
    … …  
    1819c======================================================================
    1920#include "dimensions.h"
    20 #include "dimphy.h"
    2121#include "netcdf.inc"
    2222#include "dimsoil.h"
    … …  
    7878      IF( cycle_diurne ) tab_cntrl( 7 ) = 1.
    7979      IF(   soil_model ) tab_cntrl( 8 ) = 1.
    80       IF(     new_oliq ) tab_cntrl( 9 ) = 1.
    8180      IF(     ok_orodr ) tab_cntrl(10 ) = 1.
    8281      IF(     ok_orolf ) tab_cntrl(11 ) = 1.

  • trunk/libf/phyvenus/physiq.F

    r97 r101  
    1010     .            d_u, d_v, d_t, d_qx, d_ps)
    1111
    12       USE ioipsl
    13       USE histcom
    14       USE infotrac
    15       USE control_mod
    16       IMPLICIT none
    1712c======================================================================
    1813c
    … …  
    2823cAA                  -  stockage des moyennes des champs necessaires
    2924cAA                     en mode traceur off-line
    30 c======================================================================
    31 c   CLEFS CPP POUR LES IO
    32 c   =====================
    33 c#define histhf
    34 #define histday
    35 #define histmth
    36 #define histins
    37 c======================================================================
    3825c    modif   ( P. Le Van ,  12/10/98 )
    3926c
    … …  
    6754c d_ps----output-R-tendance physique de la pression au sol
    6855c======================================================================
     56      USE ioipsl
     57      USE histcom
     58      USE infotrac
     59      USE control_mod
     60      use dimphy
     61      USE comgeomphy
     62      IMPLICIT none
     63c======================================================================
     64c   CLEFS CPP POUR LES IO
     65c   =====================
     66c#define histhf
     67#define histday
     68#define histmth
     69#define histins
     70c======================================================================
    6971#include "dimensions.h"
    7072      integer jjmp1
    7173      parameter (jjmp1=jjm+1-1/jjm)
    72 #include "dimphy.h"
    7374#include "dimsoil.h"
    7475#include "clesphys.h"
    7576#include "temps.h"
    76 #include "comgeomphy.h"
    7777#include "iniprint.h"
    78 #include "raddim.h"
    7978#include "timerad.h"
    8079#include "logic.h"
    … …  
    117116      REAL qx(klon,klev,nqmax)
    118117
    119       REAL t_ancien(klon,klev)
    120       REAL u_ancien(klon,klev)
    121       SAVE t_ancien, u_ancien
     118      REAL,save,allocatable :: t_ancien(:,:)
     119      REAL,save,allocatable :: u_ancien(:,:)
    122120      LOGICAL ancien_ok
    123121      SAVE ancien_ok
    … …  
    134132      REAL d_ps(klon)
    135133
    136       INTEGER klevp1, klevm1
    137       PARAMETER(klevp1=klev+1,klevm1=klev-1)
    138 c
    139       REAL swnet(klon,klevp1)
    140       SAVE swnet
    141 c
    142       REAL lwnet(klon,klevp1)
    143       SAVE lwnet
    144 c
    145       REAL LWdnTOA(klon)
    146       SAVE LWdnTOA
     134      REAL,save,allocatable :: swnet(:,:)
     135      REAL,save,allocatable :: lwnet(:,:)
    147136c
    148137      logical ok_hf
    … …  
    165154      SAVE radpas                 ! frequence d'appel rayonnement
    166155c
    167       REAL radsol(klon)
    168       SAVE radsol               ! bilan radiatif au sol calcule par code radiatif
    169 c
    170       REAL rlev(klon,klev+1)      ! altitude a chaque niveau (interface inferieure de la couche)
    171       SAVE rlev     
    172 ci 
    173       INTEGER itap
    174       SAVE itap                   ! compteur pour la physique
    175 c
    176       REAL ftsol(klon)
    177       SAVE ftsol                  ! temperature du sol
    178 c
    179       REAL ftsoil(klon,nsoilmx)
    180       SAVE ftsoil                 ! temperature dans le sol
    181 c
    182       REAL falbe(klon)
    183       SAVE falbe                  ! albedo
    184 
     156      REAL,save,allocatable :: radsol(:) ! bilan radiatif au sol calcule par code radiatif
     157      REAL,save,allocatable :: rlev(:,:) ! altitude a chaque niveau (interface inferieure de la couche)
     158      INTEGER,save :: itap        ! compteur pour la physique
     159      REAL,save,allocatable :: ftsol(:)    ! temperature du sol
     160      REAL,save,allocatable :: ftsoil(:,:) ! temperature dans le sol
     161      REAL,save,allocatable :: falbe(:)    ! albedo
    185162      REAL delp(klon,klev)        ! epaisseur d'une couche
    186163     
    … …  
    189166c  Parametres de l'Orographie a l'Echelle Sous-Maille (OESM):
    190167c
    191       REAL zmea(klon)
    192       SAVE zmea                   ! orographie moyenne
    193 c
    194       REAL zstd(klon)
    195       SAVE zstd                   ! deviation standard de l'OESM
    196 c
    197       REAL zsig(klon)
    198       SAVE zsig                   ! pente de l'OESM
    199 c
    200       REAL zgam(klon)
    201       save zgam                   ! anisotropie de l'OESM
    202 c
    203       REAL zthe(klon)
    204       SAVE zthe                   ! orientation de l'OESM
    205 c
    206       REAL zpic(klon)
    207       SAVE zpic                   ! Maximum de l'OESM
    208 c
    209       REAL zval(klon)
    210       SAVE zval                   ! Minimum de l'OESM
    211 c
    212       REAL rugoro(klon)
    213       SAVE rugoro                 ! longueur de rugosite de l'OESM
     168      REAL,save,allocatable :: zmea(:)   ! orographie moyenne
     169      REAL,save,allocatable :: zstd(:)   ! deviation standard de l'OESM
     170      REAL,save,allocatable :: zsig(:)   ! pente de l'OESM
     171      REAL,save,allocatable :: zgam(:)   ! anisotropie de l'OESM
     172      REAL,save,allocatable :: zthe(:)   ! orientation de l'OESM
     173      REAL,save,allocatable :: zpic(:)   ! Maximum de l'OESM
     174      REAL,save,allocatable :: zval(:)   ! Minimum de l'OESM
     175      REAL,save,allocatable :: rugoro(:) ! longueur de rugosite de l'OESM
    214176
    215177      INTEGER igwd,idx(klon),itest(klon)
    … …  
    241203
    242204      REAL zustrph(klon),zvstrph(klon)
    243 c
    244       REAL zuthe(klon),zvthe(klon)
    245       SAVE zuthe
    246       SAVE zvthe
     205      REAL,save,allocatable :: zuthe(:),zvthe(:)
    247206
    248207c Variables locales:
    … …  
    259218
    260219      REAL sens(klon), dsens(klon) ! chaleur sensible et sa derivee
    261       REAL dlw(klon)    ! derivee infra rouge
    262 cym
    263       SAVE dlw
    264 cym
    265       REAL fder(klon) ! Derive de flux (sensible et latente)
    266       save fder
     220      REAL,save,allocatable :: dlw(:)  ! derivee infra rouge
     221      REAL,save,allocatable :: fder(:) ! Derive de flux (sensible et latente)
    267222      REAL ve(klon) ! integr. verticale du transport meri. de l'energie
    268223      REAL vq(klon) ! integr. verticale du transport meri. de l'eau
    … …  
    307262      REAL flux_ec(klon,klev)   ! flux de chaleur Ec
    308263c
    309       REAL heat(klon,klev)    ! chauffage solaire
    310       REAL cool(klon,klev)    ! refroidissement infrarouge
    311       REAL dtrad(klon,klev)   ! K s-1
    312       REAL tmpout(klon,klev)  ! K s-1
    313       REAL topsw(klon), toplw(klon), solsw(klon), sollw(klon)
    314       real sollwdown(klon)    ! downward LW flux at surface
    315 
    316264c Le rayonnement n'est pas calcule tous les pas, il faut donc
    317265c                      sauvegarder les sorties du rayonnement
    318       SAVE  heat,cool,dtrad,topsw,toplw,solsw,sollw,sollwdown
    319 c
     266      REAL,save,allocatable :: heat(:,:)    ! chauffage solaire
     267      REAL,save,allocatable :: cool(:,:)    ! refroidissement infrarouge
     268      REAL,save,allocatable :: dtrad(:,:)   ! K s-1
     269      REAL,save,allocatable :: topsw(:), toplw(:)
     270      REAL,save,allocatable :: solsw(:), sollw(:)
     271      REAL,save,allocatable :: sollwdown(:) ! downward LW flux at surface
     272      REAL tmpout(klon,klev)  ! K s-1
     273
    320274      INTEGER itaprad
    321275      SAVE itaprad
    … …  
    453407c========================
    454408      IF (debut) THEN
     409         allocate(t_ancien(klon,klev),u_ancien(klon,klev))
     410         allocate(swnet(klon,klevp1),lwnet(klon,klevp1))
     411         allocate(radsol(klon),ftsol(klon),falbe(klon))
     412         allocate(rlev(klon,klevp1),ftsoil(klon,nsoilmx))
     413         allocate(zmea(klon),zstd(klon),zsig(klon),zgam(klon))
     414         allocate(zthe(klon),zpic(klon),zval(klon),rugoro(klon))
     415         allocate(zuthe(klon),zvthe(klon),dlw(klon),fder(klon))
     416         allocate(heat(klon,klev),cool(klon,klev))
     417         allocate(dtrad(klon,klev),topsw(klon),toplw(klon))
     418         allocate(solsw(klon),sollw(klon),sollwdown(klon))
    455419
    456420         CALL suphec ! initialiser constantes et parametres phys.
    … …  
    576540         ENDIF
    577541
    578          ecrit_ins = NINT(RDAY/dtime/*ecritphy)  ! Fraction de jour reglable
     542         ecrit_ins = NINT(RDAY/dtime*ecritphy)  ! Fraction de jour reglable
    579543         IF (ok_instan) THEN
    580544         WRITE(lunout,*)'La frequence de sortie instant. est de ',
    … …  
    852816
    853817      if (iflag_trac.eq.1) then
     818         DO k = 1, klev
     819         DO i = 1, klon
     820            delp(i,k) = paprs(i,k)-paprs(i,k+1)
     821         ENDDO
     822         ENDDO
    854823         DO iq=1, nqmax
    855824             CALL cltrac(dtime,ycoefh,t_seri,

  • trunk/libf/phyvenus/phytrac.F

    r97 r101  
    2525     O                    tr_seri)
    2626
     27c======================================================================
     28c Auteur(s) FH
     29c Objet: Moniteur general des tendances traceurs
     30c
     31cAA Remarques en vrac:
     32cAA--------------------
     33cAA 1/ le call phytrac se fait avec nqmax
     34c======================================================================
    2735      USE ioipsl
    2836      USE infotrac
    2937      USE control_mod
    30 
     38      use dimphy
     39      USE comgeomphy
    3140      IMPLICIT none
    32 c======================================================================
    33 c Auteur(s) FH
    34 c Objet: Moniteur general des tendances traceurs
    35 c
    36 cAA Remarques en vrac:
    37 cAA--------------------
    38 cAA 1/ le call phytrac se fait avec nqmax
    39 c======================================================================
    4041#include "YOMCST.h"
    4142#include "dimensions.h"
    42 #include "dimphy.h"
    4343#include "clesphys.h" !///utile?
    4444#include "temps.h"
    4545#include "paramet.h"
    46 #include "comgeomphy.h"
    4746c======================================================================
    4847
    … …  
    8281c      REAL flxmass_w(nlon,nlev)
    8382
    84 cAA Rem : nqmax : nombre de vrais traceurs est defini dans dimphy.h
    85 
    8683cAA ----------------------------
    8784cAA  VARIABLES LOCALES TRACEURS
    … …  
    124121c   k,i,it,tr_seri(klon,klev,nqmax),pplay(nlon,nlev),
    125122      integer nqCO_OCS
    126       real zprof(klev,nqtot)
    127123c      real pzero,gamma
    128124c      parameter (pzero=85000.)
    … …  
    130126      REAL alpha
    131127      real deltatr(klon,klev,nqtot) ! ecart au profil de ref zprof
    132       real tau(klev,nqtot)          ! temps de relaxation vers le profil (s)
    133       save zprof,tau
     128      real,save,allocatable :: zprof(:,:)
     129      real,save,allocatable :: tau(:,:) ! temps de relaxation vers le profil (s)
    134130c======================================================================
    135131c
    … …  
    151147      REAL d_tr_cv(klon,klev,nqmax) ! tendance de traceurs  conv pour chq traceur
    152148C
     149      character*20 modname
    153150      character*80 abort_message
    154151c
    … …  
    161158     s     /.true.,.false.,.false./
    162159
    163 c======================================================================
     160      modname = 'phytrac'
     161c======================================================================
     162
     163       if(first) then
     164         allocate(zprof(klev,nqtot),tau(klev,nqtot))
     165         first = .false.
     166       endif
     167
    164168         ps(:)=paprs(:,1)
    165169c TRACEURS TYPE CO ET OCS

  • trunk/libf/phyvenus/printflag.F

    r3 r101  
    1212
    1313       REAL tabcntr0( 100 )
    14        LOGICAL cycle_diurn0,soil_model0,new_oliq0,ok_orodr0
     14       LOGICAL cycle_diurn0,soil_model0,ok_orodr0
    1515       LOGICAL ok_orolf0,ok_gw_nonoro0
    1616       LOGICAL ok_journe,ok_instan
    … …  
    3232       PRINT 100
    3333
    34 c      PRINT 11, new_oliq, ok_orodr, ok_orolf, ok_gw_nonoro 
    35 c      PRINT 100
     34       PRINT 11, ok_orodr, ok_orolf, ok_gw_nonoro 
     35       PRINT 100
    3636
    3737       PRINT 12, nbapp_rad
    … …  
    4848        cycle_diurn0  = .FALSE.
    4949        soil_model0   = .FALSE.
    50         new_oliq0     = .FALSE.
    5150        ok_orodr0     = .FALSE.
    5251        ok_orolf0     = .FALSE.
    … …  
    5554        IF( tabcntr0( 7 ).EQ. 1. )   cycle_diurn0 = .TRUE.
    5655        IF( tabcntr0( 8 ).EQ. 1. )    soil_model0 = .TRUE.
    57         IF( tabcntr0( 9 ).EQ. 1. )      new_oliq0 = .TRUE.
    5856        IF( tabcntr0(10 ).EQ. 1. )      ok_orodr0 = .TRUE.
    5957        IF( tabcntr0(11 ).EQ. 1. )      ok_orolf0 = .TRUE.
    … …  
    8280     ,        soil_model )     THEN
    8381        PRINT 14, soil_model0, soil_model
    84         PRINT 100
    85        ENDIF
    86 
    87        IF( new_oliq0.AND..NOT.new_oliq.OR..NOT.new_oliq0.AND.
    88      ,        new_oliq )     THEN
    89         PRINT 16, new_oliq0, new_oliq
    9082        PRINT 100
    9183       ENDIF
    … …  
    127119
    128120
    129  11    FORMAT(2x,5(1H*),'  new_oliq = ',l3,3x,', Ok_orodr = ',
     121 11    FORMAT(2x,5(1H*),', Ok_orodr = ',
    130122     , l3,3x,', Ok_orolf = ',l3,3x,5(1H*) )
    131123
    … …  
    141133
    142134 15    FORMAT(2x,'$$$$$$$$   Attention !!      ok_orodr  different  sur',
    143      , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
    144 
    145  16    FORMAT(2x,'$$$$$$$$   Attention !!      new_oliq  different  sur',
    146135     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
    147136

  • trunk/libf/phyvenus/radlwsw.1mat

    r97 r101  
    99     .                  lwnet, swnet)
    1010c     
    11       IMPLICIT none
    1211c======================================================================
    1312c Auteur(s): Z.X. Li (LMD/CNRS) date: 19960719
    … …  
    3534     
    3635c======================================================================
     36      use dimphy
     37      IMPLICIT none
    3738#include "dimensions.h"
    38 #include "dimphy.h"
    39 #include "raddim.h"
    4039c
    4140      real rmu0(klon), fract(klon), dist

  • trunk/libf/phyvenus/radlwsw.F

    r97 r101  
    99     .                  lwnet, swnet)
    1010c     
    11       IMPLICIT none
    1211c======================================================================
    1312c Auteur(s): Z.X. Li (LMD/CNRS) date: 19960719
    … …  
    3938
    4039c======================================================================
     40      use dimphy
     41      USE comgeomphy
     42      IMPLICIT none
    4143#include "dimensions.h"
    42 #include "dimphy.h"
    43 #include "raddim.h"
    4444#include "YOMCST.h"
    4545#include "clesphys.h"
    46 #include "comgeomphy.h"
    4746#include "comcstVE.h"
    4847c
    … …  
    7170      REAL   zsollwdown
    7271cIM END
    73       real   ksive(0:kflev+1,0:kflev+1,nnuve,nbmat) ! ksi matrixes in Vincent's file
     72      real,save,allocatable :: ksive(:,:,:,:) ! ksi matrixes in Vincent's file
    7473      real    psimap(0:kflev+1,0:kflev+1,klon)
    7574      real    deltapsimap(0:kflev+1,0:kflev+1,klon)
    7675      real    psi(0:kflev+1,0:kflev+1)
    7776      real    deltapsi(0:kflev+1,0:kflev+1)
    78       real    latdeg,ztop(klon) ! in km
     77      real    latdeg
    7978      real    pt0(klon,0:kflev+1)
    80 
    81       save    ksive,ztop
     79      real,save,allocatable :: ztop(:) ! in km
    8280
    8381      logical firstcall
    … …  
    9290
    9391      if (firstcall) then
     92
     93c ---------- ksive --------------
     94        allocate(ksive(0:kflev+1,0:kflev+1,nnuve,nbmat))
    9495        call load_ksi(ksive)
    95 
    9696c ---------- ztop --------------
     97        allocate(ztop(klon))
    9798        DO i = 1, klon
    9899             ztop(i) = 70.
    99100        ENDDO !i
    100 
    101101c ztop: d'apres fit à figure 16 du papier Zavosa et al (tmp) traitant des
    102102c       donnees Venera

  • trunk/libf/phyvenus/radlwsw.NewtonCool

    r97 r101  
    99     .                  lwnet, swnet)
    1010c     
    11       IMPLICIT none
    1211c======================================================================
    1312c Auteur(s): Z.X. Li (LMD/CNRS) date: 19960719
    … …  
    3837
    3938c======================================================================
     39      use dimphy
     40      USE comgeomphy
     41      IMPLICIT none
    4042#include "dimensions.h"
    41 #include "dimphy.h"
    42 #include "raddim.h"
    4343#include "YOMCST.h"
    4444#include "clesphys.h"
    45 #include "comgeomphy.h"
    4645c
    4746      real rmu0(klon), fract(klon), dist

  • trunk/libf/phyvenus/radlwsw.multimat

    r97 r101  
    99     .                  lwnet, swnet)
    1010c     
    11       IMPLICIT none
    1211c======================================================================
    1312c Auteur(s): Z.X. Li (LMD/CNRS) date: 19960719
    … …  
    3938
    4039c======================================================================
     40      use dimphy
     41      USE comgeomphy
     42      IMPLICIT none
    4143#include "dimensions.h"
    42 #include "dimphy.h"
    43 #include "raddim.h"
    4444#include "YOMCST.h"
    4545#include "clesphys.h"
    46 #include "comgeomphy.h"
    4746#include "comcstVE.h"
    4847c

  • trunk/libf/phyvenus/readstartphy.F

    r3 r101  
    99     .           fder,radsol,
    1010     .           tabcntr0)
    11       IMPLICIT none
    1211c======================================================================
    1312c Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
    1413c Objet: Lecture de l'etat initial pour la physique
    1514c======================================================================
     15      use dimphy
     16      IMPLICIT none
    1617#include "dimensions.h"
    17 #include "dimphy.h"
    1818#include "netcdf.inc"
    1919#include "dimsoil.h"

  • trunk/libf/phyvenus/soil.F

    r3 r101  
    44      SUBROUTINE soil(ptimestep, knon, ptsrf, ptsoil,
    55     s          pcapcal, pfluxgrd)
    6       IMPLICIT NONE
    76
    87c=======================================================================
    … …  
    4746c   -------------
    4847
     48      use dimphy
     49      IMPLICIT NONE
    4950#include "dimensions.h"
    5051#include "YOMCST.h"
    51 #include "dimphy.h"
    5252#include "dimsoil.h"
    5353#include "clesphys.h"
    … …  
    7474c   ----------------------
    7575      REAL dz1(nsoilmx),dz2(nsoilmx)
    76       REAL zc(klon,nsoilmx),zd(klon,nsoilmx)
     76      REAL,allocatable :: zc(:,:),zd(:,:)
    7777      REAL lambda
    7878      SAVE dz1,dz2,zc,zd,lambda
    … …  
    9595
    9696      IF (firstcall) THEN
     97
     98      allocate(zc(klon,nsoilmx),zd(klon,nsoilmx))
    9799
    98100c-----------------------------------------------------------------------

  • trunk/libf/phyvenus/sugwd.F

    r3 r101  
    7474C        LAST:  99-07-09     (FRANCOIS LOTT,LMD)
    7575C     ------------------------------------------------------------------
     76      use dimphy
    7677      IMPLICIT NONE
    7778
    7879#include "dimensions.h"
    7980#include "paramet.h"
    80 #include "dimphy.h"
    8181
    8282#include "YOEGWD.h"

  • trunk/libf/phyvenus/sw_venus_dc.F

    r3 r101  
    44     S              PTOPSW,PSOLSW,ZFSNET)
    55     
     6      use dimphy
    67      IMPLICIT none
    78
    89#include "dimensions.h"
    9 #include "dimphy.h"
    10 #include "raddim.h"
    1110#include "YOMCST.h"
    1211C
    … …  
    7170      logical firstcall
    7271      data    firstcall/.true./
    73       REAL   zsolVE(KFLEV+1) ! net solar flux at ppb levels, fichiers VE
    74      
    75       save   solza,zsnet,presdc,tempdc,altdc,zsolVE
     72      REAL,save,allocatable :: zsolVE(:) ! net solar flux at ppb levels, fichiers VE
     73      save   solza,zsnet,presdc,tempdc,altdc
    7674      save   firstcall
    7775     
    … …  
    8179
    8280      if (firstcall) then
     81       allocate(zsolVE(klevp1))
     82
    8383       open(11,file='dataDCrisp.dat')
    8484       read(11,*) nullchar

  • trunk/libf/phyvenus/transp.F

    r3 r101  
    66     s                   vtran_e, vtran_q, utran_e, utran_q)
    77c
    8       IMPLICIT none
    98c======================================================================
    109c Auteur(s): Z.X.Li (LMD/CNRS)
    … …  
    1211c Objet: Calculer le transport de l'energie et de la vapeur d'eau
    1312c======================================================================
    14 c
     13      use dimphy
     14      IMPLICIT none
    1515#include "dimensions.h"
    16 #include "dimphy.h"
    1716#include "YOMCST.h"
    1817c

  • trunk/libf/phyvenus/ustarhb.F

    r97 r101  
    33!
    44      SUBROUTINE ustarhb(knon,u,v,cd_m, ustar)
    5       IMPLICIT none
    65c======================================================================
    76c Laurent Li (LMD/CNRS), le 30 septembre 1998
    … …  
    1716c model. J. of Climate, vol. 6, 1825-1842.
    1817c======================================================================
     18      use dimphy
     19      IMPLICIT none
    1920#include "dimensions.h"
    20 #include "dimphy.h"
    2121#include "YOMCST.h"
    2222c

  • trunk/libf/phyvenus/vdif_kcay.F

    r3 r101  
    55     s   ,zlev,zlay,u,v,teta,cd,q2,q2diag,km,kn,ustar
    66     s   ,l_mix)
     7c.......................................................................
     8      use dimphy
    79      IMPLICIT NONE
    8 c.......................................................................
    910#include "dimensions.h"
    10 #include "dimphy.h"
    1111c.......................................................................
    1212c
    … …  
    195195      PARAMETER (q2min=1.e-5)
    196196      PARAMETER (q2max=1.E+2)
    197       PARAMETER (nlay=klev)
    198       PARAMETER (nlev=klev+1)
    199197c
    200198      PARAMETER (
    … …  
    221219      save first
    222220      data first/.true./
     221
     222      nlay=klev
     223      nlev=klev+1
    223224c.......................................................................
    224225c  traitment des valeur de q2 en entree

  • trunk/libf/phyvenus/write_histmth.h

    r97 r101  
    109109       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
    110110       WRITE(str2,'(i2.2)') iq
    111        CALL histwrite(nid_mth,tnom(iq),itau_w,zx_tmp_3d,
     111       CALL histwrite(nid_mth,tname(iq),itau_w,zx_tmp_3d,
    112112     .                                   iim*jjmp1*klev,ndex3d)
    113113          ELSE

  • trunk/libf/phyvenus/writerestartphy.F

    r3 r101  
    66     .           t_ancien)
    77
    8       IMPLICIT none
    98
    109c======================================================================
    1110c Ecriture de l'etat de redemarrage pour la physique (pour newstart)
    1211c======================================================================
     12      use dimphy
     13      IMPLICIT none
    1314#include "dimensions.h"
    14 #include "dimphy.h"
    1515#include "netcdf.inc"
    1616#include "dimsoil.h"

  • trunk/libf/phyvenus/yamada.F

    r3 r101  
    55     s   ,zlev,zlay,u,v,teta,cd,q2,km,kn,ustar
    66     s   ,l_mix)
     7c.......................................................................
     8      use dimphy
    79      IMPLICIT NONE
    8 c.......................................................................
    910#include "dimensions.h"
    10 #include "dimphy.h"
    1111c.......................................................................
    1212c
    … …  
    4646      integer l_mix,ngrid
    4747
    48 
    49       integer nlay,nlev
    50       PARAMETER (nlay=klev)
    51       PARAMETER (nlev=klev+1)
    52 
    5348      logical first
    5449      save first
    5550      data first/.true./
    56 
    5751
    5852      integer ig,k
    … …  
    9892
    9993c  Correction d'un bug sauvage a verifier.
    100 c      do k=2,nlev
    101       do k=2,nlay
     94c      do k=2,klevp1
     95      do k=2,klev
    10296                                                          do ig=1,ngrid
    10397         dz(ig,k)=zlay(ig,k)-zlay(ig,k-1)

  • trunk/libf/phyvenus/yamada4.F

    r3 r101  
    66     s   ,zlev,zlay,u,v,teta,cd,km,kn,kq,ustar
    77     s   ,iflag_pbl)
     8c.......................................................................
     9      use dimphy
    810      IMPLICIT NONE
    9 c.......................................................................
    1011#include "dimensions.h"
    11 #include "dimphy.h"
    1212c.......................................................................
    1313c
    … …  
    5050      REAL teta(klon,klev)
    5151      REAL cd(klon)
    52       REAL q2(klon,klev+1),qpre
     52      REAL qpre
    5353      REAL unsdz(klon,klev)
    5454      REAL unsdzdec(klon,klev+1)
    … …  
    6565
    6666      integer nlay,nlev
    67       PARAMETER (nlay=klev)
    68       PARAMETER (nlev=klev+1)
    6967
    7068      logical first
    … …  
    8381      real dtetadz(klon,klev+1)
    8482      real m2cstat,mcstat,kmcstat
    85       real l(klon,klev+1),l0(klon)
    86       save l0
     83      real l(klon,klev+1)
     84      real,save,allocatable :: l0(:)
    8785c  ATTENTION! mis ici car j'ai enlevé q2 des arguments...
    8886c   sinon, c'est au-dessus que ça se passe...
    89       save q2
     87      REAL,save,allocatable :: q2(:,:)
    9088
    9189      real sq(klon),sqz(klon),zz(klon,klev+1)
    … …  
    9997      real fl,zzz,zl0,zq2,zn2
    10098
    101       real rino(klon,klev+1),smyam(klon,klev),styam(klon,klev)
    102      s  ,lyam(klon,klev),knyam(klon,klev)
    103      s  ,w2yam(klon,klev),t2yam(klon,klev)
    104       common/pbldiag/rino,smyam,styam,lyam,knyam,w2yam,t2yam
     99c     real rino(klon,klev+1),smyam(klon,klev),styam(klon,klev)
     100c    s  ,lyam(klon,klev),knyam(klon,klev)
     101c    s  ,w2yam(klon,klev),t2yam(klon,klev)
    105102
    106103      frif(ri)=0.6588*(ri+0.1776-sqrt(ri*ri-0.3221*ri+0.03156))
    … …  
    116113
    117114c===================================
    118 c INITIALISATIONS (surtout pour k=1, à cause diagnostiques...)
     115c INITIALISATIONS
     116      nlay=klev
     117      nlev=klev+1
     118
    119119      if (first) then
     120        allocate(l0(klon))
     121        allocate(q2(klon,klevp1))
     122
     123c (surtout pour k=1, à cause diagnostiques...)
    120124        q2 = 0.
    121125      endif
    … …  
    462466
    463467c     print*,'YAMADA4 1'
    464 c   Diagnostique pour stockage
    465 
    466       rino=rif
    467       smyam(:,1:klev)=sm(:,1:klev)
    468       styam=sm(:,1:klev)*alpha(:,1:klev)
    469       lyam(1:klon,1:klev)=l(:,1:klev)
    470       knyam(1:klon,1:klev)=kn(:,1:klev)
    471468
    472469c   Estimations de w'2 et T'2 d'apres Abdela et McFarlane
    473470
    474         if(1.eq.0)then
    475       w2yam=q2(:,1:klev)*0.24
    476      s    +lyam(:,1:klev)*5.17*kn(:,1:klev)*n2(:,1:klev)
    477      s   /sqrt(q2(:,1:klev))
    478 
    479       t2yam=9.1*kn(:,1:klev)*dtetadz(:,1:klev)**2/sqrt(q2(:,1:klev))
    480      s  *lyam(:,1:klev)
    481         endif
     471c     if(1.eq.0)then
     472c      w2yam=q2(:,1:klev)*0.24
     473c    s    +lyam(:,1:klev)*5.17*kn(:,1:klev)*n2(:,1:klev)
     474c    s   /sqrt(q2(:,1:klev))
     475c
     476c      t2yam=9.1*kn(:,1:klev)*dtetadz(:,1:klev)**2/sqrt(q2(:,1:klev))
     477c    s  *lyam(:,1:klev)
     478c     endif
    482479
    483480c     print*,'OKFIN'

  • trunk/libf/phyvenus/zenang.F

    r3 r101  
    22      SUBROUTINE zenang(longi,gmtime,pdtrad,lat,long,
    33     s                  pmu0,frac)
    4       IMPLICIT none
    54c=============================================================
    65c Auteur : O. Boucher (LMD/CNRS)
    … …  
    2524c frac-----OUTPUT: ensoleillement moyen entre gmtime et gmtime+pdtrad
    2625c================================================================
     26      use dimphy
     27      IMPLICIT none
    2728#include "dimensions.h"
    28 #include "dimphy.h"
    2929#include "YOMCST.h"
    3030c================================================================

  • trunk/makelmdz

    r66 r101  
    273273if [[ "$physique" != "nophys" ]]
    274274then
     275   CPP_KEY="$CPP_KEY CPP_PHYS"
     276fi
     277
     278if [[ "$physique" == "lmd" ]]
     279then
    275280   #Default planet type is Earth
    276281   CPP_KEY="$CPP_KEY CPP_EARTH"
Note: See TracChangeset for help on using the changeset viewer.