Context Navigation

phys_cosp.F90

Last change on this file was 3435, checked in by Laurent Fairhead, 6 years ago

"Historic" :-) commit merging the physics branch used for DYNAMICO with the LMDZ trunk.
The same physics branch can now be used seamlessly with the traditional lon-lat LMDZ
dynamical core and DYNAMICO.
Testing consisted in running a lon-lat LMDZ bucket simulation with the NPv6.1 physics package
with the original trunk sources and the merged sources. Tests were succesful in the sense that
numeric continuity was preserved in the restart files from both simulation. Further tests
included running both versions of the physics codes for one year in a LMDZOR setting in which
the restart files also came out identical.

Caution:

as the physics package now manages unstructured grids, grid information needs to be transmitted

to the surface scheme ORCHIDEE. This means that the interface defined in surf_land_orchidee_mod.F90
is only compatible with ORCHIDEE version orchidee2.1 and later versions. If previous versions of
ORCHIDEE need to be used, the CPP key ORCHIDEE_NOUNSTRUCT needs to be set at compilation time.
This is done automatically if makelmdz/makelmdz_fcm are called with the veget orchidee2.0 switch

due to a limitation in XIOS, the time at which limit conditions will be read in by DYNAMICO will be

delayed by one physic timestep with respect to the time it is read in by the lon-lat model. This is caused
by the line

IF (MOD(itime-1, lmt_pas) == 0 .OR. (jour_lu /= jour .AND. grid_type /= unstructured)) THEN ! time to read

in limit_read_mod.F90

Work still needed on COSP integration and XML files for DYNAMICO

EM, YM, LF

Property copyright set to
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory

File size: 19.2 KB

Line
1	! Simulateur COSP : Cfmip Observation Simulator Package
2
3	! ISCCP, Radar (QuickBeam), Lidar et Parasol (ACTSIM), MISR, RTTOVS
4	!Idelkadi Abderrahmane Aout-Septembre 2009 First Version
5	!Idelkadi Abderrahmane Nov 2015 version v1.4.0
6
7	subroutine phys_cosp( itap,dtime,freq_cosp, &
8	ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP, &
9	ecrit_mth,ecrit_day,ecrit_hf, ok_all_xml, missing_val, &
10	Nptslmdz,Nlevlmdz,lon,lat, presnivs,overlaplmdz,sunlit, &
11	ref_liq,ref_ice,fracTerLic,u_wind,v_wind,phis,phi,ph,p,skt,t, &
12	sh,rh,tca,cca,mr_lsliq,mr_lsice,fl_lsrainI,fl_lssnowI, &
13	fl_ccrainI,fl_ccsnowI,mr_ozone,dtau_s,dem_s)
14
15	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
16	!!!! Inputs :
17	! itap, !Increment de la physiq
18	! dtime, !Pas de temps physiq
19	! overlap, !Overlap type in SCOPS
20	! Npoints, !Nb de points de la grille physiq
21	! Nlevels, !Nb de niveaux verticaux
22	! Ncolumns, !Number of subcolumns
23	! lon,lat, !Longitudes et latitudes de la grille LMDZ
24	! ref_liq,ref_ice, !Rayons effectifs des particules liq et ice (en microm)
25	! fracTerLic, !Fraction terre a convertir en masque
26	! u_wind,v_wind, !Vents a 10m ???
27	! phi, !Geopotentiel
28	! phis, !Geopotentiel sol
29	! ph, !pression pour chaque inter-couche
30	! p, !Pression aux milieux des couches
31	! skt,t, !Temp au sol et temp 3D
32	! sh, !Humidite specifique
33	! rh, !Humidite relatif
34	! tca, !Fraction nuageuse
35	! cca !Fraction nuageuse convective
36	! mr_lsliq, !Liq Cloud water content
37	! mr_lsice, !Ice Cloud water content
38	! mr_ccliq, !Convective Cloud Liquid water content
39	! mr_ccice, !Cloud ice water content
40	! fl_lsrain, !Large scale precipitation lic
41	! fl_lssnow, !Large scale precipitation ice
42	! fl_ccrain, !Convective precipitation lic
43	! fl_ccsnow, !Convective precipitation ice
44	! mr_ozone, !Concentration ozone (Kg/Kg)
45	! dem_s !Cloud optical emissivity
46	! dtau_s !Cloud optical thickness
47	! emsfc_lw = 1. !Surface emissivity dans radlwsw.F90
48
49	!!! Outputs :
50	! calipso2D, !Lidar Low/heigh/Mean/Total-level Cloud Fraction
51	! calipso3D, !Lidar Cloud Fraction (532 nm)
52	! cfadlidar, !Lidar Scattering Ratio CFAD (532 nm)
53	! parasolrefl, !PARASOL-like mono-directional reflectance
54	! atb, !Lidar Attenuated Total Backscatter (532 nm)
55	! betamol, !Lidar Molecular Backscatter (532 nm)
56	! cfaddbze, !Radar Reflectivity Factor CFAD (94 GHz)
57	! clcalipso2, !Cloud frequency of occurrence as seen by CALIPSO but not CloudSat
58	! dbze, !Efective_reflectivity_factor
59	! cltlidarradar, !Lidar and Radar Total Cloud Fraction
60	! clMISR, !Cloud Fraction as Calculated by the MISR Simulator
61	! clisccp2, !Cloud Fraction as Calculated by the ISCCP Simulator
62	! boxtauisccp, !Optical Depth in Each Column as Calculated by the ISCCP Simulator
63	! boxptopisccp, !Cloud Top Pressure in Each Column as Calculated by the ISCCP Simulator
64	! tclisccp, !Total Cloud Fraction as Calculated by the ISCCP Simulator
65	! ctpisccp, !Mean Cloud Top Pressure as Calculated by the ISCCP Simulator
66	! tauisccp, !Mean Optical Depth as Calculated by the ISCCP Simulator
67	! albisccp, !Mean Cloud Albedo as Calculated by the ISCCP Simulator
68	! meantbisccp, !Mean all-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator
69	! meantbclrisccp !Mean clear-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator
70
71	!!! AI rajouter les nouvelles sorties
72	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
73
74	!! AI rajouter
75	#include "cosp_defs.h"
76	USE MOD_COSP_CONSTANTS
77	USE MOD_COSP_TYPES
78	USE MOD_COSP
79	USE mod_phys_lmdz_para
80	USE mod_grid_phy_lmdz
81	use ioipsl
82	use iophy
83	use cosp_output_mod
84	use cosp_output_write_mod
85	! use MOD_COSP_Modis_Simulator, only : cosp_modis
86	#ifdef CPP_XIOS
87	USE xios, ONLY: xios_field_is_active
88	#endif
89	use cosp_read_otputkeys
90
91	IMPLICIT NONE
92
93	! Local variables
94	character(len=64),PARAMETER :: cosp_input_nl='cosp_input_nl.txt'
95	character(len=64),PARAMETER :: cosp_output_nl='cosp_output_nl.txt'
96	integer, save :: isccp_topheight,isccp_topheight_direction,overlap
97	integer,save :: Ncolumns ! Number of subcolumns in SCOPS
98	integer, save :: Npoints ! Number of gridpoints
99	!$OMP THREADPRIVATE(Npoints)
100	integer, save :: Nlevels ! Number of levels
101	Integer :: Nptslmdz,Nlevlmdz ! Nb de points issus de physiq.F
102	integer, save :: Nlr ! Number of levels in statistical outputs
103	integer, save :: Npoints_it ! Max number of gridpoints to be processed in one iteration
104	integer :: i
105	type(cosp_config),save :: cfg ! Configuration options
106	!$OMP THREADPRIVATE(cfg)
107	type(cosp_gridbox) :: gbx ! Gridbox information. Input for COSP
108	type(cosp_subgrid) :: sgx ! Subgrid outputs
109	type(cosp_sgradar) :: sgradar ! Output from radar simulator
110	type(cosp_sglidar) :: sglidar ! Output from lidar simulator
111	type(cosp_isccp) :: isccp ! Output from ISCCP simulator
112	!! AI rajout modis
113	type(cosp_modis) :: modis ! Output from MODIS simulator
114	!!
115	type(cosp_misr) :: misr ! Output from MISR simulator
116	!! AI rajout rttovs
117	! type(cosp_rttov) :: rttov ! Output from RTTOV
118	!!
119	type(cosp_vgrid) :: vgrid ! Information on vertical grid of stats
120	type(cosp_radarstats) :: stradar ! Summary statistics from radar simulator
121	type(cosp_lidarstats) :: stlidar ! Summary statistics from lidar simulator
122
123	integer :: t0,t1,count_rate,count_max
124	integer :: Nlon,Nlat
125	real,save :: radar_freq,k2,ZenAng,co2,ch4,n2o,co,emsfc_lw
126	!$OMP THREADPRIVATE(emsfc_lw)
127	integer,dimension(RTTOV_MAX_CHANNELS),save :: Channels
128	real,dimension(RTTOV_MAX_CHANNELS),save :: Surfem
129	integer, save :: surface_radar,use_mie_tables,use_gas_abs,do_ray,melt_lay
130	integer, save :: Nprmts_max_hydro,Naero,Nprmts_max_aero,lidar_ice_type
131	integer, save :: platform,satellite,Instrument,Nchannels
132	logical, save :: use_vgrid,csat_vgrid,use_precipitation_fluxes,use_reff
133
134	! Declaration necessaires pour les sorties IOIPSL
135	integer :: ii
136	real :: ecrit_day,ecrit_hf,ecrit_mth, missing_val
137	logical :: ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP, ok_all_xml
138
139	logical, save :: debut_cosp=.true.
140	!$OMP THREADPRIVATE(debut_cosp)
141
142	logical, save :: first_write=.true.
143	!$OMP THREADPRIVATE(first_write)
144
145	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Input variables from LMDZ-GCM
146	integer :: overlaplmdz ! overlap type: 1=max, 2=rand, 3=max/rand ! cosp input (output lmdz)
147	real,dimension(Nptslmdz,Nlevlmdz) :: height,phi,p,ph,T,sh,rh,tca,cca,mr_lsliq,mr_lsice,mr_ccliq,mr_ccice, &
148	fl_lsrain,fl_lssnow,fl_ccrain,fl_ccsnow,fl_lsgrpl, &
149	zlev,zlev_half,mr_ozone,radliq,radice,dtau_s,dem_s,ref_liq,ref_ice
150	real,dimension(Nptslmdz,Nlevlmdz) :: fl_lsrainI,fl_lssnowI,fl_ccrainI,fl_ccsnowI
151	real,dimension(Nptslmdz) :: lon,lat,skt,fracTerLic,u_wind,v_wind,phis,sunlit
152	real,dimension(Nlevlmdz) :: presnivs
153	integer :: itap,k,ip
154	real :: dtime,freq_cosp
155	real,dimension(2) :: time_bnds
156
157	double precision :: d_dtime
158	double precision,dimension(2) :: d_time_bnds
159
160	real,dimension(2,SR_BINS) :: sratio_bounds
161	real,dimension(SR_BINS) :: sratio_ax
162
163	namelist/COSP_INPUT/overlap,isccp_topheight,isccp_topheight_direction, &
164	npoints_it,ncolumns,use_vgrid,nlr,csat_vgrid, &
165	radar_freq,surface_radar,use_mie_tables, &
166	use_gas_abs,do_ray,melt_lay,k2,Nprmts_max_hydro,Naero,Nprmts_max_aero, &
167	lidar_ice_type,use_precipitation_fluxes,use_reff, &
168	platform,satellite,Instrument,Nchannels, &
169	Channels,Surfem,ZenAng,co2,ch4,n2o,co
170
171	!---------------- End of declaration of variables --------------
172
173	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
174	! Read namelist with COSP inputs
175	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
176
177	if (debut_cosp) then
178	NPoints=Nptslmdz
179	Nlevels=Nlevlmdz
180
181	! Lecture du namelist input
182	CALL read_cosp_input
183
184	! Clefs Outputs initialisation
185	#ifdef CPP_XIOS
186	call cosp_outputkeys_init(cfg)
187	#else
188	call read_cosp_output_nl(itap,cosp_output_nl,cfg)
189	#endif
190	!!! call cosp_outputkeys_test(cfg)
191	print*,' Cles des differents simulateurs cosp a itap :',itap
192	print*,'Lradar_sim,Llidar_sim,Lisccp_sim,Lmisr_sim,Lmodis_sim,Lrttov_sim,Lstats', &
193	cfg%Lradar_sim,cfg%Llidar_sim,cfg%Lisccp_sim,cfg%Lmisr_sim,cfg%Lmodis_sim, &
194	cfg%Lrttov_sim,cfg%Lstats
195
196	if (overlaplmdz.ne.overlap) then
197	print*,'Attention overlaplmdz different de overlap lu dans namelist '
198	endif
199	print*,'Fin lecture Namelists, debut_cosp =',debut_cosp
200
201	endif ! debut_cosp
202
203	!!! Ici on modifie les cles logiques pour les outputs selon les champs actives dans les .xml
204	if ((itap.gt.1).and.(first_write))then
205	#ifdef CPP_XIOS
206	call read_xiosfieldactive(cfg)
207	#endif
208	first_write=.false.
209
210	print*,' Cles des differents simulateurs cosp a itap :',itap
211	print*,'Lradar_sim,Llidar_sim,Lisccp_sim,Lmisr_sim,Lmodis_sim,Lrttov_sim,Lstats', &
212	cfg%Lradar_sim,cfg%Llidar_sim,cfg%Lisccp_sim,cfg%Lmisr_sim,cfg%Lmodis_sim, &
213	cfg%Lrttov_sim,cfg%Lstats
214	endif
215
216	time_bnds(1) = dtime-dtime/2.
217	time_bnds(2) = dtime+dtime/2.
218
219	d_time_bnds=time_bnds
220	d_dtime=dtime
221
222	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
223	! Allocate memory for gridbox type
224	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
225	! AI mars 2017
226	! print *, 'Allocating memory for gridbox type...'
227
228	! Surafce emissivity
229	emsfc_lw = 1.
230
231	call construct_cosp_gridbox(d_dtime,d_time_bnds,radar_freq,surface_radar,use_mie_tables,use_gas_abs, &
232	do_ray,melt_lay,k2, &
233	Npoints,Nlevels,Ncolumns,N_HYDRO,Nprmts_max_hydro,Naero,Nprmts_max_aero,Npoints_it, &
234	lidar_ice_type,isccp_topheight,isccp_topheight_direction,overlap,emsfc_lw, &
235	use_precipitation_fluxes,use_reff, &
236	Platform,Satellite,Instrument,Nchannels,ZenAng, &
237	channels(1:Nchannels),surfem(1:Nchannels),co2,ch4,n2o,co,gbx)
238
239	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
240	! Here code to populate input structure
241	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
242
243	! print *, 'Populating input structure...'
244	gbx%longitude = lon
245	gbx%latitude = lat
246
247	gbx%p = p !
248	gbx%ph = ph
249	gbx%zlev = phi/9.81
250
251	zlev_half(:,1) = phis(:)/9.81
252	do k = 2, Nlevels
253	do ip = 1, Npoints
254	zlev_half(ip,k) = phi(ip,k)/9.81 + &
255	(phi(ip,k)-phi(ip,k-1))/9.81 * (ph(ip,k)-p(ip,k)) / (p(ip,k)-p(ip,k-1))
256	enddo
257	enddo
258	gbx%zlev_half = zlev_half
259
260	gbx%T = T
261	gbx%q = rh*100.
262	gbx%sh = sh
263	! On ne veut pas que cosp distingue les nuages stratiformes et convectifs
264	! on passe les contenus totaux (conv+strat)
265	gbx%cca = 0. !convective_cloud_amount (1)
266	gbx%tca = tca ! total_cloud_amount (1)
267	gbx%psfc = ph(:,1) !pression de surface
268	gbx%skt = skt !Skin temperature (K)
269
270	do ip = 1, Npoints
271	if (fracTerLic(ip).ge.0.5) then
272	gbx%land(ip) = 1.
273	else
274	gbx%land(ip) = 0.
275	endif
276	enddo
277	gbx%mr_ozone = mr_ozone !mass_fraction_of_ozone_in_air (kg/kg)
278	! A voir l equivalent LMDZ (u10m et v10m)
279	gbx%u_wind = u_wind !eastward_wind (m s-1)
280	gbx%v_wind = v_wind !northward_wind
281
282	! sunlit calcule a partir de la fraction d ensoleillement par jour
283	! do ip = 1, Npoints
284	! if (sunlit(ip).le.0.) then
285	! gbx%sunlit(ip)=0.
286	! else
287	! gbx%sunlit(ip)=1.
288	! endif
289	! enddo
290	gbx%sunlit=sunlit
291
292	! A voir l equivalent LMDZ
293	mr_ccliq = 0.0
294	mr_ccice = 0.0
295	gbx%mr_hydro(:,:,I_LSCLIQ) = mr_lsliq !mixing_ratio_large_scale_cloud_liquid (kg/kg)
296	gbx%mr_hydro(:,:,I_LSCICE) = mr_lsice !mixing_ratio_large_scale_cloud_ic
297	gbx%mr_hydro(:,:,I_CVCLIQ) = mr_ccliq !mixing_ratio_convective_cloud_liquid
298	gbx%mr_hydro(:,:,I_CVCICE) = mr_ccice !mixing_ratio_convective_cloud_ice
299	! A revoir
300	fl_lsrain = fl_lsrainI + fl_ccrainI
301	fl_lssnow = fl_lssnowI + fl_ccsnowI
302	gbx%rain_ls = fl_lsrain !flux_large_scale_cloud_rain (kg m^-2 s^-1)
303	gbx%snow_ls = fl_lssnow !flux_large_scale_cloud_snow
304	! A voir l equivalent LMDZ
305	fl_lsgrpl=0.
306	fl_ccsnow = 0.
307	fl_ccrain = 0.
308	gbx%grpl_ls = fl_lsgrpl !flux_large_scale_cloud_graupel
309	gbx%rain_cv = fl_ccrain !flux_convective_cloud_rain
310	gbx%snow_cv = fl_ccsnow !flux_convective_cloud_snow
311
312	gbx%Reff(:,:,I_LSCLIQ) = ref_liq*1e-6
313	gbx%Reff(:,:,I_LSCICE) = ref_ice*1e-6
314	!! AI A revoir
315	gbx%Reff(:,:,I_CVCLIQ) = ref_liq*1e-6
316	gbx%Reff(:,:,I_CVCICE) = ref_ice*1e-6
317
318	! ISCCP simulator
319	gbx%dtau_s = dtau_s
320	gbx%dtau_c = 0.
321	gbx%dem_s = dem_s
322	gbx%dem_c = 0.
323
324	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
325	! Define new vertical grid
326	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
327	! print *, 'Defining new vertical grid...'
328	call construct_cosp_vgrid(gbx,Nlr,use_vgrid,csat_vgrid,vgrid)
329
330	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
331	! Allocate memory for other types
332	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
333	! print *, 'Allocating memory for other types...'
334	call construct_cosp_subgrid(Npoints, Ncolumns, Nlevels, sgx)
335	call construct_cosp_sgradar(cfg,Npoints,Ncolumns,Nlevels,N_HYDRO,sgradar)
336	call construct_cosp_radarstats(cfg,Npoints,Ncolumns,vgrid%Nlvgrid,N_HYDRO,stradar)
337	call construct_cosp_sglidar(cfg,Npoints,Ncolumns,Nlevels,N_HYDRO,PARASOL_NREFL,sglidar)
338	call construct_cosp_lidarstats(cfg,Npoints,Ncolumns,vgrid%Nlvgrid,N_HYDRO,PARASOL_NREFL,stlidar)
339	call construct_cosp_isccp(cfg,Npoints,Ncolumns,Nlevels,isccp)
340	!! AI rajout
341	call construct_cosp_modis(cfg,Npoints,modis)
342	!!
343	call construct_cosp_misr(cfg,Npoints,misr)
344	! call construct_cosp_rttov(cfg,Npoints,Nchannels,rttov)
345
346	!+++++++++++++ Open output files and define output files axis !+++++++++++++
347	if (debut_cosp) then
348
349	!$OMP MASTER
350	! print *, ' Open outpts files and define axis'
351	call cosp_output_open(Nlevlmdz, Ncolumns, presnivs, dtime, freq_cosp, &
352	ok_mensuelCOSP, ok_journeCOSP, ok_hfCOSP, ok_all_xml, &
353	ecrit_mth, ecrit_day, ecrit_hf, use_vgrid, vgrid, stlidar)
354	!$OMP END MASTER
355	!$OMP BARRIER
356	endif ! debut_cosp
357	! else
358	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
359	! Call simulator
360	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
361	! print *, 'Calling simulator...'
362	!! AI
363	! call cosp(overlap,Ncolumns,cfg,vgrid,gbx,sgx,sgradar,sglidar,isccp,misr,stradar,stlidar)
364	!#ifdef RTTOV
365	! call cosp(overlap,Ncolumns,cfg,vgrid,gbx,sgx,sgradar,sglidar,isccp,misr,modis,rttov,stradar,stlidar)
366	!#else
367	if (.NOT. debut_cosp) call cosp(overlap,Ncolumns,cfg,vgrid,gbx,sgx,sgradar,sglidar,isccp,misr,modis,stradar,stlidar)
368	!#endif
369	!!
370	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
371
372	!!!!!!!!!!!!!!!!!! Ecreture des sorties Cosp !!!!!!!!!!!!!!r!!!!!!:!!!!!
373
374	! print *, 'Calling write output'
375	if (.NOT. debut_cosp) call cosp_output_write(Nlevlmdz, Npoints, Ncolumns, itap, dtime, freq_COSP, missing_val, &
376	cfg, gbx, vgrid, sglidar, sgradar, stlidar, stradar, &
377	isccp, misr, modis)
378	! endif !debut_cosp
379	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
380	! Deallocate memory in derived types
381	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
382	! print *, 'Deallocating memory...'
383	call free_cosp_gridbox(gbx)
384	call free_cosp_subgrid(sgx)
385	call free_cosp_sgradar(sgradar)
386	call free_cosp_radarstats(stradar)
387	call free_cosp_sglidar(sglidar)
388	call free_cosp_lidarstats(stlidar)
389	call free_cosp_isccp(isccp)
390	call free_cosp_misr(misr)
391	!! AI
392	call free_cosp_modis(modis)
393	! call free_cosp_rttov(rttov)
394	!!
395	call free_cosp_vgrid(vgrid)
396
397	!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
398	! Time in s. Only for testing purposes
399	! call system_clock(t1,count_rate,count_max)
400	! print ,(t1-t0)1.0/count_rate
401	if (debut_cosp) then
402	debut_cosp=.false.
403	endif
404
405	CONTAINS
406
407	SUBROUTINE read_cosp_input
408
409	IF (is_master) THEN
410	OPEN(10,file=cosp_input_nl,status='old')
411	READ(10,nml=cosp_input)
412	CLOSE(10)
413	ENDIF
414	!$OMP BARRIER
415	CALL bcast(overlap)
416	CALL bcast(isccp_topheight)
417	CALL bcast(isccp_topheight_direction)
418	CALL bcast(npoints_it)
419	CALL bcast(ncolumns)
420	CALL bcast(use_vgrid)
421	CALL bcast(nlr)
422	CALL bcast(csat_vgrid)
423	CALL bcast(radar_freq)
424	CALL bcast(surface_radar)
425	CALL bcast(use_mie_tables)
426	CALL bcast(use_gas_abs)
427	CALL bcast(do_ray)
428	CALL bcast(melt_lay)
429	CALL bcast(k2)
430	CALL bcast(Nprmts_max_hydro)
431	CALL bcast(Naero)
432	CALL bcast(Nprmts_max_aero)
433	CALL bcast(lidar_ice_type)
434	CALL bcast(use_precipitation_fluxes)
435	CALL bcast(use_reff)
436	CALL bcast(platform)
437	CALL bcast(satellite)
438	CALL bcast(Instrument)
439	CALL bcast(Nchannels)
440	CALL bcast(Channels)
441	CALL bcast(Surfem)
442	CALL bcast(ZenAng)
443	CALL bcast(co2)
444	CALL bcast(ch4)
445	CALL bcast(n2o)
446	CALL bcast(co)
447
448	END SUBROUTINE read_cosp_input
449
450	end subroutine phys_cosp

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

Download in other formats:

Original Format