Changeset 102


Ignore:
Timestamp:
Mar 24, 2011, 11:46:27 AM (14 years ago)
Author:
slebonnois
Message:

SL : corrections et modifications dans phytitan correspondant a celles
faites apres compilation Venus. Titan pas encore compile.

Location:
trunk/libf
Files:
19 deleted
36 edited
4 copied

Legend:

Unmodified
Added
Removed

  • trunk/libf/phytitan/YOMCST.h

    r97 r102  
    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/phytitan/aaam_bud.F

    r3 r102  
    22     i                   rea,rg,ome,     
    33     i                   plat,plon,phis,
    4      i                   dragu,liftu,phyu,
    5      i                   dragv,liftv,phyv,
     4     i                   dragu,liftu,clu,
     5     i                   dragv,liftv,clv,
    66     i                   p, u, v)
    77c
     8      use dimphy
    89      implicit none
    910c======================================================================
    … …  
    1718C    C               ra,rg,romega,
    1819C    C               rlat,rlon,pphis,
    19 C    C               zustrdr,zustrli,zustrph,
    20 C    C               zvstrdr,zvstrli,zvstrph,
     20C    C               zustrdr,zustrli,zustrcl,
     21C    C               zvstrdr,zvstrli,zvstrcl,
    2122C    C               paprs,u,v)
    2223C
    … …  
    3940c dragu---input-R-orodrag stress (zonal)
    4041c liftu---input-R-orolift stress (zonal)
    41 c phyu----input-R-Stress total de la physique (zonal)
     42c clu-----input-R-Boundary layer stress (zonal)
    4243c dragv---input-R-orodrag stress (Meridional)
    4344c liftv---input-R-orolift stress (Meridional)
    44 c phyv----input-R-Stress total de la physique (Meridional)
     45c clv-----input-R-Boundary layer stress (Meridional)
    4546c p-------input-R-Pressure (Pa) at model half levels
    4647c u-------input-R-Horizontal wind (m/s)
    … …  
    7879
    7980#include "dimensions.h"
    80 #include "dimphy.h"
    8181c
    8282c ARGUMENTS
    … …  
    8585      REAL rjour,rsec,rea,rg,ome
    8686      REAL plat(nlon),plon(nlon),phis(nlon)
    87       REAL dragu(nlon),liftu(nlon),phyu(nlon)             
    88       REAL dragv(nlon),liftv(nlon),phyv(nlon)             
     87      REAL dragu(nlon),liftu(nlon),clu(nlon)             
     88      REAL dragv(nlon),liftv(nlon),clv(nlon)             
    8989      REAL p(nlon,nlev+1), u(nlon,nlev), v(nlon,nlev)
    9090c
    … …  
    148148          ssou(i,1)=dragu(l)+liftu(l)
    149149          ssov(i,1)=dragv(l)+liftv(l)
    150           blsu(i,1)=phyu(l)-dragu(l)-liftu(l)
    151           blsv(i,1)=phyv(l)-dragv(l)-liftv(l)
     150          blsu(i,1)=clu(l)
     151          blsv(i,1)=clv(l)
    152152
    153153        enddo
    … …  
    162162          ssou(iim+1,j)=dragu(l+1)+liftu(l+1)
    163163          ssov(iim+1,j)=dragv(l+1)+liftv(l+1)
    164           blsu(iim+1,j)=phyu(l+1)-dragu(l+1)-liftu(l+1)
    165           blsv(iim+1,j)=phyv(l+1)-dragv(l+1)-liftv(l+1)
     164          blsu(iim+1,j)=clu(l+1)
     165          blsv(iim+1,j)=clv(l+1)
    166166      zlon(iim+1)=-plon(l+1)*xpi/180.
    167167      zlat(j)=plat(l+1)*xpi/180.
    … …  
    182182          ssou(i,j)=dragu(l)+liftu(l)
    183183          ssov(i,j)=dragv(l)+liftv(l)
    184           blsu(i,j)=phyu(l)-dragu(l)-liftu(l)
    185           blsv(i,j)=phyv(l)-dragv(l)-liftv(l)
     184          blsu(i,j)=clu(l)
     185          blsv(i,j)=clv(l)
    186186      zlon(i)=plon(l)*xpi/180.
    187187
    … …  
    214214          ssou(i,jjm+1)=dragu(l)+liftu(l)
    215215          ssov(i,jjm+1)=dragv(l)+liftv(l)
    216           blsu(i,jjm+1)=phyu(l)-dragu(l)-liftu(l)
    217           blsv(i,jjm+1)=phyv(l)-dragv(l)-liftv(l)
     216          blsu(i,jjm+1)=clu(l)
     217          blsv(i,jjm+1)=clv(l)
    218218      ub(i,jjm+1)=ub(1,jjm+1)                               
    219219      vb(i,jjm+1)=vb(1,jjm+1)                               
    … …  
    333333     c      tmou(3)/hadley,tsso(3)/hadley,tbls(3)/hadley
    334334100   format(F12.5,15(1x,F12.5))
    335 c100   format(F12.5,5(1x,F12.5))
     335c00   format(F12.5,5(1x,F12.5))
    336336
    337337      write(iam+1,*)((zs(i,j),i=1,iim),j=1,jjm+1)

  • trunk/libf/phytitan/ajsec.F

    r3 r102  
    33!
    44! ADAPTATION GCM POUR CP(T)
    5       SUBROUTINE ajsec(paprs, pplay, ppk, tfi, ufi, vfi, nq, qfi,
     5      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======================================================================
    … …  
    2325c d_vfi-----output-R-Incrementation du vent meridien
    2426c d_qfi-----output-R-Incrementation des traceurs
    25 
    2627c======================================================================
    2728#include "dimensions.h"
    28 #include "dimphy.h"
    2929#include "YOMCST.h"
    3030      REAL paprs(klon,klev+1), pplay(klon,klev)
    … …  
    3636      REAL qfi(klon,klev,nq), d_qfi(klon,klev,nq)
    3737c
    38       INTEGER limbas, limhau ! les couches a ajuster
    39 ccc      PARAMETER (limbas=klev-3, limhau=klev)
    40       PARAMETER (limbas=1, limhau=klev)
     38      INTEGER,save :: limbas, limhau ! les couches a ajuster
    4139c
    4240      REAL zh(klon,klev)
    … …  
    5149c Initialisation:
    5250c
     51      limbas=1
     52      limhau=klev
     53
    5354      DO k = 1, klev
    5455      DO i = 1, klon
    … …  
    5657         d_ufi(i,k) = 0.0
    5758         d_vfi(i,k) = 0.0
    58         d_qfi(i,k,:) = 0.0
     59        d_qfi(i,k,:) = 0.0
    5960         zu(i,k)    = ufi(i,k)
    6061         zv(i,k)    = vfi(i,k)
    61         zq(i,k,:)  = qfi(i,k,:)
     62        zq(i,k,:)  = qfi(i,k,:)
    6263      ENDDO
    6364      ENDDO
    … …  
    119120              zum=0.
    120121              zvm=0.
    121               zqm=0.
     122              zqm=0.
    122123              DO k = k1, k2
    123124                zalpha=zalpha+ABS(zh(i,k)-hm)*zdp(i,k)
    … …  
    125126                zum=zum+ufi(i,k)*zdp(i,k)
    126127                zvm=zvm+vfi(i,k)*zdp(i,k)
    127                 do iq=1,nq
    128                    zqm(iq)=zqm(iq)+qfi(i,k,iq)*zdp(i,k)
    129                 enddo
     128                do iq=1,nq
     129                  zqm(iq)=zqm(iq)+qfi(i,k,iq)*zdp(i,k)
     130                enddo
    130131              ENDDO
    131132              zalpha=zalpha/(hm*(paprs(i,k1)-paprs(i,k2+1)))
    132133              zum=zum/(paprs(i,k1)-paprs(i,k2+1))
    133134              zvm=zvm/(paprs(i,k1)-paprs(i,k2+1))
    134               do iq=1,nq
    135                  zqm(iq)=zqm(iq)/(paprs(i,k1)-paprs(i,k2+1))
    136               enddo
     135              do iq=1,nq
     136                zqm(iq)=zqm(iq)/(paprs(i,k1)-paprs(i,k2+1))
     137              enddo
    137138
    138139              IF(zalpha.GT.1.) THEN
    … …  
    153154                 zv(i,k)=vfi(i,k)+zalpha*(zvm-vfi(i,k))
    154155                 do iq=1,nq
    155                    zq(i,k,iq)=qfi(i,k,iq)+zalpha*(zqm(iq)-qfi(i,k,iq))
     156                   zq(i,k,iq)=qfi(i,k,iq)+zalpha*(zqm(iq)-qfi(i,k,iq))
    156157                 enddo
    157158              ENDDO
    … …  
    187188         do iq=1,nq
    188189           d_qfi(i,k,iq) = 0.0
    189          enddo 
     190         enddo
    190191      ENDDO
    191192      ENDDO
    … …  
    200201         do iq=1,nq
    201202           d_qfi(i,k,iq) = 0.0
    202          enddo 
     203         enddo
    203204      ENDDO
    204205      ENDDO

  • trunk/libf/phytitan/ballon.F

    r3 r102  
    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
    … …  
    255255     .                       phi,lam,ubal,vbal,abal)
    256256
     257      use dimphy
    257258      implicit none
    258259
    … …  
    283284
    284285#include "dimensions.h"
    285 #include "dimphy.h"
    286286#include "YOMCST.h"
    287287c
    … …  
    354354      subroutine pos_implem(phi,lam,ubal,vbal,dt)
    355355
     356      use dimphy
    356357      implicit none
    357358
    … …  
    370371
    371372#include "dimensions.h"
    372 #include "dimphy.h"
    373373#include "YOMCST.h"
    374374c

  • trunk/libf/phytitan/calchim.F

    r97 r102  
    22     .                   ctemp,cplay,cplev,
    33     .                   dqyc)
    4       implicit none
    54     
    65c-------------------------------------------------
    76c 
    8 c     Introduction d'une routine chimique
     7c     Introduction d une routine chimique
    98c
    109c     Auteur: S. Lebonnois,  01/2000 | 09/2003
    … …  
    1312c-------------------------------------------------
    1413c
    15 
     14      use dimphy
     15      implicit none
    1616#include "dimensions.h"
    17 #include "dimphy.h"
    1817#include "clesphys.h"
    1918#include "paramet.h"

  • trunk/libf/phytitan/clmain.F

    r3 r102  
    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"
    … …  
    477477
    478478      USE interface_surf
     479      use dimphy
    479480
    480481      IMPLICIT none
    … …  
    484485c======================================================================
    485486#include "dimensions.h"
    486 #include "dimphy.h"
    487487#include "YOMCST.h"
    488488#include "dimsoil.h"
    … …  
    755755     e                  paprs,pplay,delp,
    756756     s                  d_ven,flux_v)
     757
     758      use dimphy
    757759      IMPLICIT none
    758760c======================================================================
    … …  
    778780c======================================================================
    779781#include "dimensions.h"
    780 #include "dimphy.h"
    781782#include "iniprint.h"
    782783      INTEGER knon
    … …  
    883884     .                  ts,u,v,t,
    884885     .                  pcfm, pcfh)
     886
     887      use dimphy
    885888      IMPLICIT none
    886889c======================================================================
    … …  
    903906c======================================================================
    904907#include "dimensions.h"
    905 #include "dimphy.h"
    906908#include "YOMCST.h"
    907909#include "iniprint.h"
    … …  
    938940      PARAMETER (prandtl=0.4)
    939941      INTEGER isommet ! le sommet de la couche limite
    940 c TEST VENUS
    941       PARAMETER (isommet=klev)
    942 c     PARAMETER (isommet=5)
    943942
    944943      LOGICAL tvirtu ! calculer Ri d'une maniere plus performante
    … …  
    984983      DATA appel1er /.TRUE./
    985984c
     985      isommet=klev
     986
    986987      IF (appel1er) THEN
    987988        if (prt_level > 9) THEN
    … …  
    11661167      SUBROUTINE coefkz2(knon, paprs, pplay,t,
    11671168     .                  pcfm, pcfh)
     1169
     1170      use dimphy
    11681171      IMPLICIT none
    11691172c======================================================================
    … …  
    11821185c======================================================================
    11831186#include "dimensions.h"
    1184 #include "dimphy.h"
    11851187#include "YOMCST.h"
    11861188#include "iniprint.h"

  • trunk/libf/phytitan/cltrac.F

    r3 r102  
    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)
    … …  
    118119      ENDDO
    119120c
     121c ATTENTION SHUNTE!!!!!!
     122c     DO k = 1, klev
     123c     DO i = 1, klon
     124c        d_tr(i,k) = 0.
     125c     ENDDO
     126c     ENDDO
     127
    120128      RETURN
    121129      END

  • trunk/libf/phytitan/cooling.F

    r3 r102  
    11      SUBROUTINE COOLING(ngrid,NL,PRESS,TEMP,Z,Q0,lwnet,pfluxi)
    2       IMPLICIT NONE
    32
    43c=======================================================================
    … …  
    5251c   ------------
    5352
     53      use dimphy
     54      IMPLICIT NONE
    5455#include "dimensions.h"
    55 #include "dimphy.h"
    5656#include "YOMCST.h"
    5757      INTEGER NLAYER,NSPECI,NSPC1I

  • trunk/libf/phytitan/diagphy.F

    r3 r102  
    77     $    , d_etp_tot, d_qt_tot, d_ec_tot
    88     $    , fs_bound, fq_bound)
     9
     10! ATTENTION !! PAS DU TOUT A JOUR POUR VENUS OU TITAN...
     11
    912C======================================================================
    1013C
    … …  
    4548C======================================================================
    4649C
     50      use dimphy
    4751      implicit none
    4852
    4953#include "dimensions.h"
    50 #include "dimphy.h"
    5154#include "YOMCST.h"
    5255C
    … …  
    205208c======================================================================
    206209 
     210      use dimphy
    207211      IMPLICIT NONE
    208212C
    209213#include "dimensions.h"
    210 #include "dimphy.h"
    211214#include "YOMCST.h"
    212215C

  • trunk/libf/phytitan/heating.F

    r3 r102  
    11       SUBROUTINE heating(dist,rmu0,fract,sol_htg,swnet)
    2 
    3       IMPLICIT NONE
    42
    53
    … …  
    3028c   -------------
    3129
     30      use dimphy
     31      IMPLICIT NONE
    3232#include "dimensions.h"
    33 #include "dimphy.h"
    3433
    3534      INTEGER NLEVEL,NLAYER,NSPECV

  • trunk/libf/phytitan/hgardfou.F

    r3 r102  
    33!
    44      SUBROUTINE hgardfou (t,tsol,text)
    5       IMPLICIT none
    65c======================================================================
    76c Verifier la temperature
    87c======================================================================
     8      use dimphy
     9      IMPLICIT none
    910#include "dimensions.h"
    10 #include "dimphy.h"
    1111#include "YOMCST.h"
    1212      REAL t(klon,klev), tsol(klon)

  • trunk/libf/phytitan/inicondens.F

    r3 r102  
    99c   -------------
    1010
     11      use dimphy
    1112      IMPLICIT NONE
    1213#include "dimensions.h"
    13 #include "dimphy.h"
    1414#include "YOMCST.h"
    1515

  • trunk/libf/phytitan/iniphysiq.F

    r3 r102  
    99     $           plat,plon,parea,pcu,pcv,
    1010     $           prad,pg,pr,pcpp)
    11       IMPLICIT NONE
     11
    1212c
    1313c=======================================================================
    … …  
    4343c   -------------
    4444 
     45      use dimphy
     46      USE comgeomphy
     47      IMPLICIT NONE
    4548#include "dimensions.h"
    46 #include "dimphy.h"
    47 #include "comgeomphy.h"
    4849
    4950      REAL prad,pg,pr,pcpp,punjours

  • trunk/libf/phytitan/optci.F

    r97 r102  
    11      SUBROUTINE OPTCI(ykim,nmicro,IPRINT)
     2      use dimphy
    23#include "dimensions.h"
    3 #include "dimphy.h"
    44#include "microtab.h"
    55#include "numchimrad.h"

  • trunk/libf/phytitan/optci_1pt.F

    r3 r102  
    22     .            COSBI_1pt,DTAUI_1pt,TAUHI_1pt,TAUHID_1pt,TAUCI_1pt,
    33     .            TAUGI_1pt,TAUGID_1pt,WBARI_1pt,TAUI_1pt,IPRINT)
     4      use dimphy
    45#include "dimensions.h"
    5 #include "dimphy.h"
    66#include "microtab.h"
    77#include "numchimrad.h"

  • trunk/libf/phytitan/optcv.F

    r97 r102  
    22
    33
     4      use dimphy
    45#include "dimensions.h"
    5 #include "dimphy.h"
    66#include "microtab.h"
    77#include "clesphys.h"

  • trunk/libf/phytitan/optcv_1pt.F

    r3 r102  
    44
    55
     6      use dimphy
    67#include "dimensions.h"
    7 #include "dimphy.h"
    88#include "microtab.h"
    99#include "clesphys.h"

  • trunk/libf/phytitan/pg2.F

    r3 r102  
    4141c------------------------------------------------------
    4242
     43      use dimphy
     44      USE comgeomphy
    4345#include "dimensions.h"
    44 #include "dimphy.h"
    4546#include "microtab.h"
    4647#include "paramet.h"
    4748#include "aerprod.h"
    4849#include "clesphys.h"
    49 #include "comgeomphy.h"
    5050#include "YOMCST.h"
    5151

  • trunk/libf/phytitan/phyetat0.F

    r97 r102  
    1010     .           tabcntr0,
    1111     .           t_ancien,ancien_ok)
    12       IMPLICIT none
    1312c======================================================================
    1413c Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
    1514c Objet: Lecture de l'etat initial pour la physique
    1615c======================================================================
     16      use dimphy
     17      IMPLICIT none
    1718#include "dimensions.h"
    18 #include "dimphy.h"
    1919#include "netcdf.inc"
    2020#include "dimsoil.h"

  • trunk/libf/phytitan/phyredem.F

    r97 r102  
    1111
    1212      USE control_mod
     13      use dimphy
    1314      IMPLICIT none
    1415c======================================================================
    … …  
    1718c======================================================================
    1819#include "dimensions.h"
    19 #include "dimphy.h"
    2020#include "netcdf.inc"
    2121#include "dimsoil.h"
    … …  
    7575      IF( cycle_diurne ) tab_cntrl( 7 ) = 1.
    7676      IF(   soil_model ) tab_cntrl( 8 ) = 1.
    77       IF(     new_oliq ) tab_cntrl( 9 ) = 1.
    7877      IF(     ok_orodr ) tab_cntrl(10 ) = 1.
    7978      IF(     ok_orolf ) tab_cntrl(11 ) = 1.
     79      IF( ok_gw_nonoro ) tab_cntrl(12 ) = 1.
    8080
    8181      tab_cntrl(13) = day_end

  • trunk/libf/phytitan/physiq.F

    r97 r102  
    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
    … …  
    2924cAA                  -  stockage des moyennes des champs necessaires
    3025cAA                     en mode traceur off-line
    31 c======================================================================
    32 c   CLEFS CPP POUR LES IO
    33 c   =====================
    34 #define histmth
    35 #define histday
    36 #define histins
    37 c======================================================================
    3826c    modif   ( P. Le Van ,  12/10/98 )
    3927c
    … …  
    6755c d_ps----output-R-tendance physique de la pression au sol
    6856c======================================================================
     57      USE ioipsl
     58      USE histcom
     59      USE infotrac
     60      USE control_mod
     61      use dimphy
     62      USE comgeomphy
     63      IMPLICIT none
     64c======================================================================
     65c   CLEFS CPP POUR LES IO
     66c   =====================
     67#define histmth
     68#define histday
     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"
    … …  
    583582      do iq=1,nqmax
    584583         nom = tname(iq)
    585 c        print*,iq,"nom=",nom,"tname=",tname(iq),"tnom=",tnom(iq)
     584c        print*,iq,"nom=",nom,"tname=",tname(iq)
    586585         print*,iq,"nom=",nom
    587586         if (nom(1:1).eq."q") then

  • trunk/libf/phytitan/phytrac.F

    r97 r102  
    33     .                   pplev,pplay,delp,ptemp,pmu0,pfract,pdecli,
    44     .                   lonsol,tr_seri,d_tr_mph,d_tr_kim)
    5 
    6       USE infotrac
    7 
    8       IMPLICIT none
    95
    106c======================================================================
    … …  
    3228c d_tr_kim----output-R-tendance chimique de "qx" (kg/kg/s)
    3329c======================================================================
     30      USE infotrac
     31      use dimphy
     32      IMPLICIT none
    3433#include "dimensions.h"
    35 #include "dimphy.h"
    3634#include "clesphys.h"
    3735#include "YOMCST.h"

  • trunk/libf/phytitan/printflag.F

    r3 r102  
    1212
    1313       REAL tabcntr0( 100 )
    14        LOGICAL cycle_diurn0,soil_model0,new_oliq0,ok_orodr0
    15        LOGICAL ok_orolf0,ok_limitvr0
     14       LOGICAL cycle_diurn0,soil_model0,ok_orodr0
     15       LOGICAL ok_orolf0,ok_gw_nonoro0
    1616       LOGICAL ok_mensuel, ok_journe, ok_instan
    1717       INTEGER radpas , radpas0
    … …  
    3333       PRINT 100
    3434
    35 c      PRINT 11, new_oliq, ok_orodr, ok_orolf   
    36 c      PRINT 100
     35       PRINT 11, ok_orodr, ok_orolf, ok_gw_nonoro   
     36       PRINT 100
    3737
    3838       PRINT 12, nbapp_rad
    … …  
    5555        cycle_diurn0  = .FALSE.
    5656        soil_model0   = .FALSE.
    57         new_oliq0     = .FALSE.
    5857        ok_orodr0     = .FALSE.
    5958        ok_orolf0     = .FALSE.
    60         ok_limitvr0  = .FALSE.
     59        ok_gw_nonoro0 = .FALSE.
    6160
    6261        IF( tabcntr0( 7 ).EQ. 1. )   cycle_diurn0 = .TRUE.
    6362        IF( tabcntr0( 8 ).EQ. 1. )    soil_model0 = .TRUE.
    64         IF( tabcntr0( 9 ).EQ. 1. )      new_oliq0 = .TRUE.
    6563        IF( tabcntr0(10 ).EQ. 1. )      ok_orodr0 = .TRUE.
    6664        IF( tabcntr0(11 ).EQ. 1. )      ok_orolf0 = .TRUE.
    67         IF( tabcntr0(12 ).EQ. 1. )    ok_limitvr0 = .TRUE.
     65        IF( tabcntr0(12 ).EQ. 1. )  ok_gw_nonoro0 = .TRUE.
    6866
    6967        PRINT *,' $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
    … …  
    108106       ENDIF
    109107
    110        IF( new_oliq0.AND..NOT.new_oliq.OR..NOT.new_oliq0.AND.
    111      ,        new_oliq )     THEN
    112         PRINT 16, new_oliq0, new_oliq
    113         PRINT 100
    114        ENDIF
    115 
    116108       IF( ok_orodr0.AND..NOT.ok_orodr.OR..NOT.ok_orodr0.AND.
    117109     ,        ok_orodr )     THEN
    … …  
    126118       ENDIF
    127119
    128        IF( ok_limitvr0.AND..NOT.ok_limitvrai.OR..NOT.ok_limitvr0.
    129      ,     AND.ok_limitvrai )     THEN
    130         PRINT 18, ok_limitvr0, ok_limitvrai
     120       IF( ok_gw_nonoro0.AND..NOT.ok_gw_nonoro.OR..NOT.ok_gw_nonoro0.
     121     ,     AND.ok_gw_nonoro )     THEN
     122        PRINT 18, ok_gw_nonoro0, ok_gw_nonoro
    131123        PRINT 100
    132124       ENDIF
    … …  
    140132     . , ok_instan = ',l3,8x,5(1H*) )
    141133
    142  7     FORMAT(2x,5(1H*),15x,'      ok_limitvrai   = ',l3,16x,5(1h*) )
    143 
    144134 8     FORMAT(2x,'*****             radpas    =                      ' ,
    145135     , i4,6x,' *****')
    … …  
    149139
    150140
    151  11    FORMAT(2x,5(1H*),'  new_oliq = ',l3,3x,', Ok_orodr = ',
    152      , l3,3x,', Ok_orolf = ',l3,3x,5(1H*) )
     141 11    FORMAT(2x,5(1H*),', Ok_orodr = ',l3,3x,', Ok_orolf = ',l3,3x,
     142     ,    ' ok_gw_nonoro = ',l3,3x,5(1H*) )
    153143
    154144
    … …  
    165155     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
    166156
    167  16    FORMAT(2x,'$$$$$$$$   Attention !!      new_oliq  different  sur',
    168      , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
    169 
    170157 17    FORMAT(2x,'$$$$$$$$   Attention !!      ok_orolf  different  sur',
    171158     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
    172159
    173  18    FORMAT(2x,'$$$$$$$$   Attention !!  ok_limitvrai  different  sur',
     160 18    FORMAT(2x,'$$$$$$$$   Attention !!  ok_gw_nonoro  different  sur',
    174161     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
    175162

  • trunk/libf/phytitan/radlwsw.F

    r97 r102  
    66     .                  lwnet, swnet)
    77c     
    8       IMPLICIT none
    98c======================================================================
    109c 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"
    4746c
    4847c ARGUMENTS

  • trunk/libf/phytitan/radtitan.F

    r97 r102  
    11       SUBROUTINE RADTITAN(p,nq,nmicro,ycomp)
    2 
    3       IMPLICIT NONE
    4 
    52
    63c=======================================================================
    … …  
    3431c   -------------
    3532
     33      USE infotrac
     34      use dimphy
     35      USE comgeomphy
     36      IMPLICIT NONE
    3637#include "dimensions.h"
    37 #include "dimphy.h"
    3838#include "clesphys.h"
    3939#include "microtab.h"
    4040#include "numchimrad.h"
    41 #include "comgeomphy.h"
    4241#include "YOMCST.h"
    43 #include "advtrac.h"     !! pour noms des traceurs
    4442
    4543c Pour le CRAY, les block data doivent etre declares external
    … …  
    297295
    298296          do i=1,nq
    299              if (tnom(i).eq."CH4") then
     297             if (tname(i).eq."CH4") then
    300298                iradch4=i
    301              elseif (tnom(i).eq."C2H2") then
     299             elseif (tname(i).eq."C2H2") then
    302300                iradc2h2=i
    303              elseif (tnom(i).eq."C2H6") then
     301             elseif (tname(i).eq."C2H6") then
    304302                iradc2h6=i
    305              elseif (tnom(i).eq."HCN") then
     303             elseif (tname(i).eq."HCN") then
    306304                iradhcn=i
    307              elseif (tnom(i).eq."N2") then
     305             elseif (tname(i).eq."N2") then
    308306                iradn2=i
    309              elseif (tnom(i).eq."H2") then
     307             elseif (tname(i).eq."H2") then
    310308                iradh2=i
    311309             endif

  • trunk/libf/phytitan/readstartphy.F

    r3 r102  
    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/phytitan/sfluxv.F

    r3 r102  
    11      SUBROUTINE SFLUXV(IPRINT,IG,dist_sol)
    22
     3      use dimphy
    34      IMPLICIT NONE
    4 
    55#include "dimensions.h"
    6 #include "dimphy.h"
    76#include "comorbit.h"
    87

  • trunk/libf/phytitan/soil.F

    r3 r102  
    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/phytitan/ustarhb.F

    r3 r102  
    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/phytitan/vdif_kcay.F

    r3 r102  
    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/phytitan/writerestartphy.F

    r3 r102  
    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/phytitan/yamada.F

    r3 r102  
    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/phytitan/yamada4.F

    r3 r102  
    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/phytitan/zenang.F

    r3 r102  
    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/RDAY)
    2625c================================================================
     26      use dimphy
     27      IMPLICIT none
    2728#include "dimensions.h"
    28 #include "dimphy.h"
    2929#include "YOMCST.h"
    3030#include "comorbit.h"

  • trunk/libf/phyvenus/iniphysiq.F

    r101 r102  
    1010     $           prad,pg,pr,pcpp)
    1111
    12       use dimphy
    13       USE comgeomphy
    14       IMPLICIT NONE
    1512c
    1613c=======================================================================
    … …  
    4643c   -------------
    4744 
     45      use dimphy
     46      USE comgeomphy
     47      IMPLICIT NONE
    4848#include "dimensions.h"
    4949

  • trunk/libf/phyvenus/printflag.F

    r101 r102  
    119119
    120120
    121  11    FORMAT(2x,5(1H*),', Ok_orodr = ',
    122      , l3,3x,', Ok_orolf = ',l3,3x,5(1H*) )
     121 11    FORMAT(2x,5(1H*),', Ok_orodr = ',l3,3x,', Ok_orolf = ',l3,3x,
     122     ,    ' ok_gw_nonoro = ',l3,3x,5(1H*) )
    123123
    124124
Note: See TracChangeset for help on using the changeset viewer.