Context Navigation

suecrad.F90 @ 5503

Last change on this file since 5503 was 5160, checked in by abarral, 6 months ago
Put .h into modules
File size: 56.8 KB

Rev	Line
[5159]	1
[5133]	2	! $Id: suecrad.F90 4251 2022-09-20 00:22:43Z fhourdin $
[5159]	3
[5154]	4	SUBROUTINE SUECRAD (KULOUT, KLEV, PETAH)
[5133]	5
[5154]	6	!**** SUECRAD - INITIALIZE COMMONS YOERxx CONTROLLING RADIATION
[5133]	7
[5154]	8	! PURPOSE.
	9	! --------
	10	! INITIALIZE YOERAD, THE COMMON THAT CONTROLS THE
	11	! RADIATION OF THE MODEL, AND YOERDU THAT INCLUDES
	12	! ADJUSTABLE PARAMETERS FOR RADIATION COMPUTATIONS
[5133]	13
[5154]	14	!** INTERFACE.
	15	! ----------
	16	! CALL SUECRAD FROM SUPHEC
	17	! ------- ------
[5133]	18
[5154]	19	! EXPLICIT ARGUMENTS :
	20	! --------------------
	21	! NONE
[5133]	22
[5154]	23	! IMPLICIT ARGUMENTS :
	24	! --------------------
	25	! COMMONS YOERAD, YOERDU
[5133]	26
[5154]	27	! METHOD.
	28	! -------
	29	! SEE DOCUMENTATION
[5133]	30
[5154]	31	! EXTERNALS.
	32	! ----------
	33	! SUAER, SUAERH, SUAERV, SULW, SUSW, SUOCST, SUSAT
	34	! SUAERL, SUAERSN, SUSRTAER, SRTM_INIT, SUSRTCOP
[5133]	35
[5154]	36	! REFERENCE.
	37	! ----------
	38	! ECMWF Research Department documentation of the IFS
[5133]	39
[5154]	40	! AUTHOR.
	41	! -------
	42	! JEAN-JACQUES MORCRETTE ECMWF
[5133]	43
[5154]	44	! MODIFICATIONS.
	45	! --------------
	46	! ORIGINAL : 88-12-15
	47	! P.COURTIER AND M.HAMRUD NAME SURAD ALREADY USED
	48	! Modified 93-11-15 by Ph. Dandin : FMR scheme with MF
	49	! Modified 95-12 by PhD : Cloud overlapping hypothesis for FMR
	50	! 980317 JJMorcrette clean-up (NRAD, NFLUX)
	51	! 000118 JJMorcrette variable concentr. uniformly mixed gases
	52	! 990525 JJMorcrette GISS volcanic and new tropospheric aerosols
	53	! 990831 JJMorcrette RRTM
	54	! R. El Khatib 01-02-02 proper initialization of NFRRC moved in SUCFU
	55	! 010129 JJMorcrette clean-up LERAD1H, NLNGR1H
	56	! 011105 GMozdzynski support new radiation grid
	57	! 011005 JJMorcrette CCN --> Re Water clouds
	58	! R. El Khatib 01-02-02 LRRTM=lecmwf by default
	59	! 020909 GMozdzynski support NRADRES to specify radiation grid
	60	! 021001 GMozdzynski support on-demand radiation communications
	61	! 030422 GMozdzynski automatic min-halo
	62	! 030501 JJMorcrette new radiation grid on, new aerosols on (default)
	63	! 030513 JJMorcrette progn. O3 / radiation interactions off (default)
	64	! M.Hamrud 01-Oct-2003 CY28 Cleaning
	65	! 050315 JJMorcrette prog.aerosols v1
	66	! 041214 JJMorcrette SRTM
	67	! 050111 JJMorcrette new cloud optical properties
	68	! 050415 GMozdzynski Reduced halo support for radiation interpolation
	69	! 051004 JJMorcrette UV surface radiation processor
	70	! 051220 JJMorcrette SRTM112g+LWSCAT+UVprocessor+(bgfx:swclr, radaca)
	71	! 060510 JJMorcrette MODIS albedo (UVis, NIR)x(parallel+diffuse)
	72	! 060510 JJMorcrette MODIS albedo (UVis, NIR)x(parallel+diffuse)
	73	! JJMorcrette 20060721 PP of clear-sky PAR and TOA incident solar radiation
	74	! 060625 JJMorcrette MODIS albedo (UVis, NIR)x(parallel+diffuse)
	75	! 060726 JJMorcrette McICA default operational configuration
	76	! ------------------------------------------------------------------
[5133]	77
[5154]	78	USE PARKIND1, ONLY: JPIM, JPRB
	79	USE YOMHOOK, ONLY: LHOOK, DR_HOOK
[5133]	80
[5154]	81	USE PARDIM, ONLY: JPMXGL
	82	USE PARRRTM, ONLY: JPLAY
	83	USE PARSRTM, ONLY: JPGPT
	84	USE YOMCT0, ONLY: LOUTPUT, NPRINTLEV, LALLOPR, &
	85	& NPROC, N_REGIONS_NS, N_REGIONS_EW
	86	USE YOMDIM, ONLY: NDLON, NSMAX, NDGENL, &
	87	& NDGSAL, NDGLG, NDGSAG, NDGENG, NDSUR1, &
	88	& NDLSUR, NDGSUR, NGPBLKS, NFLEVG, NPROMA
	89	USE YOMCT0B, ONLY: LECMWF
	90	USE YOMDYN, ONLY: TSTEP
	91	! Ce qui concerne NULRAD commente par MPL le 15.04.09
	92	!USE YOMLUN , ONLY : NULNAM ,NULRAD ,NULOUT
	93	USE YOMLUN, ONLY: NULRAD, NULOUT
	94	USE YOMCST, ONLY: RDAY, RG, RCPD, RPI, RI0
	95	USE YOMPHY, ONLY: LMPHYS, LRAYFM, LRAYFM15
	96	USE YOEPHY, ONLY: LEPHYS, LERADI, LE4ALB
	97	USE YOERDI, ONLY: RCCO2, RCCH4, RCN2O, RCCFC11, RCCFC12, RSOLINC
	98	USE YOERAD, ONLY: NAER, NOZOCL, &
	99	& NRADFR, NRADPFR, NRADPLA, NRINT, &
	100	& NRADNFR, NRADSFR, NOVLP, NRPROMA, &
	101	!& NLW ,NSW ,NTSW ,NCSRADF ,&
	102	! NSW mis dans .def MPL 20140211
	103	& NLW, NTSW, NCSRADF, &
	104	& NMODE, NLNGR1H, NSWNL, NSWTL, NUV, &
	105	& LERAD1H, LERADHS, LEPO3RA, LRADLB, LONEWSW, &
	106	& LCCNL, LCCNO, &
	107	& LECSRAD, LHVOLCA, LNEWAER, LRRTM, LSRTM, LDIFFC, &
	108	& NRADINT, NRADRES, CRTABLEDIR, CRTABLEFIL, &
	109	& NICEOPT, NLIQOPT, NRADIP, NRADLP, NINHOM, NLAYINH, &
	110	& LRAYL, LOPTRPROMA, &
	111	& RCCNLND, RCCNSEA, RLWINHF, RSWINHF, RRe2De, &
	112	& RPERTOZ, NPERTOZ, NMCICA, &
	113	& LNOTROAER, NPERTAER, LECO2VAR, LHGHG, NHINCSOL, NSCEN, &
	114	& LEDBUG
	115	USE YOERDU, ONLY: NUAER, NTRAER, RCDAY, R10E, &
	116	& REPLOG, REPSC, REPSCO, REPSCQ, REPSCT, &
	117	& REPSCW, DIFF
	118	USE YOEAERD, ONLY: CVDAES, CVDAEL, CVDAEU, CVDAED, &
	119	& RCAEOPS, RCAEOPL, RCAEOPU, RCAEOPD, RCTRBGA, &
	120	& RCVOBGA, RCSTBGA, RCTRPT, RCAEADM, RCAEROS, &
	121	& RCAEADK
	122	USE YOE_UVRAD, ONLY: JUVLAM, LUVPROC, LUVTDEP, LUVDBG, NRADUV, NUVTIM, RUVLAM, RMUZUV
[5133]	123
[5154]	124	USE YOMMP, ONLY: MYPROC, NPRCIDS, LSPLIT, NAPSETS, &
	125	& NPTRFLOFF, NFRSTLOFF, MYFRSTACTLAT, MYLSTACTLAT, &
	126	& NSTA, NONL, NPTRFRSTLAT, NFRSTLAT, NLSTLAT, &
	127	& MY_REGION_NS, MY_REGION_EW, NGLOBALINDEX, &
	128	& NRISTA, NRIONL, NRIOFF, NRIEXT, NRICORE, &
	129	& NRISENDPOS, NRIRECVPOS, NRISENDPTR, NRIRECVPTR, &
	130	& NARIB1, NRIPROCS, NRIMPBUFSZ, NRISPT, NRIRPT, &
	131	& NRICOMM, &
	132	& NROSTA, NROONL, NROOFF, NROEXT, NROCORE, &
	133	& NROSENDPOS, NRORECVPOS, NROSENDPTR, NRORECVPTR, &
	134	& NAROB1, NROPROCS, NROMPBUFSZ, NROSPT, NRORPT, &
	135	& NROCOMM
	136	USE YOMGC, ONLY: GELAT, GELAM
	137	USE YOMLEG, ONLY: RMU, RSQM2
	138	USE YOMSC2, ONLY: &
	139	& NRIWIDEN, NRIWIDES, NRIWIDEW, NRIWIDEE, &
	140	& NROWIDEN, NROWIDES, NROWIDEW, NROWIDEE
	141	USE YOMGEM, ONLY: NGPTOT, NGPTOTG, NGPTOTMX, NLOENG
	142	USE YOMTAG, ONLY: MTAGRAD
	143	USE YOMPRAD, ONLY: LODBGRADI, LODBGRADL, RADGRID, &
	144	& LRADONDEM
	145	USE YOMRADF, ONLY: EMTD, TRSW, EMTC, TRSC, &
	146	& SRSWD, SRLWD, SRSWDCS, SRLWDCS, SRSWDV, &
	147	& SRSWDUV, EDRO, SRSWPAR, SRSWUVB, SRSWPARC, SRSWTINC, &
	148	& EMTU, RMOON
	149	! Commente par MPL 26.11.08
	150	!USE YOPHNC , ONLY : LERADN2
	151	! MPLefebvre 6-11-08 commente tout ce qui concerne MPL_MODULE
	152	!USE MPL_MODULE , ONLY : MPL_BROADCAST, MPL_SEND, MPL_RECV
	153	USE YOM_YGFL, ONLY: YO3
	154	!!!!! A REVOIR (MPL) NDLNPR devrait etre initialise dans sudyn.F90
	155	USE YOMDYN, ONLY: NDLNPR
	156	USE lmdz_clesphys
[5133]	157
[5154]	158	IMPLICIT NONE
[5133]	159
[5154]	160	INTEGER(KIND = JPIM), INTENT(IN) :: KLEV
	161	INTEGER(KIND = JPIM), INTENT(IN) :: KULOUT
	162	REAL(KIND = JPRB), INTENT(IN) :: PETAH(KLEV + 1)
	163	! LOCAL ARRAYS FOR THE PURPOSE OF READING NAMRGRI (RADIATION GRID)
	164	INTEGER(KIND = JPIM) :: NRGRI(JPMXGL)
[5133]	165
[5154]	166	INTEGER(KIND = JPIM) :: IDGL, INBLW, IRADFR, IST1HR, ISTNHR, IDIR, IFIL
	167	INTEGER(KIND = JPIM) :: IRIRPTSUR, IRISPTSUR, IRIMAPLEN
	168	INTEGER(KIND = JPIM) :: JLON, JGLAT, JGL, JGLSUR, IDLSUR, IOFF, ILAT, ISTLON, IENDLON
	169	INTEGER(KIND = JPIM) :: IRORPTSUR, IROSPTSUR, IROMAPLEN
	170	INTEGER(KIND = JPIM) :: ILBRLATI, IUBRLATI, IGLGLO, IDUM, IU
	171	INTEGER(KIND = JPIM) :: J, JROC, IGPTOT
	172	INTEGER(KIND = JPIM) :: IROWIDEMAXN, IROWIDEMAXS, IROWIDEMAXW, IROWIDEMAXE
	173	INTEGER(KIND = JPIM) :: IRIWIDEMAXN, IRIWIDEMAXS, IRIWIDEMAXW, IRIWIDEMAXE
	174	INTEGER(KIND = JPIM) :: IARIB1MAX, IAROB1MAX
	175	INTEGER(KIND = JPIM) :: IWIDE(10)
	176	INTEGER(KIND = JPIM) :: ILATS_DIFF_F, ILATS_DIFF_C
	177	INTEGER(KIND = JPIM), PARAMETER :: JP_MIN_HALO = 5
	178	INTEGER(KIND = JPIM) :: ISW, JUV, IDAYUV
[5133]	179
[5154]	180	LOGICAL :: LLINEAR_GRID
	181	LOGICAL :: LLDEBUG, LLP
[5133]	182
[5154]	183	REAL(KIND = JPRB) :: ZSTPHR, ZTSTEP, ZGEMU, ZLON, ZD1, ZD2, ZD3, ZD4, ZD5, ZD6
	184	REAL(KIND = JPRB) :: ZMINRADLAT, ZMAXRADLAT, ZMINRADLON, ZMAXRADLON
	185	REAL(KIND = JPRB) :: ZMINMDLLAT, ZMAXMDLLAT, ZMINMDLLON, ZMAXMDLLON
	186	REAL(KIND = JPRB) :: ZLAT
	187	!REAL(KIND=JPRB) :: RLATVOL, RLONVOL
[5133]	188
[5154]	189	CHARACTER (LEN = 300) :: CLFN
	190	INTEGER(KIND = JPIM), PARAMETER :: JPIOMASTER = 1
[5133]	191
[5154]	192	INTEGER(KIND = JPIM), ALLOCATABLE :: IRISENDPOS(:)
	193	INTEGER(KIND = JPIM), ALLOCATABLE :: IRIRECVPOS(:)
	194	INTEGER(KIND = JPIM), ALLOCATABLE :: IRISENDPTR(:)
	195	INTEGER(KIND = JPIM), ALLOCATABLE :: IRIRECVPTR(:)
	196	INTEGER(KIND = JPIM), ALLOCATABLE :: IRICOMM(:)
	197	INTEGER(KIND = JPIM), ALLOCATABLE :: IRIMAP(:, :)
	198	INTEGER(KIND = JPIM), ALLOCATABLE :: IROSENDPOS(:)
	199	INTEGER(KIND = JPIM), ALLOCATABLE :: IRORECVPOS(:)
	200	INTEGER(KIND = JPIM), ALLOCATABLE :: IROSENDPTR(:)
	201	INTEGER(KIND = JPIM), ALLOCATABLE :: IRORECVPTR(:)
	202	INTEGER(KIND = JPIM), ALLOCATABLE :: IROCOMM(:)
	203	INTEGER(KIND = JPIM), ALLOCATABLE :: IROMAP(:, :)
	204	INTEGER(KIND = JPIM), ALLOCATABLE :: IGLOBALINDEX(:)
[5133]	205
[5154]	206	REAL(KIND = JPRB), ALLOCATABLE :: ZLATX(:)
	207	REAL(KIND = JPRB), ALLOCATABLE :: ZLONX(:)
	208	REAL(KIND = JPRB) :: ZHOOK_HANDLE
[5133]	209
[5154]	210	INTERFACE
[5133]	211	#include "setup_trans.h"
	212	#include "trans_inq.h"
[5154]	213	END INTERFACE
[5133]	214
	215	#include "abor1.intfb.h"
	216	#include "posnam.intfb.h"
	217	#include "rrtm_init_140gp.intfb.h"
	218
	219	#include "rdcset.intfb.h"
	220	#include "suaerh.intfb.h"
	221	#include "suaerl.intfb.h"
	222	#include "suaersn.intfb.h"
	223	#include "suaerv.intfb.h"
	224	#include "suclopn.intfb.h"
	225	#include "suecradi.intfb.h"
	226	#include "suecradl.intfb.h"
	227	#include "sulwn.intfb.h"
	228	#include "sulwneur.intfb.h"
	229	#include "suovlp.intfb.h"
	230	#include "surdi.intfb.h"
	231	#include "surrtab.intfb.h"
	232	#include "surrtftr.intfb.h"
	233	#include "surrtpk.intfb.h"
	234	#include "surrtrf.intfb.h"
	235	#include "susat.intfb.h"
	236	#include "suswn.intfb.h"
	237	#include "susrtaer.intfb.h"
	238	#include "srtm_init.intfb.h"
	239	#include "susrtcop.intfb.h"
	240	#include "su_aerw.intfb.h"
	241	#include "su_uvrad.intfb.h"
	242	#include "su_mcica.intfb.h"
	243
[5154]	244	! ----------------------------------------------------------------
[5133]	245
	246	#include "naerad.h"
	247	#include "namrgri.h"
[5154]	248	!MPL/IM 20160915 on prend GES de phylmd
[5133]	249
[5154]	250	!* 1. INITIALIZE NEUROFLUX LONGWAVE RADIATION
	251	! ---------------------------------------
[5133]	252
[5154]	253	IF (LHOOK) CALL DR_HOOK('SUECRAD', 0, ZHOOK_HANDLE)
	254	!CALL GSTATS(1818,0) MPL 2.12.08
	255	!IF (LERADN2) THEN
	256	! CALL SULWNEUR(KLEV)
	257	!ENDIF
[5133]	258
[5154]	259	!* 2. SET DEFAULT VALUES.
	260	! -------------------
[5133]	261
[5154]	262	!* 2.1 PRESET INDICES IN YOERAD
	263	! --------------------------
[5133]	264
[5154]	265	LERAD1H = .FALSE.
	266	NLNGR1H = 6
[5133]	267
[5154]	268	LERADHS = .TRUE.
	269	LONEWSW = .TRUE.
	270	LECSRAD = .FALSE.
[5133]	271
[5154]	272	!LE4ALB=.FALSE.
	273	!This is read from SU0PHY in NAEPHY and put in YOEPHY
[5133]	274
[5154]	275	!- default setting of cloud optical properties
	276	! liquid water cloud 0: Fouquart (SW), Smith-Shi (LW)
	277	! 1: Slingo (SW), Savijarvi (LW)
	278	! 2: Slingo (SW), Lindner-Li (LW)
	279	! ice water cloud 0: Ebert-Curry (SW), Smith-Shi (LW)
	280	! 1: Ebert-Curry (SW), Ebert-Curry (LW)
	281	! 2: Fu-Liou'93 (SW), Fu-Liou'93 (LW)
	282	! 3: Fu'96 (SW), Fu et al'98 (LW)
	283	NLIQOPT = 2 ! before 3?R1 default=0 2
	284	NICEOPT = 3 ! before 3?R1 default=1 3
[5133]	285
[5154]	286	!- default setting of cloud effective radius/diameter
	287	! liquid water cloud 0: f(P) 10 to 45
	288	! 1: 13: ocean; 10: land
	289	! 2: Martin et al. CCN 50 over ocean, 900 over land
	290	! ice water cloud 0: 40 microns
	291	! 1: f(T) 40 to 130 microns
	292	! 2: f(T) 30 to 60
	293	! 3: f(T,IWC) Sun'01: 22.5 to 175 microns
	294	! conversion factor between effective radius and particle size for ice
[5158]	295	NRADIP = 3 ! before 3?R1 default=2 3
	296	NRADLP = 2 ! before 3?R1 default=2 2
[5160]	297	PRINT *, 'SUECRAD: NRADLP, NRADIP=', NRADLP, NRADIP
[5154]	298	RRe2De = 0.64952_JPRB ! before 3?R1 default=0.5_JPRB
[5133]	299
[5154]	300	!- RRTM as LW scheme
	301	LRRTM = .FALSE.
	302	LECMWF = .FALSE.
	303	IF (iflag_rrtm.EQ.1) THEN
	304	LRRTM = .TRUE.
	305	LECMWF = .TRUE.
	306	! LRRTM = .FALSE. ! Utiliser pour faire tourner le "vieux" rayonnement
	307	! LECMWF = .FALSE.
	308	ENDIF
[5133]	309
[5154]	310	!LRRTM = .FALSE.
[5133]	311
[5154]	312	!- SRTM as SW scheme
	313	!!!!! A REVOIR (MPL) verifier signification de LSRTM
	314	LSRTM = .FALSE. ! before 3?R1 default was .FALSE. true
[5133]	315
[5154]	316	! -- McICA treatment of cloud-radiation interactions
	317	! - 1 is maximum-random, 2 is generalized cloud overlap (before 31R1 default=0 no McICA)
	318	NMcICA = 2 ! 2 for generalized overlap
[5133]	319
[5154]	320	!- Inhomogeneity factors in LW and SW (0=F, 1=0.7 in both, 2=Barker's, 3=Cairns)
	321	NINHOM = 0 ! before 3?R1 default=1
	322	NLAYINH = 0
	323	RLWINHF = 1.0_JPRB ! before 3?R1 default=0.7
	324	RSWINHF = 1.0_JPRB ! before 3?R1 default=0.7
	325	!- Diffusivity correction a la Savijarvi
	326	LDIFFC = .FALSE. ! before 31R1 default=.FALSE.
[5133]	327
[5154]	328	!- history of volcanic aerosols
	329	LHVOLCA = .FALSE.
	330	!- monthly climatol. of tropospheric aerosols from Tegen et al. (1997)
	331	LNEWAER = .TRUE.
	332	!!! cpl LNOTROAER=.FALSE.
	333	LNOTROAER = .TRUE.
	334	NPERTAER = 0
[5133]	335
[5154]	336	!- New Rayleigh formulation
	337	LRAYL = .TRUE.
[5133]	338
[5154]	339	!- Number concentration of aerosols if specified
	340	LCCNL = .TRUE. ! before 3?R1 default=.FALSE. true
	341	LCCNO = .TRUE. ! before 3?R1 default=.FALSE. true
	342	RCCNLND = 900._JPRB ! before 3?R1 default=900. now irrelevant
	343	RCCNSEA = 50._JPRB ! before 3?R1 default=50. now irrelevant
[5133]	344
[5154]	345	!- interaction radiation / prognostic O3 off by default
	346	LEPO3RA = .FALSE.
[5160]	347	PRINT *, 'SUECRAD-0'
[5154]	348	IF (.NOT.YO3%LGP) THEN
	349	LEPO3RA = .FALSE.
	350	ENDIF
	351	RPERTOZ = 0._JPRB
	352	NPERTOZ = 0
[5133]	353
[5154]	354	!NAER: CONFIGURATION INDEX FOR AEROSOLS
	355	!!!!! A REVOIR (MPL) a mettre dans un fichier .def
	356	NAER = 1
	357	NMODE = 0
	358	NOZOCL = 1
	359	NRADFR = -3
	360	IF (NSMAX >= 511) NRADFR = -1
	361	NRADPFR = 0
	362	NRADPLA = 15
[5133]	363
[5154]	364	! -- UV diagnostic of surface fluxes over the 280-400 nm interval
	365	! with up-to 24 values (5 nm wide spectral intervals)
	366	LUVPROC = .FALSE.
	367	LUVTDEP = .TRUE.
	368	LUVDBG = .FALSE.
	369	NRADUV = -3
	370	NUVTIM = 0
	371	NUV = 24
	372	RMUZUV = 1.E-01_JPRB
	373	DO JUV = 1, NUV
	374	RUVLAM(JUV) = 280._JPRB + (JUV - 1) * 5._JPRB
	375	ENDDO
[5133]	376
[5154]	377	!- radiation interpolation (George M's grid on by default)
	378	LLDEBUG = .TRUE.
	379	LEDBUG = .FALSE.
	380	NRADINT = 3
	381	NRADRES = 0
[5133]	382
[5154]	383	NRINT = 4
[5133]	384
[5154]	385	LRADLB = .TRUE.
	386	CRTABLEDIR = './'
	387	CRTABLEFIL = 'not set'
	388	LRADONDEM = .TRUE.
	389	!GM Temporary as per trans/external/setup_trans.F90
	390	LLINEAR_GRID = NSMAX > (NDLON + 3) / 3
	391	IF(LLDEBUG)THEN
	392	WRITE(NULOUT, '("SUECRAD: NSMAX=",I6)')NSMAX
	393	WRITE(NULOUT, '("SUECRAD: NDLON=",I6)')NDLON
	394	WRITE(NULOUT, '("SUECRAD: LLINEAR_GRID=",L5)')LLINEAR_GRID
	395	ENDIF
[5133]	396
[5154]	397	NUAER = 24
	398	NTRAER = 15
	399	! 1: max-random, 2: max, 3: random (5,6,7,8 pour meso-NH)
	400	! le CASE qui suit car les conventions sont differentes dans ARP et LMDZ (MPL 20100415)
	401	SELECT CASE (overlap)
[5133]	402	CASE (:1)
[5154]	403	NOVLP = 2
[5133]	404	CASE (2)
[5154]	405	NOVLP = 3
[5133]	406	CASE (3:)
[5154]	407	NOVLP = 1
[5133]	408	END SELECT
[5160]	409	PRINT *, 'SUECRAD: NOVLP=', NOVLP
[5154]	410	NLW = 16
	411	NTSW = 14
	412	!NSW = 6 !!!!! Maintenant dans config.def (MPL 20140213)
	413	NSWNL = 6
	414	NSWTL = 2
	415	NCSRADF = 1
	416	IF(NSMAX >= 106) THEN
	417	NRPROMA = 80
	418	ELSEIF(NSMAX == 63) THEN
	419	NRPROMA = 48
	420	ELSE
	421	NRPROMA = 64
	422	ENDIF
[5133]	423
[5154]	424	!* 2.3 SET SECURITY PARAMETERS
	425	! -----------------------
[5133]	426
[5154]	427	REPSC = 1.E-04_JPRB
	428	REPSCO = 1.E-12_JPRB
	429	REPSCQ = 1.E-12_JPRB
	430	REPSCT = 1.E-12_JPRB
	431	REPSCW = 1.E-12_JPRB
	432	REPLOG = 1.E-12_JPRB
[5133]	433
	434
[5154]	435	!* 2.4 BACKGROUND GAS CONCENTRATIONS (IPCC/SACC, 1990)
	436	! -----------------------------------------------
[5133]	437
[5154]	438	LECO2VAR = .FALSE.
	439	LHGHG = .FALSE.
	440	NHINCSOL = 0
	441	NSCEN = 1
	442	RSOLINC = RI0
[5133]	443
[5154]	444	! Valeurs d origine MPL 18052010
	445	!RCCO2 = 353.E-06_JPRB
	446	!RCCH4 = 1.72E-06_JPRB
	447	!RCN2O = 310.E-09_JPRB
	448	!RCCFC11 = 280.E-12_JPRB
	449	!RCCFC12 = 484.E-12_JPRB
[5133]	450
[5154]	451	! Valeurs LMDZ (physiq.def) MPL 18052010
	452	!RCCO2 = 348.E-06_JPRB
	453	!RCCH4 = 1.65E-06_JPRB
	454	!RCN2O = 306.E-09_JPRB
	455	!RCCFC11 = 280.E-12_JPRB
	456	!RCCFC12 = 484.E-12_JPRB
[5133]	457
[5154]	458	!MPL/IM 20160915 on prend GES de phylmd
	459	RCCO2 = CO2_ppm * 1.0e-06
	460	RCCH4 = CH4_ppb * 1.0e-09
	461	RCN2O = N2O_ppb * 1.0e-09
	462	RCCFC11 = CFC11_ppt * 1.0e-12
	463	RCCFC12 = CFC12_ppt * 1.0e-12
[5160]	464	!PRINT *,'LMDZSUECRAD-1 RCCO2=',RCCO2
	465	!PRINT *,'LMDZSUECRAD-1 RCCH4=',RCCH4
	466	!PRINT *,'LMDZSUECRAD-1 RCN2O=',RCN2O
	467	!PRINT *,'LMDZSUECRAD-1 RCCFC11=',RCCFC11
	468	!PRINT *,'LMDZSUECRAD-1 RCCFC12=',RCCFC12
[5154]	469	! ------------------------------------------------------------------
[5133]	470
[5154]	471	!* 3. READ VALUES OF RADIATION CONFIGURATION
	472	! --------------------------------------
[5133]	473
[5154]	474	!CALL POSNAM(NULNAM,'NAERAD')
	475	!READ (NULNAM,NAERAD)
[5160]	476	PRINT *, 'SUECRAD-2'
[5133]	477
[5154]	478	!CALL POSNAM(NULNAM,'NAEAER')
	479	!READ (NULNAM,NAEAER)
[5133]	480
[5154]	481	!IF (NTYPAER(9) /= 0) THEN
	482	! RGEMUV=(RLATVOL+90._JPRB)*RPI/180._JPRB
	483	! RGELAV=RLONVOL*RPI/180._JPRB
	484	! RCLONV=COS(RGELAV)
	485	! RSLONV=SIN(RGELAV)
	486	! DO J=1,NGPTOT-1
	487	! IF (RGELAV > GELAM(J) .AND. RGELAV <= GELAM(J+1) .AND. &
	488	! & RGEMUV < RMU(JL) .AND. RGEMUV >= RMU(JL+1) ) THEN
	489	! RDGMUV=ABS( RMU(J+1) - RMU(J))
	490	! RDGLAV=ABS( GELAM(J+1)-GELAM(J) )
	491	! RDSLONV=ABS( SIN(GELAM(JL+1))-SIN(GELAM(JL)) )
	492	! RDCLONV=ABS( COS(GELAM(JL+1))-COS(GELAM(JL)) )
	493	! END IF
	494	! END DO
	495	!END IF
[5133]	496
[5154]	497	!- reset some parameters if SW6 is used (revert to pre-CY3?R1 operational configuration)
	498	IF (.NOT.LSRTM) THEN
	499	NMcICA = 0
	500	LCCNL = .FALSE.
	501	LCCNO = .FALSE.
	502	LDIFFC = .FALSE.
	503	NICEOPT = 1
	504	NLIQOPT = 0
	505	NRADIP = 4
	506	NRADLP = 3
	507	RRe2De = 0.5_JPRB
	508	NINHOM = 1
	509	RLWINHF = 0.7_JPRB
	510	RSWINHF = 0.7_JPRB
	511	ENDIF
[5160]	512	PRINT *, 'SUECRAD-3'
[5133]	513
[5154]	514	!- for McICA computations, make sure these parameters are as follows ...
	515	IF (NMCICA /= 0) THEN
	516	NINHOM = 0
	517	RLWINHF = 1.0_JPRB
	518	RSWINHF = 1.0_JPRB
	519	!-- read the XCW values for Raisanen-Cole-Barker cloud generator
	520	CALL SU_McICA
	521	ENDIF
[5160]	522	PRINT *, 'SUECRAD-4'
[5133]	523
[5154]	524	IF(LLDEBUG)THEN
	525	WRITE(NULOUT, '("SUECRAD: NRADINT=",I2)')NRADINT
	526	WRITE(NULOUT, '("SUECRAD: NRADRES=",I4)')NRADRES
	527	ENDIF
[5133]	528
[5154]	529	! DETERMINE WHETHER NRPROMA IS NEGATIVE AND SET LOPTRPROMA
[5133]	530
[5154]	531	LOPTRPROMA = NRPROMA > 0
	532	NRPROMA = ABS(NRPROMA)
[5133]	533
[5154]	534	IF(NRADINT > 0 .AND. NRADRES == NSMAX)THEN
	535	WRITE(NULOUT, '("SUECRAD: NRADINT > 0 .AND. NRADRES = NSMAX, NRADINT RESET TO 0")')
	536	NRADINT = 0
	537	ENDIF
[5133]	538
[5154]	539	IF(NRADINT > 0 .AND. LRAYFM .AND. NAER /= 0 .AND. .NOT.LHVOLCA)THEN
	540	! This combination is not supported as aerosol data would be
	541	! required to be interpolated (see radintg)
	542	WRITE(NULOUT, '("SUECRAD: NRADINT>0, LRAYFM=T NAER /= 0 .AND. LHVOLCA=F,",&
	543	& " NRADRES RESET TO NSMAX (NO INTERPOLATION)")')
	544	NRADRES = NSMAX
	545	ENDIF
	546	!CALL GSTATS(1818,1) MPL 2.12.08
[5133]	547
[5154]	548	100 CONTINUE
[5133]	549
[5154]	550	IF(LERADI)THEN ! START OF LERADI BLOCK
[5133]	551
[5154]	552	IF(NRADINT == -1)THEN
[5133]	553
[5154]	554	! INITIALISE DATA STRUCTURES REQUIRED FOR RADIATION INTERPOLATION
[5133]	555
[5154]	556	LODBGRADI = .FALSE.
	557	CALL SUECRADI
[5133]	558
[5154]	559	! INITIALISE DATA STRUCTURES REQUIRED FOR RADIATION COURSE GRID
	560	! LOAD BALANCING
[5133]	561
[5154]	562	LODBGRADL = .FALSE.
	563	! CALL SUECRADL ! MPL 1.12.08
	564	CALL ABOR1('JUSTE APRES CALL SUECRADL COMMENTE')
[5133]	565
[5154]	566	ELSEIF(NRADINT == 0)THEN
[5133]	567
[5154]	568	IF(NRADRES /= NSMAX)THEN
	569	WRITE(NULOUT, '("SUECRAD: NRADINT=0 REQUESTED, NRADRES RESET TO NSMAX")')
	570	NRADRES = NSMAX
	571	ENDIF
	572	RADGRID%NGPTOT = NGPTOT
[5133]	573
[5154]	574	NARIB1 = 0
	575	NAROB1 = 0
[5133]	576
[5154]	577	ELSEIF(NRADINT >=1 .AND. NRADINT <= 3)THEN
[5133]	578
[5154]	579	NARIB1 = 0
	580	NAROB1 = 0
[5133]	581
[5154]	582	! set the default radiation grid resolution for the current model resolution
	583	! if not already specified
	584	IF(NRADRES == 0)THEN
	585	IF(LLINEAR_GRID)THEN ! RATIO OF GRID-POINTS (MODEL/RAD)
	586	IF(NSMAX == 63)THEN
	587	NRADRES = 21 ! 3.62
	588	LLINEAR_GRID = .FALSE.
	589	ENDIF
	590	IF(NSMAX == 95) NRADRES = 95 ! 1.00
	591	IF(NSMAX == 159) NRADRES = 63 ! 5.84
	592	IF(NSMAX == 255) NRADRES = 95 ! 6.69
	593	IF(NSMAX == 319) NRADRES = 159 ! 3.87
	594	IF(NSMAX == 399) NRADRES = 159 ! 5.99
	595	IF(NSMAX == 511) NRADRES = 255 ! 3.92
	596	IF(NSMAX == 639) NRADRES = 319 ! 3.92
	597	IF(NSMAX == 799) NRADRES = 399 ! 3.94
	598	IF(NSMAX == 1023) NRADRES = 511 ! 3.94
	599	IF(NSMAX == 1279) NRADRES = 639 !
	600	IF(NSMAX == 2047) NRADRES = 1023 !
	601	ELSE ! NOT LINEAR GRID
	602	IF(NSMAX == 21) NRADRES = 21 ! 1.00
	603	IF(NSMAX == 42) NRADRES = 21 ! 3.62
	604	IF(NSMAX == 63) NRADRES = 42 ! 2.17
	605	IF(NSMAX == 106) NRADRES = 63 ! 2.69
	606	IF(NSMAX == 170) NRADRES = 63 ! 6.69
	607	IF(NSMAX == 213) NRADRES = 106 ! 3.87
	608	IF(NSMAX == 266) NRADRES = 106 ! 5.99
	609	IF(NSMAX == 341) NRADRES = 170 ! 3.92
	610	IF(NSMAX == 426) NRADRES = 213 ! 3.92
	611	IF(NSMAX == 533) NRADRES = 266 ! 3.94
	612	IF(NSMAX == 682) NRADRES = 341 ! 3.94
[5133]	613	ENDIF
	614	ENDIF
[5160]	615	PRINT *, 'SUECRAD-5'
[5133]	616
[5154]	617	! test if radiation grid resolution has been set
	618	IF(NRADRES == 0)THEN
	619	WRITE(NULOUT, '("SUECRAD: NRADRES NOT SET OR DEFAULT FOUND,NSMAX=",I4)')NSMAX
	620	CALL ABOR1('SUECRAD: NRADRES NOT SET OR DEFAULT FOUND')
	621	ENDIF
[5133]	622
[5154]	623	! test if no interpolation is required
	624	IF(NRADINT > 0 .AND. NRADRES == NSMAX)THEN
	625	WRITE(NULOUT, '("SUECRAD: NRADINT > 0 .AND. NRADRES = NSMAX, NRADINT RESET TO 0")')
	626	NRADINT = 0
	627	GOTO 100
	628	ENDIF
[5133]	629
[5154]	630	! CALL GSTATS(1818,0) MPL 2.12.08
	631	IF(CRTABLEFIL == 'not set')THEN
	632	IF(LLINEAR_GRID)THEN
	633	IF(NRADRES < 1000)THEN
	634	WRITE(CRTABLEFIL, '("rtablel_2",I3.3)')NRADRES
	635	ELSE
	636	WRITE(CRTABLEFIL, '("rtablel_2",I4.4)')NRADRES
	637	ENDIF
[5133]	638	ELSE
[5154]	639	IF(NRADRES < 1000)THEN
	640	WRITE(CRTABLEFIL, '("rtable_2" ,I3.3)')NRADRES
	641	ELSE
	642	WRITE(CRTABLEFIL, '("rtable_2" ,I4.4)')NRADRES
	643	ENDIF
[5133]	644	ENDIF
	645	ENDIF
[5154]	646	! CALL GSTATS(1818,1) MPL 2.12.08
[5133]	647
[5154]	648	RADGRID%NSMAX = NRADRES
[5133]	649
[5154]	650	IF(MYPROC == JPIOMASTER)THEN
	651	IDIR = LEN_TRIM(CRTABLEDIR)
	652	IFIL = LEN_TRIM(CRTABLEFIL)
	653	CLFN = CRTABLEDIR(1:IDIR) // CRTABLEFIL(1:IFIL)
	654	! Ce qui concerne NULRAD commente par MPL le 15.04.09
	655	! OPEN(NULRAD,FILE=CLFN,ACTION="READ",ERR=999)
	656	! GOTO 1000
	657	! 999 CONTINUE
	658	! WRITE(NULOUT,'("SUECRAD: UNABLE TO OPEN FILE ",A)')CLFN
	659	! CALL ABOR1('SUECRAD: UNABLE TO OPEN RADIATION GRID RTABLE FILE')
	660	! 1000 CONTINUE
	661	NRGRI(:) = 0
	662	! Ce qui concerne NAMRGRI commente par MPL le 15.04.09
	663	! CALL POSNAM(NULRAD,'NAMRGRI')
	664	! READ (NULRAD,NAMRGRI)
	665	IDGL = 1
	666	DO WHILE(NRGRI(IDGL)>0)
	667	IF(LLDEBUG)THEN
	668	WRITE(NULOUT, '("SUECRAD: NRGRI(",I4,")=",I4)')IDGL, NRGRI(IDGL)
	669	ENDIF
	670	IDGL = IDGL + 1
	671	ENDDO
	672	IDGL = IDGL - 1
	673	RADGRID%NDGLG = IDGL
	674	IF(LLDEBUG)THEN
	675	WRITE(NULOUT, '("SUECRAD: RADGRID%NDGLG=",I4)')RADGRID%NDGLG
[5133]	676	ENDIF
[5154]	677	! CLOSE(NULRAD)
[5133]	678	ENDIF
[5154]	679	! CALL GSTATS(667,0) MPL 2.12.08
	680	IF(NPROC > 1)THEN
	681	stop 'Pas pret pour proc > 1'
	682	! CALL MPL_BROADCAST (RADGRID%NDGLG,MTAGRAD,JPIOMASTER,CDSTRING='SUECRAD:')
	683	ENDIF
	684	ALLOCATE(RADGRID%NRGRI(RADGRID%NDGLG))
	685	IF(MYPROC == JPIOMASTER)THEN
	686	RADGRID%NRGRI(1:RADGRID%NDGLG) = NRGRI(1:RADGRID%NDGLG)
	687	ENDIF
	688	IF(NPROC > 1)THEN
	689	stop 'Pas pret pour proc > 1'
	690	! CALL MPL_BROADCAST (RADGRID%NRGRI(1:RADGRID%NDGLG),MTAGRAD,JPIOMASTER,CDSTRING='SUECRAD:')
	691	ENDIF
	692	! CALL GSTATS(667,1) MPL 2.12.08
[5133]	693
[5154]	694	! CALL GSTATS(1818,0) MPL 2.12.08
	695	IF (NRADINT == 1)THEN
	696	WRITE(NULOUT, '("SUECRAD: INTERPOLATION METHOD - SPECTRAL TRANSFORM")')
	697	RADGRID%NDGSUR = 0
	698	NRIWIDEN = 0
	699	NRIWIDES = 0
	700	NRIWIDEW = 0
	701	NRIWIDEE = 0
	702	NROWIDEN = 0
	703	NROWIDES = 0
	704	NROWIDEW = 0
	705	NROWIDEE = 0
	706	ELSEIF(NRADINT == 2)THEN
	707	WRITE(NULOUT, '("SUECRAD: INTERPOLATION METHOD - 4 POINT")')
	708	RADGRID%NDGSUR = 2
	709	ELSEIF(NRADINT == 3)THEN
	710	WRITE(NULOUT, '("SUECRAD: INTERPOLATION METHOD - 12 POINT")')
	711	RADGRID%NDGSUR = 2
	712	ENDIF
	713	WRITE(NULOUT, '("SUECRAD: RADGRID%NDGSUR =",I8)')RADGRID%NDGSUR
[5133]	714
[5154]	715	RADGRID%NDGSAG = 1 - RADGRID%NDGSUR
	716	RADGRID%NDGENG = RADGRID%NDGLG + RADGRID%NDGSUR
	717	RADGRID%NDLON = RADGRID%NRGRI(RADGRID%NDGLG / 2)
	718	WRITE(NULOUT, '("SUECRAD: RADGRID%NDGSAG =",I8)')RADGRID%NDGSAG
	719	WRITE(NULOUT, '("SUECRAD: RADGRID%NDGENG =",I8)')RADGRID%NDGENG
	720	WRITE(NULOUT, '("SUECRAD: RADGRID%NDGLG =",I8)')RADGRID%NDGLG
	721	WRITE(NULOUT, '("SUECRAD: RADGRID%NDLON =",I8)')RADGRID%NDLON
	722	CALL FLUSH(NULOUT)
[5133]	723
[5154]	724	ALLOCATE(RADGRID%NLOENG(RADGRID%NDGSAG:RADGRID%NDGENG))
	725	RADGRID%NLOENG(1:RADGRID%NDGLG) = RADGRID%NRGRI(1:RADGRID%NDGLG)
	726	IF(RADGRID%NDGSUR >= 1)THEN
	727	DO JGLSUR = 1, RADGRID%NDGSUR
	728	RADGRID%NLOENG(1 - JGLSUR) = RADGRID%NLOENG(JGLSUR)
	729	ENDDO
	730	DO JGLSUR = 1, RADGRID%NDGSUR
	731	RADGRID%NLOENG(RADGRID%NDGLG + JGLSUR) = RADGRID%NLOENG(RADGRID%NDGLG + 1 - JGLSUR)
	732	ENDDO
	733	ENDIF
	734	! CALL GSTATS(1818,1) MPL 2.12.08
[5133]	735
[5154]	736	! Setup the transform package for the radiation grid
	737	CALL SETUP_TRANS (KSMAX = RADGRID%NSMAX, &
	738	& KDGL = RADGRID%NDGLG, &
	739	& KLOEN = RADGRID%NLOENG(1:RADGRID%NDGLG), &
	740	& LDLINEAR_GRID = LLINEAR_GRID, &
	741	& LDSPLIT = LSPLIT, &
	742	& KAPSETS = NAPSETS, &
	743	& KRESOL = RADGRID%NRESOL_ID)
[5133]	744
[5154]	745	ALLOCATE(RADGRID%NSTA(RADGRID%NDGSAG:RADGRID%NDGENG + N_REGIONS_NS - 1, N_REGIONS_EW))
	746	ALLOCATE(RADGRID%NONL(RADGRID%NDGSAG:RADGRID%NDGENG + N_REGIONS_NS - 1, N_REGIONS_EW))
	747	ALLOCATE(RADGRID%NPTRFRSTLAT(N_REGIONS_NS))
	748	ALLOCATE(RADGRID%NFRSTLAT(N_REGIONS_NS))
	749	ALLOCATE(RADGRID%NLSTLAT(N_REGIONS_NS))
	750	ALLOCATE(RADGRID%RMU(RADGRID%NDGSAG:RADGRID%NDGENG))
	751	ALLOCATE(RADGRID%RSQM2(RADGRID%NDGSAG:RADGRID%NDGENG))
	752	ALLOCATE(RADGRID%RLATIG(RADGRID%NDGSAG:RADGRID%NDGENG))
[5133]	753
[5154]	754	! Interrogate the transform package for the radiation grid
	755	! CALL GSTATS(1818,0) MPL 2.12.08
	756	CALL TRANS_INQ (KRESOL = RADGRID%NRESOL_ID, &
	757	& KSPEC2 = RADGRID%NSPEC2, &
	758	& KNUMP = RADGRID%NUMP, &
	759	& KGPTOT = RADGRID%NGPTOT, &
	760	& KGPTOTG = RADGRID%NGPTOTG, &
	761	& KGPTOTMX = RADGRID%NGPTOTMX, &
	762	& KPTRFRSTLAT = RADGRID%NPTRFRSTLAT, &
	763	& KFRSTLAT = RADGRID%NFRSTLAT, &
	764	& KLSTLAT = RADGRID%NLSTLAT, &
	765	& KFRSTLOFF = RADGRID%NFRSTLOFF, &
	766	& KSTA = RADGRID%NSTA(1:RADGRID%NDGLG + N_REGIONS_NS - 1, :), &
	767	& KONL = RADGRID%NONL(1:RADGRID%NDGLG + N_REGIONS_NS - 1, :), &
	768	& KPTRFLOFF = RADGRID%NPTRFLOFF, &
	769	& PMU = RADGRID%RMU(1:))
[5133]	770
[5154]	771	IF(NRADINT == 2 .OR. NRADINT == 3)THEN
	772	DO JGL = 1, RADGRID%NDGLG
	773	RADGRID%RSQM2(JGL) = SQRT(1.0_JPRB - RADGRID%RMU(JGL) * RADGRID%RMU(JGL))
	774	RADGRID%RLATIG(JGL) = ASIN(RADGRID%RMU(JGL))
	775	! WRITE(NULOUT,'("SUECRAD: JGL=",I6," RADGRID%RLATIG=",F10.3)')&
	776	! & JGL,RADGRID%RLATIG(JGL)
[5133]	777	ENDDO
[5154]	778	IF(RADGRID%NDGSUR >= 1)THEN
	779	DO JGLSUR = 1, RADGRID%NDGSUR
	780	RADGRID%RMU(1 - JGLSUR) = RADGRID%RMU(JGLSUR)
	781	RADGRID%RSQM2(1 - JGLSUR) = RADGRID%RSQM2(JGLSUR)
	782	RADGRID%RLATIG(1 - JGLSUR) = RPI - RADGRID%RLATIG(JGLSUR)
	783	ENDDO
	784	DO JGLSUR = 1, RADGRID%NDGSUR
	785	RADGRID%RMU(RADGRID%NDGLG + JGLSUR) = RADGRID%RMU(RADGRID%NDGLG + 1 - JGLSUR)
	786	RADGRID%RSQM2(RADGRID%NDGLG + JGLSUR) = RADGRID%RSQM2(RADGRID%NDGLG + 1 - JGLSUR)
	787	RADGRID%RLATIG(RADGRID%NDGLG + JGLSUR) = -RPI - RADGRID%RLATIG(RADGRID%NDGLG + 1 - JGLSUR)
	788	ENDDO
	789	ENDIF
[5133]	790	ENDIF
	791
[5154]	792	RADGRID%NDGSAL = 1
	793	RADGRID%NDGENL = RADGRID%NLSTLAT(MY_REGION_NS) - RADGRID%NFRSTLOFF
	794	RADGRID%NDSUR1 = 3 - MOD(RADGRID%NDLON, 2)
	795	IDLSUR = MAX(RADGRID%NDLON, 2 * RADGRID%NSMAX + 1)
	796	RADGRID%NDLSUR = IDLSUR + RADGRID%NDSUR1
	797	RADGRID%MYFRSTACTLAT = RADGRID%NFRSTLAT(MY_REGION_NS)
	798	RADGRID%MYLSTACTLAT = RADGRID%NLSTLAT(MY_REGION_NS)
[5133]	799
[5154]	800	WRITE(NULOUT, '("SUECRAD: RADGRID%NRESOL_ID =",I8)')RADGRID%NRESOL_ID
	801	WRITE(NULOUT, '("SUECRAD: RADGRID%NSMAX =",I8)')RADGRID%NSMAX
	802	WRITE(NULOUT, '("SUECRAD: RADGRID%NSPEC2 =",I8)')RADGRID%NSPEC2
	803	WRITE(NULOUT, '("SUECRAD: RADGRID%NGPTOT =",I8)')RADGRID%NGPTOT
	804	WRITE(NULOUT, '("SUECRAD: RADGRID%NGPTOTG =",I8)')RADGRID%NGPTOTG
	805	WRITE(NULOUT, '("SUECRAD: RADGRID%NDGSAL =",I8)')RADGRID%NDGSAL
	806	WRITE(NULOUT, '("SUECRAD: RADGRID%NDGENL =",I8)')RADGRID%NDGENL
	807	WRITE(NULOUT, '("SUECRAD: RADGRID%NDSUR1 =",I8)')RADGRID%NDSUR1
	808	WRITE(NULOUT, '("SUECRAD: RADGRID%NDLSUR =",I8)')RADGRID%NDLSUR
	809	WRITE(NULOUT, '("SUECRAD: RADGRID%MYFRSTACTLAT =",I8)')RADGRID%MYFRSTACTLAT
	810	WRITE(NULOUT, '("SUECRAD: RADGRID%MYLSTACTLAT =",I8)')RADGRID%MYLSTACTLAT
	811	CALL FLUSH(NULOUT)
[5133]	812
[5154]	813	ALLOCATE(RADGRID%NASM0(0:RADGRID%NSPEC2))
	814	ALLOCATE(RADGRID%MYMS(RADGRID%NUMP))
	815	CALL TRANS_INQ (KRESOL = RADGRID%NRESOL_ID, &
	816	& KASM0 = RADGRID%NASM0, &
	817	& KMYMS = RADGRID%MYMS)
[5133]	818
[5154]	819	ALLOCATE(RADGRID%GELAM(RADGRID%NGPTOT))
	820	ALLOCATE(RADGRID%GELAT(RADGRID%NGPTOT))
	821	ALLOCATE(RADGRID%GESLO(RADGRID%NGPTOT))
	822	ALLOCATE(RADGRID%GECLO(RADGRID%NGPTOT))
	823	ALLOCATE(RADGRID%GEMU (RADGRID%NGPTOT))
[5133]	824
[5154]	825	IOFF = 0
	826	ILAT = RADGRID%NPTRFLOFF
	827	DO JGLAT = RADGRID%NFRSTLAT(MY_REGION_NS), &
	828	& RADGRID%NLSTLAT(MY_REGION_NS)
	829	ZGEMU = RADGRID%RMU(JGLAT)
	830	ILAT = ILAT + 1
	831	ISTLON = RADGRID%NSTA(ILAT, MY_REGION_EW)
	832	IENDLON = ISTLON - 1 + RADGRID%NONL(ILAT, MY_REGION_EW)
[5133]	833
[5154]	834	DO JLON = ISTLON, IENDLON
	835	ZLON = REAL(JLON - 1, JPRB) * 2.0_JPRB * RPI &
	836	& / REAL(RADGRID%NLOENG(JGLAT), JPRB)
	837	IOFF = IOFF + 1
	838	RADGRID%GELAM(IOFF) = ZLON
	839	RADGRID%GELAT(IOFF) = ASIN(ZGEMU)
	840	RADGRID%GESLO(IOFF) = SIN(ZLON)
	841	RADGRID%GECLO(IOFF) = COS(ZLON)
	842	RADGRID%GEMU (IOFF) = ZGEMU
	843	ENDDO
[5133]	844	ENDDO
	845
[5154]	846	IF(NRADINT == 2 .OR. NRADINT == 3)THEN
[5133]	847
[5154]	848	! For grid point interpolations we need to calculate the halo size
	849	! required by each processor
[5133]	850
[5154]	851	ALLOCATE(ZLATX(RADGRID%NGPTOTMX))
	852	ALLOCATE(ZLONX(RADGRID%NGPTOTMX))
	853	DO J = 1, RADGRID%NGPTOT
	854	ZLATX(J) = RADGRID%GELAT(J) / RPI * 2.0_JPRB * 90.0
	855	ZLONX(J) = (RADGRID%GELAM(J) - RPI) / RPI * 180.0
	856	ENDDO
	857	ZMINRADLAT = MINVAL(ZLATX(1:RADGRID%NGPTOT))
	858	ZMAXRADLAT = MAXVAL(ZLATX(1:RADGRID%NGPTOT))
	859	ZMINRADLON = MINVAL(ZLONX(1:RADGRID%NGPTOT))
	860	ZMAXRADLON = MAXVAL(ZLONX(1:RADGRID%NGPTOT))
	861	IF(LLDEBUG)THEN
	862	WRITE(NULOUT, '("RADGRID,BEGIN")')
	863	IF(MYPROC /= 1)THEN
	864	stop 'Pas pret pour proc > 1'
	865	! CALL MPL_SEND(RADGRID%NGPTOT,KDEST=NPRCIDS(1),KTAG=1,CDSTRING='SUECRAD.R')
	866	! CALL MPL_SEND(ZLATX(1:RADGRID%NGPTOT),KDEST=NPRCIDS(1),KTAG=2,CDSTRING='SUECRAD.R')
	867	! CALL MPL_SEND(ZLONX(1:RADGRID%NGPTOT),KDEST=NPRCIDS(1),KTAG=3,CDSTRING='SUECRAD.R')
	868	ENDIF
	869	IF(MYPROC == 1)THEN
	870	DO JROC = 1, NPROC
	871	IF(JROC == MYPROC)THEN
	872	DO J = 1, RADGRID%NGPTOT
	873	WRITE(NULOUT, '(F7.2,2X,F7.2,2X,I6)')ZLATX(J), ZLONX(J), MYPROC
	874	ENDDO
	875	ELSE
	876	stop 'Pas pret pour proc > 1'
	877	! CALL MPL_RECV(IGPTOT,KSOURCE=NPRCIDS(JROC),KTAG=1,CDSTRING='SUECRAD.M')
	878	! CALL MPL_RECV(ZLATX(1:IGPTOT),KSOURCE=NPRCIDS(JROC),KTAG=2,CDSTRING='SUECRAD.M')
	879	! CALL MPL_RECV(ZLONX(1:IGPTOT),KSOURCE=NPRCIDS(JROC),KTAG=3,CDSTRING='SUECRAD.M')
	880	DO J = 1, IGPTOT
	881	WRITE(NULOUT, '(F7.2,2X,F7.2,2X,I6)')ZLATX(J), ZLONX(J), JROC
	882	ENDDO
	883	ENDIF
	884	ENDDO
	885	ENDIF
	886	WRITE(NULOUT, '("RADGRID,END")')
[5133]	887	ENDIF
[5154]	888	DEALLOCATE(ZLATX)
	889	DEALLOCATE(ZLONX)
	890
	891	ALLOCATE(ZLATX(NGPTOTMX))
	892	ALLOCATE(ZLONX(NGPTOTMX))
	893	DO J = 1, NGPTOT
	894	ZLATX(J) = GELAT(J) / RPI * 2.0_JPRB * 90.0
	895	ZLONX(J) = (GELAM(J) - RPI) / RPI * 180.0
	896	ENDDO
	897	ZMINMDLLAT = MINVAL(ZLATX(1:NGPTOT))
	898	ZMAXMDLLAT = MAXVAL(ZLATX(1:NGPTOT))
	899	ZMINMDLLON = MINVAL(ZLONX(1:NGPTOT))
	900	ZMAXMDLLON = MAXVAL(ZLONX(1:NGPTOT))
	901	IF(LLDEBUG)THEN
	902	WRITE(NULOUT, '("MODELGRID,BEGIN")')
	903	IF(MYPROC /= 1)THEN
	904	stop 'Pas pret pour proc > 1'
	905	! CALL MPL_SEND(NGPTOT,KDEST=NPRCIDS(1),KTAG=1,CDSTRING='SUECRAD')
	906	! CALL MPL_SEND(ZLATX(1:NGPTOT),KDEST=NPRCIDS(1),KTAG=2,CDSTRING='SUECRAD')
	907	! CALL MPL_SEND(ZLONX(1:NGPTOT),KDEST=NPRCIDS(1),KTAG=3,CDSTRING='SUECRAD')
	908	! CALL MPL_SEND(NGLOBALINDEX(1:NGPTOT),KDEST=NPRCIDS(1),KTAG=4,CDSTRING='SUECRAD')
	909	ENDIF
	910	IF(MYPROC == 1)THEN
	911	DO JROC = 1, NPROC
	912	IF(JROC == MYPROC)THEN
	913	DO J = 1, NGPTOT
	914	WRITE(NULOUT, '(F7.2,2X,F7.2,2X,I6,2X,I12)')ZLATX(J), ZLONX(J), MYPROC, NGLOBALINDEX(J)
	915	ENDDO
	916	ELSE
	917	stop 'Pas pret pour proc > 1'
	918	! CALL MPL_RECV(IGPTOT,KSOURCE=NPRCIDS(JROC),KTAG=1,CDSTRING='SUECRAD')
	919	! CALL MPL_RECV(ZLATX(1:IGPTOT),KSOURCE=NPRCIDS(JROC),KTAG=2,CDSTRING='SUECRAD')
	920	! CALL MPL_RECV(ZLONX(1:IGPTOT),KSOURCE=NPRCIDS(JROC),KTAG=3,CDSTRING='SUECRAD')
	921	ALLOCATE(IGLOBALINDEX(1:IGPTOT))
	922	! CALL MPL_RECV(IGLOBALINDEX(1:IGPTOT),KSOURCE=NPRCIDS(JROC),KTAG=4,CDSTRING='SUECRAD')
	923	DO J = 1, IGPTOT
	924	WRITE(NULOUT, '(F7.2,2X,F7.2,2X,I6,2X,I12)')ZLATX(J), ZLONX(J), JROC, IGLOBALINDEX(J)
	925	ENDDO
	926	DEALLOCATE(IGLOBALINDEX)
	927	ENDIF
	928	ENDDO
	929	ENDIF
	930	WRITE(NULOUT, '("MODELGRID,END")')
[5133]	931	ENDIF
[5154]	932	DEALLOCATE(ZLATX)
	933	DEALLOCATE(ZLONX)
	934
	935	IF(LLDEBUG)THEN
	936	WRITE(NULOUT, '("ZMINRADLAT=",F10.2)')ZMINRADLAT
	937	WRITE(NULOUT, '("ZMINMDLLAT=",F10.2)')ZMINMDLLAT
	938	WRITE(NULOUT, '("ZMAXRADLAT=",F10.2)')ZMAXRADLAT
	939	WRITE(NULOUT, '("ZMAXMDLLAT=",F10.2)')ZMAXMDLLAT
	940	WRITE(NULOUT, '("ZMINRADLON=",F10.2)')ZMINRADLON
	941	WRITE(NULOUT, '("ZMINMDLLON=",F10.2)')ZMINMDLLON
	942	WRITE(NULOUT, '("ZMAXRADLON=",F10.2)')ZMAXRADLON
	943	WRITE(NULOUT, '("ZMAXMDLLON=",F10.2)')ZMAXMDLLON
[5133]	944	ENDIF
[5154]	945
	946	ZLAT = NDGLG / 180.
	947	ILATS_DIFF_C = CEILING(ABS(ZMINRADLAT - ZMINMDLLAT) * ZLAT)
	948	ILATS_DIFF_F = FLOOR (ABS(ZMINRADLAT - ZMINMDLLAT) * ZLAT)
	949	IF(ZMINRADLAT < ZMINMDLLAT)THEN
	950	NRIWIDES = JP_MIN_HALO + ILATS_DIFF_C
	951	ELSE
	952	NRIWIDES = MAX(0, JP_MIN_HALO - ILATS_DIFF_F)
[5133]	953	ENDIF
[5154]	954	ILATS_DIFF_C = CEILING(ABS(ZMAXRADLAT - ZMAXMDLLAT) * ZLAT)
	955	ILATS_DIFF_F = FLOOR (ABS(ZMAXRADLAT - ZMAXMDLLAT) * ZLAT)
	956	IF(ZMAXRADLAT < ZMAXMDLLAT)THEN
	957	NRIWIDEN = MAX(0, JP_MIN_HALO - ILATS_DIFF_F)
	958	ELSE
	959	NRIWIDEN = JP_MIN_HALO + ILATS_DIFF_C
	960	ENDIF
	961	ILATS_DIFF_C = CEILING(ABS(ZMINRADLON - ZMINMDLLON) * ZLAT)
	962	ILATS_DIFF_F = FLOOR (ABS(ZMINRADLON - ZMINMDLLON) * ZLAT)
	963	IF(ZMINRADLON < ZMINMDLLON)THEN
	964	NRIWIDEW = JP_MIN_HALO + ILATS_DIFF_C
	965	ELSE
	966	NRIWIDEW = MAX(0, JP_MIN_HALO - ILATS_DIFF_F)
	967	ENDIF
	968	ILATS_DIFF_C = CEILING(ABS(ZMAXRADLON - ZMAXMDLLON) * ZLAT)
	969	ILATS_DIFF_F = FLOOR (ABS(ZMAXRADLON - ZMAXMDLLON) * ZLAT)
	970	IF(ZMAXRADLON < ZMAXMDLLON)THEN
	971	NRIWIDEE = MAX(0, JP_MIN_HALO - ILATS_DIFF_F)
	972	ELSE
	973	NRIWIDEE = JP_MIN_HALO + ILATS_DIFF_C
	974	ENDIF
[5133]	975
[5154]	976	ZLAT = RADGRID%NDGLG / 180.
	977	ILATS_DIFF_C = CEILING(ABS(ZMINRADLAT - ZMINMDLLAT) * ZLAT)
	978	ILATS_DIFF_F = FLOOR (ABS(ZMINRADLAT - ZMINMDLLAT) * ZLAT)
	979	IF(ZMINMDLLAT < ZMINRADLAT)THEN
	980	NROWIDES = JP_MIN_HALO + ILATS_DIFF_C
	981	ELSE
	982	NROWIDES = MAX(0, JP_MIN_HALO - ILATS_DIFF_F)
	983	ENDIF
	984	ILATS_DIFF_C = CEILING(ABS(ZMAXRADLAT - ZMAXMDLLAT) * ZLAT)
	985	ILATS_DIFF_F = FLOOR (ABS(ZMAXRADLAT - ZMAXMDLLAT) * ZLAT)
	986	IF(ZMAXMDLLAT < ZMAXRADLAT)THEN
	987	NROWIDEN = MAX(0, JP_MIN_HALO - ILATS_DIFF_F)
	988	ELSE
	989	NROWIDEN = JP_MIN_HALO + ILATS_DIFF_C
	990	ENDIF
	991	ILATS_DIFF_C = CEILING(ABS(ZMINRADLON - ZMINMDLLON) * ZLAT)
	992	ILATS_DIFF_F = FLOOR (ABS(ZMINRADLON - ZMINMDLLON) * ZLAT)
	993	IF(ZMINMDLLON < ZMINRADLON)THEN
	994	NROWIDEW = JP_MIN_HALO + ILATS_DIFF_C
	995	ELSE
	996	NROWIDEW = MAX(0, JP_MIN_HALO - ILATS_DIFF_F)
	997	ENDIF
	998	ILATS_DIFF_C = CEILING(ABS(ZMAXRADLON - ZMAXMDLLON) * ZLAT)
	999	ILATS_DIFF_F = FLOOR (ABS(ZMAXRADLON - ZMAXMDLLON) * ZLAT)
	1000	IF(ZMAXMDLLON < ZMAXRADLON)THEN
	1001	NROWIDEE = MAX(0, JP_MIN_HALO - ILATS_DIFF_F)
	1002	ELSE
	1003	NROWIDEE = JP_MIN_HALO + ILATS_DIFF_C
	1004	ENDIF
[5133]	1005
	1006	ENDIF
	1007
[5154]	1008	RADGRID%NDGSAH = MAX(RADGRID%NDGSAG, &
	1009	& RADGRID%NDGSAL + RADGRID%NFRSTLOFF - NROWIDEN) - RADGRID%NFRSTLOFF
	1010	RADGRID%NDGENH = MIN(RADGRID%NDGENG, &
	1011	& RADGRID%NDGENL + RADGRID%NFRSTLOFF + NROWIDES) - RADGRID%NFRSTLOFF
	1012	WRITE(NULOUT, '("SUECRAD: RADGRID%NDGSAH =",I8)')RADGRID%NDGSAH
	1013	WRITE(NULOUT, '("SUECRAD: RADGRID%NDGENH =",I8)')RADGRID%NDGENH
[5133]	1014
[5154]	1015	IF(NRADINT == 2 .OR. NRADINT == 3)THEN
[5133]	1016
[5154]	1017	ILBRLATI = MAX(RADGRID%NDGSAG, &
	1018	& RADGRID%NDGSAL + RADGRID%NFRSTLOFF - NROWIDEN) - RADGRID%NFRSTLOFF
	1019	IUBRLATI = MIN(RADGRID%NDGENG, &
	1020	& RADGRID%NDGENL + RADGRID%NFRSTLOFF + NROWIDES) - RADGRID%NFRSTLOFF
	1021	ALLOCATE(RADGRID%RLATI(ILBRLATI:IUBRLATI))
	1022	ALLOCATE(RADGRID%RIPI0(ILBRLATI:IUBRLATI))
	1023	ALLOCATE(RADGRID%RIPI1(ILBRLATI:IUBRLATI))
	1024	ALLOCATE(RADGRID%RIPI2(ILBRLATI:IUBRLATI))
[5133]	1025
[5154]	1026	DO JGL = ILBRLATI, IUBRLATI
	1027	IGLGLO = JGL + RADGRID%NFRSTLOFF
	1028	IF(IGLGLO >= 0.AND.IGLGLO <= RADGRID%NDGLG) THEN
	1029	ZD1 = RADGRID%RLATIG(IGLGLO - 1) - RADGRID%RLATIG(IGLGLO)
	1030	ZD2 = RADGRID%RLATIG(IGLGLO - 1) - RADGRID%RLATIG(IGLGLO + 1)
	1031	ZD3 = RADGRID%RLATIG(IGLGLO - 1) - RADGRID%RLATIG(IGLGLO + 2)
	1032	ZD4 = RADGRID%RLATIG(IGLGLO) - RADGRID%RLATIG(IGLGLO + 1)
	1033	ZD5 = RADGRID%RLATIG(IGLGLO) - RADGRID%RLATIG(IGLGLO + 2)
	1034	ZD6 = RADGRID%RLATIG(IGLGLO + 1) - RADGRID%RLATIG(IGLGLO + 2)
	1035	RADGRID%RIPI0(JGL) = -1.0_JPRB / (ZD1 * ZD4 * ZD5)
	1036	RADGRID%RIPI1(JGL) = 1.0_JPRB / (ZD2 * ZD4 * ZD6)
	1037	RADGRID%RIPI2(JGL) = -1.0_JPRB / (ZD3 * ZD5 * ZD6)
	1038	ENDIF
	1039	RADGRID%RLATI(JGL) = RADGRID%RLATIG(IGLGLO)
	1040	ENDDO
	1041
	1042	IF(NPROC > 1)THEN
	1043	IRIRPTSUR = NGPTOTG
	1044	IRISPTSUR = 2 * NGPTOTG
	1045	ELSE
	1046	IRIRPTSUR = 0
	1047	IRISPTSUR = 0
[5133]	1048	ENDIF
	1049
[5154]	1050	ALLOCATE(NRISTA(NDGSAL - NRIWIDEN:NDGENL + NRIWIDES))
	1051	ALLOCATE(NRIONL(NDGSAL - NRIWIDEN:NDGENL + NRIWIDES))
	1052	ALLOCATE(NRIOFF(NDGSAL - NRIWIDEN:NDGENL + NRIWIDES))
	1053	ALLOCATE(NRIEXT(1 - NDLON:NDLON + NDLON, 1 - NRIWIDEN:NDGENL + NRIWIDES))
	1054	ALLOCATE(NRICORE(NGPTOT))
	1055	ALLOCATE(IRISENDPOS(IRISPTSUR))
	1056	ALLOCATE(IRIRECVPOS(IRIRPTSUR))
	1057	ALLOCATE(IRISENDPTR(NPROC + 1))
	1058	ALLOCATE(IRIRECVPTR(NPROC + 1))
	1059	ALLOCATE(IRICOMM(NPROC))
	1060	ALLOCATE(IRIMAP(4, NDGLG))
	1061	! MPL 1.12.08
	1062	! CALL RDCSET('RI',NRIWIDEN,NRIWIDES,NRIWIDEW,NRIWIDEE,&
	1063	! & IRIRPTSUR,IRISPTSUR,&
	1064	! & NDGLG,NDLON,NDGSAG,NDGENG,IDUM,IDUM,NDGSAL,NDGENL,&
	1065	! & NDSUR1,NDLSUR,NDGSUR,NGPTOT,IDUM,&
	1066	! & NPTRFLOFF,NFRSTLOFF,MYFRSTACTLAT,MYLSTACTLAT,&
	1067	! & NSTA,NONL,NLOENG,NPTRFRSTLAT,NFRSTLAT,NLSTLAT,&
	1068	! & RMU,RSQM2,&
	1069	! & NRISTA,NRIONL,NRIOFF,NRIEXT,NRICORE,NARIB1,&
	1070	! & NRIPROCS,NRIMPBUFSZ,NRIRPT,NRISPT,&
	1071	! & IRISENDPOS,IRIRECVPOS,IRISENDPTR,IRIRECVPTR,IRICOMM,IRIMAP,IRIMAPLEN)
	1072	CALL ABOR1('JUSTE APRES CALL RDCSET COMMENTE')
	1073	WRITE(NULOUT, '("SUECRAD: NARIB1=",I12)')NARIB1
	1074	ALLOCATE(NRISENDPOS(NRISPT))
	1075	ALLOCATE(NRIRECVPOS(NRIRPT))
	1076	ALLOCATE(NRISENDPTR(NRIPROCS + 1))
	1077	ALLOCATE(NRIRECVPTR(NRIPROCS + 1))
	1078	ALLOCATE(NRICOMM(NRIPROCS))
	1079	NRISENDPOS(1:NRISPT) = IRISENDPOS(1:NRISPT)
	1080	NRIRECVPOS(1:NRIRPT) = IRIRECVPOS(1:NRIRPT)
	1081	NRISENDPTR(1:NRIPROCS + 1) = IRISENDPTR(1:NRIPROCS + 1)
	1082	NRIRECVPTR(1:NRIPROCS + 1) = IRIRECVPTR(1:NRIPROCS + 1)
	1083	NRICOMM(1:NRIPROCS) = IRICOMM(1:NRIPROCS)
	1084	DEALLOCATE(IRISENDPOS)
	1085	DEALLOCATE(IRIRECVPOS)
	1086	DEALLOCATE(IRISENDPTR)
	1087	DEALLOCATE(IRIRECVPTR)
	1088	DEALLOCATE(IRICOMM)
	1089	DEALLOCATE(IRIMAP)
[5133]	1090
[5154]	1091	IF(NPROC > 1)THEN
	1092	IRORPTSUR = RADGRID%NGPTOTG
	1093	IROSPTSUR = 2 * RADGRID%NGPTOTG
	1094	ELSE
	1095	IRORPTSUR = 0
	1096	IROSPTSUR = 0
	1097	ENDIF
[5133]	1098
[5154]	1099	ALLOCATE(NROSTA(RADGRID%NDGSAL - NROWIDEN:RADGRID%NDGENL + NROWIDES))
	1100	ALLOCATE(NROONL(RADGRID%NDGSAL - NROWIDEN:RADGRID%NDGENL + NROWIDES))
	1101	ALLOCATE(NROOFF(RADGRID%NDGSAL - NROWIDEN:RADGRID%NDGENL + NROWIDES))
	1102	ALLOCATE(NROEXT(1 - RADGRID%NDLON:RADGRID%NDLON + RADGRID%NDLON, &
	1103	& 1 - NROWIDEN:RADGRID%NDGENL + NROWIDES))
	1104	ALLOCATE(NROCORE(RADGRID%NGPTOT))
	1105	ALLOCATE(IROSENDPOS(IROSPTSUR))
	1106	ALLOCATE(IRORECVPOS(IRORPTSUR))
	1107	ALLOCATE(IROSENDPTR(NPROC + 1))
	1108	ALLOCATE(IRORECVPTR(NPROC + 1))
	1109	ALLOCATE(IROCOMM(NPROC))
	1110	ALLOCATE(IROMAP(4, RADGRID%NDGLG))
	1111	! MPL 1.12.08
	1112	! CALL RDCSET('RO',NROWIDEN,NROWIDES,NROWIDEW,NROWIDEE,&
	1113	! & IRORPTSUR,IROSPTSUR,&
	1114	! & RADGRID%NDGLG,RADGRID%NDLON,RADGRID%NDGSAG,&
	1115	! & RADGRID%NDGENG,IDUM,IDUM,RADGRID%NDGSAL,RADGRID%NDGENL,&
	1116	! & RADGRID%NDSUR1,RADGRID%NDLSUR,RADGRID%NDGSUR,RADGRID%NGPTOT,IDUM,&
	1117	! & RADGRID%NPTRFLOFF,RADGRID%NFRSTLOFF,RADGRID%MYFRSTACTLAT,RADGRID%MYLSTACTLAT,&
	1118	! & RADGRID%NSTA,RADGRID%NONL,RADGRID%NLOENG,RADGRID%NPTRFRSTLAT,&
	1119	! & RADGRID%NFRSTLAT,RADGRID%NLSTLAT,&
	1120	! & RADGRID%RMU,RADGRID%RSQM2,&
	1121	! & NROSTA,NROONL,NROOFF,NROEXT,NROCORE,NAROB1,&
	1122	! & NROPROCS,NROMPBUFSZ,NRORPT,NROSPT,&
	1123	! & IROSENDPOS,IRORECVPOS,IROSENDPTR,IRORECVPTR,IROCOMM,IROMAP,IROMAPLEN)
	1124	CALL ABOR1('JUSTE APRES CALL RDCSET COMMENTE')
	1125	WRITE(NULOUT, '("SUECRAD: NAROB1=",I12)')NAROB1
	1126	ALLOCATE(NROSENDPOS(NROSPT))
	1127	ALLOCATE(NRORECVPOS(NRORPT))
	1128	ALLOCATE(NROSENDPTR(NROPROCS + 1))
	1129	ALLOCATE(NRORECVPTR(NROPROCS + 1))
	1130	ALLOCATE(NROCOMM(NROPROCS))
	1131	NROSENDPOS(1:NROSPT) = IROSENDPOS(1:NROSPT)
	1132	NRORECVPOS(1:NRORPT) = IRORECVPOS(1:NRORPT)
	1133	NROSENDPTR(1:NROPROCS + 1) = IROSENDPTR(1:NROPROCS + 1)
	1134	NRORECVPTR(1:NROPROCS + 1) = IRORECVPTR(1:NROPROCS + 1)
	1135	NROCOMM(1:NROPROCS) = IROCOMM(1:NROPROCS)
	1136	DEALLOCATE(IROSENDPOS)
	1137	DEALLOCATE(IRORECVPOS)
	1138	DEALLOCATE(IROSENDPTR)
	1139	DEALLOCATE(IRORECVPTR)
	1140	DEALLOCATE(IROCOMM)
	1141	DEALLOCATE(IROMAP)
[5133]	1142
[5154]	1143	IF(LLDEBUG)THEN
	1144	WRITE(NULOUT, '("")')
	1145	IRIWIDEMAXN = 0
	1146	IRIWIDEMAXS = 0
	1147	IRIWIDEMAXW = 0
	1148	IRIWIDEMAXE = 0
	1149	IROWIDEMAXN = 0
	1150	IROWIDEMAXS = 0
	1151	IROWIDEMAXW = 0
	1152	IROWIDEMAXE = 0
	1153	IARIB1MAX = 0
	1154	IAROB1MAX = 0
	1155	IWIDE(1) = NRIWIDEN
	1156	IWIDE(2) = NRIWIDES
	1157	IWIDE(3) = NRIWIDEW
	1158	IWIDE(4) = NRIWIDEE
	1159	IWIDE(5) = NROWIDEN
	1160	IWIDE(6) = NROWIDES
	1161	IWIDE(7) = NROWIDEW
	1162	IWIDE(8) = NROWIDEE
	1163	IWIDE(9) = NARIB1
	1164	IWIDE(10) = NAROB1
	1165	IF(MYPROC /= 1)THEN
	1166	stop 'Pas pret pour proc > 1'
	1167	! CALL MPL_SEND(IWIDE(1:10),KDEST=NPRCIDS(1),KTAG=1,CDSTRING='SUECRAD.W')
	1168	ENDIF
	1169	IF(MYPROC == 1)THEN
	1170	DO JROC = 1, NPROC
	1171	IF(JROC /= MYPROC)THEN
	1172	stop 'Pas pret pour proc > 1'
	1173	! CALL MPL_RECV(IWIDE(1:10),KSOURCE=NPRCIDS(JROC),KTAG=1,CDSTRING='SUECRAD.W')
	1174	ENDIF
	1175	WRITE(NULOUT, '("SUECRAD: PROC=",I5,2X,"NRIWIDEN=",I3,2X,"NROWIDEN=",I3 )')&
	1176	& JROC, IWIDE(1), IWIDE(5)
	1177	WRITE(NULOUT, '("SUECRAD: PROC=",I5,2X,"NRIWIDES=",I3,2X,"NROWIDES=",I3 )')&
	1178	& JROC, IWIDE(2), IWIDE(6)
	1179	WRITE(NULOUT, '("SUECRAD: PROC=",I5,2X,"NRIWIDEW=",I3,2X,"NROWIDEW=",I3 )')&
	1180	& JROC, IWIDE(3), IWIDE(7)
	1181	WRITE(NULOUT, '("SUECRAD: PROC=",I5,2X,"NRIWIDEE=",I3,2X,"NROWIDEE=",I3 )')&
	1182	& JROC, IWIDE(4), IWIDE(8)
	1183	WRITE(NULOUT, '("SUECRAD: PROC=",I5,2X,"NARIB1=",I10,2X,"NAROB1=",I10 )')&
	1184	& JROC, IWIDE(9), IWIDE(10)
	1185	WRITE(NULOUT, '("")')
	1186	IF(IWIDE(1) > IRIWIDEMAXN) IRIWIDEMAXN = IWIDE(1)
	1187	IF(IWIDE(2) > IRIWIDEMAXS) IRIWIDEMAXS = IWIDE(2)
	1188	IF(IWIDE(3) > IRIWIDEMAXW) IRIWIDEMAXW = IWIDE(3)
	1189	IF(IWIDE(4) > IRIWIDEMAXE) IRIWIDEMAXE = IWIDE(4)
	1190	IF(IWIDE(5) > IROWIDEMAXN) IROWIDEMAXN = IWIDE(5)
	1191	IF(IWIDE(6) > IROWIDEMAXS) IROWIDEMAXS = IWIDE(6)
	1192	IF(IWIDE(7) > IROWIDEMAXW) IROWIDEMAXW = IWIDE(7)
	1193	IF(IWIDE(8) > IROWIDEMAXE) IROWIDEMAXE = IWIDE(8)
	1194	IF(IWIDE(9) > IARIB1MAX) IARIB1MAX = IWIDE(9)
	1195	IF(IWIDE(10) > IAROB1MAX) IAROB1MAX = IWIDE(10)
	1196	ENDDO
	1197	WRITE(NULOUT, '("")')
	1198	WRITE(NULOUT, '("SUECRAD: NRIWIDEN(MAX) =",I8)')IRIWIDEMAXN
	1199	WRITE(NULOUT, '("SUECRAD: NRIWIDES(MAX) =",I8)')IRIWIDEMAXS
	1200	WRITE(NULOUT, '("SUECRAD: NRIWIDEW(MAX) =",I8)')IRIWIDEMAXW
	1201	WRITE(NULOUT, '("SUECRAD: NRIWIDEE(MAX) =",I8)')IRIWIDEMAXE
	1202	WRITE(NULOUT, '("SUECRAD: NROWIDEN(MAX) =",I8)')IROWIDEMAXN
	1203	WRITE(NULOUT, '("SUECRAD: NROWIDES(MAX) =",I8)')IROWIDEMAXS
	1204	WRITE(NULOUT, '("SUECRAD: NROWIDEW(MAX) =",I8)')IROWIDEMAXW
	1205	WRITE(NULOUT, '("SUECRAD: NROWIDEE(MAX) =",I8)')IROWIDEMAXE
	1206	WRITE(NULOUT, '("SUECRAD: NARIB1(MAX) =",I10)')IARIB1MAX
	1207	WRITE(NULOUT, '("SUECRAD: NAROB1(MAX) =",I10)')IAROB1MAX
	1208	WRITE(NULOUT, '("")')
	1209	ENDIF
	1210	CALL FLUSH(NULOUT)
	1211	ENDIF
[5133]	1212
	1213	ENDIF
[5154]	1214	! CALL GSTATS(1818,1) MPL 2.12.08
[5133]	1215
[5154]	1216	ELSE
	1217
	1218	WRITE(NULOUT, '("SUECRAD: INVALID VALUE FOR NRADINT=",I6)')NRADINT
	1219	CALL ABOR1('SUECRAD: NRADINT INVALID')
	1220
[5133]	1221	ENDIF
	1222
[5154]	1223	ENDIF ! END OF LERADI BLOCK
[5133]	1224
[5154]	1225	! ----------------------------------------------------------------
[5133]	1226
[5154]	1227	!* 4. INITIALIZE RADIATION COEFFICIENTS.
	1228	! ----------------------------------
[5133]	1229
[5154]	1230	RCDAY = RDAY * RG / RCPD
	1231	DIFF = 1.66_JPRB
	1232	R10E = 0.4342945_JPRB
[5133]	1233
[5154]	1234	! CALL GSTATS(1818,0) MPL 2.12.08
	1235	CALL SURDI
[5133]	1236
[5154]	1237	IF (NINHOM == 0) THEN
	1238	RLWINHF = 1._JPRB
	1239	RSWINHF = 1._JPRB
	1240	ENDIF
[5133]	1241
[5154]	1242	! ----------------------------------------------------------------
[5133]	1243
[5154]	1244	!* 5. INITIALIZE RADIATION ABSORPTION COEFFICIENTS
	1245	! --------------------------------------------
[5133]	1246
[5154]	1247	!* 5.1. Initialization routine for RRTM
	1248	! -------------------------------
[5133]	1249
[5154]	1250	CALL SURRTAB
	1251	CALL SURRTPK
	1252	CALL SURRTRF
	1253	CALL SURRTFTR
[5133]	1254
[5154]	1255	IF (LRRTM) THEN
	1256	IF (KLEV > JPLAY) THEN
	1257	WRITE(UNIT = KULOUT, &
	1258	& FMT = '('' RRTM MAXIMUM NUMBER OF LAYERS IS REACHED'',&
	1259	& '' CALL ABORT'')')
	1260	CALL ABOR1(' ABOR1 CALLED SUECRAD')
	1261	ENDIF
[5133]	1262
[5154]	1263	! Read the absorption coefficient data and reduce from 256 to 140 g-points
[5133]	1264
[5154]	1265	CALL RRTM_INIT_140GP
[5133]	1266
[5154]	1267	INBLW = 16
[5133]	1268
[5154]	1269	ELSE
	1270	INBLW = 6
[5133]	1271
[5154]	1272	ENDIF
[5133]	1273
[5154]	1274	CALL SULWN
	1275	CALL SUSWN (NTSW, NSW)
	1276	CALL SUCLOPN (NTSW, NSW, KLEV)
[5133]	1277
[5154]	1278	!-- routines specific to SRTM
	1279	IF (LSRTM) THEN
	1280	NTSW = 14
	1281	ISW = 14
	1282	CALL SRTM_INIT
	1283	CALL SUSRTAER
	1284	CALL SUSRTCOP
	1285	WRITE(UNIT = KULOUT, FMT = '(''SRTM Configuration'',L8,3I4)')LSRTM, NTSW, ISW, JPGPT
[5133]	1286
[5154]	1287	ELSE
	1288	IF (.NOT.LONEWSW .OR. ((NSW /= 2).AND.(NSW /= 4).AND.(NSW /= 6))) THEN
	1289	WRITE(UNIT = KULOUT, FMT = '(''Wrong SW Configuration'',L8,I3)')LONEWSW, NSW
	1290	ENDIF
[5133]	1291
[5154]	1292	CALL SUSWN (NTSW, NSW)
	1293	CALL SUAERSN (NTSW, NSW)
[5133]	1294	ENDIF
[5154]	1295	WRITE(UNIT = KULOUT, FMT = '('' NLW,NTSW,NSW SET EQUAL TO:'',3I3)') INBLW, NTSW, NSW
[5133]	1296
	1297
[5154]	1298	!-- routine specific to the UV processor
	1299	IF (LUVPROC) THEN
	1300	NUVTIM = NUVTIM * 86400
	1301	CALL SU_UVRAD (NUV)
	1302	ENDIF
[5133]	1303
[5154]	1304	! ----------------------------------------------------------------
[5133]	1305
[5154]	1306	!* 6. INITIALIZE AEROSOL OPTICAL PARAMETERS AND DISTRIBUTION
	1307	! ------------------------------------------------------
[5133]	1308
[5154]	1309	!- LW optical properties
	1310	CALL SUAERL
	1311	!- SW optical properties moved above
	1312	!CALL SUAERSN (NTSW,NSW)
[5133]	1313
[5154]	1314	!- horizontal distribution
	1315	CALL SUAERH
[5133]	1316
[5154]	1317	!- vertical distribution
	1318	CALL SUAERV (KLEV, PETAH, &
	1319	& CVDAES, CVDAEL, CVDAEU, CVDAED, &
	1320	& RCTRBGA, RCVOBGA, RCSTBGA, RCAEOPS, RCAEOPL, RCAEOPU, &
	1321	& RCAEOPD, RCTRPT, RCAEADK, RCAEADM, RCAEROS &
	1322	&)
[5133]	1323
[5154]	1324	!-- Overlap function (only used if NOVLP=4)
	1325	! Appel supprime par MPL (30042010) car NOVLP=4 pas utilise
	1326	! sinon il faudrait calculer le geopotentiel STZ
	1327	!CALL SUOVLP ( KLEV )
[5133]	1328
[5154]	1329	!-- parameters for prognostic aerosols
	1330	CALL SU_AERW
[5133]	1331
[5154]	1332	! ----------------------------------------------------------------
[5133]	1333
[5154]	1334	!* 7. INITIALIZE SATELLITE GEOMETRICAL/RADIOMETRIC PARAMETERS
	1335	! -------------------------------------------------------
[5133]	1336
[5154]	1337	IF (LEPHYS .AND. NMODE > 1) THEN
	1338	CALL SUSAT
	1339	ENDIF
	1340	!CALL GSTATS(1818,1) MPL 2.12.08
[5133]	1341
[5154]	1342	! ----------------------------------------------------------------
[5133]	1343
[5154]	1344	!* 8. INITIALIZE CLIMATOLOGICAL OZONE DISTRIBUTION
	1345	! --------------------------------------------
	1346	! (not done here!!! called from APLPAR as it depends
	1347	! on model pressure levels!)
[5133]	1348
[5154]	1349	! ----------------------------------------------------------------
[5133]	1350
[5154]	1351	!* 9. SET UP MODEL CONFIGURATION FOR TIME-SPACE INTERPOLATION
	1352	! -------------------------------------------------------
[5133]	1353
[5154]	1354	ZTSTEP = MAX(TSTEP, 1.0_JPRB)
	1355	ZSTPHR = 3600._JPRB / ZTSTEP
	1356	IRADFR = NRADFR
	1357	IF(NRADFR < 0) THEN
	1358	NRADFR = -NRADFR * ZSTPHR + 0.5_JPRB
	1359	ENDIF
	1360	NRADPFR = NRADPFR * NRADFR
	1361	IF (MOD(NRADPLA, 2) == 0.AND. NRADPLA /= 0) THEN
	1362	NRADPLA = NRADPLA + 1
	1363	ENDIF
[5133]	1364
[5154]	1365	IF(NRADUV < 0) THEN
	1366	NRADUV = -NRADUV * ZSTPHR + 0.5_JPRB
	1367	ENDIF
[5133]	1368
[5154]	1369	IST1HR = ZSTPHR + 0.05_JPRB
	1370	ISTNHR = NLNGR1H * ZSTPHR + 0.05_JPRB
	1371	IF (MOD(3600._JPRB, ZTSTEP) > 0.1_JPRB) THEN
	1372	801 CONTINUE
	1373	IST1HR = IST1HR + 1
	1374	IF (MOD(ISTNHR, IST1HR) /= 0) GO TO 801
	1375	ENDIF
	1376	IF (NRADFR == 1) THEN
	1377	NRADSFR = NRADFR
	1378	ELSE
	1379	NRADSFR = IST1HR
	1380	ENDIF
	1381	NRADNFR = NRADFR
[5133]	1382
[5154]	1383	IF(LRAYFM) THEN
	1384	NRPROMA = NDLON + 6 + (1 - MOD(NDLON, 2))
	1385	ENDIF
[5133]	1386
[5154]	1387	! ----------------------------------------------------------------
[5133]	1388
[5154]	1389	!* 10. ALLOCATE WORK ARRAYS
	1390	! --------------------
[5133]	1391
[5154]	1392	IU = NULOUT
	1393	LLP = NPRINTLEV >= 1.OR. LALLOPR
[5133]	1394
[5154]	1395	IF (LEPHYS) THEN
	1396	ALLOCATE(EMTD(NPROMA, NFLEVG + 1, NGPBLKS))
	1397	IF(LLP)WRITE(IU, 9) 'EMTD ', SIZE(EMTD), SHAPE(EMTD)
	1398	ALLOCATE(TRSW(NPROMA, NFLEVG + 1, NGPBLKS))
	1399	IF(LLP)WRITE(IU, 9) 'TRSW ', SIZE(TRSW), SHAPE(TRSW)
	1400	ALLOCATE(EMTC(NPROMA, NFLEVG + 1, NGPBLKS))
	1401	IF(LLP)WRITE(IU, 9) 'EMTC ', SIZE(EMTC), SHAPE(EMTC)
	1402	ALLOCATE(TRSC(NPROMA, NFLEVG + 1, NGPBLKS))
	1403	IF(LLP)WRITE(IU, 9) 'TRSC ', SIZE(TRSC), SHAPE(TRSC)
	1404	ALLOCATE(SRSWD(NPROMA, NGPBLKS))
	1405	IF(LLP)WRITE(IU, 9) 'SRSWD ', SIZE(SRSWD), SHAPE(SRSWD)
	1406	ALLOCATE(SRLWD(NPROMA, NGPBLKS))
	1407	IF(LLP)WRITE(IU, 9) 'SRLWD ', SIZE(SRLWD), SHAPE(SRLWD)
	1408	ALLOCATE(SRSWDCS(NPROMA, NGPBLKS))
	1409	IF(LLP)WRITE(IU, 9) 'SRSWDCS ', SIZE(SRSWDCS), SHAPE(SRSWDCS)
	1410	ALLOCATE(SRLWDCS(NPROMA, NGPBLKS))
	1411	IF(LLP)WRITE(IU, 9) 'SRLWDCS ', SIZE(SRLWDCS), SHAPE(SRLWDCS)
	1412	ALLOCATE(SRSWDV(NPROMA, NGPBLKS))
	1413	IF(LLP)WRITE(IU, 9) 'SRSWDV ', SIZE(SRSWDV), SHAPE(SRSWDV)
	1414	ALLOCATE(SRSWDUV(NPROMA, NGPBLKS))
	1415	IF(LLP)WRITE(IU, 9) 'SRSWDUV ', SIZE(SRSWDUV), SHAPE(SRSWDUV)
	1416	ALLOCATE(EDRO(NPROMA, NGPBLKS))
	1417	IF(LLP)WRITE(IU, 9) 'EDRO ', SIZE(EDRO), SHAPE(EDRO)
	1418	ALLOCATE(SRSWPAR(NPROMA, NGPBLKS))
	1419	IF(LLP)WRITE(IU, 9) 'SRSWPAR ', SIZE(SRSWPAR), SHAPE(SRSWPAR)
	1420	ALLOCATE(SRSWUVB(NPROMA, NGPBLKS))
	1421	IF(LLP)WRITE(IU, 9) 'SRSWUVB ', SIZE(SRSWUVB), SHAPE(SRSWUVB)
[5133]	1422
[5154]	1423	ELSEIF(LMPHYS .AND. (LRAYFM.OR.LRAYFM15)) THEN
	1424	ALLOCATE(EMTD(NPROMA, NFLEVG + 1, NGPBLKS))
	1425	IF(LLP)WRITE(IU, 9) 'EMTD ', SIZE(EMTD), SHAPE(EMTD)
	1426	ALLOCATE(TRSW(NPROMA, NFLEVG + 1, NGPBLKS))
	1427	IF(LLP)WRITE(IU, 9) 'TRSW ', SIZE(TRSW), SHAPE(TRSW)
	1428	ALLOCATE(EMTU(NPROMA, NFLEVG + 1, NGPBLKS))
	1429	IF(LLP)WRITE(IU, 9) 'EMTC ', SIZE(EMTU), SHAPE(EMTU)
	1430	ALLOCATE(RMOON(NPROMA, NGPBLKS))
	1431	IF(LLP)WRITE(IU, 9) 'RMOON ', SIZE(RMOON), SHAPE(RMOON)
	1432	ENDIF
	1433	ALLOCATE(SRSWPARC(NPROMA, NGPBLKS))
	1434	IF(LLP)WRITE(IU, 9) 'SRSWPARC ', SIZE(SRSWPARC), SHAPE(SRSWPARC)
	1435	ALLOCATE(SRSWTINC(NPROMA, NGPBLKS))
	1436	IF(LLP)WRITE(IU, 9) 'SRSWTINC ', SIZE(SRSWTINC), SHAPE(SRSWTINC)
[5133]	1437
[5154]	1438	9 FORMAT(1X, 'ARRAY ', A10, ' ALLOCATED ', 8I8)
[5133]	1439
[5154]	1440	! ----------------------------------------------------------------
[5133]	1441
[5154]	1442	!* 10. PRINT FINAL VALUES.
	1443	! -------------------
[5133]	1444
[5154]	1445	IF (LOUTPUT) THEN
	1446	WRITE(UNIT = KULOUT, FMT = '('' COMMON YOERAD '')')
	1447	WRITE(UNIT = KULOUT, FMT = '('' LERADI = '',L5 &
	1448	& ,'' LERAD1H = '',L5,'' LECO2VAR= '',L5,'' LHGHG = '',L5 &
	1449	& ,'' NLNGR1H = '',I2,'' NRADSFR = '',I2)')&
	1450	& LERADI, LERAD1H, LECO2VAR, LHGHG, NLNGR1H, NRADSFR
	1451	WRITE(UNIT = KULOUT, FMT = '('' LEPO3RA = '',L5,'' YO3%LGP = '',L5 )') LEPO3RA, YO3%LGP
	1452	WRITE(UNIT = KULOUT, FMT = '('' NRADFR = '',I2 &
	1453	& ,'' NRADPFR = '',I3 &
	1454	& ,'' NRADPLA = '',I2 &
	1455	& ,'' NRINT = '',I1 &
	1456	& ,'' NRPROMA = '',I5 &
	1457	& )')&
	1458	& NRADFR, NRADPFR, NRADPLA, NRINT, NRPROMA
	1459	WRITE(UNIT = KULOUT, FMT = '('' LERADHS= '',L5 &
	1460	& ,'' LRRTM = '',L5 &
	1461	& ,'' LSRTM = '',L5 &
	1462	& ,'' NMODE = '',I1 &
	1463	& ,'' NOZOCL= '',I1 &
	1464	& ,'' NAER = '',I1 &
	1465	& ,'' NHINCSOL='',I2 &
	1466	& )')&
	1467	& LERADHS, LRRTM, LSRTM, NMODE, NOZOCL, NAER, NHINCSOL
	1468	IF (.NOT.LHGHG .AND. .NOT.LECO2VAR) WRITE(UNIT = KULOUT, FMT = '('' RCCO2= '',E10.3 &
	1469	&,'' RCCH4= '',E10.3,'' RCN2O= '',E10.3,'' RCCFC11= '',E10.3,'' RCFC12= '',E10.3 &
	1470	&)')&
	1471	& RCCO2, RCCH4, RCN2O, RCCFC11, RCCFC12
	1472	WRITE(UNIT = KULOUT, FMT = '('' NINHOM = '',I1 &
	1473	& ,'' NLAYINH='',I1 &
	1474	& ,'' RLWINHF='',F4.2 &
	1475	& ,'' RSWINHF='',F4.2 &
	1476	& )')&
	1477	& NINHOM, NLAYINH, RLWINHF, RSWINHF
	1478	IF (NPERTAER /= 0 .OR. NPERTOZ /= 0) THEN
	1479	WRITE(UNIT = KULOUT, FMT = '('' NPERTAER= '',I2 &
	1480	& ,'' LNOTROAER='',L5 &
	1481	& ,'' NPERTOZ = '',I1 &
	1482	& ,'' RPERTOZ = '',F5.0 &
	1483	& )')&
	1484	& NPERTAER, LNOTROAER, NPERTOZ, RPERTOZ
	1485	ENDIF
	1486	WRITE(UNIT = KULOUT, FMT = '('' NRADINT = '',I2)')NRADINT
	1487	WRITE(UNIT = KULOUT, FMT = '('' NRADRES = '',I4)')NRADRES
	1488	WRITE(UNIT = KULOUT, FMT = '('' LRADONDEM = '',L5)')LRADONDEM
	1489	IF(NRADINT > 0)THEN
	1490	IDIR = LEN_TRIM(CRTABLEDIR)
	1491	IFIL = LEN_TRIM(CRTABLEFIL)
	1492	WRITE(UNIT = KULOUT, FMT = '('' CRTABLEDIR = '',A,'' CRTABLEFIL = '',A)')&
	1493	& CRTABLEDIR(1:IDIR), CRTABLEFIL(1:IFIL)
	1494	ENDIF
	1495	WRITE(UNIT = KULOUT, FMT = '('' LCCNL = '',L5 &
	1496	& ,'' LCCNO = '',L5 &
	1497	& ,'' RCCNLND= '',F5.0 &
	1498	& ,'' RCCNSEA= '',F5.0 &
	1499	& ,'' LE4ALB = '',L5 &
	1500	&)')&
	1501	& LCCNL, LCCNO, RCCNLND, RCCNSEA, LE4ALB
	1502	IF (LHVOLCA) THEN
	1503	WRITE(UNIT = KULOUT, FMT = '('' HISTORY OF VOLCANIC AEROSOLS= '',L5)')LHVOLCA
	1504	ENDIF
	1505	WRITE(UNIT = KULOUT, FMT = '('' LONEWSW= '',L5 &
	1506	& ,'' NRADIP = '',I1 &
	1507	& ,'' NRADLP = '',I1 &
	1508	& ,'' NICEOPT= '',I1 &
	1509	& ,'' NLIQOPT= '',I1 &
	1510	& ,'' LDIFFC = '',L5 &
	1511	& )')&
	1512	& LONEWSW, NRADIP, NRADLP, NICEOPT, NLIQOPT, LDIFFC
	1513	WRITE(UNIT = KULOUT, FMT = '('' WARNING! CLOUD OVERLAP ASSUMPT. IS''&
	1514	& ,'' NOVLP = '',I2 &
	1515	& )')&
	1516	& NOVLP
	1517	IF (LUVPROC) THEN
	1518	IDAYUV = NUVTIM / 86400
	1519	WRITE(UNIT = KULOUT, FMT = '('' LUVPROC = '',L5 &
	1520	& ,'' LUVTDEP= '',L5 &
	1521	& ,'' NRADUV = '',I2 &
	1522	& ,'' NUV = '',I2 &
	1523	& ,'' NDAYUV = '',I5 &
	1524	& ,'' RMUZUV = '',E9.3 &
	1525	& )')&
	1526	& LUVPROC, LUVTDEP, NRADUV, NUV, IDAYUV, RMUZUV
	1527	WRITE(UNIT = KULOUT, FMT = '('' RUVLAM = '',24F6.1)') (RUVLAM(JUV), JUV = 1, NUV)
	1528	WRITE(UNIT = KULOUT, FMT = '('' JUVLAM = '',24(3X,I1,2X))') (JUVLAM(JUV), JUV = 1, NUV)
	1529	ENDIF
	1530	WRITE(UNIT = KULOUT, FMT = '('' NMCICA= '',I2 &
	1531	& )')&
	1532	& NMCICA
[5133]	1533	ENDIF
	1534
[5154]	1535	! ------------------------------------------------------------------
[5133]	1536
[5154]	1537	IF (LHOOK) CALL DR_HOOK('SUECRAD', 1, ZHOOK_HANDLE)
[5133]	1538	END SUBROUTINE SUECRAD

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

Download in other formats:

Original Format