source: LMDZ4/trunk/libf/dyn3d/create_etat0_limit.F @ 1279

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

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.9 KB
Line 
1!
2! $Id: create_etat0_limit.F 1279 2009-12-10 09:02:56Z fairhead $
3!
4       PROGRAM create_etat0_limit
5#ifdef CPP_EARTH
6! This prog. is designed to work for Earth
7       USE dimphy
8       USE comgeomphy
9       USE infotrac
10#ifdef CPP_IOIPSL
11       use ioipsl, only: ioconf_calendar
12#endif
13       IMPLICIT NONE
14c
15c
16c     Programme d'appel a etat0, creation des etats initiaux et limit_netcdf
17c   
18c
19c     interbar = .T . si appel a  interpol. barycentrique inter_barxy
20c
21c     extrap   = .T . si on fait une extrapolation de donnees , comme pour
22c       les  SST  lorsque  le fichier ne contient pas uniquement  des points
23c     oceaniques .
24c
25c     oldice   = .T. si l'on veut garder les anciennes glaces , obtenues
26c     par  grille_m  ( grid_atob ) .
27c
28c     on cree le masque dans etat0 que l'on passe ensuite dans limit pour
29c     garder les coherences
30
31      LOGICAL interbar, extrap , oldice
32      PARAMETER ( interbar = .true. , extrap = .FALSE. , oldice=.false.)
33#include "dimensions.h"
34#include "paramet.h"
35#include "indicesol.h"
36#include  "control.h"
37      REAL :: masque(iip1,jjp1)
38!      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)
39
40      IF (config_inca /= 'none') THEN
41#ifdef INCA
42         call init_const_lmdz(
43     $        nbtr,anneeref,dayref,
44     $        iphysiq, day_step,nday)
45#endif
46         print *, 'nbtr =' , nbtr
47      END IF
48
49      CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
50      PRINT *,'---> klon=',klon
51      call InitComgeomphy
52
53#ifdef CPP_IOIPSL
54      call ioconf_calendar('360d')
55#endif
56
57      WRITE(6,*) '  *********************  '
58      WRITE(6,*) ' interbar = ',interbar
59      CALL etat0_netcdf ( interbar, masque )
60c
61      WRITE(6,1)
62      WRITE(6,*) '  *********************  '
63      WRITE(6,*) '  ***  Limit_netcdf ***  '
64      WRITE(6,*) '  *********************  '
65      WRITE(6,1)
66     
67c     
68      CALL  limit_netcdf ( interbar, extrap , oldice, masque)
69
701     FORMAT(//)
71
72#endif
73! of #ifdef CPP_EARTH
74      STOP
75      END
Note: See TracBrowser for help on using the repository browser.