source: LMDZ6/branches/DYNAMICO-conv/libf/phylmd/inifis_mod.F90 @ 5445

Last change on this file since 5445 was 3049, checked in by Laurent Fairhead, 7 years ago

Replacing iniprint

  • 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.8 KB
RevLine 
[1403]1! $Id: inifis_mod.F90 3049 2017-10-29 10:13:59Z fhourdin $
[2311]2MODULE inifis_mod
[524]3
[2311]4CONTAINS
5
6  SUBROUTINE inifis(punjours, prad, pg, pr, pcpp)
7  ! Initialize some physical constants and settings
[3049]8  USE init_print_control_mod, ONLY : init_print_control
9  USE print_control_mod, ONLY: lunout
[1992]10  IMPLICIT NONE
[1403]11
[2311]12  include "YOMCST.h"
13  REAL,INTENT(IN) :: prad, pg, pr, pcpp, punjours
[524]14
[1992]15  CHARACTER (LEN=20) :: modname = 'inifis'
16  CHARACTER (LEN=80) :: abort_message
17
[2311]18  ! Initialize flags lunout, prt_level, debug
19  CALL init_print_control
[1992]20
[2311]21  ! suphel => initialize some physical constants (orbital parameters,
22  !           geoid, gravity, thermodynamical constants, etc.) in the
23  !           physics
24  CALL suphel
[1992]25
[2311]26  ! check that physical constants set in 'suphel' are coherent
27  ! with values set in the dynamics:
28  IF (rday/=punjours) THEN
29    WRITE (lunout, *) 'inifis: length of day discrepancy!!!'
30    WRITE (lunout, *) '  in the dynamics punjours=', punjours
31    WRITE (lunout, *) '   but in the physics RDAY=', rday
32    IF (abs(rday-punjours)>0.01*punjours) THEN
33        ! stop here if the relative difference is more than 1%
34      abort_message = 'length of day discrepancy'
35      CALL abort_physic(modname, abort_message, 1)
36    END IF
[1992]37  END IF
[2311]38  IF (rg/=pg) THEN
39    WRITE (lunout, *) 'inifis: gravity discrepancy !!!'
40    WRITE (lunout, *) '     in the dynamics pg=', pg
41    WRITE (lunout, *) '  but in the physics RG=', rg
42    IF (abs(rg-pg)>0.01*pg) THEN
43        ! stop here if the relative difference is more than 1%
44      abort_message = 'gravity discrepancy'
45      CALL abort_physic(modname, abort_message, 1)
46    END IF
47  END IF
48  IF (ra/=prad) THEN
49    WRITE (lunout, *) 'inifis: planet radius discrepancy !!!'
50    WRITE (lunout, *) '   in the dynamics prad=', prad
51    WRITE (lunout, *) '  but in the physics RA=', ra
52    IF (abs(ra-prad)>0.01*prad) THEN
53        ! stop here if the relative difference is more than 1%
54      abort_message = 'planet radius discrepancy'
55      CALL abort_physic(modname, abort_message, 1)
56    END IF
57  END IF
58  IF (rd/=pr) THEN
59    WRITE (lunout, *) 'inifis: reduced gas constant discrepancy !!!'
60    WRITE (lunout, *) '     in the dynamics pr=', pr
61    WRITE (lunout, *) '  but in the physics RD=', rd
62    IF (abs(rd-pr)>0.01*pr) THEN
63        ! stop here if the relative difference is more than 1%
64      abort_message = 'reduced gas constant discrepancy'
65      CALL abort_physic(modname, abort_message, 1)
66    END IF
67  END IF
68  IF (rcpd/=pcpp) THEN
69    WRITE (lunout, *) 'inifis: specific heat discrepancy !!!'
70    WRITE (lunout, *) '     in the dynamics pcpp=', pcpp
71    WRITE (lunout, *) '  but in the physics RCPD=', rcpd
72    IF (abs(rcpd-pcpp)>0.01*pcpp) THEN
73        ! stop here if the relative difference is more than 1%
74      abort_message = 'specific heat discrepancy'
75      CALL abort_physic(modname, abort_message, 1)
76    END IF
77  END IF
[1992]78
[2311]79  END SUBROUTINE inifis
80 
81END MODULE inifis_mod
Note: See TracBrowser for help on using the repository browser.