[1989] | 1 | SUBROUTINE SUINIT(klon,klev) |
---|
| 2 | #ifdef DOC |
---|
| 3 | |
---|
| 4 | ! **** *SUINIT* - SCM initialization. |
---|
| 5 | |
---|
| 6 | ! Purpose. |
---|
| 7 | ! -------- |
---|
| 8 | |
---|
| 9 | ! ** Interface. |
---|
| 10 | ! ---------- |
---|
| 11 | |
---|
| 12 | ! Explicit arguments : None. |
---|
| 13 | ! -------------------- |
---|
| 14 | |
---|
| 15 | ! Implicit arguments : None. |
---|
| 16 | ! -------------------- |
---|
| 17 | |
---|
| 18 | ! Method. |
---|
| 19 | ! ------- |
---|
| 20 | |
---|
| 21 | ! Externals. None. |
---|
| 22 | ! ---------- |
---|
| 23 | |
---|
| 24 | ! Reference. |
---|
| 25 | ! ---------- |
---|
| 26 | |
---|
| 27 | ! Author. |
---|
| 28 | ! ------- |
---|
| 29 | ! Eric Bazile, Francois Bouyssel et Jean-Marcel Piriou |
---|
| 30 | |
---|
| 31 | ! Modifications. |
---|
| 32 | ! -------------- |
---|
| 33 | ! Original :97-02-01 |
---|
| 34 | ! Jozef Vivoda, SHMI: calling sequence as in 3D model |
---|
| 35 | ! and ECMWF setup |
---|
| 36 | ! 2001-11-27 P. Marquet : several printout on listing (NULOUT=15) |
---|
| 37 | |
---|
| 38 | ! ------------------------------------------------------------------ |
---|
| 39 | #endif |
---|
| 40 | |
---|
| 41 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
---|
| 42 | !#include "tsmbkind.h" |
---|
| 43 | |
---|
| 44 | USE PARDIM, ONLY : JPMXLE |
---|
| 45 | USE YOMCT0B , ONLY : LECMWF |
---|
| 46 | USE YOMRIP , ONLY : NINDAT ,NSSSSS |
---|
| 47 | USE YOMDIM |
---|
| 48 | USE YOMDPHY |
---|
| 49 | ! MPL 29042010: NDLNPR,RHYDR0 non initialises et pour ne pas mettre tout sudyn.F90 |
---|
| 50 | USE YOMDYN , ONLY : TSTEP , NDLNPR , RHYDR0 ! MPL 29042010 |
---|
| 51 | !USE YOMEVOL , ONLY : TECH ,FREQFS ,FREQFE , FREQDDH |
---|
| 52 | !USE YOMCT0 , ONLY : LFROG |
---|
| 53 | ! quelques ajouts qui viennent de suallo |
---|
| 54 | USE YOMGEM , ONLY : VDELA , VDELB ,VC ,NLOEN ,NLOENG ,NGPTOT |
---|
| 55 | USE YOMSTA , ONLY : STZ ,STPREH ,STPRE ,STPHI ,STTEM ,STDEN |
---|
| 56 | USE YOEAERD , ONLY : CVDAES ,CVDAEL ,CVDAEU ,CVDAED |
---|
| 57 | USE YOEOVLP , ONLY : RA1OVLP |
---|
| 58 | USE YOECLD , ONLY : CETA |
---|
| 59 | USE YOECND , ONLY : CEVAPCU |
---|
| 60 | USE YOMTOPH , ONLY : RMESOU ,RMESOT ,RMESOQ |
---|
| 61 | USE YOMGC , ONLY : GEMU ,GELAM ,GELAT ,GECLO ,GESLO ,GM ,GAW |
---|
| 62 | |
---|
| 63 | |
---|
| 64 | IMPLICIT NONE |
---|
| 65 | LOGICAL LLTRACE, LLDEBUG |
---|
| 66 | integer klon,klev |
---|
| 67 | CHARACTER*200 CFICP |
---|
| 68 | CHARACTER*200 CFLUX |
---|
| 69 | CHARACTER*200 CLIST |
---|
| 70 | CHARACTER*200 CFDDH |
---|
| 71 | CHARACTER*80 CNMEXP |
---|
| 72 | |
---|
| 73 | |
---|
| 74 | LLTRACE=.TRUE. |
---|
| 75 | LLDEBUG=.TRUE. |
---|
| 76 | |
---|
| 77 | ! ------------------------ |
---|
| 78 | ! * READ NAMELISTS. |
---|
| 79 | ! ------------------------ |
---|
| 80 | |
---|
| 81 | !---------------------------------------------------------------- |
---|
| 82 | ! Elements indispensables de SUNAM pour faire tourner RRTM dans LMDZ |
---|
| 83 | !------------------------------------------------------------------- |
---|
| 84 | CFICP='Profile' |
---|
| 85 | CFLUX='Output' |
---|
| 86 | CLIST='Listing' |
---|
| 87 | CFDDH='DHFDL' |
---|
| 88 | CNMEXP='SCM' |
---|
| 89 | TSTEP=450 |
---|
| 90 | ! MPL 29042010 - RHYDR0 - upper boundary contition for hydrostatic |
---|
| 91 | RHYDR0=1._JPRB |
---|
| 92 | ! MPL 29042010 |
---|
| 93 | ! NDLNPR : NDLNPR=0: conventional formulation of delta, i.e. ln(P(l)/P(l-1)). |
---|
| 94 | ! NDLNPR=1: formulation of delta used in non hydrostatic model, |
---|
| 95 | NDLNPR=0 |
---|
[5160] | 96 | PRINT *,'SUINIT: RHYDR0 NDLNPR',RHYDR0,NDLNPR |
---|
[1989] | 97 | |
---|
| 98 | !---------------------------------------------------------------- |
---|
| 99 | ! Elements indispensables de SUDIM pour faire tourner RRTM dans LMDZ |
---|
| 100 | !------------------------------------------------------------------- |
---|
| 101 | NDLON=klon |
---|
| 102 | NFLEVG=klev |
---|
| 103 | NPROMA=klon |
---|
| 104 | |
---|
| 105 | !------------------------------------------------------------------- |
---|
| 106 | !JV Initialize constants |
---|
| 107 | ! --------------------- |
---|
| 108 | !JV |
---|
| 109 | IF (LLTRACE) WRITE(*,*) " coucou SUINIT : avant SUCST" |
---|
| 110 | WRITE(*,FMT='('' ---------------- '')') |
---|
| 111 | WRITE(*,FMT='('' SUCST : '')') |
---|
| 112 | WRITE(*,FMT='('' ---------------- '')') |
---|
| 113 | NINDAT=20090408 !!!!! A REVOIR (MPL) |
---|
| 114 | NSSSSS=0 ! LMDZ demarre tjrs a 00h -- MPL 15.04.09 |
---|
| 115 | CALL SUCST(6,NINDAT,NSSSSS,1) |
---|
[5160] | 116 | PRINT *,'SUINIT: NINDAT, NSSSSS',NINDAT, NSSSSS |
---|
[1989] | 117 | |
---|
| 118 | IF (LLDEBUG) THEN |
---|
| 119 | WRITE(*,FMT='('' SUINIT / apres : SUCST '')') |
---|
| 120 | ENDIF |
---|
| 121 | |
---|
| 122 | |
---|
| 123 | ! ------------------------ |
---|
| 124 | ! * ALLOCATES RECUPERES DE SUALLO |
---|
| 125 | ! ------------------------ |
---|
| 126 | ALLOCATE(VDELA (MAX(JPMXLE,NFLEVG))) |
---|
| 127 | ALLOCATE(VDELB (MAX(JPMXLE,NFLEVG))) |
---|
[3435] | 128 | VDELB = 0 !ym missing init |
---|
[1989] | 129 | ALLOCATE( VC (NFLEVG) ) |
---|
[3435] | 130 | VC = 0 !ym missing init |
---|
[1989] | 131 | ALLOCATE( NLOEN (NPROMA) ) |
---|
| 132 | ALLOCATE( NLOENG (NPROMA) ) |
---|
| 133 | ALLOCATE( STZ (NFLEVG) ) |
---|
| 134 | ALLOCATE( CVDAES (NFLEVG+1)) |
---|
| 135 | ALLOCATE( CVDAEL (NFLEVG+1)) |
---|
| 136 | ALLOCATE( CVDAEU (NFLEVG+1)) |
---|
| 137 | ALLOCATE( CVDAED (NFLEVG+1)) |
---|
| 138 | ALLOCATE(RA1OVLP(NFLEVG)) |
---|
| 139 | |
---|
| 140 | ALLOCATE(STPREH(0:NFLEVG)) ! Nouvel ajout MPL 22062010 |
---|
| 141 | ALLOCATE(STPRE(NFLEVG)) |
---|
| 142 | ALLOCATE(STPHI(NFLEVG)) |
---|
| 143 | ALLOCATE(STTEM(NFLEVG)) |
---|
| 144 | ALLOCATE(STDEN(NFLEVG)) |
---|
| 145 | |
---|
| 146 | ALLOCATE(CETA(NFLEVG)) ! Nouvel ajout MPL 28062010 |
---|
| 147 | ALLOCATE(CEVAPCU(NFLEVG)) |
---|
| 148 | ALLOCATE(RMESOU(NFLEVG)) |
---|
| 149 | ALLOCATE(RMESOT(NFLEVG)) |
---|
| 150 | ALLOCATE(RMESOQ(NFLEVG)) |
---|
| 151 | |
---|
| 152 | ! ------------------------ |
---|
| 153 | ! * ALLOCATES RECUPERES DE SUGEM2 |
---|
| 154 | ! ------------------------ |
---|
| 155 | |
---|
| 156 | ALLOCATE(GEMU (NGPTOT)) ! Nouvel ajout MPL 28062010 |
---|
| 157 | ALLOCATE(GELAM (NGPTOT)) |
---|
| 158 | ALLOCATE(GELAT (NGPTOT)) |
---|
| 159 | ALLOCATE(GECLO (NGPTOT)) |
---|
| 160 | ALLOCATE(GESLO (NGPTOT)) |
---|
| 161 | ALLOCATE(GM (NGPTOT)) |
---|
| 162 | ALLOCATE(GAW (NGPTOT)) |
---|
| 163 | ! |
---|
| 164 | ! ------------------------------------------------------------------ |
---|
| 165 | |
---|
| 166 | END SUBROUTINE SUINIT |
---|