Changeset 985 for LMDZ4/trunk/libf/dyn3dpar/conf_gcm.F

Timestamp:

Jul 30, 2008, 5:50:03 PM (16 years ago)

Author:

Laurent Fairhead

Message:

Mise a jour de dyn3dpar par rapport a dyn3d, inclusion OpenMP et filtre FFT YM
LF

File:

• LMDZ4/trunk/libf/dyn3dpar/conf_gcm.F (modified) (8 diffs)

LMDZ4/trunk/libf/dyn3dpar/conf_gcm.F

@@ -8 +8 @@
       use IOIPSL
       use misc_mod
+      use mod_filtre_fft, ONLY : use_filtre_fft
+      use mod_hallo, ONLY : use_mpi_alloc
+      use parallel, ONLY : omp_chunk
       IMPLICIT NONE
 c-----------------------------------------------------------------------
@@ -34 +37 @@
 #include "serre.h"
 #include "comdissnew.h"
-#include "clesphys.h"
+!#include "clesphys.h"
 #include "iniprint.h"
+
+! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
 c
 c
@@ -298 +303 @@
 c
 
-!Config  Key  = ecritphy
-!Config  Desc = Frequence d'ecriture de la physique
-!Config  Def  = 1
-!Config  Help = frequence  de l'ecriture du fichier histphy
-!Config         en jours.
-       ecritphy = 1
-       CALL getin('ecritphy',ecritphy)
-
-!Config  Key  = cycle_diurne
-!Config  Desc = Cycle ddiurne
-!Config  Def  = y
-!Config  Help = Cette option permet d'eteidre le cycle diurne.
-!Config         Peut etre util pour accelerer le code !
-       cycle_diurne = .TRUE.
-       CALL getin('cycle_diurne',cycle_diurne)
-
-!Config  Key  = soil_model
-!Config  Desc = Modele de sol
-!Config  Def  = y
-!Config  Help = Choix du modele de sol (Thermique ?)
-!Config         Option qui pourait un string afin de pouvoir
-!Config         plus de choix ! Ou meme une liste d'options !
-       soil_model = .TRUE.
-       CALL getin('soil_model',soil_model)
-
-!Config  Key  = new_oliq
-!Config  Desc = Nouvelle eau liquide
-!Config  Def  = y
-!Config  Help = Permet de mettre en route la
-!Config         nouvelle parametrisation de l'eau liquide !
-       new_oliq = .TRUE.
-       CALL getin('new_oliq',new_oliq)
-
-!Config  Key  = ok_orodr
-!Config  Desc = Orodr ???
-!Config  Def  = y
-!Config  Help = Y en a pas comprendre !
-!Config         
-       ok_orodr = .TRUE.
-       CALL getin('ok_orodr',ok_orodr)
-
-!Config  Key  =  ok_orolf
-!Config  Desc = Orolf ??
-!Config  Def  = y
-!Config  Help = Connais pas !
-       ok_orolf = .TRUE.
-       CALL getin('ok_orolf', ok_orolf)
-
-!Config  Key  = ok_limitvrai
-!Config  Desc = Force la lecture de la bonne annee
-!Config  Def  = n
-!Config  Help = On peut forcer le modele a lire le
-!Config         fichier SST de la bonne annee. C'est une tres bonne
-!Config         idee, pourquoi ne pas mettre toujours a y ???
-       ok_limitvrai = .FALSE.
-       CALL getin('ok_limitvrai',ok_limitvrai)
-
-!Config  Key  = nbapp_rad
-!Config  Desc = Frequence d'appel au rayonnement
-!Config  Def  = 12
-!Config  Help = Nombre  d'appels des routines de rayonnements
-!Config         par jour.
-       nbapp_rad = 12
-       CALL getin('nbapp_rad',nbapp_rad)
-
-!Config  Key  = iflag_con
-!Config  Desc = Flag de convection
-!Config  Def  = 2
-!Config  Help = Flag  pour la convection les options suivantes existent :
-!Config         1 pour LMD,
-!Config         2 pour Tiedtke,
-!Config         3 pour CCM(NCAR)  
-       iflag_con = 2
-       CALL getin('iflag_con',iflag_con)
-!
+
 !Config  Key  = ip_ebil_dyn
 !Config  Desc = PRINT level for energy conserv. diag.
@@ -384 +315 @@
        CALL getin('ip_ebil_dyn',ip_ebil_dyn)
 !
-
-      DO i = 1, longcles
-       clesphy0(i) = 0.
-      ENDDO
-                          clesphy0(1) = FLOAT( iflag_con )
-                          clesphy0(2) = FLOAT( nbapp_rad )
-
-       IF( cycle_diurne  ) clesphy0(3) =  1.
-       IF(   soil_model  ) clesphy0(4) =  1.
-       IF(     new_oliq  ) clesphy0(5) =  1.
-       IF(     ok_orodr  ) clesphy0(6) =  1.
-       IF(     ok_orolf  ) clesphy0(7) =  1.
-       IF(  ok_limitvrai ) clesphy0(8) =  1.
 
 
@@ -633 +551 @@
       write(lunout,*)' purmats = ', purmats
       write(lunout,*)' iflag_phys = ', iflag_phys
-      write(lunout,*)' iphysiq = ', iphysiq
-      write(lunout,*)' ecritphy = ', ecritphy
-      write(lunout,*)' cycle_diurne = ', cycle_diurne
-      write(lunout,*)' soil_model = ', soil_model
-      write(lunout,*)' new_oliq = ', new_oliq
-      write(lunout,*)' ok_orodr = ', ok_orodr
-      write(lunout,*)' ok_orolf = ', ok_orolf
-      write(lunout,*)' ok_limitvrai = ', ok_limitvrai
-      write(lunout,*)' nbapp_rad = ', nbapp_rad
-      write(lunout,*)' iflag_con = ', iflag_con
       write(lunout,*)' clonn = ', clonn 
       write(lunout,*)' clatt = ', clatt
@@ -777 +685 @@
       CALL getin('config_inca',config_inca)
 
+!Config  Key  = use_filtre_fft
+!Config  Desc = flag d'activation des FFT pour le filtre
+!Config  Def  = false
+!Config  Help = permet d'activer l'utilisation des FFT pour effectuer
+!Config         le filtrage aux poles. 
+      use_filtre_fft=.FALSE.
+      CALL getin('use_filtre_fft',use_filtre_fft)
+      
+!Config  Key  = use_mpi_alloc
+!Config  Desc = Utilise un buffer MPI en mémoire globale
+!Config  Def  = false
+!Config  Help = permet d'activer l'utilisation d'un buffer MPI
+!Config         en mémoire globale a l'aide de la fonction MPI_ALLOC.
+!Config         Cela peut améliorer la bande passante des transferts MPI
+!Config         d'un facteur 2  
+      use_mpi_alloc=.FALSE.
+      CALL getin('use_mpi_alloc',use_mpi_alloc)
+
+!Config  Key  = omp_chunk
+!Config  Desc = taille des blocs openmp
+!Config  Def  = 1
+!Config  Help = defini la taille des packets d'itération openmp
+!Config         distribuée à chaque tâche lors de l'entrée dans une
+!Config         boucle parallélisée
+  
+      omp_chunk=1
+      CALL getin('omp_chunk',omp_chunk)
 
       write(lunout,*)' #########################################'
@@ -799 +734 @@
       write(lunout,*)' purmats = ', purmats
       write(lunout,*)' iflag_phys = ', iflag_phys
-      write(lunout,*)' iphysiq = ', iphysiq
-      write(lunout,*)' ecritphy = ', ecritphy
-      write(lunout,*)' cycle_diurne = ', cycle_diurne
-      write(lunout,*)' soil_model = ', soil_model
-      write(lunout,*)' new_oliq = ', new_oliq
-      write(lunout,*)' ok_orodr = ', ok_orodr
-      write(lunout,*)' ok_orolf = ', ok_orolf
-      write(lunout,*)' ok_limitvrai = ', ok_limitvrai
-      write(lunout,*)' nbapp_rad = ', nbapp_rad
-      write(lunout,*)' iflag_con = ', iflag_con
-      write(lunout,*)' clonn = ', clonn 
-      write(lunout,*)' clatt = ', clatt
+      write(lunout,*)' clon = ', clon 
+      write(lunout,*)' clat = ', clat
       write(lunout,*)' grossismx = ', grossismx
       write(lunout,*)' grossismy = ', grossismy
@@ -820 +745 @@
       write(lunout,*)' offline = ', offline
       write(lunout,*)' config_inca = ', config_inca
+      write(lunout,*)' use_filtre_fft = ', use_filtre_fft
+      write(lunout,*)' use_mpi_alloc = ', use_mpi_alloc
+      write(lunout,*)' omp_chunk = ', omp_chunk
 c
       RETURN