source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/sucape.F90 @ 5308

Last change on this file since 5308 was 3331, checked in by acozic, 6 years ago

Add modification for isotopes

  • Property svn:executable set to *
File size: 2.5 KB
Line 
1SUBROUTINE SUCAPE(KULOUT)
2
3!**** *SUCAPE  * - ROUTINE TO INITIALIZE THE VARIABLES FOR CAPE
4!                  AND CIN COMPUTATION
5
6!     PURPOSE.
7!     --------
8!        SET DEFAULT VALUES, THEN READS NAMELIST NAMCAPE
9
10!**   INTERFACE.
11!     ----------
12!        *CALL* *SUCAPE(KULOUT)*
13
14!         EXPLICIT ARGUMENTS :  KULOUT
15!         --------------------
16
17!         IMPLICIT ARGUMENTS :
18!         --------------------
19!            COMMON  YOMCAPE
20!            COMMON  YOMPHY
21!            COMMON  YOMLUN
22
23!     METHOD.
24!     -------
25
26!     EXTERNALS.
27!     ----------
28
29!     REFERENCE.
30!     ----------
31
32!     AUTHOR.
33!     -------
34!        N.Pristov 03/2001
35
36!     MODIFICATIONS.
37!     --------------
38!        ORIGINAL : 03/2001
39!        M.Hamrud      01-Oct-2003 CY28 Cleaning
40!-----------------------------------------------------------------------
41
42USE PARKIND1  ,ONLY : JPIM     ,JPRB
43USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
44
45USE YOMCAPE  , ONLY :  NCAPEITER   ,NETAPES  ,GCAPERET   ,GCAPEPSD
46USE YOMPHY   , ONLY :  NBITER
47! Ce qui concerne NULNAM commente par MPL le 15.04.09
48!USE YOMLUN   , ONLY :  NULNAM
49
50IMPLICIT NONE
51
52INTEGER(KIND=JPIM),INTENT(IN)    :: KULOUT
53REAL(KIND=JPRB) :: ZHOOK_HANDLE
54
55#include "abor1.intfb.h"
56#include "posnam.intfb.h"
57
58#include "namcape.h"
59
60!      ----------------------------------------------------------------
61!*       1.    SET DEFAULT VALUES.
62!              -------------------
63IF (LHOOK) CALL DR_HOOK('SUCAPE',0,ZHOOK_HANDLE)
64NCAPEITER=NBITER
65NETAPES=2
66GCAPERET=0._JPRB
67GCAPEPSD=30000._JPRB
68
69!      ----------------------------------------------------------------
70!*       2.    MODIFIES DEFAULT VALUES. READ NAMELIST.
71!              ------------------------
72
73! Ce qui concerne NAMCAPE commente par MPL le 15.04.09
74!CALL POSNAM(NULNAM,'NAMCAPE')
75!READ(NULNAM,NAMCAPE)
76
77IF (NCAPEITER <= 0 ) CALL ABOR1('SUCAPE:  INVALID VALUE FOR NCAPEITER')
78IF (NETAPES <= 0 ) CALL ABOR1('SUCAPE:  INVALID VALUE FOR NETAPES')
79IF ((GCAPERET < 0.0_JPRB ) .OR. (GCAPERET > 1.0_JPRB ))  &
80 & CALL ABOR1('SUCAPE:  INVALID VALUE FOR GCAPERET') 
81IF (GCAPEPSD <= 0.0_JPRB ) CALL ABOR1('SUCAPE:  INVALID VALUE FOR GCAPEPSD')
82
83!      -----------------------------------------------------------
84!*       3.    PRINT FINAL VALUES.
85!              -------------------
86
87WRITE(KULOUT,'('' NCAPEITER = '',I2,'' NETAPES = '',I2,'' GCAPERET = '',&
88 & E13.6,'' GCAPEPSD = '',E13.6)') NCAPEITER, NETAPES, GCAPERET, GCAPEPSD 
89 
90IF (LHOOK) CALL DR_HOOK('SUCAPE',1,ZHOOK_HANDLE)
91END SUBROUTINE SUCAPE
Note: See TracBrowser for help on using the repository browser.