source: LMDZ4/trunk/libf/dyn3dpar/create_etat0_limit.F @ 1076

Last change on this file since 1076 was 1017, checked in by lsce, 16 years ago
  • Copure des lignes trop longues

JG

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.3 KB
Line 
1!
2! $Header$
3!
4       PROGRAM create_etat0_limit
5       USE dimphy
6       USE comgeomphy
7       USE mod_phys_lmdz_para
8       USE mod_const_mpi
9       USE phys_state_var_mod     
10       IMPLICIT NONE
11c
12c
13c     Programme d'appel a etat0, creation des etats initiaux et limit_netcdf
14c   
15c
16c     interbar = .T . si appel a  interpol. barycentrique inter_barxy
17c
18c     extrap   = .T . si on fait une extrapolation de donnees , comme pour
19c       les  SST  lorsque  le fichier ne contient pas uniquement  des points
20c     oceaniques .
21c
22c     oldice   = .T. si l'on veut garder les anciennes glaces , obtenues
23c     par  grille_m  ( grid_atob ) .
24c
25c     on cree le masque dans etat0 que l'on passe ensuite dans limit pour
26c     garder les coherences
27
28      LOGICAL interbar, extrap , oldice
29      PARAMETER ( interbar = .true. , extrap = .FALSE. , oldice=.false.)
30#include "dimensions.h"
31#include "paramet.h"
32#include "indicesol.h"
33#include "advtrac.h"
34#include  "control.h"
35#include "clesphys.h"
36      REAL :: masque(iip1,jjp1)
37!      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)
38
39c initialisation traceurs
40      hadv_flg(:) = 0.
41      vadv_flg(:) = 0.
42      conv_flg(:) = 0.
43      pbl_flg(:)  = 0.
44      tracnam(:)  = '        '
45      nprath = 1
46      nbtrac = 0
47      mmt_adj(:,:,:,:) = 1
48
49      IF (config_inca /= 'none') THEN
50#ifdef INCA
51         call init_const_lmdz(
52     $        nbtrac,anneeref,dayref,
53     $        iphysiq,day_step,nday)
54#endif
55         print *, 'nbtrac =' , nbtrac
56      END IF
57
58      CALL init_mpi
59
60
61      CALL Init_Phys_lmdz(iim,jjp1,llm,nqmx-2,1,(/(jjm-1)*iim+2/))
62      PRINT *,'---> klon=',klon
63
64      IF (mpi_size>1 .OR. omp_size>1) THEN
65        CALL abort_gcm('create_etat0_limit','In parallel mode,
66     &                 create_etat0_limit must be called only
67     &                 for 1 process and 1 task')
68      ENDIF
69      CALL phys_state_var_init
70      call InitComgeomphy
71     
72
73      WRITE(6,*) '  *********************  '
74      WRITE(6,*) ' interbar = ',interbar
75      CALL etat0_netcdf ( interbar, masque )
76c
77      WRITE(6,1)
78      WRITE(6,*) '  *********************  '
79      WRITE(6,*) '  ***  Limit_netcdf ***  '
80      WRITE(6,*) '  *********************  '
81      WRITE(6,1)
82     
83c     
84      CALL  limit_netcdf ( interbar, extrap , oldice, masque)
85
861     FORMAT(//)
87
88      STOP
89      END
Note: See TracBrowser for help on using the repository browser.