source: LMDZ5/trunk/libf/cosp/phys_cosp.F90 @ 1671

Last change on this file since 1671 was 1539, checked in by musat, 13 years ago

Ajouts CFMIP2/CMIP5

  • 6eme fichier de sortie "stations" histstn.nc qui necessite 2 fichiers (voir DefLists?): npCFMIP_param.data(_*) contenant le nombre de points (120 pour simulations AMIP, 73 pour aqua) pointlocations.txt(_*) contenant le numero, les coordonnees (lon,lat) et le nom de chaque station
  • flag LOGICAL dans tous les appels histwrite_phy pour pouvoir sortir le fichier histstn.nc

NB: 1) les flags de type phys_ que l'on met dans le physiq.def_L* pour ajouter plus de sorties

necessitent dorenavant 6 valeurs, la 6eme correspondant au fichier histstn.nc

2) par defaut le fichier histstn.nc ne sort pas; pour le sortir ajouter les lignes suivantes

dans physiq.def_L*

### Type de fichier : global (n) ou stations (y)
phys_out_filestations = n n n n n y

  • introduction de 2 jeux de flags pour les taux des GES; taux actuels avec suffixes _act, taux futurs avec "_per" avec 2 appels au rayonnement si taux "_per" different des taux "_act" (utiles pour diags. CFMIP 4CO2)
  • flags "betaCRF" pour calculs CRF pour experiences sensibilite proprietes optiques eau liquide nuageuse avec initialisations par defaut; sinon besoin de fichier beta_crf.data

Ajout flag LOGICAL lCOSP necessaire pour sortir un fichier stations
IM

File size: 19.5 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
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
505end subroutine phys_cosp
Note: See TracBrowser for help on using the repository browser.