Changeset 119


Ignore:
Timestamp:
May 17, 2011, 10:56:57 AM (14 years ago)
Author:
slebonnois
Message:

Sebastien Lebonnois: apres validation des versions Venus et Titan,
correction d'un certain nombre de bugs.

Location:
trunk
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • trunk/chantiers/commit_importants.log

    r117 r119  
    874874Petites MAJ routines graphiques
    875875
     876*********************
     877**** commit_v119 ****
     878*********************
     879
     880Sebastien Lebonnois: apres validation des versions Venus et Titan,
     881correction d'un certain nombre de bugs.
     882

  • trunk/chantiers/compilation.txt

    r37 r119  
    1616   (qui me permettent de travailler sur Bellonzi avec IOIPSL  et NETCDF
    1717    déjà là)
     18
     19Seb (02/05/2011):
     20L'outil makelmdz marche, adapte egalement pour Titan.
     21Compils ok.
     22

  • trunk/deftanks/titan/z2sig.def

    r5 r119  
    11 1.00000      H: atmospheric scale height (km) (used as a reference only)
    2   0.00288765924   Typical log(sigma) for 1st layer (z=H*log(sigma))
    3   0.0143499468   Typical log(sigma) for 2nd layer, etc...
    4   0.0368630514
    5   0.0749709532
    6   0.132121697
    7   0.210579455
     2  0.00100000     Typical log(sigma) for 1st layer (z=H*log(sigma))
     3  0.00600000     Typical log(sigma) for 2nd layer, etc...
     4  0.0160000000
     5  0.0400000000
     6  0.090000000
     7  0.180000000
    88  0.311430246
    99  0.434752703

  • trunk/libf/dyn3d/conf_gcm.F

    r108 r119  
    376376       iphysiq = 5
    377377       CALL getin('iphysiq', iphysiq)
     378
     379!Config  Key  = iflag_trac
     380!Config  Desc = traceurs presents ou non
     381!Config  Def  = 1
     382!Config  Help = Permet de faire tourner le modele sans traceurs
     383!Config         
     384       iflag_trac = 1
     385       CALL getin('iflag_trac',iflag_trac)
    378386
    379387!Config  Key  = ip_ebil_dyn
    … …  
    657665      write(lunout,*)' iflag_phys = ', iflag_phys
    658666      write(lunout,*)' iphysiq = ', iphysiq
     667      write(lunout,*)' iflag_trac = ', iflag_trac
    659668      write(lunout,*)' clonn = ', clonn
    660669      write(lunout,*)' clatt = ', clatt
    … …  
    894903      write(lunout,*)' iflag_phys = ', iflag_phys
    895904      write(lunout,*)' iphysiq = ', iphysiq
     905      write(lunout,*)' iflag_trac = ', iflag_trac
    896906      write(lunout,*)' clon = ', clon
    897907      write(lunout,*)' clat = ', clat

  • trunk/libf/dyn3d/leapfrog.F

    r110 r119  
    448448         IF (first) THEN
    449449#include "ini_paramLMDZ_dyn.h"
     450          first=.false.
    450451         ENDIF
    451452c
    … …  
    863864      END IF ! of IF(.not.purmats)
    864865
    865       first=.false.
    866 
    867866      STOP
    868867      END

  • trunk/libf/dyn3d/logic.h

    r101 r119  
    1515     & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus                      &
    1616     &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
    17      &  ,ok_limit,ok_etat0,iflag_trac
     17     &  ,ok_limit,ok_etat0
    1818
    19       INTEGER iflag_phys
     19      INTEGER iflag_phys,iflag_trac
    2020!-----------------------------------------------------------------------

  • trunk/libf/dyn3dpar/conf_gcm.F

    r108 r119  
    388388       iphysiq = 5
    389389       CALL getin('iphysiq', iphysiq)
     390
     391!Config  Key  = iflag_trac
     392!Config  Desc = traceurs presents ou non
     393!Config  Def  = 1
     394!Config  Help = Permet de faire tourner le modele sans traceurs
     395!Config         
     396       iflag_trac = 1
     397       CALL getin('iflag_trac',iflag_trac)
    390398
    391399!Config  Key  = ip_ebil_dyn
    … …  
    679687      write(lunout,*)' iflag_phys = ', iflag_phys
    680688      write(lunout,*)' iphysiq = ', iphysiq
     689      write(lunout,*)' iflag_trac = ', iflag_trac
    681690      write(lunout,*)' clonn = ', clonn
    682691      write(lunout,*)' clatt = ', clatt
    … …  
    947956      write(lunout,*)' iflag_phys = ', iflag_phys
    948957      write(lunout,*)' iphysiq = ', iphysiq
     958      write(lunout,*)' iflag_trac = ', iflag_trac
    949959      write(lunout,*)' clon = ', clon
    950960      write(lunout,*)' clat = ', clat

  • trunk/libf/dyn3dpar/logic.h

    r101 r119  
    1515     & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus                      &
    1616     &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
    17      &  ,ok_limit,ok_etat0,iflag_trac
     17     &  ,ok_limit,ok_etat0
    1818
    19       INTEGER iflag_phys
     19      INTEGER iflag_phys,iflag_trac
    2020!$OMP THREADPRIVATE(/logic/)
    2121!-----------------------------------------------------------------------

  • trunk/libf/phytitan/clcdrag.F90

    r97 r119  
    3939! Variables locales :
    4040      INTEGER :: i
    41       REAL :: zdu2, ztsolv, ztvd, zscf
    42       REAL :: zucf, zcr
    43       REAL :: friv, frih
     41      REAL :: zscf
     42      REAL :: zucf
     43      REAL :: FRIV,FRIH
    4444      REAL, dimension(klon) :: zcfm1, zcfm2
    4545      REAL, dimension(klon) :: zcfh1, zcfh2
    … …  
    9393            pcfh(i) = zcdn(i)* fins(zri(i))
    9494          ENDIF
    95           zcr = (0.0016/(zcdn(i)*SQRT(zdu2)))*ABS(ztvd-ztsolv)**(1./3.)
    9695        ENDIF
    9796      END DO

  • trunk/libf/phytitan/init_phys_lmdz.F90

    r102 r119  
    1515
    1616    CALL init_grid_phy_lmdz(iim,jjp1,llm)
    17     CALL init_dimphy(klon,nbp_lev)
     17    CALL init_dimphy(klon_glo,nbp_lev)
    1818 
    1919END SUBROUTINE init_phys_lmdz 

  • trunk/libf/phytitan/optci.F

    r104 r119  
    11      SUBROUTINE OPTCI(ykim,qaer,nmicro,IPRINT)
    22      use dimphy
     3      use infotrac
    34#include "dimensions.h"
    45#include "microtab.h"
    … …  
    910c   ---------
    1011      REAL    ykim(klon,klev,nqtot)
    11       real    qaer(klon,nlayer,nqtot)
     12      real    qaer(klon,klev,nqtot)
    1213      integer nmicro
    1314c   ---------
    … …  
    198199        if (microfi.eq.1) then
    199200           do iq=1,nrad
    200               do j=1,NLAYER
    201                  zqaer_1pt(j,iq)=qaer(ig,j,iq)
    202               enddo
    203            enddo
     201             do j=1,NLAYER
     202                 zqaer_1pt(j,iq)=qaer(ig,j,iq)
     203             enddo
     204           enddo
    204205        else
    205206         if (ig.eq.1)  then
    … …  
    260261      RETURN
    261262      END
     263

  • trunk/libf/phytitan/optcv.F

    r104 r119  
    11      SUBROUTINE OPTCV(qaer,nmicro,IPRINT)
    22
    3 
    43      use dimphy
     4      use infotrac
    55#include "dimensions.h"
    66#include "microtab.h"
    … …  
    99c   Argument:
    1010c   ---------
    11       REAL    qaer(klon,nlayer,nqtot)
     11      REAL    qaer(klon,klev,nqtot)
    1212      integer nmicro
    1313c   ---------

  • trunk/libf/phytitan/suphec.F

    r106 r119  
    5858c     RSIDAY=1.37889e6   ! 15.96 j
    5959      RDAY=1.37889e6   ! 15.96 j
    60       ROMEGA=2.*RPI/RSIDAY
    6160c 1/(duree du jour) = 1/(periode rotation) - 1/(periode revolution)
    6261c     RDAY=RSIDAY/(1.-RSIDAY/RSIYEA)
    6362      RSIDAY=RDAY/(1.+RDAY/RSIYEA)
     63      ROMEGA=2.*RPI/RSIDAY
    6464      REA=1.5e12
    6565      REPSM=0.  ! 0. veut dire qu'on commence au point vernal

  • trunk/libf/phytitan/write_histmth.h

    r110 r119  
    156156      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    157157      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    158       CALL histwrite(nid_day,"thv"//str1,itau_w,zx_tmp_3d,
     158      CALL histwrite(nid_mth,"thv"//str1,itau_w,zx_tmp_3d,
    159159     .                                   iim*jjmp1*klev,ndex3d)
    160160       enddo      ! fin boucle NSPECV
    … …  
    176176      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    177177      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    178       CALL histwrite(nid_day,"khv"//str1,itau_w,zx_tmp_3d,
     178      CALL histwrite(nid_mth,"khv"//str1,itau_w,zx_tmp_3d,
    179179     .                                   iim*jjmp1*klev,ndex3d)
    180180       enddo      ! fin boucle NSPECV
    … …  
    189189      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    190190      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    191       CALL histwrite(nid_day,"tgv"//str1,itau_w,zx_tmp_3d,
     191      CALL histwrite(nid_mth,"tgv"//str1,itau_w,zx_tmp_3d,
    192192     .                                   iim*jjmp1*klev,ndex3d)
    193193       enddo      ! fin boucle NSPECV
    … …  
    209209      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    210210      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    211       CALL histwrite(nid_day,"kgv"//str1,itau_w,zx_tmp_3d,
     211      CALL histwrite(nid_mth,"kgv"//str1,itau_w,zx_tmp_3d,
    212212     .                                   iim*jjmp1*klev,ndex3d)
    213213       enddo      ! fin boucle NSPECV
    … …  
    222222      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    223223      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    224       CALL histwrite(nid_day,"thi"//str1,itau_w,zx_tmp_3d,
     224      CALL histwrite(nid_mth,"thi"//str1,itau_w,zx_tmp_3d,
    225225     .                                   iim*jjmp1*klev,ndex3d)
    226226       enddo      ! fin boucle NSPECI
    … …  
    242242      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    243243      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    244       CALL histwrite(nid_day,"khi"//str1,itau_w,zx_tmp_3d,
     244      CALL histwrite(nid_mth,"khi"//str1,itau_w,zx_tmp_3d,
    245245     .                                   iim*jjmp1*klev,ndex3d)
    246246       enddo      ! fin boucle NSPECI
    … …  
    255255      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    256256      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    257       CALL histwrite(nid_day,"tgi"//str1,itau_w,zx_tmp_3d,
     257      CALL histwrite(nid_mth,"tgi"//str1,itau_w,zx_tmp_3d,
    258258     .                                   iim*jjmp1*klev,ndex3d)
    259259       enddo      ! fin boucle NSPECI
    … …  
    275275      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    276276      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    277       CALL histwrite(nid_day,"kgi"//str1,itau_w,zx_tmp_3d,
     277      CALL histwrite(nid_mth,"kgi"//str1,itau_w,zx_tmp_3d,
    278278     .                                   iim*jjmp1*klev,ndex3d)
    279279       enddo      ! fin boucle NSPECI

  • trunk/libf/phytitan/zenang.F

    r102 r119  
    150150      ENDIF   !---comparaison omega1 et omega2
    151151c
     152c Petit test rajoute pour les cas pathologiques aux poles
     153      if (pmu0(i).lt.0.) pmu0(i) = -1.*pmu0(i)
     154
    152155      ENDDO
    153156c

  • trunk/libf/phyvenus/init_phys_lmdz.F90

    r101 r119  
    1515
    1616    CALL init_grid_phy_lmdz(iim,jjp1,llm)
    17     CALL init_dimphy(klon,nbp_lev)
     17    CALL init_dimphy(klon_glo,nbp_lev)
    1818 
    1919END SUBROUTINE init_phys_lmdz 

  • trunk/libf/phyvenus/physiq.F

    r105 r119  
    390390      ok_sync=.TRUE.
    391391
    392       bilansmc = 1
     392      bilansmc = 0
    393393      ballons  = 0
    394394
Note: See TracChangeset for help on using the changeset viewer.