Changeset 634 for LMDZ4/branches/LMDZ4_par_0/libf/phylmd/phytrac.F

Timestamp:

May 4, 2005, 5:11:29 PM (19 years ago)

Author:

Laurent Fairhead

Message:

Modifications faites à la physique pour la rendre parallele YM
Une branche de travail LMDZ4_par_0 a été créée provisoirement afin de tester
les modifs pleinement avant leurs inclusions dans le tronc principal
LF

File:

• LMDZ4/branches/LMDZ4_par_0/libf/phylmd/phytrac.F (modified) (12 diffs)

LMDZ4/branches/LMDZ4_par_0/libf/phylmd/phytrac.F

@@ -64 +64 @@
 
       USE ioipsl
-
+      USE dimphy
+      USE comgeomphy
+      USE iophy
+#ifdef CPP_PARALLEL
+      USE vampir
+#endif
 #ifdef INCA
       USE sflx
@@ -76 +81 @@
 #ifdef INCA_AER
       USE AEROSOL_MOD, only : ntr,trmx,trnx
+cym      USE AEROSOL_DIAG, only : cla,las,tausum,angst,aload,scon
+cym     .                            ,scavcoef_st,scavcoef_cv
+cym      USE AEROSOL_PROGNOS, ONLY : md
       USE AEROSOL_DIAG,only : cla,las,tausum,angst,aload,cload,totaerh2o,tau,
      $  emiss20,sconc,scavcoef_st,scavcoef_cv
@@ -103 +111 @@
 #include "YOMCST.h"
 #include "dimensions.h"
-#include "dimphy.h"
+cym#include "dimphy.h"
 #include "indicesol.h"
 #include "temps.h"
 #include "paramet.h"
 #include "control.h"
-#include "comgeomphy.h"
+cym#include "comgeomphy.h"
 #include "advtrac.h"
 c======================================================================
@@ -225 +233 @@
 cAA ................................................
  
-      REAL trs(klon,nbtr)    ! Conc. radon ds le sol
-      SAVE trs
-
-      REAL masktr(klon,nbtr) ! Masque reservoir de sol traceur
+      REAL,save,allocatable :: trs(:,:)    ! Conc. radon ds le sol
+cym      SAVE trs
+      REAL :: trs_tmp(klon2)
+      
+      REAL,save,allocatable :: masktr(:,:) ! Masque reservoir de sol traceur
 c                            Masque de l'echange avec la surface
 c                           (1 = reservoir) ou (possible => 1 )
-      SAVE masktr
-      REAL fshtr(klon,nbtr)  ! Flux surfacique dans le reservoir de sol
-      SAVE fshtr
-      REAL hsoltr(nbtr)      ! Epaisseur equivalente du reservoir de sol
-      SAVE hsoltr
-      REAL tautr(nbtr)       ! Constante de decroissance radioactive
-      SAVE tautr
-      REAL vdeptr(nbtr)      ! Vitesse de depot sec dans la couche Brownienne
-      SAVE vdeptr
-      REAL scavtr(nbtr)      ! Coefficient de lessivage
-      SAVE scavtr
+cym      SAVE masktr
+      REAL,save,allocatable :: fshtr(:,:)  ! Flux surfacique dans le reservoir de sol
+cym      SAVE fshtr
+      REAL,save,allocatable :: hsoltr(:)      ! Epaisseur equivalente du reservoir de sol
+cym      SAVE hsoltr
+      REAL,save,allocatable :: tautr(:)       ! Constante de decroissance radioactive
+cym      SAVE tautr
+      REAL,save,allocatable :: vdeptr(:)      ! Vitesse de depot sec dans la couche Brownienne
+cym      SAVE vdeptr
+      REAL,save,allocatable :: scavtr(:)      ! Coefficient de lessivage
+cym      SAVE scavtr
 cAA
       CHARACTER*2 itn
@@ -273 +282 @@
 C nature du traceur
 c
-      logical aerosol(nbtr)  ! Nature du traceur
+      logical,save,allocatable :: aerosol(:)  ! Nature du traceur
 c                            ! aerosol(it) = true  => aerosol 
 c                            ! aerosol(it) = false => gaz 
 c                            ! nat_trac(it) = 1. aerosol
-      logical clsol(nbtr)    ! clsol(it) = true => CL sol calculee
-      logical radio(nbtr)    ! radio(it)=true => decroisssance radioactive
-      save aerosol,clsol,radio
+      logical,save,allocatable :: clsol(:)    ! clsol(it) = true => CL sol calculee
+      logical,save,allocatable :: radio(:)    ! radio(it)=true => decroisssance radioactive
+cym      save aerosol,clsol,radio
 C
 c======================================================================
@@ -298 +307 @@
       REAL d_tr(klon,klev), d_trs(klon) ! tendances de traceurs 
       REAL d_tr_cl(klon,klev,nbtr) ! tendance de traceurs  couche limite
-      REAL d_tr_cv(klon,klev,nbtr) ! tendance de traceurs  conv pour chq traceur
+      REAL d_tr_cli(klon,klev,nbtr) ! tendance de traceurs  CL pour chq traceur
+      REAL d_tr_cv(klon,klev,nbtr) ! tendance de traceurs  convection
+      REAL d_tr_cvi(klon,klev,nbtr) ! tendance de traceurs  conv pour chq traceur
       REAL d_tr_th(klon,klev,nbtr) ! la tendance des thermiques
       REAL d_tr_dec(klon,klev,nbtr) ! la tendance de la decroissance 
@@ -367 +378 @@
 #endif
      .                                         /)
-
-
-      REAL, PARAMETER   :: dry_mass = 28.966
-      REAL, POINTER     :: hbuf(:)
-      REAL, ALLOCATABLE :: obuf(:)
-      REAL              :: calday
-      REAL              :: pdel(klon,klev)
-      REAL              :: dummy(klon,klev) = 0.
-#endif
-#ifdef INCA_AER
-      integer la
+      dummy(:,:)=0.
 #endif
 c
@@ -386 +387 @@
 
          if (debutphy) then
+           allocate( trs(klon,nbtr) )  
+           allocate( masktr(klon,nbtr)) 
+           allocate( fshtr(klon,nbtr) )
+           allocate( hsoltr(nbtr))
+           allocate( tautr(nbtr))
+           allocate( vdeptr(nbtr))
+           allocate( scavtr(nbtr))
+           allocate( aerosol(nbtr))
+           allocate( clsol(nbtr))
+           allocate( radio(nbtr))
 
           ecrit_tra = NINT(86400./pdtphys *ecritphy) 
@@ -416 +427 @@
 c        print*,'valeur de debut dans phytrac :',debutphy
          trs(:,:) = 0.
-
-         open (99,file='starttrac',status='old',
+         
+       if (phy_rank==0) then
+         trs_tmp(:)=0.
+         open (99,file='starttrac',status='old',
      .         err=999,form='formatted')
-         read(99,*) (trs(i,1),i=1,klon)
+         read(99,*) (trs_tmp(i),i=1,klon2)
 999      close(99)
+       endif
+       call ScatterField(trs_tmp,trs(:,1),1)
 c         print*, 'apres starttrac'
 
@@ -451 +466 @@
       endif
 #ifdef INCA
+      call VTe(VTphysiq)
+      call VTb(VTinca)
 !======================================================================
 !     Chimie
@@ -539 +556 @@
       END DO
 #endif
+      call VTe(VTinca)
+      call VTb(VTphysiq)
 #else
 
@@ -861 +880 @@
       if (lafin) then
          print*, 'c est la fin de la physique'
-         open (99,file='restarttrac',  form='formatted')
-         do i=1,klon
-             write(99,*) trs(i,1)
-         enddo
-         PRINT*, 'Ecriture du fichier restarttrac'
-         close(99)
+         
+         call GatherField(trs(:,1),trs_tmp,1)
+         if (phy_rank==0) then
+         
+           open (99,file='restarttrac',  form='formatted')
+           do i=1,klon
+               write(99,*) trs_tmp(i)
+           enddo
+           PRINT*, 'Ecriture du fichier restarttrac'
+           close(99)
+         endif
       else
          print*, 'physique pas fini'