Changeset 762 for LMDZ4/trunk


Ignore:
Timestamp:
Jun 4, 2007, 3:34:26 PM (17 years ago)
Author:
Laurent Fairhead
Message:

Merge entre la version V3_conv et le HEAD
YM, JG, LF

Location:
LMDZ4/trunk/libf/dyn3d
Files:
11 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/dyn3d/advtrac.F

    r616 r762  
    88     *                   p,  masse,q,iapptrac,teta,
    99     *                  flxw,
    10      *                  pk,
    11      *                  mmt_adj,
    12      *                  hadv_flg)
     10     *                  pk)
    1311#else
    1412      SUBROUTINE advtrac(pbaru,pbarv ,
     
    5149      REAL pk(ip1jmp1,llm)
    5250#ifdef INCA
    53       INTEGER            :: hadv_flg(nqmx)
    54       REAL               :: mmt_adj(ip1jmp1,llm,1)
    5551      REAL               :: flxw(ip1jmp1,llm)
    5652#endif
     
    215211#ifdef INCA
    216212       do iiq = iq+1, iq+3
    217          q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
     213c         q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
     214          q(:,:,iiq)=q(:,:,iiq)*1
    218215       enddo
    219216#endif
     
    233230#ifdef INCA
    234231       do iiq = iq+1, iq+9
    235          q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
     232c         q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
     233         q(:,:,iiq)=q(:,:,iiq)*1
    236234       enddo
    237235#endif
  • LMDZ4/trunk/libf/dyn3d/advtrac.h

    r524 r762  
    55c INCLUDE 'advtrac.h'
    66
    7       COMMON/advtr/iadv,hadv,vadv,tnom,tname,ttext,niadv
     7      COMMON/advtr/iadv,hadv,vadv,tnom,tname,ttext,niadv,
     8     &     nbtrac, nprath, mmt_adj, hadv_flg, vadv_flg, conv_flg,
     9     &     pbl_flg, tracnam
    810      INTEGER iadv(nqmx) ! indice schema de transport
    911      INTEGER hadv(nqmx) ! indice schema transport horizontal
     
    1315      character*10 tname(nqmx) ! nom du traceur pour restart
    1416      character*13 ttext(nqmx) ! nom long du traceur pour sorties
     17
     18
     19      integer nbtrac
     20      integer nprath
     21      real    mmt_adj(iim+1,jjm+1,llm, 1)
     22      integer hadv_flg(nqmx)
     23      integer vadv_flg(nqmx)
     24      integer conv_flg(nqmx-2)
     25      integer pbl_flg(nqmx-2)
     26      character*8 tracnam(nqmx-2)
    1527c-----------------------------------------------------------------------
  • LMDZ4/trunk/libf/dyn3d/caladvtrac.F

    r616 r762  
    88     *                   p ,masse, dq ,  teta,
    99     *                   flxw,
    10      *                   pk,
    11      *                   mmt_adj,
    12      *                   hadv_flg)
     10     *                   pk)
    1311#else
    1412            SUBROUTINE caladvtrac(q,pbaru,pbarv ,
     
    4240      REAL teta( ip1jmp1,llm),pk( ip1jmp1,llm)
    4341#ifdef INCA
    44       INTEGER            :: hadv_flg(nqmx)
    45       REAL               :: mmt_adj(iip1,jjp1,llm,1)
    4642      REAL               :: flxw(ip1jmp1,llm)
    4743#endif
     
    7773     *             p,  masse,q,iapptrac, teta,
    7874     .             flxw,
    79      .             pk,
    80      .             mmt_adj,
    81      .             hadv_flg)
     75     .             pk)
    8276#else
    8377      CALL advtrac( pbaru,pbarv,
  • LMDZ4/trunk/libf/dyn3d/create_etat0_limit.F

    r524 r762  
    33!
    44       PROGRAM create_etat0_limit
     5       USE dimphy
     6       USE comgeomphy
     7
    58c
    69c
     
    1821c
    1922c     on cree le masque dans etat0 que l'on passe ensuite dans limit pour
    20 c     garder les cohérences
     23c     garder les coherences
    2124
    2225      LOGICAL interbar, extrap , oldice
     
    2528#include "paramet.h"
    2629#include "indicesol.h"
    27 #include "dimphy.h"
     30cym#include "dimphy.h"
    2831      REAL :: masque(iip1,jjp1)
    29       REAL :: pctsrf(klon, nbsrf)
     32      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)
    3033
     34     
     35      call InitDimphy
     36      call init_phys_openmp
     37      call InitComgeomphy
     38     
    3139      WRITE(6,*) '  *********************  '
    3240      WRITE(6,*) ' interbar = ',interbar
     
    3846      WRITE(6,*) '  *********************  '
    3947      WRITE(6,1)
    40 c
     48     
     49c     
    4150      CALL  limit_netcdf ( interbar, extrap , oldice, masque, pctsrf )
    4251
  • LMDZ4/trunk/libf/dyn3d/etat0_netcdf.F

    r724 r762  
    88      USE startvar
    99      USE ioipsl
     10      USE dimphy
    1011      !
    1112      IMPLICIT NONE
     
    2324#include "comconst.h"
    2425#include "indicesol.h"
    25 #include "dimphy.h"
     26c#include "dimphy.h"
    2627#include "dimsoil.h"
    2728#include "temps.h"
     
    510511C verif que somme des sous surface = 1
    511512C
    512       ji=count(abs(sum(pctsrf(1 : klon, 1 : nbsrf), dim = 2) - 1.0 )
     513      ji=count( (abs( sum(pctsrf(1 : klon, 1 : nbsrf),dim=2))-1.0)
    513514     $    .GT. EPSFRA)
    514515      IF (ji .NE. 0) THEN
     
    671672      rnebcon = 0.0
    672673      ratqs = 0.0
    673       run_off_lic_0 = 0.0
     674      run_off_lic_0 = 0.0 
    674675
    675676cIM   call phyredem("startphy.nc",phystep,radpas, co2_ppm, solaire,
     677   
    676678      call phyredem("startphy.nc",phystep,radpas,
    677679     $    latfi, lonfi, pctsrf, tsolsrf, tsoil, tslab, seaice,
     
    683685     $    t_ancien, q_ancien, rnebcon, ratqs, clwcon,
    684686     $    run_off_lic_0)
    685       print*,'sortie phyredem'
    686687
    687688C     Sortie Visu pour les champs dynamiques
  • LMDZ4/trunk/libf/dyn3d/gcm.F

    r691 r762  
    99      USE IOIPSL
    1010#endif
    11 
     11      USE dimphy
     12      USE comgeomphy
    1213      IMPLICIT NONE
    1314
     
    5758#include "iniprint.h"
    5859#include "tracstoke.h"
    59 
     60#include "advtrac.h"
    6061
    6162      INTEGER         longcles
     
    144145      dynhistave_file = 'dyn_hist_ave.nc'
    145146
     147c initialisation Anne
     148      hadv_flg(:) = 0.
     149      vadv_flg(:) = 0.
     150      conv_flg(:) = 0.
     151      pbl_flg(:)  = 0.
     152      tracnam(:)  = '        '
     153      nprath = 1
     154      nbtrac = 0
     155      mmt_adj(:,:,:,:) = 1
     156
     157
    146158c--------------------------------------------------------------------------
    147159c   Iflag_phys controle l'appel a la physique :
     
    183195      CALL defrun( 99, .TRUE. , clesphy0 )
    184196#endif
     197
     198      call InitDimphy
     199      call init_phys_openmp
     200      call InitComgeomphy
     201
     202#ifdef INCA
     203      call init_const_lmdz(nbtrac,anneeref,dayref,iphysiq,day_step,nday)
     204      call init_inca_para(iim,jjm+1,klon2,phy_size,klon_para_nb)
     205#endif
     206
     207c
     208c
     209c------------------------------------
     210c   Initialisation partie parallele
     211c------------------------------------
     212
    185213c
    186214c
     
    211239      endif
    212240
     241#ifdef INCA
     242      call init_inca_dim(klon,llm,iim,jjm,
     243     $     rlonu,rlatu,rlonv,rlatv)
     244#endif
    213245
    214246
  • LMDZ4/trunk/libf/dyn3d/guide.F

    r644 r762  
    367367                if (first.and.ini_anal) vcov(ij,l)=a
    368368            enddo
    369             if (first.and.ini_anal) vcov(ij,l)=a
    370369         enddo
    371370      endif
  • LMDZ4/trunk/libf/dyn3d/iniadvtrac.F

    r616 r762  
    66      subroutine iniadvtrac(nq)
    77      USE ioipsl
    8 #ifdef INCA
    9       USE transport_controls, only : hadv_flg, vadv_flg
    10       USE species_names
    11       USE chemshut
    12 #endif
    138      IMPLICIT NONE
    149c=======================================================================
     
    5954      descrq(20)='SLP'
    6055      descrq(30)='PRA'
     56
     57#ifdef INCA
     58
     59      CALL init_transport(
     60     $     hadv_flg,
     61     $     vadv_flg,
     62     $     conv_flg,
     63     $     pbl_flg,
     64     $     tracnam)
     65#endif
    6166
    6267c-----------------------------------------------------------------------
  • LMDZ4/trunk/libf/dyn3d/leapfrog.F

    r692 r762  
    77     &                    time_0)
    88
    9 #ifdef INCA
    10       USE transport_controls, ONLY : hadv_flg, mmt_adj
    11 #endif
    129
    1310cIM : pour sortir les param. du modele dans un fis. netcdf 110106
     
    6158#include "com_io_dyn.h"
    6259#include "iniprint.h"
    63 
     60#include "advtrac.h"
    6461c#include "tracstoke.h"
    6562
     
    297294     *                      p, masse, dq,  teta,
    298295     .             flxw,
    299      .             pk,
    300      .             mmt_adj,
    301      .             hadv_flg)
     296     .             pk)
    302297#else
    303298             CALL caladvtrac(q,pbaru,pbarv,
     
    515510
    516511            IF( itau. EQ. itaufinp1 ) then 
    517 c$$$       write(79,*) 'ucov',ucov
    518 c$$$       write(80,*) 'vcov',vcov
    519 c$$$       write(81,*) 'teta',teta
    520 c$$$       write(82,*) 'ps',ps
    521 c$$$       write(83,*) 'q',q
    522 c$$$       WRITE(85,*) 'q1 = ',q(:,:,1)
    523 c$$$       WRITE(86,*) 'q3 = ',q(:,:,3)
     512              write(79,*) 'ucov',ucov
     513              write(80,*) 'vcov',vcov
     514              write(81,*) 'teta',teta
     515              write(82,*) 'ps',ps
     516              write(83,*) 'q',q
     517              WRITE(85,*) 'q1 = ',q(:,:,1)
     518              WRITE(86,*) 'q3 = ',q(:,:,3)
    524519
    525520              abort_message = 'Simulation finished'
  • LMDZ4/trunk/libf/dyn3d/limit_netcdf.F

    r677 r762  
    55C
    66      SUBROUTINE limit_netcdf(interbar, extrap, oldice, masque, pctsrf)
    7 c
     7      USE dimphy
    88      IMPLICIT none
    99c
     
    3030#include "comgeom2.h"
    3131#include "comconst.h"
    32 #include "dimphy.h"
     32cy#include "dimphy.h"
    3333#include "indicesol.h"
    3434c
  • LMDZ4/trunk/libf/dyn3d/write_paramLMDZ_dyn.h

    r721 r762  
    120120     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    121121c
    122       zx_tmp_2d(1:iip1,1:jjp1)=ecritphy
     122      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(ecritphy)
    123123      CALL histwrite(nid_ctesGCM, "ecritphy", itau_dyn,
    124124     .               zx_tmp_2d,iip1*jjp1,ndex2d)
Note: See TracChangeset for help on using the changeset viewer.