Changeset 270 for trunk/LMDZ.COMMON/libf/dyn3dpar/ce0l.F90

Timestamp:

Aug 18, 2011, 12:09:27 PM (13 years ago)

Author:

emillour

Message:

Ehouarn: Mise a jour des dynamiques (seq et ) pour suivre
les changements et developpements de LMDZ5 terrestre
(mise a niveau avec LMDZ5 trunk, rev 1560). Ce qui ne devrais pas changer grand chose au fonctionnement de base du GCM).
Modification importante: correction du bug sur le cumul des flux de masse pour le transport des traceurs (mais dans les faits semble avoir peu d'impact).
(voir "commit_importants.log" pour les details des ajouts et modifications).

File:

• trunk/LMDZ.COMMON/libf/dyn3dpar/ce0l.F90 (modified) (6 diffs)

trunk/LMDZ.COMMON/libf/dyn3dpar/ce0l.F90

@@ -1 +1 @@
 !
-! $Id: ce0l.F90 1482 2011-02-09 15:03:09Z jghattas $
+! $Id: ce0l.F90 1511 2011-04-28 15:21:47Z jghattas $
 !
 !-------------------------------------------------------------------------------
@@ -19 +19 @@
   USE dimphy
   USE comgeomphy
-  USE mod_phys_lmdz_para
-  USE mod_const_mpi
   USE infotrac
-  USE parallel, ONLY: finalize_parallel
 
 #ifdef CPP_IOIPSL
@@ -30 +27 @@
 #endif
   IMPLICIT NONE
-#include "iniprint.h"
 #ifndef CPP_EARTH
   WRITE(lunout,*)'limit_netcdf: Earth-specific routine, needs Earth physics'
@@ -40 +36 @@
 #include "paramet.h"
 #include "indicesol.h"
+#include "iniprint.h"
 #include "temps.h"
 #include "logic.h"
+  INTEGER, PARAMETER            :: longcles=20
+  REAL,    DIMENSION(longcles)  :: clesphy0
   REAL,    DIMENSION(iip1,jjp1) :: masque
   CHARACTER(LEN=15)             :: calnd
+  REAL,    DIMENSION(iip1,jjp1) :: phis ! geopotentiel au sol
 !-------------------------------------------------------------------------------
-  CALL conf_gcm( 99, .TRUE. )
-
-  CALL init_mpi
+  CALL conf_gcm( 99, .TRUE. , clesphy0 )
 
   CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
   WRITE(lunout,*)'---> klon=',klon
-  IF (mpi_size>1 .OR. omp_size>1) THEN
-       CALL abort_gcm('ce0l','In parallel mode,                         &
- &                 ce0l must be called only                             &
- &                 for 1 process and 1 task',1)
-  ENDIF
-
   CALL InitComgeomphy
 
@@ -89 +81 @@
   WRITE(lunout,'(//)')
   WRITE(lunout,*) ' interbar = ',interbar
-  CALL etat0_netcdf(interbar,masque,ok_etat0)
+  CALL etat0_netcdf(interbar,masque,phis,ok_etat0)
 
   IF(ok_limit) THEN
@@ -100 +92 @@
   END IF
 
-!$OMP MASTER
-  CALL finalize_parallel
-!$OMP END MASTER
-
+  IF (grilles_gcm_netcdf) THEN
+     WRITE(lunout,'(//)')
+     WRITE(lunout,*) '  ***************************  '
+     WRITE(lunout,*) '  ***  grilles_gcm_netcdf ***  '
+     WRITE(lunout,*) '  ***************************  '
+     WRITE(lunout,'(//)')
+     CALL grilles_gcm_netcdf_sub(masque,phis)
+  END IF
 #endif
 ! of #ifndef CPP_EARTH #else