Ignore:
Timestamp:
Nov 6, 2006, 4:51:16 PM (18 years ago)
Author:
Laurent Fairhead
Message:

Modifications pour rendre INCA plus independant de LMDZ ACo
LF

Location:
LMDZ4/branches/V3_test/libf/dyn3dpar
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/V3_test/libf/dyn3dpar/advtrac.h

    r709 r726  
    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      integer nbtrac
     19      integer nprath
     20      real    mmt_adj(iim+1,jjm+1,llm, 1)
     21      integer hadv_flg(nqmx)
     22      integer vadv_flg(nqmx)
     23      integer conv_flg(nqmx-2)
     24      integer pbl_flg(nqmx-2)
     25      character*8 tracnam(nqmx-2)
    1526c-----------------------------------------------------------------------

  • LMDZ4/branches/V3_test/libf/dyn3dpar/caladvtrac_p.F

    r709 r726  
    99     *                   flxw,
    1010     *                   pk,
    11      *                   mmt_adj,
    12      *                   hadv_flg,iapptrac)
     11     *                   iapptrac)
    1312#else
    1413            SUBROUTINE caladvtrac_p(q,pbaru,pbarv ,
    … …  
    4140      REAL p( ip1jmp1,llmp1),q( ip1jmp1,llm,nqmx),dq( ip1jmp1,llm,2 )
    4241      REAL teta( ip1jmp1,llm),pk( ip1jmp1,llm)
    43 #ifdef INCA_CH4
    44 cym      INTEGER            :: hadv_flg(nq)
    45       INTEGER            :: hadv_flg(nqmx)
    46       REAL               :: mmt_adj(iip1,jjp1,llm)
     42#ifdef INCA
    4743      REAL               :: flxw(ip1jmp1,llm)
    4844#endif
    … …  
    7874c      print *,'appel a advtrac'
    7975
    80 #ifdef INCA_CH4
     76#ifdef INCA
    8177      CALL advtrac_p( pbaru,pbarv,
    8278     *             p,  masse,q,iapptrac, teta,
    8379     .             flxw,
    84      .             pk,
    85      .             mmt_adj,
    86      .             hadv_flg)
     80     .             pk)
    8781#else
    8882      CALL advtrac_p( pbaru,pbarv,

  • LMDZ4/branches/V3_test/libf/dyn3dpar/gcm.F

    r709 r726  
    1414      USE mod_hallo
    1515      USE Bands
    16 #ifdef INCA
    17       USE inca_dim
    18 #endif
    1916      IMPLICIT NONE
    2017
    … …  
    6461#include "iniprint.h"
    6562#include "tracstoke.h"
    66 
     63#include "advtrac.h"
    6764
    6865      INTEGER         longcles
    … …  
    157154      dynhist_file = 'dyn_hist'
    158155      dynhistave_file = 'dyn_hist_ave'
     156
     157
     158c initialisation Anne
     159      hadv_flg(:) = 0.
     160      vadv_flg(:) = 0.
     161      conv_flg(:) = 0.
     162      pbl_flg(:)  = 0.
     163      tracnam(:)  = '        '
     164      nprath = 1
     165      nbtrac = 0
     166      mmt_adj(:,:,:,:) = 1
     167
    159168
    160169c--------------------------------------------------------------------------
    … …  
    221230c$OMP END PARALLEL
    222231#ifdef INCA
    223       call init_inca_dim
     232      call init_const_lmdz(nbtrac,anneeref,dayref,iphysiq,day_step,nday)
     233      call init_inca_para(iim,jjm+1,klon2,phy_size,klon_para_nb)
    224234#endif
    225235
    … …  
    250260      endif
    251261
     262#ifdef INCA
     263      call init_inca_dim(klon,llm,iim,jjm,
     264     $     rlonu,rlatu,rlonv,rlatv)
     265#endif
    252266
    253267

  • LMDZ4/branches/V3_test/libf/dyn3dpar/iniadvtrac.F

    r709 r726  
    66      subroutine iniadvtrac(nq)
    77      USE ioipsl
    8 #ifdef INCA
    9       USE transport_controls, only : hadv_flg, vadv_flg
    10 cym      USE chemshut
    11       USE species_names
    12 #endif
    138      IMPLICIT NONE
    149c=======================================================================
    … …  
    6055      descrq(30)='PRA'
    6156
     57#ifdef INCA
     58
     59      CALL init_transport(
     60     $     hadv_flg,
     61     $     vadv_flg,
     62     $     conv_flg,
     63     $     pbl_flg,
     64     $     tracnam)
     65#endif
     66
    6267c-----------------------------------------------------------------------
    6368c        Choix  des schemas d'advection pour l'eau et les traceurs
    … …  
    110115      tnom(2)='H2Ol'
    111116      nq=nbtrac+2
     117     
    112118       if (nq.gt.nqmx) then
    113        print*,'nombre de traceurs incompatible INCA/LMDZT'
     119       print*,'nombre de traceurs incompatible INCA/LMDZT', nq, nbtrac
    114120       stop
    115121       endif

  • LMDZ4/branches/V3_test/libf/dyn3dpar/leapfrog_p.F

    r709 r726  
    1919       USE Write_Field_p
    2020       USE vampir
    21        
    22 #ifdef INCA
    23       USE transport_controls, ONLY : hadv_flg, mmt_adj
    24 #endif
    2521
    2622      IMPLICIT NONE
    … …  
    7672#include "academic.h"
    7773#include "clesphys.h"
    78 
     74#include "advtrac.h"
    7975     
    8076      include 'mpif.h'
    … …  
    148144
    149145      REAL alpha(ip1jmp1,llm),beta(ip1jmp1,llm)
    150 #ifdef INCA_CH4
     146#ifdef INCA
    151147      REAL :: flxw(ip1jmp1,llm)
    152148#endif
    … …  
    528524c$OMP PARALLEL DEFAULT(SHARED)
    529525c
    530 #ifdef INCA_CH4
     526#ifdef INCA
    531527             CALL caladvtrac_p(q,pbaru,pbarv,
    532528     *                      p, masse, dq,  teta,
    533529     .             flxw,
    534      .             pk,
    535      .             mmt_adj,
    536      .             hadv_flg,iapptrac)
     530     .             pk,
     531     .             iapptrac)
    537532#else
    538533             CALL caladvtrac_p(q,pbaru,pbarv,
    … …  
    680675     *                               jj_Nb_physic,2,2,Request_physic)
    681676        enddo
    682 #ifdef INCA_CH4
     677#ifdef INCA
    683678        call Register_SwapFieldHallo(flxw,flxw,ip1jmp1,llm,
    684679     *                               jj_Nb_physic,2,2,Request_physic)
    … …  
    698693     $               ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
    699694     $               du,dv,dteta,dq,w,
    700 #ifdef INCA_CH4
     695#ifdef INCA
    701696     $               flxw,
    702697#endif

  • LMDZ4/branches/V3_test/libf/dyn3dpar/mod_hallo.F90

    r718 r726  
    33implicit none
    44!  include 'mpif.h'
    5   integer, parameter :: MaxRequest=80
     5  integer, parameter :: MaxRequest=200
    66  integer, parameter :: MaxProc=80
    77  integer, parameter :: MaxBufferSize=1024*1024*16
    … …  
    430430      type(request_SR),pointer :: Req
    431431      type(Hallo),pointer :: PtrHallo
    432       integer, dimension(2*mpi_size) :: TabRequest
    433       integer, dimension(MPI_STATUS_SIZE,2*mpi_size) :: TabStatus
     432      integer, dimension(4) :: TabRequest
     433      integer, dimension(MPI_STATUS_SIZE,4) :: TabStatus
    434434      integer :: NbRequest
    435435      integer :: i,rank,pos,ij,l,ierr
    … …  
    512512      type(request_SR),pointer :: Req
    513513      type(Hallo),pointer :: PtrHallo
    514       integer, dimension(mpi_size) :: TabRequest
    515       integer, dimension(MPI_STATUS_SIZE,mpi_size) :: TabStatus
     514      integer, dimension(4) :: TabRequest
     515      integer, dimension(MPI_STATUS_SIZE,4) :: TabStatus
    516516      integer :: NbRequest
    517517      integer :: i,rank,pos,ij,l,ierr
    … …  
    553553      type(request_SR),pointer :: Req
    554554      type(Hallo),pointer :: PtrHallo
    555       integer, dimension(mpi_size) :: TabRequest
    556       integer, dimension(MPI_STATUS_SIZE,mpi_size) :: TabStatus
     555      integer, dimension(4) :: TabRequest
     556      integer, dimension(MPI_STATUS_SIZE,4) :: TabStatus
    557557      integer :: NbRequest
    558558      integer :: i,rank,pos,ij,l,ierr
Note: See TracChangeset for help on using the changeset viewer.