source: LMDZ4/trunk/libf/dyn3d/ce0l.F90 @ 1378

Last change on this file since 1378 was 1323, checked in by Laurent Fairhead, 15 years ago

Changes made in r1293 are integrated into the trunk
Start files are identical between r1293 and this version


Les modifications de la r1293 sont intégrées à la trunk
Les fichiers start et startphy sont identiques entre la version 1293 et celle-ci

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