Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

Atm_RT_INI.F90 @ 3642

Last change on this file since 3642 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR

Integration of MAR physics

File size: 19.4 KB

Rev	Line
[2089]	1	subroutine Atm_RT_INI
	2
	3	!--------------------------------------------------------------------------+
	4	! Mon 24-Jun-2013 MAR \|
	5	! subroutine Atm_RT_INI initializes MAR PHYsics interface \|
	6	! with ECMWF Solar/Infrared Radiative Transfer Scheme \|
	7	! \|
	8	! version 3.p.4.1 created by H. Gallee, Tue 2-Apr-2013 \|
	9	! Last Modification by H. Gallee, Mon 24-Jun-2013 \|
	10	! \|
	11	!--------------------------------------------------------------------------+
	12	! \|
	13	! Content: CALL of - ECMWF Code initializing the Radiation Transfert \|
	14	! - ECMWF Radiation Transfert \|
	15	! \|
	16	! ECMWF Code Source: J.-J. Morcrette, 28 nov 2002 \|
	17	! \|
	18	!--------------------------------------------------------------------------+
	19
	20	use Mod_Real
	21	use Mod_PHY_RT_dat
	22	use Mod_PHY____grd
	23	use Mod_PHY_S0_grd
	24	use Mod_PHY_RT_grd
	25	use Mod_PHY_DY_kkl
	26	use Mod_SISVAT_gpt
	27
	28
	29
	30	#include "tsmbkind.h"
	31
	32
	33
	34
	35	! Global Variables (ECMWF)
	36	! ========================
	37
	38	USE PARRTM1D , ONLY : JP_LON ,JP_IDIA ,JP_FDIA ,JP_TDIA ,&
	39	& JP_LEV ,JP_LW ,JP_SW ,JP_NUA ,JP_MODE ,&
	40	& JP_AER ,JP_LEVP1
	41	USE YOMCST , ONLY : RD ,RG ,RTT ,RSIGMA ,&
	42	& RCPD ,RPI ,RDAY ,REA ,RI0 ,&
	43	& REPSM ,RMD ,RKBOL ,RNAVO ,R ,&
	44	& RLVTT ,RLSTT
	45	USE YOERAD , ONLY : NSW ,NTSW ,NRADFR ,&
	46	& LRRTM ,LINHOM ,LOIFUEC ,LTEMPDS ,LOWASYF ,&
	47	& LOWHSSS,LONEWSW ,LNEWAER ,LHVOLCA ,&
	48	& NRADIP ,NRADLP ,NOZOCL ,&
	49	& NICEOPT,NLIQOPT ,NOVLP ,NHOWINH ,RMINICE
	50	USE YOEAERD , ONLY : CVDAES ,CVDAEL ,CVDAEU ,CVDAED ,&
	51	& RCTRBGA,RCVOBGA ,RCSTBGA ,RCAEOPS ,RCAEOPL ,RCAEOPU ,&
	52	& RCAEOPD,RCTRPT ,RCAEADK ,RCAEADM ,RCAEROS
	53	USE YOERDI , ONLY : RCH4 ,RN2O ,RO3 ,RCFC11 ,&
	54	& RCFC12
	55	USE YOERDU , ONLY : RCDAY ,R10E ,DIFF ,REPLOG ,&
	56	& REPSC ,REPSCO ,REPSCQ ,REPSCT ,REPSCW ,&
	57	& NTRAER
	58	USE YOETHF , ONLY : R2ES ,R3LES ,R3IES ,R4LES ,&
	59	& R4IES ,R5LES ,R5IES ,R5ALVCP ,R5ALSCP ,&
	60	& RALVDCP,RALSDCP ,RTWAT ,RTICE ,RTICECU
	61
	62
	63
	64
	65	IMPLICIT NONE
	66
	67
	68
	69
	70	! INTERNAL VARIABLES
	71	! ==================
	72
	73	integer :: i ,j ,ikl
	74
	75	! For Use in radlsw
	76	! ^^^^^^^^^^^^^^^^^
	77	INTEGER_M :: KIDIA ,KFDIA ,KTDIA ,KLON ,KLEV
	78	INTEGER_M :: KMODE ,KAER ,KSW
	79
	80	INTEGER_M :: KBOX ,NBOX
	81	INTEGER_M :: NDUMP ,ILWRAD
	82	!
	83	REAL_B :: PRII05
	84
	85	REAL_B :: PAER5(JP_LON,JP_AER,JP_LEV) ! Aerosol Optical Depth
	86	REAL_B :: PALBD5(JP_LON,JP_SW)
	87	REAL_B :: PALBP5(JP_LON,JP_SW)
	88	!
	89	REAL_B :: PAPH5(JP_LON,JP_LEVP1)
	90	REAL_B :: PAP5(JP_LON,JP_LEV)
	91	!
	92	REAL_B :: PCCO25
	93	REAL_B :: PCLFR5(JP_LON,JP_LEV)
	94	REAL_B :: PDP5(JP_LON,JP_LEV)
	95	REAL_B :: PEMIS5(JP_LON)
	96	REAL_B :: PEMIW5(JP_LON)
	97	REAL_B :: PLSM5(JP_LON)
	98	REAL_B :: PMU05(JP_LON)
	99	REAL_B :: POZON5(JP_LON,JP_LEV)
	100	REAL_B :: PQ5(JP_LON,JP_LEV)
	101	!
	102	REAL_B :: PQIWP5(JP_LON,JP_LEV)
	103	REAL_B :: PQLWP5(JP_LON,JP_LEV) ! Dropplets Concentration
	104	REAL_B :: PSQIW5(JP_LON,JP_LEV) ! Ice Crystals Concentration
	105	REAL_B :: PSQLW5(JP_LON,JP_LEV) !
	106	REAL_B :: PQS5(JP_LON,JP_LEV)
	107	REAL_B :: PQRAIN5(JP_LON,JP_LEV)
	108	REAL_B :: PRAINT5(JP_LON,JP_LEV)
	109	REAL_B :: PRLVRI5(JP_LON,JP_LEV)
	110	REAL_B :: PRLVRL5(JP_LON,JP_LEV)
	111	REAL_B :: PTH5(JP_LON,JP_LEVP1)
	112	REAL_B :: PT5(JP_LON,JP_LEV)
	113	REAL_B :: PTS5(JP_LON)
	114	REAL_B :: PNBAS5(JP_LON)
	115	REAL_B :: PNTOP5(JP_LON)
	116
	117	! For Use in SUCLD
	118	! ^^^^^^^^^^^^^^^^
	119	REAL_B :: ZETA(JP_LEV)
	120	REAL_B :: ZETAH(JP_LEVP1)
	121
	122	! For Use in SUOVLP
	123	! ^^^^^^^^^^^^^^^^^
	124	REAL_B :: ZTVIR
	125	REAL_B :: ZFACT
	126	REAL_B :: ZAZ(JP_LEV)
	127	REAL_B :: ZAZH(JP_LEVP1)
	128
	129
	130
	131	! Local Variables
	132	! ----------------
	133
	134	REAL_B :: RTIMTR
	135
	136	INTEGER_M :: JL , JK , JAER ,JNU
	137	INTEGER_M :: KULOUT , NINDAT , NSSSSS ,KPRTLEV
	138	INTEGER_M :: IYR , MONTH , IDAY ,IMINUT
	139	INTEGER_M :: KPRINT , SKSHIFT
	140
	141
	142
	143
	144	! Load External Functions
	145	! =======================
	146
	147	#include "fctast.h"
	148	#include "fcttim.h"
	149	#include "fcttre.h"
	150
	151
	152	KPRTLEV= 1
	153	KPRINT = 1
	154	SKSHIFT= 0
	155
	156
	157
	158
	159	! Time Base
	160	! =========
	161
	162	NINDAT = min(yearTU,2004)10000+mon_TU100+Day_TU ! Date in the form yyyyMMdd
	163	NSSSSS = HourTU * 3600+minuTU* 60+sec_TU ! Nb of second since day Begin
	164	IYR = NAA(NINDAT)
	165	MONTH = NMM(NINDAT)
	166	IDAY = NDD(NINDAT)
	167	RTIMTR = RTIME(IYR,MONTH,IDAY,NSSSSS)
	168	IMINUT = INT(FLOAT(NSSSSS)/60.)
	169
	170
	171
	172
	173	! Basic Initialization
	174	! ====================
	175
	176	! Dimensions (auxiliary variables)
	177	! --------------------------------
	178
	179	KIDIA = 1 ! DO'NT CHANGE
	180	KFDIA = JP_LON ! Nb Columns
	181	KTDIA = 1 !
	182	KLON = JP_LON ! Nb Columns
	183	KLEV = JP_LEV ! Nb Levels
	184	KMODE = JP_MODE ! Used in Planck Fcts Specification
	185	KAER = JP_AER !
	186
	187	! Nb of Solar Spectral Intervals
	188	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	189	KSW = JP_SW ! SW Nb of Spectral Intervals (max is JP_SW=6)
	190	NSW = JP_SW ! SW Nb of Spectral Intervals (max is JP_SW=6)
	191	NTSW = JP_SW ! SW Nb of Spectral Intervals (max is JP_SW=6)
	192
	193	KBOX = 0 ! \VER
	194	NBOX = 1 ! \VER
	195
	196	ILWRAD = 1 ! 0: Morcrette, 1991 operational before 20000627
	197	! 1: Mlawer et al., 1997 now ECMWF-operational
	198	! 2: Morcrette, 1991 original as in ERA-15'
	199
	200	NDUMP = 3 ! No Print
	201	! NDUMP = 2 ! 1D Results
	202	! NDUMP = 1 ! Debug
	203	! NDUMP = 0 ! ALL
	204	KULOUT = 6 ! Output Device for SUCST
	205
	206
	207	! Verification of the Dimensions
	208	! ------------------------------
	209
	210	write(6,*) 'INITIALISATION OF ECMWF RADIATIVE TRANSFERT: BEGIN'
	211
	212	write(6,*) 'CONTROL'
	213	write(6,*) 'JP_LON=',JP_LON
	214	write(6,*) 'kcolp',kcolp
	215	IF (kcolp .GT.JP_LON) THEN
	216	write(6,6001) kcolp ,JP_LON
	217	6001 format(' @!#& BAD dimensions SET-UP (kcolp > JP_LON); = ('&
	218	& ,i3,',',i3,')')
	219	STOP
	220	END IF
	221	IF (mzp .NE.JP_LEV) THEN
	222	write(6,6002) mzp ,JP_LEV
	223	6002 format(' @!#& BAD dimensions SET-UP (mzp /= JP_LEV); = ('&
	224	& ,i3,',',i3,')')
	225	STOP
	226	END IF
	227	IF (naero .NE.JP_AER) THEN
	228	write(6,6003) naero ,JP_AER
	229	6003 format(' @!#& BAD dimensions SET-UP (naero/= JP_AER); = ('&
	230	& ,i3,',',i3,')')
	231	STOP
	232	END IF
	233
	234
	235	! Mathematical Constants
	236	! ----------------------
	237
	238	REPLOG = 1.E-12 ! Minimum Logarithm Argument
	239	REPSC = 1.E-12
	240	REPSCO = 1.E-12
	241	REPSCQ = 1.E-12
	242	REPSCT = 1.E-12
	243	REPSCW = 1.E-12
	244
	245
	246	! Switches (general)
	247	! ------------------
	248
	249	LONEWSW= .TRUE. ! .TRUE. SWSN radiative routine is used
	250	IF (ILWRAD.EQ.1) THEN
	251	LRRTM = .TRUE. ! .TRUE. RRTM radiative routine is used
	252	ELSE
	253	LRRTM = .FALSE. ! .FALSE.RRTM radiative routine is NOT used
	254	END IF
	255	LTEMPDS= .FALSE. ! .TRUE. ALLOWS FOR SURF. T DISCONTIN. IN RAD.COMPUT.
	256
	257	NTRAER = 19 ! NUMBER OF TRANSMISSION FUNCTIONS W OR W/O AEROSOLS
	258
	259
	260	! Switches (Clouds Optical Properties)
	261	! ------------------------------------
	262
	263	LINHOM = .FALSE. ! Tiedke (1995) correct. factor (0.7) of tau not used
	264	NHOWINH= 2 ! Tau correction factor: exp(-(sig/tau)^2)
	265	! (used if LINHOM = .TRUE.)
	266	LOIFUEC= .FALSE. ! .FALSE. IF ICE CLOUDS AS EBERT-CURRY (LW & SW)
	267	LOWASYF= .FALSE. ! .FALSE. IF WATER CLOUDS AS FOUQUART (SW)
	268	LOWHSSS= .FALSE. ! .FALSE. IF WATER CLOUDS AS SMITH-SHI (LW)
	269	NRADIP = 3 ! Ice effective Radius:
	270	! 0 fixed 40 microns
	271	! 1 f(T) 40 - 130 microns
	272	! 2 f(T) 30 - 60 microns Jakob-Klein
	273	! 3 f(T,IWC) Sun-Rikus, 1999
	274	RMINICE= 15 ! Minimum Diameter for Ice Particles (micronm)
	275	! Needed only if NRADIP = 3
	276	! (see von Walden et al., 2003 (Oct) Tab. 2 p.1393)
	277	NRADLP = 0 ! Liquid effective Radius: f(Pressure)
	278	! 0 effective radius - liquid as f(Pressure)
	279	! 1 fixed 10 microns over land, 13 over ocean
	280	! 2 computed from LWC Martin et al, 1994
	281	NLIQOPT= 1 ! Cloud Optical Properties (Water): 1=ECMWF Operat.
	282	! 0 LW: Smith-Shi, 1992; SW: Fouquart, 1987
	283	! 1 LW: Savijarvi, 1997; SW: Slingo , 1989
	284	! 2 LW: Lindner,Li, 2000; SW: Slingo , 1989
	285	NICEOPT= 1 ! Cloud Optical Properties (Ice) : 1=ECMWF Operat.
	286	! 0 LW: Smith,Shi , 1992; SW: Ebert-Curry, 1992
	287	! 1 LW: Ebert,Curry, 1992; SW: Ebert-Curry, 1992
	288	! 2 LW: Fu,Liou , 1993; SW: Fu,Liou , 1993
	289	! 3 LW: Fu et al. , 1998; SW: Fu , 1996
	290	NOVLP = 2 ! CLOUD OVERLAP CONFIGURATION:
	291	! 1=MRN, 2=MAX, 3=RAN, 4=Hogan
	292
	293
	294	! Switches (Aerosols/O3)
	295	! ----------------------
	296
	297	LNEWAER= .TRUE. ! Climatology of Aerosols: TEGEN ET AL. 1997 / GADS
	298	LHVOLCA= .TRUE. ! .TRUE. IF GISS HISTORY OF VOLCANIC AEROSOLS IS ON
	299	! NOZOCL = -1 ! TESTS the vertical quadrature (NO absorber)
	300	! NOZOCL = 0 ! whatever is read for O3 as input profile
	301	! NOZOCL = 1 ! OLD ECMWF O3 climatology and aerosols
	302	NOZOCL = 2 ! Fortuin-Langematz O3 climatology and aerosols
	303	! NOZOCL = 3 ! OLD ECMWF O3 climatology and NO aerosols
	304	! NOZOCL = 4 ! Fortuin-Langematz O3 climatology and NO aerosols
	305
	306
	307	! BASIC CONSTANTS
	308	! ---------------
	309
	310	! *****
	311	CALL SUCST (KULOUT, NINDAT, NSSSSS, KPRTLEV) ! Initialize common YOMCST
	312	! ***** ! (Basic Constants)
	313	! Initialize common YOMRIP
	314	! (only date and time)
	315
	316	! YOENCST - THERMODYNAMIC TRANSITION OF PHASE
	317	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	318	RTT=273.16 !
	319
	320	! YOETHF - DERIVED CONSTANTS SPECIFIC TO ECMWF THERMODYNAMICS
	321	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	322	RTWAT=RTT !
	323	RTICE=RTT-23. !
	324
	325	! YOERDU - CONTROL, PARAMETERS AND SECURITY IN RADIATION
	326	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	327	RCDAY = RDAY * RG / RCPD !
	328	R10E = 0.4342945 ! DECIMAL / NATURAL
	329	! LOG. FACTOR
	330	DIFF = 1.66 ! DIFFUSIVITY FACTOR
	331
	332
	333	! Space/Time Independant Coefficients
	334	! -----------------------------------
	335
	336	CALL SURDI ! ECMWF Surface Albedo, Emissivity
	337	CALL SULWN ! Initialize common YOELW (new LW Coeff.)
	338	CALL SUOLW ! Initialize common YOELW (old LW Coeff.)
	339
	340	! Initialization routine for RRTM
	341	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	342	CALL SURRTAB ! AER'S RRTM LW RADIATION
	343	CALL SURRTPK ! Initialize common YOERRTWN
	344	! (k-coefficients in spectral intervals)
	345	CALL SURRTRF ! Initialize common YOERRTRF
	346	! (RRTM Reference Atmosphere)
	347	CALL SURRTFTR ! Initialize common YOERRTRF
	348
	349	! RRTM routine ! BAND [cm-1] ! low ! high
	350	! ----------------------+---------------+---------------+---------
	351	CALL RRTM_KGB1 ! 1: 10- 250 ! H2O ! H2O
	352	CALL RRTM_KGB2 ! 2: 250- 500 ! H2O ! H2O
	353	CALL RRTM_KGB3 ! 3: 500- 630 ! H2O,CO2 ! H2O,CO2
	354	CALL RRTM_KGB4 ! 4: 630- 700 ! H2O,CO2 ! O3,CO2
	355	CALL RRTM_KGB5 ! 5: 700- 820 ! H2O,CO2 ! O3,CO2
	356	CALL RRTM_KGB6 ! 6: 820- 980 ! H2O ! nothing
	357	CALL RRTM_KGB7 ! 7: 980-1080 ! H2O,O3 ! O3
	358	CALL RRTM_KGB8 ! 8: 1080-1180 ! (i.e.>~300mb) ! O3
	359	! ! H2O !
	360	CALL RRTM_KGB9 ! 9: 1180-1390 ! H2O,CH4 ! CH4
	361	CALL RRTM_KGB10 ! 10: 1390-1480 ! H2O ! H2O
	362	CALL RRTM_KGB11 ! 11: 1480-1800 ! H2O ! H2O
	363	CALL RRTM_KGB12 ! 12: 1800-2080 ! H2O,CO2 ! nothing
	364	CALL RRTM_KGB13 ! 13: 2080-2250 ! H2O,N2O ! nothing
	365	CALL RRTM_KGB14 ! 14: 2250-2380 ! CO2 ! CO2
	366	CALL RRTM_KGB15 ! 15: 2380-2600 ! N2O,CO2 ! nothing
	367	CALL RRTM_KGB16 ! 16: 2600-3000 ! H2O,CH4 ! nothing
	368
	369	! Reduce absorption coefficient data from 256 to 140 g-points
	370	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	371	CALL RRTM_INIT_140GP
	372
	373	! Initialization routine for SW (6 spectral interval resolution)
	374	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	375	CALL SUSWN ( NTSW , KSW ) ! Initialize common YOESW
	376
	377
	378	write(6,*) 'INITIALISATION OF ECMWF RADIATIVE TRANSFERT: END '
	379
	380
	381
	382
	383	! Radiation: Global (Time dependant) Parameters
	384	! =============================================
	385
	386	PRII05 = RI0 ! INSOLATION
	387	! PRII05 = RI0/(dST_UA*dST_UA) ! INSOLATION (dST_UA: Distance Soleil-Terre [UA])
	388	PCCO25 = 360.E-06*44./29. ! CONCENTRATION IN CO2 [Pa/Pa]
	389
	390
	391
	392
	393	! Surface Properties
	394	! ==================
	395
	396	! Land/sea Mask
	397	! -------------
	398
	399	! Martin CONTROL
	400	!PRINT*, 'CONTROL dans Atm_RT_INI'
	401	!PRINT*,'SIZE(psa_DY)=)',SIZE(psa_DY)
	402	!PRINT*, 'psa_DY=',psa_DY
	403	! Martin CONTROL
	404
	405	DO jl=1,KLON
	406	ikl = min(jl,kcolp)
	407	i = ii__AP(ikl)
	408	j = jj__AP(ikl)
	409	PLSM5 (JL) = 1 - MaskSV_gpt(ikl)
	410
	411
	412
	413
	414	! Atmospheric Thermodynamics (Time and Space dependant)
	415	! =====================================================
	416
	417	! Pressure Distribution
	418	! --------------------------
	419
	420	JK=1+KLEV
	421	PAPH5 (JL,JK) = (psa_DY(ikl) + pt__DY) * 1000. ! Pressure (Layer Interface)[Pa]
	422	DO JK=1,KLEV
	423	PAPH5 (JL,JK) = (psa_DY(ikl) * sigma(JK) + pt__DY) * 1000. ! Pressure (Layer) [Pa]
	424	PAP5 (JL,JK) = (psa_DY(ikl) * sigmi(JK) + pt__DY) * 1000. ! Pressure (Layer Interface)[Pa]
	425	PDP5 (JL,JK) = psa_DY(ikl) dsigmi(JK) 1000. ! Pressure (Layer Thickness)[Pa]
	426
	427
	428	! Temperature Distribution
	429	! --------------------------
	430
	431	PT5 (JL,JK) = Ta__DY(ikl,JK)
	432
	433
	434	! Water Species Distribution
	435	! --------------------------
	436
	437	PQ5 (JL,JK) = qv__DY(ikl,JK) ! Water Vapor Concentr. [kg/kg]
	438
	439	ENDDO
	440	ENDDO
	441
	442
	443
	444
	445	! Initialization (Climatologies, Time independant)
	446	! ================================================
	447
	448	! Aerosols Radiative Characteristics (YOEAER)
	449	! ----------------------------------
	450
	451	! *******
	452	CALL SUAERL ! Aerosols LW Radiative Charact.
	453	CALL SUAERSN ( NTSW , KSW ) ! Aerosols SW Radiative Charact.
	454	! *******
	455
	456
	457	! Aerosols Optical Thickness Horizontal Distribution (model grid
	458	! -------------------------------------------------- independant)
	459
	460	! ********
	461	CALL SUAERH !
	462	CALL SUECAEBC ! BLACK CARBON (URBAN/FOREST FIRE ORIGIN)
	463	CALL SUECAEOR ! ORGANIC-TYPE
	464	CALL SUECAESD ! SOIL-DUST ORIGIN
	465	CALL SUECAESS ! SEA -SALT ORIGIN
	466	CALL SUECAESU ! SULFATE-TYPE
	467	! ********
	468
	469
	470	! Clouds (YOECLD)
	471	! ---------------
	472
	473	DO jk=1,klev
	474	ZETA(jk) =PAP5(1,jk) /PAPH5(1,klev+1)
	475	ENDDO
	476	DO jk=1,klev+1
	477	ZETAH(jk)=PAPH5(1,jk)/PAPH5(1,klev+1)
	478	ENDDO
	479	! Martin Control
	480	!PRINT*, 'PAPH5(1,:)=',PAPH5(1,:)
	481	!PRINT*, 'ZETAH=',ZETAH
	482	! Martin Control
	483
	484	! *****
	485	CALL SUCLD ( KLEV , ZETA )
	486	! *****
	487
	488
	489	! Cloud Optical Parameters SW/LW (all parameterizations)
	490	! ------------------------------------------------------
	491
	492	! *******
	493	CALL SUCLOPN ( NTSW , KSW , KLEV ) ! Initialize YOECLOP
	494	! *******
	495
	496
	497	! Radar Reflectivity
	498	! ------------------
	499
	500	ZAZH(KLEV+1)= 0.
	501	ZAZ(1) =100000.
	502	JL=KIDIA
	503	DO jk=KLEV,2,-1
	504	ZTVIR = PT5(JL,jk) /(1.-0.608*PQ5(jl,jk))
	505	ZFACT = LOG(PAPH5(jl,jk+1))- LOG(PAPH5(jl,jk))
	506	ZAZH(jk) = ZAZH(jk+1) + R * ZTVIR/(RMDRG)ZFACT
	507	ZAZ(jk) = 0.5(ZAZH(jk+1)+ZAZH(jk))1000.
	508	END DO
	509
	510	! ******
	511	CALL SUOVLP ( KLEV , ZAZ ) ! Initialize ALPHA1 (%radar refl.)
	512	! ****** ! (Hogan & Illingsworth, 1999)
	513
	514
	515	! NO Absorber
	516	! -----------
	517
	518	IF (NOZOCL.EQ.-1) THEN
	519	RCH4 = 1.E-18
	520	RN2O = 1.E-18
	521	RO3 = 1.E-18
	522	RCFC11 = 1.E-18
	523	RCFC12 = 1.E-18
	524	PCCO25 = 1.E-18
	525	DO jk=1,klev
	526	DO jl=KIDIA,KFDIA
	527	POZON5(JL,JK) = 0.
	528	ENDDO
	529	ENDDO
	530	DO JK=1,KLEV
	531	DO JAER=1,KAER
	532	DO JL=KIDIA,KFDIA
	533	PAER5(JL,JAER,JK)= ZEPAER
	534	ENDDO
	535	ENDDO
	536	ENDDO
	537	END IF
	538
	539
	540
	541	return
	542	end subroutine Atm_RT_INI

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format