My Project
 All Classes Files Functions Variables Macros
phys_cosp.F90
Go to the documentation of this file.
1 ! Simulateur COSP : Cfmip Observation Simulator Package
2 
3 ! ISCCP, Radar (QuickBeam), Lidar et Parasol (ACTSIM), MISR, RTTOVS
4 !Idelkadi Abderrahmane Aout-Septembre 2009
5 
6  subroutine phys_cosp( itap,dtime,freq_cosp, &
7  ok_mensuelcosp,ok_journecosp,ok_hfcosp, &
8  ecrit_mth,ecrit_day,ecrit_hf, &
9  nptslmdz,nlevlmdz,lon,lat, presnivs,overlaplmdz, &
10  ref_liq,ref_ice,fracterlic,u_wind,v_wind,phis,phi,ph,p,skt,t, &
11  sh,rh,tca,cca,mr_lsliq,mr_lsice,fl_lsraini,fl_lssnowi, &
12  fl_ccraini,fl_ccsnowi,mr_ozone,dtau_s,dem_s)
13 
14 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
15 !!!! Inputs :
16 ! itap, !Increment de la physiq
17 ! dtime, !Pas de temps physiq
18 ! overlap, !Overlap type in SCOPS
19 ! Npoints, !Nb de points de la grille physiq
20 ! Nlevels, !Nb de niveaux verticaux
21 ! Ncolumns, !Number of subcolumns
22 ! lon,lat, !Longitudes et latitudes de la grille LMDZ
23 ! ref_liq,ref_ice, !Rayons effectifs des particules liq et ice (en microm)
24 ! fracTerLic, !Fraction terre a convertir en masque
25 ! u_wind,v_wind, !Vents a 10m ???
26 ! phi, !Geopotentiel
27 ! phis, !Geopotentiel sol
28 ! ph, !pression pour chaque inter-couche
29 ! p, !Pression aux milieux des couches
30 ! skt,t, !Temp au sol et temp 3D
31 ! sh, !Humidite specifique
32 ! rh, !Humidite relatif
33 ! tca, !Fraction nuageuse
34 ! cca !Fraction nuageuse convective
35 ! mr_lsliq, !Liq Cloud water content
36 ! mr_lsice, !Ice Cloud water content
37 ! mr_ccliq, !Convective Cloud Liquid water content
38 ! mr_ccice, !Cloud ice water content
39 ! fl_lsrain, !Large scale precipitation lic
40 ! fl_lssnow, !Large scale precipitation ice
41 ! fl_ccrain, !Convective precipitation lic
42 ! fl_ccsnow, !Convective precipitation ice
43 ! mr_ozone, !Concentration ozone (Kg/Kg)
44 ! dem_s !Cloud optical emissivity
45 ! dtau_s !Cloud optical thickness
46 ! emsfc_lw = 1. !Surface emissivity dans radlwsw.F90
47 
48 !!! Outputs :
49 ! calipso2D, !Lidar Low/heigh/Mean/Total-level Cloud Fraction
50 ! calipso3D, !Lidar Cloud Fraction (532 nm)
51 ! cfadlidar, !Lidar Scattering Ratio CFAD (532 nm)
52 ! parasolrefl, !PARASOL-like mono-directional reflectance
53 ! atb, !Lidar Attenuated Total Backscatter (532 nm)
54 ! betamol, !Lidar Molecular Backscatter (532 nm)
55 ! cfaddbze, !Radar Reflectivity Factor CFAD (94 GHz)
56 ! clcalipso2, !Cloud frequency of occurrence as seen by CALIPSO but not CloudSat
57 ! dbze, !Efective_reflectivity_factor
58 ! cltlidarradar, !Lidar and Radar Total Cloud Fraction
59 ! clMISR, !Cloud Fraction as Calculated by the MISR Simulator
60 ! clisccp2, !Cloud Fraction as Calculated by the ISCCP Simulator
61 ! boxtauisccp, !Optical Depth in Each Column as Calculated by the ISCCP Simulator
62 ! boxptopisccp, !Cloud Top Pressure in Each Column as Calculated by the ISCCP Simulator
63 ! tclisccp, !Total Cloud Fraction as Calculated by the ISCCP Simulator
64 ! ctpisccp, !Mean Cloud Top Pressure as Calculated by the ISCCP Simulator
65 ! tauisccp, !Mean Optical Depth as Calculated by the ISCCP Simulator
66 ! albisccp, !Mean Cloud Albedo as Calculated by the ISCCP Simulator
67 ! meantbisccp, !Mean all-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator
68 ! meantbclrisccp !Mean clear-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator
69 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
70 
71  USE mod_cosp_constants
72  USE mod_cosp_types
73  USE mod_cosp
74  USE mod_phys_lmdz_para
75  USE mod_grid_phy_lmdz
76  use ioipsl
77  use iophy
78 
79  IMPLICIT NONE
80 
81  ! Local variables
82  character(len=64),PARAMETER :: cosp_input_nl='cosp_input_nl.txt'
83  character(len=64),PARAMETER :: cosp_output_nl='cosp_output_nl.txt'
84  character(len=512), save :: finput ! Input file name
85  character(len=512), save :: cmor_nl
86  integer, save :: isccp_topheight,isccp_topheight_direction,overlap
87  integer,save :: ncolumns ! Number of subcolumns in SCOPS
88 ! integer,parameter :: Ncollmdz=20
89  integer,parameter :: ncolmax=100
90  integer, save :: npoints ! Number of gridpoints
91 !$OMP THREADPRIVATE(Npoints)
92  integer, save :: nlevels ! Number of levels
93  Integer :: nptslmdz,nlevlmdz ! Nb de points issus de physiq.F
94  integer, save :: nlr ! Number of levels in statistical outputs
95  integer, save :: npoints_it ! Max number of gridpoints to be processed in one iteration
96  integer :: i
97  type(cosp_config),save :: cfg ! Configuration options
98 !$OMP THREADPRIVATE(cfg)
99  type(cosp_gridbox) :: gbx ! Gridbox information. Input for COSP
100  type(cosp_subgrid) :: sgx ! Subgrid outputs
101  type(cosp_sgradar) :: sgradar ! Output from radar simulator
102  type(cosp_sglidar) :: sglidar ! Output from lidar simulator
103  type(cosp_isccp) :: isccp ! Output from ISCCP simulator
104  type(cosp_misr) :: misr ! Output from MISR simulator
105  type(cosp_vgrid) :: vgrid ! Information on vertical grid of stats
106  type(cosp_radarstats) :: stradar ! Summary statistics from radar simulator
107  type(cosp_lidarstats) :: stlidar ! Summary statistics from lidar simulator
108 
109  integer :: t0,t1,count_rate,count_max
110  integer :: nlon,nlat,geomode
111  real,save :: radar_freq,k2,zenang,co2,ch4,n2o,co,emsfc_lw
112 !$OMP THREADPRIVATE(emsfc_lw)
113  integer,dimension(RTTOV_MAX_CHANNELS),save :: channels
114  real,dimension(RTTOV_MAX_CHANNELS),save :: surfem
115  integer, save :: surface_radar,use_mie_tables,use_gas_abs,do_ray,melt_lay
116  integer, save :: nprmts_max_hydro,naero,nprmts_max_aero,lidar_ice_type
117  integer, save :: platform,satellite,instrument,nchannels
118  logical, save :: use_vgrid,csat_vgrid,use_precipitation_fluxes,use_reff
119 
120 ! Declaration necessaires pour les sorties IOIPSL
121  integer :: ii,idayref
122  real :: zjulian,zstoday,zstomth,zstohf,zout,ecrit_day,ecrit_hf,ecrit_mth
123  logical :: ok_mensuelcosp,ok_journecosp,ok_hfcosp
124  integer :: nhori,nvert,nvertp,nvertisccp,nvertm,nvertcol
125  integer, save :: nid_day_cosp,nid_mth_cosp,nid_hf_cosp
126 !$OMP THREADPRIVATE(nid_day_cosp,nid_mth_cosp,nid_hf_cosp)
127  logical, save :: debut_cosp=.true.
128 !$OMP THREADPRIVATE(debut_cosp)
129  integer :: itau_wcosp
130  character(len=2) :: str2
131  real,dimension(Ncolmax) :: column_ax
132  character(len=10),save,dimension(Ncolmax) :: chcol
133 
134  integer, save :: nlevout
135 !$OMP THREADPRIVATE(Nlevout)
136 
137  include "dimensions.h"
138  include "temps.h"
139 
140 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Input variables from LMDZ-GCM
141  integer :: overlaplmdz ! overlap type: 1=max, 2=rand, 3=max/rand ! cosp input (output lmdz)
142  real,dimension(Nptslmdz,Nlevlmdz) :: height,phi,p,ph,t,sh,rh,tca,cca,mr_lsliq,mr_lsice,mr_ccliq,mr_ccice, &
143  fl_lsrain,fl_lssnow,fl_ccrain,fl_ccsnow,fl_lsgrpl, &
144  zlev,zlev_half,mr_ozone,radliq,radice,dtau_s,dem_s,ref_liq,ref_ice
145  real,dimension(Nptslmdz,Nlevlmdz) :: fl_lsraini,fl_lssnowi,fl_ccraini,fl_ccsnowi
146  real,dimension(Nptslmdz) :: lon,lat,skt,fracterlic,u_wind,v_wind,phis
147  real,dimension(Nlevlmdz) :: presnivs
148  integer :: itap,k,ip
149  real :: dtime,freq_cosp
150  logical, parameter :: lcosp=.false.
151 
152 !
153  namelist/cosp_input/cmor_nl,overlap,isccp_topheight,isccp_topheight_direction, &
154  npoints_it,ncolumns,nlevels,use_vgrid,nlr,csat_vgrid,finput, &
155  radar_freq,surface_radar,use_mie_tables, &
156  use_gas_abs,do_ray,melt_lay,k2,nprmts_max_hydro,naero,nprmts_max_aero, &
157  lidar_ice_type,use_precipitation_fluxes,use_reff, &
158  platform,satellite,instrument,nchannels, &
159  channels,surfem,zenang,co2,ch4,n2o,co
160 
161 !---------------- End of declaration of variables --------------
162 
163 
164 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
165 ! Read namelist with COSP inputs
166 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
167 
168  if (debut_cosp) then
169  npoints=nptslmdz
170  nlevels=nlevlmdz
171 
172 ! Lecture du namelist input
173  CALL read_cosp_input
174 
175  do ii=1,ncolumns
176  write(str2,'(i2.2)')ii
177  chcol(ii)="c"//str2
178  column_ax(ii) = real(ii)
179  enddo
180 
181 ! Clefs Outputs
182  call read_cosp_output_nl(cosp_output_nl,cfg)
183 
184  if (overlaplmdz.ne.overlap) then
185  print*,'Attention overlaplmdz different de overlap lu dans namelist '
186  endif
187  print*,'Fin lecture Namelists, debut_cosp =',debut_cosp
188 
189  print*,' Cles sorties cosp :'
190  print*,' Lradar_sim,Llidar_sim,Lisccp_sim,Lmisr_sim,Lrttov_sim', &
191  cfg%Lradar_sim,cfg%Llidar_sim,cfg%Lisccp_sim,cfg%Lmisr_sim,cfg%Lrttov_sim
192 
193  endif ! debut_cosp
194 
195  print*,'Debut phys_cosp itap,dtime,freq_cosp,ecrit_mth,ecrit_day,ecrit_hf ', &
196  itap,dtime,freq_cosp,ecrit_mth,ecrit_day,ecrit_hf
197 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
198 ! Allocate local arrays
199 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
200 ! call system_clock(t0,count_rate,count_max) !!! Only for testing purposes
201 
202 
203 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
204 ! Allocate memory for gridbox type
205 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
206  print *, 'Allocating memory for gridbox type...'
207 
208  call construct_cosp_gridbox(dble(itap),radar_freq,surface_radar,use_mie_tables,use_gas_abs,do_ray,melt_lay,k2, &
209  npoints,nlevels,ncolumns,n_hydro,nprmts_max_hydro,naero,nprmts_max_aero,npoints_it, &
210  lidar_ice_type,isccp_topheight,isccp_topheight_direction,overlap,emsfc_lw, &
211  use_precipitation_fluxes,use_reff, &
212  platform,satellite,instrument,nchannels,zenang, &
213  channels(1:nchannels),surfem(1:nchannels),co2,ch4,n2o,co,gbx)
214 
215 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
216 ! Here code to populate input structure
217 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
218 
219  print *, 'Populating input structure...'
220  gbx%longitude = lon
221  gbx%latitude = lat
222 
223  gbx%p = p !
224  gbx%ph = ph
225  gbx%zlev = phi/9.81
226 
227  zlev_half(:,1) = phis(:)/9.81
228  do k = 2, nlevels
229  do ip = 1, npoints
230  zlev_half(ip,k) = phi(ip,k)/9.81 + &
231  (phi(ip,k)-phi(ip,k-1))/9.81 * (ph(ip,k)-p(ip,k)) / (p(ip,k)-p(ip,k-1))
232  enddo
233  enddo
234  gbx%zlev_half = zlev_half
235 
236  gbx%T = t
237  gbx%q = rh*100.
238  gbx%sh = sh
239  gbx%cca = cca !convective_cloud_amount (1)
240  gbx%tca = tca ! total_cloud_amount (1)
241  gbx%psfc = ph(:,1) !pression de surface
242  gbx%skt = skt !Skin temperature (K)
243 
244  do ip = 1, npoints
245  if (fracterlic(ip).ge.0.5) then
246  gbx%land(ip) = 1.
247  else
248  gbx%land(ip) = 0.
249  endif
250  enddo
251  gbx%mr_ozone = mr_ozone !mass_fraction_of_ozone_in_air (kg/kg)
252 ! A voir l equivalent LMDZ (u10m et v10m)
253  gbx%u_wind = u_wind !eastward_wind (m s-1)
254  gbx%v_wind = v_wind !northward_wind
255 ! Attention
256  gbx%sunlit = 1
257 
258 ! A voir l equivalent LMDZ
259  mr_ccliq = 0.0
260  mr_ccice = 0.0
261  gbx%mr_hydro(:,:,i_lscliq) = mr_lsliq !mixing_ratio_large_scale_cloud_liquid (kg/kg)
262  gbx%mr_hydro(:,:,i_lscice) = mr_lsice !mixing_ratio_large_scale_cloud_ic
263  gbx%mr_hydro(:,:,i_cvcliq) = mr_ccliq !mixing_ratio_convective_cloud_liquid
264  gbx%mr_hydro(:,:,i_cvcice) = mr_ccice !mixing_ratio_convective_cloud_ice
265 ! A revoir
266  fl_lsrain = fl_lsraini + fl_ccraini
267  fl_lssnow = fl_lssnowi + fl_ccsnowi
268  gbx%rain_ls = fl_lsrain !flux_large_scale_cloud_rain (kg m^-2 s^-1)
269  gbx%snow_ls = fl_lssnow !flux_large_scale_cloud_snow
270 ! A voir l equivalent LMDZ
271  fl_lsgrpl=0.
272  fl_ccsnow = 0.
273  fl_ccrain = 0.
274  gbx%grpl_ls = fl_lsgrpl !flux_large_scale_cloud_graupel
275  gbx%rain_cv = fl_ccrain !flux_convective_cloud_rain
276  gbx%snow_cv = fl_ccsnow !flux_convective_cloud_snow
277 
278  gbx%Reff(:,:,i_lscliq) = ref_liq*1e-6
279  gbx%Reff(:,:,i_lscice) = ref_ice*1e-6
280  gbx%Reff(:,:,i_cvcliq) = ref_liq*1e-6
281  gbx%Reff(:,:,i_cvcice) = ref_ice*1e-6
282 
283  ! ISCCP simulator
284  gbx%dtau_s = dtau_s
285  gbx%dtau_c = 0.
286  gbx%dem_s = dem_s
287  gbx%dem_c = 0.
288 
289 ! Surafce emissivity
290  emsfc_lw = 1.
291 
292 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
293  ! Define new vertical grid
294 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
295  print *, 'Defining new vertical grid...'
296  call construct_cosp_vgrid(gbx,nlr,use_vgrid,csat_vgrid,vgrid)
297 
298  if (debut_cosp) then
299 ! Creer le fichier de sorie, definir les variable de sortie
300  ! Axe verticale (Pa ou Km)
301  nlevout = vgrid%Nlvgrid
302 
303  do ii=1,ncolumns
304  column_ax(ii) = real(ii)
305  enddo
306 
307  if (ok_mensuelcosp) then
308  include "ini_histmthCOSP.h"
309  endif
310  if (ok_journecosp) then
311  include "ini_histdayCOSP.h"
312  endif
313  if (ok_hfcosp) then
314  include "ini_histhfCOSP.h"
315  endif
316 
317  debut_cosp=.false.
318  endif ! debut_cosp
319 
320 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
321  ! Allocate memory for other types
322 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
323  print *, 'Allocating memory for other types...'
324  call construct_cosp_subgrid(npoints, ncolumns, nlevels, sgx)
325  call construct_cosp_sgradar(cfg,npoints,ncolumns,nlevels,n_hydro,sgradar)
326  call construct_cosp_radarstats(cfg,npoints,ncolumns,vgrid%Nlvgrid,n_hydro,stradar)
327  call construct_cosp_sglidar(cfg,npoints,ncolumns,nlevels,n_hydro,parasol_nrefl,sglidar)
328  call construct_cosp_lidarstats(cfg,npoints,ncolumns,vgrid%Nlvgrid,n_hydro,parasol_nrefl,stlidar)
329  call construct_cosp_isccp(cfg,npoints,ncolumns,nlevels,isccp)
330  call construct_cosp_misr(cfg,npoints,misr)
331 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
332  ! Call simulator
333 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
334  print *, 'Calling simulator...'
335  call cosp(overlap,ncolumns,cfg,vgrid,gbx,sgx,sgradar,sglidar,isccp,misr,stradar,stlidar)
336 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
337  ! Write outputs to CMOR-compliant NetCDF
338 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
339 
340 ! A traiter le cas ou l on a des valeurs indefinies
341 ! Attention teste
342 
343 ! if(1.eq.0)then
344 
345 
346  do k = 1,nlevout
347  do ip = 1,npoints
348  if(stlidar%lidarcld(ip,k).eq.r_undef)then
349  stlidar%lidarcld(ip,k)=0.
350  endif
351  enddo
352 
353 
354  do ii= 1,sr_bins
355  do ip = 1,npoints
356  if(stlidar%cfad_sr(ip,ii,k).eq.r_undef)then
357  stlidar%cfad_sr(ip,ii,k)=0.
358  endif
359  enddo
360  enddo
361  enddo
362 
363  do ip = 1,npoints
364  do k = 1,nlevlmdz
365  if(sglidar%beta_mol(ip,k).eq.r_undef)then
366  sglidar%beta_mol(ip,k)=0.
367  endif
368 
369  do ii= 1,ncolumns
370  if(sglidar%beta_tot(ip,ii,k).eq.r_undef)then
371  sglidar%beta_tot(ip,ii,k)=0.
372  endif
373  enddo
374 
375  enddo !k = 1,Nlevlmdz
376  enddo !ip = 1,Npoints
377 
378  do k = 1,lidar_ncat
379  do ip = 1,npoints
380  if(stlidar%cldlayer(ip,k).eq.r_undef)then
381  stlidar%cldlayer(ip,k)=0.
382  endif
383  enddo
384  enddo
385 
386 ! endif
387 
388  do ip = 1,npoints
389  if(isccp%totalcldarea(ip).eq.-1.e+30)then
390  isccp%totalcldarea(ip)=0.
391  endif
392  if(isccp%meanptop(ip).eq.-1.e+30)then
393  isccp%meanptop(ip)=0.
394  endif
395  if(isccp%meantaucld(ip).eq.-1.e+30)then
396  isccp%meantaucld(ip)=0.
397  endif
398  if(isccp%meanalbedocld(ip).eq.-1.e+30)then
399  isccp%meanalbedocld(ip)=0.
400  endif
401  if(isccp%meantb(ip).eq.-1.e+30)then
402  isccp%meantb(ip)=0.
403  endif
404  if(isccp%meantbclr(ip).eq.-1.e+30)then
405  isccp%meantbclr(ip)=0.
406  endif
407 
408  do k=1,7
409  do ii=1,7
410  if(isccp%fq_isccp(ip,ii,k).eq.-1.e+30)then
411  isccp%fq_isccp(ip,ii,k)=0.
412  endif
413  enddo
414  enddo
415 
416  do ii=1,ncolumns
417  if(isccp%boxtau(ip,ii).eq.-1.e+30)then
418  isccp%boxtau(ip,ii)=0.
419  endif
420  enddo
421 
422  do ii=1,ncolumns
423  if(isccp%boxptop(ip,ii).eq.-1.e+30)then
424  isccp%boxptop(ip,ii)=0.
425  endif
426  enddo
427  enddo
428 
429  if (ok_mensuelcosp) then
430  include "write_histmthCOSP.h"
431  endif
432  if (ok_journecosp) then
433  include "write_histdayCOSP.h"
434  endif
435  if (ok_hfcosp ) then
436  include "write_histhfCOSP.h"
437  endif
438 
439 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
440  ! Deallocate memory in derived types
441 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
442  print *, 'Deallocating memory...'
443  call free_cosp_gridbox(gbx)
444  call free_cosp_subgrid(sgx)
445  call free_cosp_sgradar(sgradar)
446  call free_cosp_radarstats(stradar)
447  call free_cosp_sglidar(sglidar)
448  call free_cosp_lidarstats(stlidar)
449  call free_cosp_isccp(isccp)
450  call free_cosp_misr(misr)
451  call free_cosp_vgrid(vgrid)
452 
453 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
454  ! Time in s. Only for testing purposes
455 ! call system_clock(t1,count_rate,count_max)
456 ! print *,(t1-t0)*1.0/count_rate
457 
458  CONTAINS
459 
460  SUBROUTINE read_cosp_input
461 
462  IF (is_master) THEN
463  OPEN(10,file=cosp_input_nl,status='old')
464  READ(10,nml=cosp_input)
465  CLOSE(10)
466  ENDIF
467  CALL bcast(cmor_nl)
468  CALL bcast(overlap)
469  CALL bcast(isccp_topheight)
470  CALL bcast(isccp_topheight_direction)
471  CALL bcast(npoints_it)
472  CALL bcast(ncolumns)
473  CALL bcast(nlevels)
474  CALL bcast(use_vgrid)
475  CALL bcast(nlr)
476  CALL bcast(csat_vgrid)
477  CALL bcast(finput)
478  CALL bcast(radar_freq)
479  CALL bcast(surface_radar)
480  CALL bcast(use_mie_tables)
481  CALL bcast(use_gas_abs)
482  CALL bcast(do_ray)
483  CALL bcast(melt_lay)
484  CALL bcast(k2)
485  CALL bcast(nprmts_max_hydro)
486  CALL bcast(naero)
487  CALL bcast(nprmts_max_aero)
488  CALL bcast(lidar_ice_type)
489  CALL bcast(use_precipitation_fluxes)
490  CALL bcast(use_reff)
491  CALL bcast(platform)
492  CALL bcast(satellite)
493  CALL bcast(instrument)
494  CALL bcast(nchannels)
495  CALL bcast(channels)
496  CALL bcast(surfem)
497  CALL bcast(zenang)
498  CALL bcast(co2)
499  CALL bcast(ch4)
500  CALL bcast(n2o)
501  CALL bcast(co)
502 !$OMP BARRIER
503  END SUBROUTINE read_cosp_input
504 
505 end subroutine phys_cosp