source: LMDZ5/trunk/libf/dyn3d_common/iniconst.F90 @ 3411

Last change on this file since 3411 was 2601, checked in by Ehouarn Millour, 8 years ago

Cleanup in the dynamics: turn temps.h into module temps_mod.F90
EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.1 KB
RevLine 
[524]1!
[1403]2! $Id: iniconst.F90 2601 2016-07-24 09:51:55Z fairhead $
[524]3!
[1625]4SUBROUTINE iniconst
[524]5
[1625]6  USE control_mod
[1520]7#ifdef CPP_IOIPSL
[1625]8  use IOIPSL
[1520]9#else
[1625]10  ! if not using IOIPSL, we still need to use (a local version of) getin
11  use ioipsl_getincom
[1520]12#endif
[2597]13  USE comconst_mod, ONLY: im, imp1, jm, jmp1, lllm, lllmm1, lllmp1, &
14                          unsim, pi, r, kappa, cpp, dtvr, dtphys
[2600]15  USE comvert_mod, ONLY: disvert_type, pressure_exner
[2597]16 
[1625]17  IMPLICIT NONE
18  !
19  !      P. Le Van
20  !
21  !   Declarations:
22  !   -------------
23  !
24  include "dimensions.h"
25  include "paramet.h"
26  include "iniprint.h"
[524]27
[1625]28  character(len=*),parameter :: modname="iniconst"
29  character(len=80) :: abort_message
30  !
31  !
32  !
33  !-----------------------------------------------------------------------
34  !   dimension des boucles:
35  !   ----------------------
[524]36
[1625]37  im      = iim
38  jm      = jjm
39  lllm    = llm
40  imp1    = iim
41  jmp1    = jjm + 1
42  lllmm1  = llm - 1
43  lllmp1  = llm + 1
[524]44
[1625]45  !-----------------------------------------------------------------------
[524]46
[1625]47  dtphys  = iphysiq * dtvr
48  unsim   = 1./iim
49  pi      = 2.*ASIN( 1. )
[524]50
[1625]51  !-----------------------------------------------------------------------
52  !
[524]53
[1625]54  r       = cpp * kappa
[524]55
[1625]56  write(lunout,*) trim(modname),': R  CP  Kappa ',r,cpp,kappa
57  !
58  !-----------------------------------------------------------------------
[524]59
[1625]60  ! vertical discretization: default behavior depends on planet_type flag
61  if (planet_type=="earth") then
62     disvert_type=1
63  else
64     disvert_type=2
65  endif
66  ! but user can also specify using one or the other in run.def:
67  call getin('disvert_type',disvert_type)
68  write(lunout,*) trim(modname),': disvert_type=',disvert_type
[1520]69
[1625]70  pressure_exner = disvert_type == 1 ! default value
71  call getin('pressure_exner', pressure_exner)
72
73  if (disvert_type==1) then
74     ! standard case for Earth (automatic generation of levels)
[2040]75     call disvert()
[1625]76  else if (disvert_type==2) then
77     ! standard case for planets (levels generated using z2sig.def file)
78     call disvert_noterre
79  else
80     write(abort_message,*) "Wrong value for disvert_type: ", disvert_type
81     call abort_gcm(modname,abort_message,0)
82  endif
83
84END SUBROUTINE iniconst
Note: See TracBrowser for help on using the repository browser.