Context Navigation

conf_phys.F @ 2467

Last change on this file since 2467 was 2447, checked in by cmathe, 4 years ago
CO2 microphysics: correction r_sedim (co2cloud.F); add: co2 cloud radiatively active
File size: 41.0 KB

Rev	Line
[1246]	1	SUBROUTINE conf_phys(ngrid,nlayer,nq)
[1226]	2
[42]	3	!=======================================================================
	4	!
	5	! purpose:
	6	! -------
	7	!
	8	! Initialisation for the physical parametrisations of the LMD
	9	! martian atmospheric general circulation modele.
	10	!
	11	! author: Frederic Hourdin 15 / 10 /93
	12	! -------
	13	! modified: Sebastien Lebonnois 11/06/2003 (new callphys.def)
	14	! Ehouarn Millour (oct. 2008) tracers are now identified
	15	! by their names and may not be contiguously
	16	! stored in the q(:,:,:,:) array
	17	! E.M. (june 2009) use getin routine to load parameters
[234]	18	! adapted to the mesoscale use - Aymeric Spiga - 01/2007-07/2011
[1224]	19	! separated inifis into conf_phys and phys_state_var_init (A. Spiga)
[42]	20	!
	21	!
	22	! arguments:
	23	! ----------
	24	!
	25	! input:
	26	! ------
	27	!
[1224]	28	! nq Number of tracers
[42]	29	!
	30	!=======================================================================
	31	!
	32	!-----------------------------------------------------------------------
	33	! declarations:
	34	! -------------
[2281]	35	USE ioipsl_getin_p_mod, ONLY : getin_p
[1617]	36	use tracer_mod, only : nuice_sed, ccn_factor, nuiceco2_sed,
[1720]	37	& nuice_ref,nuiceco2_ref
[1224]	38	use surfdat_h, only: albedo_h2o_ice, inert_h2o_ice,
[1047]	39	& frost_albedo_threshold
[1525]	40	use time_phylmdz_mod, only: ecritphy,day_step,iphysiq,ecritstart,
	41	& daysec,dtphys
[1246]	42	use dimradmars_mod, only: naerkind, name_iaer,
	43	& ini_scatterers,tauvis
[1918]	44	use datafile_mod, only: datadir
[2164]	45	use calchim_mod, only: ichemistry
[2184]	46	use co2condens_mod, only: scavco2cond
[2409]	47	use dust_param_mod, only: dustbin, doubleq, submicron, active,
[2413]	48	& lifting, freedust, callddevil,
[2417]	49	& dustscaling_mode
[2409]	50	use aeropacity_mod, only: iddist, topdustref
[42]	51	IMPLICIT NONE
[1918]	52	include "callkeys.h"
	53	include "microphys.h"
[1112]	54
[1246]	55	INTEGER,INTENT(IN) :: ngrid,nlayer,nq
[1525]	56	INTEGER ierr,j
[2304]	57	character(len=20),parameter :: modname="conf_phys"
[42]	58
	59	CHARACTER ch1*12
[1579]	60	#ifndef MESOSCALE
[1525]	61	! read in some parameters from "run.def" for physics,
	62	! or shared between dynamics and physics.
	63	ecritphy=240 ! default value
[2304]	64	call getin_p("ecritphy",ecritphy) ! frequency of outputs in physics,
[1525]	65	! in dynamical steps
	66	day_step=960 ! default value
[2304]	67	call getin_p("day_step",day_step) ! number of dynamical steps per day
[1525]	68	iphysiq=20 ! default value
[2304]	69	call getin_p("iphysiq",iphysiq) ! call physics every iphysiq dyn step
[1525]	70	ecritstart=0 ! default value
[2304]	71	call getin_p("ecritstart",ecritstart) ! write a restart every ecristart steps
[1579]	72	#endif
	73
[42]	74	! --------------------------------------------------------------
	75	! Reading the "callphys.def" file controlling some key options
	76	! --------------------------------------------------------------
	77
	78	! check that 'callphys.def' file is around
	79	OPEN(99,file='callphys.def',status='old',form='formatted'
	80	& ,iostat=ierr)
	81	CLOSE(99)
	82
	83	IF(ierr.EQ.0) THEN
	84	PRINT*
	85	PRINT*
	86	PRINT*,'--------------------------------------------'
[1224]	87	PRINT*,' conf_phys: Parameters for the physics (callphys.def)'
[42]	88	PRINT*,'--------------------------------------------'
	89
	90	write(,) "Directory where external input files are:"
[1918]	91	! default path is set in datafile_mod
[2304]	92	call getin_p("datadir",datadir)
[1918]	93	write(,) " datadir = ",trim(datadir)
[42]	94
[2281]	95	write(,) "Initialize physics with startfi.nc file ?"
	96	startphy_file=.true.
	97	call getin_p("startphy_file",startphy_file)
	98	write(,) "startphy_file", startphy_file
	99
[42]	100	write(,) "Run with or without tracer transport ?"
	101	tracer=.false. ! default value
[2304]	102	call getin_p("tracer",tracer)
[42]	103	write(,) " tracer = ",tracer
	104
	105	write(,) "Diurnal cycle ?"
	106	write(,) "(if diurnal=False, diurnal averaged solar heating)"
	107	diurnal=.true. ! default value
[2304]	108	call getin_p("diurnal",diurnal)
[42]	109	write(,) " diurnal = ",diurnal
	110
	111	write(,) "Seasonal cycle ?"
	112	write(,) "(if season=False, Ls stays constant, to value ",
	113	& "set in 'start'"
	114	season=.true. ! default value
[2304]	115	call getin_p("season",season)
[42]	116	write(,) " season = ",season
	117
	118	write(,) "Write some extra output to the screen ?"
	119	lwrite=.false. ! default value
[2304]	120	call getin_p("lwrite",lwrite)
[42]	121	write(,) " lwrite = ",lwrite
	122
	123	write(,) "Save statistics in file stats.nc ?"
[226]	124	#ifdef MESOSCALE
[42]	125	callstats=.false. ! default value
[226]	126	#else
	127	callstats=.true. ! default value
	128	#endif
[2304]	129	call getin_p("callstats",callstats)
[42]	130	write(,) " callstats = ",callstats
	131
	132	write(,) "Save EOF profiles in file 'profiles' for ",
	133	& "Climate Database?"
	134	calleofdump=.false. ! default value
[2304]	135	call getin_p("calleofdump",calleofdump)
[42]	136	write(,) " calleofdump = ",calleofdump
	137
	138	write(,) "Dust scenario: 1=constant dust (read from startfi",
	139	& " or set as tauvis); 2=Viking scenario; =3 MGS scenario,",
[677]	140	& "=6 cold (low dust) scenario; =7 warm (high dust) scenario ",
	141	& "=24,25 ... 30 :Mars Year 24, ... or 30 from TES assimilation"
[42]	142	iaervar=3 ! default value
[2304]	143	call getin_p("iaervar",iaervar)
[42]	144	write(,) " iaervar = ",iaervar
	145
[627]	146	write(,) "Reference (visible) dust opacity at 610 Pa ",
[42]	147	& "(matters only if iaervar=1)"
	148	! NB: default value of tauvis is set/read in startfi.nc file
[2304]	149	call getin_p("tauvis",tauvis)
[42]	150	write(,) " tauvis = ",tauvis
	151
	152	write(,) "Dust vertical distribution:"
	153	write(,) "(=1 top set by topdustref parameter;",
	154	& " =2 Viking scenario; =3 MGS scenario)"
	155	iddist=3 ! default value
[2304]	156	call getin_p("iddist",iddist)
[42]	157	write(,) " iddist = ",iddist
	158
	159	write(,) "Dust top altitude (km). (Matters only if iddist=1)"
	160	topdustref= 90.0 ! default value
[2304]	161	call getin_p("topdustref",topdustref)
[42]	162	write(,) " topdustref = ",topdustref
	163
[544]	164	write(,) "Prescribed surface thermal flux (H/(rho*cp),K m/s)"
	165	tke_heat_flux=0. ! default value
[2304]	166	call getin_p("tke_heat_flux",tke_heat_flux)
[544]	167	write(,) " tke_heat_flux = ",tke_heat_flux
	168	write(,) " 0 means the usual schemes are computing"
	169
[42]	170	write(,) "call radiative transfer ?"
	171	callrad=.true. ! default value
[2304]	172	call getin_p("callrad",callrad)
[42]	173	write(,) " callrad = ",callrad
	174
[234]	175	write(,) "call slope insolation scheme ?",
	176	& "(matters only if callrad=T)"
	177	#ifdef MESOSCALE
	178	callslope=.true. ! default value
	179	#else
	180	callslope=.false. ! default value (not supported yet)
	181	#endif
[2304]	182	call getin_p("callslope",callslope)
[234]	183	write(,) " callslope = ",callslope
	184
[42]	185	write(,) "call NLTE radiative schemes ?",
	186	& "(matters only if callrad=T)"
	187	callnlte=.false. ! default value
[2304]	188	call getin_p("callnlte",callnlte)
[42]	189	write(,) " callnlte = ",callnlte
	190
[414]	191	nltemodel=0 !default value
	192	write(,) "NLTE model?"
	193	write(,) "0 -> old model, static O"
	194	write(,) "1 -> old model, dynamic O"
	195	write(,) "2 -> new model"
	196	write(,) "(matters only if callnlte=T)"
[2304]	197	call getin_p("nltemodel",nltemodel)
[414]	198	write(,) " nltemodel = ",nltemodel
	199
[42]	200	write(,) "call CO2 NIR absorption ?",
	201	& "(matters only if callrad=T)"
	202	callnirco2=.false. ! default value
[2304]	203	call getin_p("callnirco2",callnirco2)
[42]	204	write(,) " callnirco2 = ",callnirco2
[414]	205
	206	write(,) "New NIR NLTE correction ?",
	207	$ "0-> old model (no correction)",
	208	$ "1-> new correction",
	209	$ "(matters only if callnirco2=T)"
[577]	210	#ifdef MESOSCALE
	211	nircorr=0 !default value. this is OK below 60 km.
	212	#else
[633]	213	nircorr=0 !default value
[577]	214	#endif
[2304]	215	call getin_p("nircorr",nircorr)
[414]	216	write(,) " nircorr = ",nircorr
	217
[42]	218	write(,) "call turbulent vertical diffusion ?"
	219	calldifv=.true. ! default value
[2304]	220	call getin_p("calldifv",calldifv)
[42]	221	write(,) " calldifv = ",calldifv
	222
[162]	223	write(,) "call thermals ?"
	224	calltherm=.false. ! default value
[2304]	225	call getin_p("calltherm",calltherm)
[162]	226	write(,) " calltherm = ",calltherm
	227
[42]	228	write(,) "call convective adjustment ?"
	229	calladj=.true. ! default value
[2304]	230	call getin_p("calladj",calladj)
[42]	231	write(,) " calladj = ",calladj
[1239]	232
	233	if (calltherm .and. calladj) then
	234	print*,'!!! PLEASE NOTE !!!'
	235	print*,'convective adjustment is on'
	236	print*,'but since thermal plume model is on'
	237	print*,'convadj is only activated above the PBL'
[162]	238	endif
[1239]	239
	240	write(,) "used latest version of yamada scheme?"
[1240]	241	callyamada4=.true. ! default value
[2304]	242	call getin_p("callyamada4",callyamada4)
[1240]	243	write(,) " callyamada4 = ",callyamada4
[42]	244
[1240]	245	if (calltherm .and. .not.callyamada4) then
[1239]	246	print*,'!!!! WARNING WARNING WARNING !!!!'
	247	print*,'if calltherm=T we strongly advise that '
[1240]	248	print*,'you set the flag callyamada4 to T '
[1239]	249	print*,'!!!! WARNING WARNING WARNING !!!!'
	250	endif
	251
[284]	252	write(,) "call Richardson-based surface layer ?"
	253	callrichsl=.false. ! default value
[2304]	254	call getin_p("callrichsl",callrichsl)
[284]	255	write(,) " callrichsl = ",callrichsl
	256
[288]	257	if (calltherm .and. .not.callrichsl) then
	258	print*,'WARNING WARNING WARNING'
	259	print*,'if calltherm=T we strongly advise that '
	260	print*,'you use the new surface layer scheme '
	261	print*,'by setting callrichsl=T '
	262	endif
	263
[553]	264	if (calladj .and. callrichsl .and. (.not. calltherm)) then
[551]	265	print*,'You should not be calling the convective adjustment
[553]	266	& scheme with the Richardson surface-layer and without the thermals
	267	&. This approach is not
[551]	268	& physically consistent and can lead to unrealistic friction
	269	& values.'
	270	print*,'If you want to use the Ri. surface-layer, either
	271	& activate thermals OR de-activate the convective adjustment.'
[2304]	272	call abort_physic(modname,
	273	& "Richardson layer must be used with thermals",1)
[551]	274	endif
[566]	275
[42]	276	write(,) "call CO2 condensation ?"
	277	callcond=.true. ! default value
[2304]	278	call getin_p("callcond",callcond)
[42]	279	write(,) " callcond = ",callcond
	280
	281	write(,)"call thermal conduction in the soil ?"
	282	callsoil=.true. ! default value
[2304]	283	call getin_p("callsoil",callsoil)
[42]	284	write(,) " callsoil = ",callsoil
	285
	286
	287	write(,)"call Lott's gravity wave/subgrid topography ",
	288	& "scheme ?"
	289	calllott=.true. ! default value
[2304]	290	call getin_p("calllott",calllott)
[42]	291	write(,)" calllott = ",calllott
	292
[2149]	293	write(,)"call Lott's non-oro GWs parameterisation ",
	294	& "scheme ?"
	295	calllott_nonoro=.false. ! default value
[2304]	296	call getin_p("calllott_nonoro",calllott_nonoro)
[2149]	297	write(,)" calllott_nonoro = ",calllott_nonoro
	298
[1974]	299	! rocket dust storm injection scheme
[2179]	300	write(,)"call rocket dust storm parametrization"
[1974]	301	rdstorm=.false. ! default value
[2304]	302	call getin_p("rdstorm",rdstorm)
[1974]	303	write(,)" rdstorm = ",rdstorm
[2160]	304	! rocket dust storm detrainment coefficient
[2411]	305	coeff_detrainment=0.05 ! default value
[2304]	306	call getin_p("coeff_detrainment",coeff_detrainment)
[2160]	307	write(,)" coeff_detrainment = ",coeff_detrainment
[42]	308
[2199]	309	! entrainment by slope wind scheme
	310	write(,)"call slope wind lifting parametrization"
	311	slpwind=.false. ! default value
[2304]	312	call getin_p("slpwind",slpwind)
[2199]	313	write(,)" slpwind = ",slpwind
	314
[2179]	315	! latent heat release from ground water ice sublimation/condensation
	316	write(,)"latent heat release during sublimation",
	317	& " /condensation of ground water ice"
[2218]	318	latentheat_surfwater=.true. ! default value
[2304]	319	call getin_p("latentheat_surfwater",latentheat_surfwater)
[2218]	320	write(,)" latentheat_surfwater = ",latentheat_surfwater
[2179]	321
[42]	322	write(,)"rad.transfer is computed every iradia",
	323	& " physical timestep"
	324	iradia=1 ! default value
[2304]	325	call getin_p("iradia",iradia)
[42]	326	write(,)" iradia = ",iradia
	327
	328
	329	write(,)"Output of the exchange coefficient mattrix ?",
	330	& "(for diagnostics only)"
	331	callg2d=.false. ! default value
[2304]	332	call getin_p("callg2d",callg2d)
[42]	333	write(,)" callg2d = ",callg2d
	334
	335	write(,)"Rayleigh scattering : (should be .false. for now)"
	336	rayleigh=.false.
[2304]	337	call getin_p("rayleigh",rayleigh)
[42]	338	write(,)" rayleigh = ",rayleigh
	339
	340
	341	! TRACERS:
	342
[358]	343	! dustbin
[42]	344	write(,)"Transported dust ? (if >0, use 'dustbin' dust bins)"
	345	dustbin=0 ! default value
[2304]	346	call getin_p("dustbin",dustbin)
[42]	347	write(,)" dustbin = ",dustbin
[358]	348	! active
[42]	349	write(,)"Radiatively active dust ? (matters if dustbin>0)"
	350	active=.false. ! default value
[2304]	351	call getin_p("active",active)
[42]	352	write(,)" active = ",active
	353
	354	! Test of incompatibility:
	355	! if active is used, then dustbin should be > 0
	356
	357	if (active.and.(dustbin.lt.1)) then
	358	print*,'if active is used, then dustbin should > 0'
[2304]	359	call abort_physic(modname,
	360	& "active option requires dustbin < 0",1)
[42]	361	endif
[358]	362	! doubleq
[42]	363	write(,)"use mass and number mixing ratios to predict",
	364	& " dust size ?"
	365	doubleq=.false. ! default value
[2304]	366	call getin_p("doubleq",doubleq)
[42]	367	write(,)" doubleq = ",doubleq
[358]	368	! submicron
[42]	369	submicron=.false. ! default value
[2304]	370	call getin_p("submicron",submicron)
[42]	371	write(,)" submicron = ",submicron
	372
	373	! Test of incompatibility:
	374	! if doubleq is used, then dustbin should be 2
	375
	376	if (doubleq.and.(dustbin.ne.2)) then
	377	print*,'if doubleq is used, then dustbin should be 2'
[2304]	378	call abort_physic(modname,
	379	& "doubleq option requires dustbin = 2",1)
[42]	380	endif
[1036]	381	if (doubleq.and.submicron.and.(nq.LT.3)) then
[42]	382	print*,'If doubleq is used with a submicron tracer,'
	383	print*,' then the number of tracers has to be'
	384	print*,' larger than 3.'
[2304]	385	call abort_physic(modname,
	386	& "submicron option requires dustbin > 2",1)
[42]	387	endif
[1088]	388
[358]	389	! lifting
[42]	390	write(,)"dust lifted by GCM surface winds ?"
	391	lifting=.false. ! default value
[2304]	392	call getin_p("lifting",lifting)
[42]	393	write(,)" lifting = ",lifting
	394
	395	! Test of incompatibility:
	396	! if lifting is used, then dustbin should be > 0
	397
	398	if (lifting.and.(dustbin.lt.1)) then
	399	print*,'if lifting is used, then dustbin should > 0'
[2304]	400	call abort_physic(modname,
	401	& "lifting option requires dustbin > 0",1)
[42]	402	endif
[1088]	403
[1974]	404	! dust injection scheme
	405	dustinjection=0 ! default: no injection scheme
[2304]	406	call getin_p("dustinjection",dustinjection)
[1974]	407	write(,)" dustinjection = ",dustinjection
	408	! dust injection scheme coefficient
	409	coeff_injection=1. ! default value
[2304]	410	call getin_p("coeff_injection",coeff_injection)
[1974]	411	write(,)" coeff_in,jection = ",coeff_injection
[2160]	412	! timing for dust injection
	413	ti_injection=10. ! default value
	414	tf_injection=12. ! default value
[2304]	415	call getin_p("ti_injection",ti_injection)
[2160]	416	write(,)" ti_injection = ",ti_injection
[2304]	417	call getin_p("tf_injection",tf_injection)
[2160]	418	write(,)" tf_injection = ",tf_injection
[1974]	419
[1088]	420	! free evolving dust
	421	! freedust=true just says that there is no lifting and no dust opacity scaling.
	422	write(,)"dust lifted by GCM surface winds ?"
	423	freedust=.false. ! default value
[2304]	424	call getin_p("freedust",freedust)
[1088]	425	write(,)" freedust = ",freedust
	426	if (freedust.and..not.doubleq) then
	427	print*,'freedust should be used with doubleq !'
[2304]	428	call abort_physic(modname,
	429	& "freedust option requires doubleq",1)
[1088]	430	endif
[2413]	431
	432	! dust rescaling mode (if any)
	433	if (freedust) then
[2417]	434	dustscaling_mode=0
[2413]	435	else
[2417]	436	dustscaling_mode=1 ! GCMv5.3 style
[2413]	437	endif
[2417]	438	call getin_p("dustscaling_mode",dustscaling_mode)
	439	write(,) "dustscaling_mode=",dustscaling_mode
[2413]	440
[1264]	441	#ifndef MESOSCALE
	442	! this test is valid in GCM case
	443	! ... not in mesoscale case, for we want to activate mesoscale lifting
[1974]	444	if (freedust.and.dustinjection.eq.0)then
	445	if(lifting) then
	446	print*,'if freedust is used and dustinjection = 0,
	447	& then lifting should not be used'
[2304]	448	call abort_physic(modname,
	449	& "freedust option with dustinjection = 0"//
[2411]	450	& " requires lifting to be false",1)
[1974]	451	endif
[1088]	452	endif
[1264]	453	#endif
[1974]	454	if (dustinjection.eq.1)then
	455	if(.not.lifting) then
	456	print*,"if dustinjection=1, then lifting should be true"
[2304]	457	call abort_physic(modname,
	458	& "dustinjection=1 requires lifting",1)
[1974]	459	endif
	460	if(.not.freedust) then
	461	print*,"if dustinjection=1, then freedust should be true"
[2304]	462	call abort_physic(modname,
	463	& "dustinjection=1 requires freedust",1)
[1974]	464	endif
	465	endif
[2199]	466	! rocket dust storm and entrainment by slope wind
[1974]	467	! Test of incompatibility:
[2199]	468	! if rdstorm or slpwind is used, then doubleq should be true
	469	if ((rdstorm.or.slpwind).and..not.doubleq) then
	470	print*,'if rdstorm or slpwind is used, then doubleq
	471	& should be used !'
[2304]	472	call abort_physic(modname,
	473	& "rdstorm or slpwind requires doubleq",1)
[1974]	474	endif
[2199]	475	if ((rdstorm.or.slpwind).and..not.active) then
	476	print*,'if rdstorm or slpwind is used, then active
	477	& should be used !'
[2304]	478	call abort_physic(modname,
	479	& "rdstorm or slpwind requires activ",1)
[1974]	480	endif
	481	if (rdstorm.and..not.lifting) then
[2199]	482	print*,'if rdstorm is used, then lifting
	483	& should be used !'
[2304]	484	call abort_physic(modname,
	485	& "rdstorm requires lifting",1)
[1974]	486	endif
[2199]	487	if ((rdstorm.or.slpwind).and..not.freedust) then
	488	print*,'if rdstorm or slpwind is used, then freedust
	489	& should be used !'
[2304]	490	call abort_physic(modname,
	491	& "rdstorm or slpwind requires freedust",1)
[1974]	492	endif
	493	if (rdstorm.and.(dustinjection.eq.0)) then
[2199]	494	print*,'if rdstorm is used, then dustinjection
	495	& should be used !'
[2304]	496	call abort_physic(modname,
	497	& "rdstorm requires dustinjection",1)
[1974]	498	endif
[1353]	499	! Dust IR opacity
	500	write(,)" Wavelength for infrared opacity of dust ?"
	501	write(,)" Choices are:"
	502	write(,)" tes --- > 9.3 microns [default]"
	503	write(,)" mcs --- > 21.6 microns"
	504	!
	505	! WARNING WARNING WARNING WARNING WARNING WARNING
	506	!
	507	! BEFORE ADDING A NEW VALUE, BE SURE THAT THE
	508	! CORRESPONDING WAVELENGTH IS IN THE LOOKUP TABLE,
	509	! OR AT LEAST NO TO FAR, TO AVOID FALLACIOUS INTERPOLATIONS.
	510	!
	511	dustiropacity="tes" !default value - is expected to shift to mcs one day
[2304]	512	call getin_p("dustiropacity",dustiropacity)
[1353]	513	write(,)" dustiropacity = ",trim(dustiropacity)
	514	select case (trim(dustiropacity))
	515	case ("tes")
	516	dustrefir = 9.3E-6
	517	case ("mcs")
	518	dustrefir = 21.6E-6
	519	case default
	520	write(,) trim(dustiropacity),
	521	& " is not a valid option for dustiropacity"
[2304]	522	call abort_physic(modname,
	523	& "invalid dustiropacity option value",1)
[1353]	524	end select
	525
[358]	526	! callddevil
[42]	527	write(,)" dust lifted by dust devils ?"
	528	callddevil=.false. !default value
[2304]	529	call getin_p("callddevil",callddevil)
[42]	530	write(,)" callddevil = ",callddevil
	531
	532	! Test of incompatibility:
	533	! if dustdevil is used, then dustbin should be > 0
	534
	535	if (callddevil.and.(dustbin.lt.1)) then
	536	print*,'if dustdevil is used, then dustbin should > 0'
[2304]	537	call abort_physic(modname,
	538	& "callddevil requires dustbin > 0",1)
[42]	539	endif
[358]	540	! sedimentation
[42]	541	write(,) "Gravitationnal sedimentation ?"
	542	sedimentation=.true. ! default value
[2304]	543	call getin_p("sedimentation",sedimentation)
[42]	544	write(,) " sedimentation = ",sedimentation
[358]	545	! activice
[42]	546	write(,) "Radiatively active transported atmospheric ",
	547	& "water ice ?"
	548	activice=.false. ! default value
[2304]	549	call getin_p("activice",activice)
[42]	550	write(,) " activice = ",activice
[358]	551	! water
[42]	552	write(,) "Compute water cycle ?"
	553	water=.false. ! default value
[2304]	554	call getin_p("water",water)
[42]	555	write(,) " water = ",water
[2312]	556	! hdo
	557	write(,) "Compute hdo cycle ?"
	558	hdo=.false. ! default value
	559	call getin_p("hdo",hdo)
	560	write(,) " hdo = ",hdo
[42]	561
[2312]	562	write(,) "Use fractionation for hdo?"
	563	hdofrac=.true. ! default value
	564	call getin_p("hdofrac",hdofrac)
	565	write(,) " hdofrac = ",hdofrac
	566
[2447]	567	! Activeco2ice
	568	write(,) "Radiatively active transported atmospheric ",
	569	& "Co2 ice ?"
	570	activeco2ice=.false. ! default value
	571	call getin_p("activeco2ice",activeco2ice)
	572	write(,) " activeco2ice = ",activeco2ice
[2312]	573
[1711]	574	! sub-grid cloud fraction: fixed
	575	write(,) "Fixed cloud fraction?"
	576	CLFfixval=1.0 ! default value
[2304]	577	call getin_p("CLFfixval",CLFfixval)
[1711]	578	write(,) "CLFfixval=",CLFfixval
	579	! sub-grid cloud fraction: varying
	580	write(,) "Use partial nebulosity?"
	581	CLFvarying=.false. ! default value
[2304]	582	call getin_p("CLFvarying",CLFvarying)
[1711]	583	write(,)"CLFvarying=",CLFvarying
	584
[1617]	585	!CO2 clouds scheme?
[1720]	586	write(,) "Compute CO2 clouds (implies microphysical scheme)?"
[1617]	587	co2clouds=.false. ! default value
[2304]	588	call getin_p("co2clouds",co2clouds)
[1617]	589	write(,) " co2clouds = ",co2clouds
[1720]	590	!Can water ice particles serve as CCN for CO2clouds
	591	write(,) "Use water ice as CO2 clouds CCN ?"
	592	co2useh2o=.false. ! default value
[2304]	593	call getin_p("co2useh2o",co2useh2o)
[1720]	594	write(,) " co2useh2o = ",co2useh2o
	595	!Do we allow a supply of meteoritic paricles to serve as CO2 ice CCN?
	596	write(,) "Supply meteoritic particle for CO2 clouds ?"
	597	meteo_flux=.false. !Default value
[2304]	598	call getin_p("meteo_flux",meteo_flux)
[1720]	599	write(,) " meteo_flux = ",meteo_flux
	600	!Do we allow a sub-grid temperature distribution for the CO2 microphysics
	601	write(,) "sub-grid temperature distribution for CO2 clouds?"
	602	CLFvaryingCO2=.false. !Default value
[2304]	603	call getin_p("CLFvaryingCO2",CLFvaryingCO2)
[1720]	604	write(,) " CLFvaryingCO2 = ",CLFvaryingCO2
	605	!Amplitude of the sub-grid temperature distribution for the CO2 microphysics
	606	write(,) "sub-grid temperature amplitude for CO2 clouds?"
	607	spantCO2=0 !Default value
[2304]	608	call getin_p("spantCO2",spantCO2)
[1720]	609	write(,) " spantCO2 = ",spantCO2
[1818]	610	!Do you want to filter the sub-grid T distribution by a Saturation index?
	611	write(,) "filter sub-grid temperature by Saturation index?"
	612	satindexco2=.true.
[2304]	613	call getin_p("satindexco2",satindexco2)
[1818]	614	write(,) " satindexco2 = ",satindexco2
[1720]	615
[1818]	616
[833]	617	! thermal inertia feedback
	618	write(,) "Activate the thermal inertia feedback ?"
	619	tifeedback=.false. ! default value
[2304]	620	call getin_p("tifeedback",tifeedback)
[833]	621	write(,) " tifeedback = ",tifeedback
	622
[42]	623	! Test of incompatibility:
	624
[833]	625	if (tifeedback.and..not.water) then
	626	print*,'if tifeedback is used,'
	627	print*,'water should be used too'
[2304]	628	call abort_physic(modname,
	629	& "tifeedback requires water",1)
[833]	630	endif
	631
	632	if (tifeedback.and..not.callsoil) then
	633	print*,'if tifeedback is used,'
	634	print*,'callsoil should be used too'
[2304]	635	call abort_physic(modname,
	636	& "tifeedback requires callsoil",1)
[833]	637	endif
	638
[42]	639	if (activice.and..not.water) then
	640	print*,'if activice is used, water should be used too'
[2304]	641	call abort_physic(modname,
	642	& "activeice requires water",1)
[42]	643	endif
	644
	645	if (water.and..not.tracer) then
	646	print*,'if water is used, tracer should be used too'
[2304]	647	call abort_physic(modname,
	648	& "water requires tracer",1)
[42]	649	endif
[2312]	650
	651	if (hdo.and..not.water) then
	652	print*,'if hdo is used, water should be used too'
[2398]	653	call abort_physic(modname,
	654	& "hd2 requires tracer",1)
[2312]	655	endif
	656
[420]	657
[2447]	658	if (activeco2ice.and..not.co2clouds) then
	659	print*,'if activeco2ice is used, co2clouds should be used too'
	660	call abort_physic(modname,
	661	& "activeco2ice requires co2clouds",1)
	662	endif
	663
[455]	664	! water ice clouds effective variance distribution for sedimentaion
[1617]	665	write(,) "Sed effective variance for water ice clouds ?"
[455]	666	nuice_sed=0.45
[2304]	667	call getin_p("nuice_sed",nuice_sed)
[455]	668	write(,) "water_param nueff Sedimentation:", nuice_sed
[1617]	669
	670	write(,) "Sed effective variance for CO2 clouds ?"
	671	nuiceco2_sed=0.45
[2304]	672	call getin_p("nuiceco2_sed",nuiceco2_sed)
[1617]	673	write(,) "CO2 nueff Sedimentation:", nuiceco2_sed
	674
	675	write(,) "REF effective variance for CO2 clouds ?"
	676	nuiceco2_ref=0.45
[2304]	677	call getin_p("nuiceco2_ref",nuiceco2_ref)
[1617]	678	write(,) "CO2 nueff Sedimentation:", nuiceco2_ref
	679
	680	write(,) "REF effective variance for water clouds ?"
	681	nuice_ref=0.45
[2304]	682	call getin_p("nuice_ref",nuice_ref)
[1617]	683	write(,) "CO2 nueff Sedimentation:", nuice_ref
	684
	685
[420]	686	! ccn factor if no scavenging
[1617]	687	write(,) "water param CCN reduc. factor ?"
[420]	688	ccn_factor = 4.5
[2304]	689	call getin_p("ccn_factor",ccn_factor)
[420]	690	write(,)" ccn_factor = ",ccn_factor
	691	write(,)"Careful: only used when microphys=F, otherwise"
	692	write(,)"the contact parameter is used instead;"
	693
[1720]	694	! microphys
	695	write(,)"Microphysical scheme for water-ice clouds?"
	696	microphys=.false. ! default value
[2304]	697	call getin_p("microphys",microphys)
[1720]	698	write(,)" microphys = ",microphys
[1655]	699
[1720]	700	! supersat
	701	write(,)"Allow super-saturation of water vapor?"
	702	supersat=.true. ! default value
[2304]	703	call getin_p("supersat",supersat)
[1720]	704	write(,)"supersat = ",supersat
[42]	705
[358]	706	! microphysical parameter contact
[1720]	707	write(,) "water contact parameter ?"
	708	mteta = 0.95
[2304]	709	call getin_p("mteta",mteta)
[1720]	710	write(,) "mteta = ", mteta
	711
[358]	712	! scavenging
[1720]	713	write(,)"Dust scavenging by H2O/CO2 snowfall ?"
	714	scavenging=.false. ! default value
[2304]	715	call getin_p("scavenging",scavenging)
[1720]	716	write(,)" scavenging = ",scavenging
[358]	717
	718
[42]	719	! Test of incompatibility:
[358]	720	! if scavenging is used, then dustbin should be > 0
[42]	721
[1720]	722	if ((microphys.and..not.doubleq).or.
[358]	723	& (microphys.and..not.water)) then
[1720]	724	print*,'if microphys is used, then doubleq,'
	725	print*,'and water must be used!'
[2304]	726	call abort_physic(modname,
	727	& "microphys requires water and doubleq",1)
[1720]	728	endif
	729	if (microphys.and..not.scavenging) then
	730	print*,''
	731	print*,'----------------WARNING-----------------'
	732	print*,'microphys is used without scavenging !!!'
	733	print*,'----------------WARNING-----------------'
	734	print*,''
	735	endif
	736
	737	if ((scavenging.and..not.microphys).or.
[1617]	738	& (scavenging.and.(dustbin.lt.1)))then
[1720]	739	print*,'if scavenging is used, then microphys'
	740	print*,'must be used!'
[2304]	741	call abort_physic(modname,
	742	& "scavenging requires microphys",1)
[1720]	743	endif
[740]	744
[2184]	745	! Instantaneous scavenging by CO2
	746	! -> expected to be replaced by scavenging with microphysics (flag scavenging) one day
	747	write(,)"Dust scavenging by instantaneous CO2 snowfall ?"
	748	scavco2cond=.false. ! default value
[2304]	749	call getin_p("scavco2cond",scavco2cond)
[2184]	750	write(,)" scavco2cond = ",scavco2cond
[358]	751	! Test of incompatibility:
[2184]	752	! if scavco2cond is used, then dustbin should be > 0
	753	if (scavco2cond.and.(dustbin.lt.1))then
	754	print*,'if scavco2cond is used, then dustbin should be > 0'
[2304]	755	call abort_physic(modname,
	756	& "scavco2cond requires dustbin > 0",1)
[2184]	757	endif
	758	! if co2clouds is used, then there is no need for scavco2cond
	759	if (co2clouds.and.scavco2cond) then
	760	print*,''
	761	print*,'----------------WARNING-----------------'
	762	print*,' microphys scavenging is used so '
	763	print*,' no need for scavco2cond !!! '
	764	print*,'----------------WARNING-----------------'
	765	print*,''
[2304]	766	call abort_physic(modname,
	767	& "incompatible co2cloud and scavco2cond options",1)
[2184]	768	endif
	769
	770	! Test of incompatibility:
[358]	771
[42]	772	write(,) "Permanent water caps at poles ?",
	773	& " .true. is RECOMMENDED"
	774	write(,) "(with .true., North cap is a source of water ",
	775	& "and South pole is a cold trap)"
	776	caps=.true. ! default value
[2304]	777	call getin_p("caps",caps)
[42]	778	write(,) " caps = ",caps
	779
[358]	780	! albedo_h2o_ice
[283]	781	write(,) "water ice albedo ?"
	782	albedo_h2o_ice=0.45
[2304]	783	call getin_p("albedo_h2o_ice",albedo_h2o_ice)
[283]	784	write(,) " albedo_h2o_ice = ",albedo_h2o_ice
[358]	785	! inert_h2o_ice
[283]	786	write(,) "water ice thermal inertia ?"
	787	inert_h2o_ice=2400 ! (J.m^-2.K^-1.s^-1/2)
[2304]	788	call getin_p("inert_h2o_ice",inert_h2o_ice)
[283]	789	write(,) " inert_h2o_ice = ",inert_h2o_ice
[358]	790	! frost_albedo_threshold
[283]	791	write(,) "frost thickness threshold for albedo ?"
[285]	792	frost_albedo_threshold=0.005 ! 5.4 mic (i.e 0.005 kg.m-2)
[2304]	793	call getin_p("frost_albedo_threshold",
[283]	794	& frost_albedo_threshold)
	795	write(,) " frost_albedo_threshold = ",
	796	& frost_albedo_threshold
	797
[485]	798	! call Titus crocus line -- DEFAULT IS NONE
	799	write(,) "Titus crocus line ?"
	800	tituscap=.false. ! default value
[2304]	801	call getin_p("tituscap",tituscap)
[485]	802	write(,) "tituscap",tituscap
[633]	803
[2164]	804	! Chemistry:
[42]	805	write(,) "photochemistry: include chemical species"
	806	photochem=.false. ! default value
[2304]	807	call getin_p("photochem",photochem)
[42]	808	write(,) " photochem = ",photochem
[2164]	809
	810	write(,) "Compute chemistry (if photochem is .true.)",
	811	& "every ichemistry physics step (default: ichemistry=1)"
	812	ichemistry=1
[2304]	813	call getin_p("ichemistry",ichemistry)
[2164]	814	write(,) " ichemistry = ",ichemistry
[42]	815
	816
[1246]	817	! SCATTERERS
	818	write(,) "how many scatterers?"
	819	naerkind=1 ! default value
[2304]	820	call getin_p("naerkind",naerkind)
[1246]	821	write(,)" naerkind = ",naerkind
	822
	823	! Test of incompatibility
	824	c Logical tests for radiatively active water-ice clouds:
	825	IF ( (activice.AND.(.NOT.water)).OR.
	826	& (activice.AND.(naerkind.LT.2)) ) THEN
	827	WRITE(,) 'If activice is TRUE, water has to be set'
	828	WRITE(,) 'to TRUE, and "naerkind" must be at least'
	829	WRITE(,) 'equal to 2.'
[2304]	830	call abort_physic(modname,
	831	& "radiatively active dust and water"//
	832	& " require naerkind > 1",1)
[1246]	833	ENDIF
	834
	835	!------------------------------------------
	836	!------------------------------------------
	837	! once naerkind is known allocate arrays
	838	! -- we do it here and not in phys_var_init
	839	! -- because we need to know naerkind first
	840	CALL ini_scatterers(ngrid,nlayer)
	841	!------------------------------------------
	842	!------------------------------------------
	843
	844
	845	c Please name the different scatterers here ----------------
	846	name_iaer(1) = "dust_conrath" !! default choice is good old Conrath profile
	847	IF (doubleq.AND.active) name_iaer(1) = "dust_doubleq" !! two-moment scheme
[2447]	848
[1246]	849	if (nq.gt.1) then
[2447]	850	! trick to avoid problems compiling with 1 tracer
	851	! and picky compilers who know name_iaer(2) is out of bounds
	852	j=2
	853	IF (rdstorm.AND..NOT.activice.AND..NOT.slpwind) then
	854	name_iaer(j) = "stormdust_doubleq" !! storm dust two-moment scheme
	855	j = j+1
	856	END IF
	857
	858	IF (rdstorm.AND.water.AND.activice.AND..NOT.slpwind) then
	859	name_iaer(j) = "stormdust_doubleq"
	860	j = j+1
	861	END IF
	862
	863	IF (slpwind.AND..NOT.activice.AND..NOT.rdstorm) then
	864	name_iaer(j) = "topdust_doubleq" !! storm dust two-moment scheme
	865	j = j+1
	866	END IF
	867
	868	IF (slpwind.AND.water.AND.activice.AND..NOT.rdstorm) then
	869	name_iaer(j) = "topdust_doubleq"
	870	j = j+1
	871	END IF
	872
	873	IF (rdstorm.AND.slpwind.AND..NOT.activice) THEN
	874	name_iaer(j) = "stormdust_doubleq"
	875	name_iaer(j+1) = "topdust_doubleq"
	876	j = j+2
	877	ENDIF
	878
	879	IF (rdstorm.AND.slpwind.AND.water.AND.activice) THEN
	880	name_iaer(j) = "stormdust_doubleq"
	881	name_iaer(j+1) = "topdust_doubleq"
	882	j = j+2
	883	ENDIF
	884
	885	IF (water.AND.activice) then
	886	name_iaer(j) = "h2o_ice" !! radiatively-active clouds
	887	j = j+1
	888	END IF
	889
	890	IF (co2clouds.AND.activeco2ice) then
	891	name_iaer(j) = "co2_ice" !! radiatively-active co2 clouds
	892	j = j+1
	893	ENDIF
	894
	895	IF (submicron.AND.active) then
	896	name_iaer(j) = "dust_submicron" !! JBM experimental stuff
	897	j = j+1
	898	ENDIF
[1246]	899	endif ! of if (nq.gt.1)
	900	c ----------------------------------------------------------
	901
[42]	902	! THERMOSPHERE
	903
	904	write(,) "call thermosphere ?"
	905	callthermos=.false. ! default value
[2304]	906	call getin_p("callthermos",callthermos)
[42]	907	write(,) " callthermos = ",callthermos
	908
	909
	910	write(,) " water included without cycle ",
	911	& "(only if water=.false.)"
	912	thermoswater=.false. ! default value
[2304]	913	call getin_p("thermoswater",thermoswater)
[42]	914	write(,) " thermoswater = ",thermoswater
	915
	916	write(,) "call thermal conduction ?",
	917	& " (only if callthermos=.true.)"
	918	callconduct=.false. ! default value
[2304]	919	call getin_p("callconduct",callconduct)
[42]	920	write(,) " callconduct = ",callconduct
	921
	922	write(,) "call EUV heating ?",
	923	& " (only if callthermos=.true.)"
	924	calleuv=.false. ! default value
[2304]	925	call getin_p("calleuv",calleuv)
[42]	926	write(,) " calleuv = ",calleuv
	927
	928	write(,) "call molecular viscosity ?",
	929	& " (only if callthermos=.true.)"
	930	callmolvis=.false. ! default value
[2304]	931	call getin_p("callmolvis",callmolvis)
[42]	932	write(,) " callmolvis = ",callmolvis
	933
	934	write(,) "call molecular diffusion ?",
	935	& " (only if callthermos=.true.)"
	936	callmoldiff=.false. ! default value
[2304]	937	call getin_p("callmoldiff",callmoldiff)
[42]	938	write(,) " callmoldiff = ",callmoldiff
	939
	940
	941	write(,) "call thermospheric photochemistry ?",
	942	& " (only if callthermos=.true.)"
	943	thermochem=.false. ! default value
[2304]	944	call getin_p("thermochem",thermochem)
[42]	945	write(,) " thermochem = ",thermochem
	946
[705]	947	write(,) "Method to include solar variability"
[1684]	948	write(,) "0-> fixed value of E10.7 (fixed_euv_value); ",
	949	& "1-> daily evolution of E10.7 (for given solvaryear)"
[705]	950	solvarmod=1
[2304]	951	call getin_p("solvarmod",solvarmod)
[705]	952	write(,) " solvarmod = ",solvarmod
	953
[1684]	954	write(,) "Fixed euv (for solvarmod==0) 10.7 value?"
	955	write(,) " (min=80 , ave=140, max=320)"
	956	fixed_euv_value=140 ! default value
[2304]	957	call getin_p("fixed_euv_value",fixed_euv_value)
[1684]	958	write(,) " fixed_euv_value = ",fixed_euv_value
[42]	959
[705]	960	write(,) "Solar variability as observed for MY: "
	961	write(,) "Only if solvarmod=1"
	962	solvaryear=24
[2304]	963	call getin_p("solvaryear",solvaryear)
[705]	964	write(,) " solvaryear = ",solvaryear
	965
[552]	966	write(,) "UV heating efficiency:",
	967	& "measured values between 0.19 and 0.23 (Fox et al. 1996)",
	968	& "lower values may be used to compensate low 15 um cooling"
	969	euveff=0.21 !default value
[2304]	970	call getin_p("euveff",euveff)
[552]	971	write(,) " euveff = ", euveff
[42]	972
[1651]	973
[42]	974	if (.not.callthermos) then
	975	if (thermoswater) then
	976	print*,'if thermoswater is set, callthermos must be true'
[2304]	977	call abort_physic(modname,
	978	& "thermoswater requires callthermos",1)
[42]	979	endif
	980	if (callconduct) then
	981	print*,'if callconduct is set, callthermos must be true'
[2304]	982	call abort_physic(modname,
	983	& "callconduct requires callthermos",1)
[42]	984	endif
	985	if (calleuv) then
	986	print*,'if calleuv is set, callthermos must be true'
[2304]	987	call abort_physic(modname,
	988	& "calleuv requires callthermos",1)
[42]	989	endif
	990	if (callmolvis) then
	991	print*,'if callmolvis is set, callthermos must be true'
[2304]	992	call abort_physic(modname,
	993	& "callmolvis requires callthermos",1)
[42]	994	endif
	995	if (callmoldiff) then
	996	print*,'if callmoldiff is set, callthermos must be true'
[2304]	997	call abort_physic(modname,
	998	& "callmoldiff requires callthermos",1)
[42]	999	endif
	1000	if (thermochem) then
	1001	print*,'if thermochem is set, callthermos must be true'
[2304]	1002	call abort_physic(modname,
	1003	& "thermochem requires callthermos",1)
[42]	1004	endif
	1005	endif
	1006
	1007	! Test of incompatibility:
	1008	! if photochem is used, then water should be used too
	1009
	1010	if (photochem.and..not.water) then
	1011	print*,'if photochem is used, water should be used too'
[2304]	1012	call abort_physic(modname,
	1013	& "photochem requires water",1)
[42]	1014	endif
	1015
	1016	! if callthermos is used, then thermoswater should be used too
	1017	! (if water not used already)
	1018
	1019	if (callthermos .and. .not.water) then
	1020	if (callthermos .and. .not.thermoswater) then
	1021	print*,'if callthermos is used, water or thermoswater
	1022	& should be used too'
[2304]	1023	call abort_physic(modname,
	1024	& "callthermos requires water or thermoswater",1)
[42]	1025	endif
	1026	endif
	1027
	1028	PRINT*,'--------------------------------------------'
	1029	PRINT*
	1030	PRINT*
	1031	ELSE
	1032	write(,)
	1033	write(,) 'Cannot read file callphys.def. Is it here ?'
[2304]	1034	call abort_physic(modname,
	1035	& "missing callphys.def file",1)
[42]	1036	ENDIF
	1037
	1038	8000 FORMAT(t5,a12,l8)
	1039	8001 FORMAT(t5,a12,i8)
	1040
	1041	PRINT*
[1224]	1042	PRINT*,'conf_phys: daysec',daysec
[42]	1043	PRINT*
[1224]	1044	PRINT*,'conf_phys: The radiative transfer is computed:'
[42]	1045	PRINT*,' each ',iradia,' physical time-step'
	1046	PRINT,' or each ',iradiadtphys,' seconds'
	1047	PRINT*
	1048	! --------------------------------------------------------------
	1049	! Managing the Longwave radiative transfer
	1050	! --------------------------------------------------------------
	1051
	1052	! In most cases, the run just use the following values :
	1053	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
	1054	callemis=.true.
	1055	! ilwd=10*int(daysec/dtphys) ! bug before 22/10/01
	1056	ilwd=1
	1057	ilwn=1 !2
	1058	ilwb=1 !2
	1059	linear=.true.
	1060	ncouche=3
	1061	alphan=0.4
	1062	semi=0
	1063
	1064	! BUT people working hard on the LW may want to read them in 'radia.def'
	1065	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
	1066
	1067	OPEN(99,file='radia.def',status='old',form='formatted'
	1068	. ,iostat=ierr)
	1069	IF(ierr.EQ.0) THEN
[1224]	1070	write(,) 'conf_phys: Reading radia.def !!!'
[42]	1071	READ(99,fmt='(a)') ch1
	1072	READ(99,*) callemis
	1073	WRITE(*,8000) ch1,callemis
	1074
	1075	READ(99,fmt='(a)') ch1
	1076	READ(99,*) iradia
	1077	WRITE(*,8001) ch1,iradia
	1078
	1079	READ(99,fmt='(a)') ch1
	1080	READ(99,*) ilwd
	1081	WRITE(*,8001) ch1,ilwd
	1082
	1083	READ(99,fmt='(a)') ch1
	1084	READ(99,*) ilwn
	1085	WRITE(*,8001) ch1,ilwn
	1086
	1087	READ(99,fmt='(a)') ch1
	1088	READ(99,*) linear
	1089	WRITE(*,8000) ch1,linear
	1090
	1091	READ(99,fmt='(a)') ch1
	1092	READ(99,*) ncouche
	1093	WRITE(*,8001) ch1,ncouche
	1094
	1095	READ(99,fmt='(a)') ch1
	1096	READ(99,*) alphan
	1097	WRITE(,) ch1,alphan
	1098
	1099	READ(99,fmt='(a)') ch1
	1100	READ(99,*) ilwb
	1101	WRITE(*,8001) ch1,ilwb
	1102
	1103
	1104	READ(99,fmt='(a)') ch1
	1105	READ(99,'(l1)') callg2d
	1106	WRITE(*,8000) ch1,callg2d
	1107
	1108	READ(99,fmt='(a)') ch1
	1109	READ(99,*) semi
	1110	WRITE(,) ch1,semi
	1111	end if
	1112	CLOSE(99)
	1113
	1114	END

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

Download in other formats:

Original Format