source: LMDZ4/trunk/libf/dyn3dpar/ce0l.F90 @ 1403

Last change on this file since 1403 was 1403, checked in by Laurent Fairhead, 14 years ago

Merged LMDZ4V5.0-dev branch changes r1292:r1399 to trunk.

Validation:
Validation consisted in compiling the HEAD revision of the trunk,
LMDZ4V5.0-dev branch and the merged sources and running different
configurations on local and SX8 machines comparing results.

Local machine: bench configuration, 32x24x11, gfortran

  • IPSLCM5A configuration (comparison between trunk and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent
  • MH07 configuration, new physics package (comparison between LMDZ4V5.0-dev branch and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent

SX8 machine (brodie), 96x95x39 on 4 processors:

  • IPSLCM5A configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent
  • MH07 configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent

Changes to the makegcm and create_make_gcm scripts to take into account
main programs in F90 files

Fusion de la branche LMDZ4V5.0-dev (r1292:r1399) au tronc principal

Validation:
La validation a consisté à compiler la HEAD de le trunk et de la banche
LMDZ4V5.0-dev et les sources fusionnées et de faire tourner le modéle selon
différentes configurations en local et sur SX8 et de comparer les résultats

En local: 32x24x11, config bench/gfortran

  • pour une config IPSLCM5A (comparaison tronc/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux (à part sur RN et Pb)
    • fichiers histoire égaux
  • pour une config nlle physique (MH07) (comparaison LMDZ4v5.0-dev/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux
    • fichiers histoire équivalents

Sur brodie, 96x95x39 sur 4 proc:

  • pour une config IPSLCM5A:
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc
  • pour une config MH07
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc

Changement sur makegcm et create_make-gcm pour pouvoir prendre en compte des
programmes principaux en *F90
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.6 KB
Line 
1!
2! $Id: ce0l.F90 1403 2010-07-01 09:02:53Z fairhead $
3!
4!-------------------------------------------------------------------------------
5!
6PROGRAM ce0l
7!
8!-------------------------------------------------------------------------------
9! Purpose: Calls etat0, creates initial states and limit_netcdf
10!
11!     interbar=.T. for barycentric interpolation inter_barxy
12!     extrap  =.T. for data extrapolation, like for the SSTs when file does not
13!                  contain ocean points only.
14!     oldice  =.T. for old-style ice, obtained using grille_m (grid_atob).
15!     masque is created in etat0, passed to limit to ensure consistancy.
16!-------------------------------------------------------------------------------
17  USE control_mod
18#ifdef CPP_EARTH
19! This prog. is designed to work for Earth
20  USE dimphy
21  USE comgeomphy
22  USE mod_phys_lmdz_para
23  USE mod_const_mpi
24  USE infotrac
25
26#ifdef CPP_IOIPSL
27  USE ioipsl, ONLY: ioconf_calendar
28#endif
29
30#endif
31  IMPLICIT NONE
32#ifndef CPP_EARTH
33  WRITE(lunout,*)'limit_netcdf: Earth-specific routine, needs Earth physics'
34#else
35!-------------------------------------------------------------------------------
36! Local variables:
37  LOGICAL, PARAMETER :: interbar=.TRUE., extrap=.FALSE., oldice=.FALSE.
38#include "dimensions.h"
39#include "paramet.h"
40#include "indicesol.h"
41#include "iniprint.h"
42#include "temps.h"
43#include "logic.h"
44  INTEGER, PARAMETER            :: longcles=20
45  REAL,    DIMENSION(longcles)  :: clesphy0
46  REAL,    DIMENSION(iip1,jjp1) :: masque
47  CHARACTER(LEN=15)             :: calnd
48!-------------------------------------------------------------------------------
49  CALL conf_gcm( 99, .TRUE. , clesphy0 )
50
51  CALL init_mpi
52
53  CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
54  WRITE(lunout,*)'---> klon=',klon
55  IF (mpi_size>1 .OR. omp_size>1) THEN
56       CALL abort_gcm('ce0l','In parallel mode,                         &
57 &                 ce0l must be called only                             &
58 &                 for 1 process and 1 task')
59  ENDIF
60
61  CALL InitComgeomphy
62
63#ifdef CPP_IOIPSL
64  SELECT CASE(calend)
65    CASE('earth_360d');CALL ioconf_calendar('360d');      calnd='a 360 jours/an'
66    CASE('earth_365d');CALL ioconf_calendar('noleap');    calnd='a 365 jours/an'
67    CASE('earth_366d');CALL ioconf_calendar('366d');      calnd='bissextile'
68    CASE('gregorian'); CALL ioconf_calendar('gregorian'); calnd='gregorien'
69    CASE('standard');  CALL ioconf_calendar('gregorian'); calnd='gregorien'
70    CASE('julian');    CALL ioconf_calendar('julian');    calnd='julien'
71    CASE('proleptic_gregorian'); CALL ioconf_calendar('gregorian')
72  !--- DC Bof...  => IOIPSL a mettre a jour: proleptic_gregorian /= gregorian
73    CASE DEFAULT
74      CALL abort_gcm('ce0l','Mauvais choix de calendrier',1)
75  END SELECT
76  WRITE(lunout,*)'CALENDRIER CHOISI: Terrestre '//TRIM(calnd)
77#endif
78
79  IF (config_inca /= 'none') THEN
80#ifdef INCA
81    CALL init_const_lmdz(nbtr,anneeref,dayref,iphysiq,day_step,nday)
82    CALL init_inca_para(iim,jjm+1,klon,1,klon_mpi_para_nb,0)
83    WRITE(lunout,*)'nbtr =' , nbtr
84#endif
85  END IF
86
87  WRITE(lunout,'(//)')
88  WRITE(lunout,*) '  *********************  '
89  WRITE(lunout,*) '  ***  etat0_netcdf ***  '
90  WRITE(lunout,*) '  *********************  '
91  WRITE(lunout,'(//)')
92  WRITE(lunout,*) ' interbar = ',interbar
93  CALL etat0_netcdf(interbar,masque,ok_etat0)
94
95  IF(ok_limit) THEN
96  WRITE(lunout,'(//)')
97  WRITE(lunout,*) '  *********************  '
98  WRITE(lunout,*) '  ***  Limit_netcdf ***  '
99  WRITE(lunout,*) '  *********************  '
100  WRITE(lunout,'(//)')
101  CALL limit_netcdf(interbar,extrap,oldice,masque)
102  END IF
103
104#endif
105! of #ifndef CPP_EARTH #else
106  STOP
107
108END PROGRAM ce0l
109!
110!-------------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.