source: LMDZ4/branches/LMDZ4_V3_patches/libf/dyn3d/create_etat0_limit.F @ 5023

Last change on this file since 5023 was 845, checked in by Laurent Fairhead, 17 years ago

Création de la branche LMDZ4_V3_patches: à partir de la version de référence LMDZ4_V3,

on corrige différents bugs qui sont sur la branche de développement HEAD pour établir
une version de référence LMDZ4_V3 pour les utilisateurs non-développeurs

Cette branche ne servira que pour les corrections de bugs et les version y seront tagguées

LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.8 KB
Line 
1!
2! $Header$
3!
4       PROGRAM create_etat0_limit
5       USE dimphy
6       USE comgeomphy
7
8c
9c
10c     Programme d'appel a etat0, creation des etats initiaux et limit_netcdf
11c   
12c
13c     interbar = .T . si appel a  interpol. barycentrique inter_barxy
14c
15c     extrap   = .T . si on fait une extrapolation de donnees , comme pour
16c       les  SST  lorsque  le fichier ne contient pas uniquement  des points
17c     oceaniques .
18c
19c     oldice   = .T. si l'on veut garder les anciennes glaces , obtenues
20c     par  grille_m  ( grid_atob ) .
21c
22c     on cree le masque dans etat0 que l'on passe ensuite dans limit pour
23c     garder les coherences
24
25      LOGICAL interbar, extrap , oldice
26      PARAMETER ( interbar = .true. , extrap = .FALSE. , oldice=.false.)
27#include "dimensions.h"
28#include "paramet.h"
29#include "indicesol.h"
30#include "advtrac.h"
31cym#include "dimphy.h"
32      REAL :: masque(iip1,jjp1)
33      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)
34
35c initialisation traceurs
36      hadv_flg(:) = 0.
37      vadv_flg(:) = 0.
38      conv_flg(:) = 0.
39      pbl_flg(:)  = 0.
40      tracnam(:)  = '        '
41      nprath = 1
42      nbtrac = 0
43      mmt_adj(:,:,:,:) = 1
44
45#ifdef INCA
46      call init_const_lmdz(nbtrac,anneeref,dayref,iphysiq,day_step,nday)
47      print *, 'nbtrac =' , nbtrac
48#endif
49
50      CALL Init_Phys_lmdz(iim,jjp1,llm,nqmx-2,1,(jjm-1)*iim+2)
51      call InitComgeomphy
52
53      WRITE(6,*) '  *********************  '
54      WRITE(6,*) ' interbar = ',interbar
55      CALL etat0_netcdf ( interbar, masque, pctsrf )
56c
57      WRITE(6,1)
58      WRITE(6,*) '  *********************  '
59      WRITE(6,*) '  ***  Limit_netcdf ***  '
60      WRITE(6,*) '  *********************  '
61      WRITE(6,1)
62     
63c     
64      CALL  limit_netcdf ( interbar, extrap , oldice, masque, pctsrf )
65
661     FORMAT(//)
67
68      STOP
69      END
Note: See TracBrowser for help on using the repository browser.