read_cosp_output_nl.F90

Go to the documentation of this file.

!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
!--------------- SUBROUTINE READ_COSP_OUTPUT_NL -------------------------
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 SUBROUTINE read_cosp_output_nl(cosp_nl,cfg)
  USE mod_cosp_constants
  USE mod_cosp_types
  USE mod_phys_lmdz_para
  character(len=*),intent(in) :: cosp_nl
  type(cosp_config),intent(out) :: cfg
  ! Local variables
  integer :: i

  logical, save ::   Lradar_sim,Llidar_sim,Lisccp_sim,Lmisr_sim,Lrttov_sim, &
             Lalbisccp,Latb532,Lboxptopisccp,Lboxtauisccp,Lcfad_dbze94, &
             Lcfad_lidarsr532,Lclcalipso2,Lclcalipso,Lclhcalipso,Lclisccp2,Lcllcalipso, &
             Lclmcalipso,Lcltcalipso,Lcltlidarradar,Lctpisccp,Ldbze94,Ltauisccp,Ltclisccp, &
             Llongitude,Llatitude,Lparasol_refl,LclMISR,Lmeantbisccp,Lmeantbclrisccp, &
             Lfrac_out,Lbeta_mol532,Ltbrttov

  namelist/cosp_output/lradar_sim,llidar_sim,lisccp_sim,lmisr_sim,lrttov_sim, &
             lalbisccp,latb532,lboxptopisccp,lboxtauisccp,lcfad_dbze94, &
             lcfad_lidarsr532,lclcalipso2,lclcalipso,lclhcalipso,lclisccp2, &
             lcllcalipso,lclmcalipso,lcltcalipso,lcltlidarradar,lctpisccp,ldbze94,ltauisccp, &
             ltclisccp,llongitude,llatitude,lparasol_refl,lclmisr,lmeantbisccp,lmeantbclrisccp, &
             lfrac_out,lbeta_mol532,ltbrttov

  do i=1,n_out_list
    cfg%out_list(i)=''
  enddo
  
  IF (is_master) THEN
    open(10,file=cosp_nl,status='old')
    read(10,nml=cosp_output)
    close(10)
  ENDIF
  
  CALL bcast(lradar_sim)
  CALL bcast(llidar_sim)
  CALL bcast(lisccp_sim)
  CALL bcast(lmisr_sim)
  CALL bcast(lrttov_sim)
  CALL bcast(lalbisccp)
  CALL bcast(latb532)
  CALL bcast(lboxptopisccp)
  CALL bcast(lboxtauisccp)
  CALL bcast(lcfad_dbze94)
  CALL bcast(lcfad_lidarsr532)
  CALL bcast(lclcalipso2)
  CALL bcast(lclcalipso)
  CALL bcast(lclhcalipso)
  CALL bcast(lclisccp2)
  CALL bcast(lcllcalipso)
  CALL bcast(lclmcalipso)
  CALL bcast(lcltcalipso)
  CALL bcast(lcltlidarradar)
  CALL bcast(lctpisccp)
  CALL bcast(ldbze94)
  CALL bcast(ltauisccp)
  CALL bcast(ltclisccp)
  CALL bcast(llongitude)
  CALL bcast(llatitude)
  CALL bcast(lparasol_refl)
  CALL bcast(lclmisr)
  CALL bcast(lmeantbisccp)
  CALL bcast(lmeantbclrisccp)
  CALL bcast(lfrac_out)
  CALL bcast(lbeta_mol532)
  CALL bcast(ltbrttov)
!$OMP BARRIER

!  print*,' Cles sorties cosp :'
!  print*,' Lradar_sim,Llidar_sim,Lisccp_sim,Lmisr_sim,Lrttov_sim', &
!           Lradar_sim,Llidar_sim,Lisccp_sim,Lmisr_sim,Lrttov_sim

  ! Deal with dependencies
  if (.not.lradar_sim) then
    lcfad_dbze94   = .false.
    lclcalipso2    = .false.
    lcltlidarradar = .false.
    ldbze94        = .false.
  endif
  if (.not.llidar_sim) then
    latb532 = .false.
    lcfad_lidarsr532 = .false.
    lclcalipso2      = .false.
    lclcalipso       = .false.
    lclhcalipso      = .false.
    lcllcalipso      = .false.
    lclmcalipso      = .false.
    lcltcalipso      = .false.
    lcltlidarradar   = .false.
    lparasol_refl    = .false.
    lbeta_mol532     = .false.
  endif
  if (.not.lisccp_sim) then
    lalbisccp       = .false.
    lboxptopisccp   = .false.
    lboxtauisccp    = .false.
    lclisccp2       = .false.
    lctpisccp       = .false.
    ltauisccp       = .false.
    ltclisccp       = .false.
    lmeantbisccp    = .false.
    lmeantbclrisccp = .false.
  endif
  if (.not.lmisr_sim) then
    lclmisr = .false.
  endif
  if (.not.lrttov_sim) then
    ltbrttov = .false.
  endif
  if ((.not.lradar_sim).and.(.not.llidar_sim).and. &
      (.not.lisccp_sim).and.(.not.lmisr_sim)) then
    lfrac_out = .false.
  endif

  ! Diagnostics that use Radar and Lidar
  if (((lclcalipso2).or.(lcltlidarradar)).and.((lradar_sim).or.(llidar_sim))) then
    lclcalipso2    = .true.
    lcltlidarradar = .true.
    llidar_sim     = .true.
    lradar_sim     = .true.
  endif

  cfg%Lstats = .false.
  if ((lradar_sim).or.(llidar_sim).or.(lisccp_sim)) cfg%Lstats = .true.

  ! Copy instrument flags to cfg structure
  cfg%Lradar_sim = lradar_sim
  cfg%Llidar_sim = llidar_sim
  cfg%Lisccp_sim = lisccp_sim
  cfg%Lmisr_sim  = lmisr_sim
  cfg%Lrttov_sim = lrttov_sim

  ! Flag to control output to file
  cfg%Lwrite_output = .false.
  if (cfg%Lstats.or.cfg%Lmisr_sim.or.cfg%Lrttov_sim) then
    cfg%Lwrite_output = .true.
  endif

  ! Output diagnostics
  i = 1
  if (lalbisccp)        cfg%out_list(i) = 'albisccp'
  i = i+1
  if (latb532)          cfg%out_list(i) = 'atb532'
  i = i+1
  if (lboxptopisccp)    cfg%out_list(i) = 'boxptopisccp'
  i = i+1
  if (lboxtauisccp)     cfg%out_list(i) = 'boxtauisccp'
  i = i+1
  if (lcfad_dbze94)     cfg%out_list(i) = 'cfad_dbze94'
  i = i+1
  if (lcfad_lidarsr532) cfg%out_list(i) = 'cfad_lidarsr532'
  i = i+1
  if (lclcalipso2)      cfg%out_list(i) = 'clcalipso2'
  i = i+1
  if (lclcalipso)       cfg%out_list(i) = 'clcalipso'
  i = i+1
  if (lclhcalipso)      cfg%out_list(i) = 'clhcalipso'
  i = i+1
  if (lclisccp2)        cfg%out_list(i) = 'clisccp2'
  i = i+1
  if (lcllcalipso)      cfg%out_list(i) = 'cllcalipso'
  i = i+1
  if (lclmcalipso)      cfg%out_list(i) = 'clmcalipso'
  i = i+1
  if (lcltcalipso)      cfg%out_list(i) = 'cltcalipso'
  i = i+1
  if (lcltlidarradar)   cfg%out_list(i) = 'cltlidarradar'
  i = i+1
  if (lctpisccp)        cfg%out_list(i) = 'ctpisccp'
  i = i+1
  if (ldbze94)          cfg%out_list(i) = 'dbze94'
  i = i+1
  if (ltauisccp)        cfg%out_list(i) = 'tauisccp'
  i = i+1
  if (ltclisccp)        cfg%out_list(i) = 'tclisccp'
  i = i+1
  if (llongitude)       cfg%out_list(i) = 'lon'
  i = i+1
  if (llatitude)        cfg%out_list(i) = 'lat'
  i = i+1
  if (lparasol_refl)    cfg%out_list(i) = 'parasol_refl'
  i = i+1
  if (lclmisr)          cfg%out_list(i) = 'clMISR'
  i = i+1
  if (lmeantbisccp)     cfg%out_list(i) = 'meantbisccp'
  i = i+1
  if (lmeantbclrisccp)  cfg%out_list(i) = 'meantbclrisccp'
  i = i+1
  if (lfrac_out)        cfg%out_list(i) = 'frac_out'
  i = i+1
  if (lbeta_mol532)     cfg%out_list(i) = 'beta_mol532'
  i = i+1
  if (ltbrttov)         cfg%out_list(i) = 'tbrttov'

  if (i /= n_out_list) then
     print *, 'COSP_IO: wrong number of output diagnostics'
     stop
  endif

  ! Copy diagnostic flags to cfg structure
  cfg%Lalbisccp = lalbisccp
  cfg%Latb532 = latb532
  cfg%Lboxptopisccp = lboxptopisccp
  cfg%Lboxtauisccp = lboxtauisccp
  cfg%Lcfad_dbze94 = lcfad_dbze94
  cfg%Lcfad_lidarsr532 = lcfad_lidarsr532
  cfg%Lclcalipso2 = lclcalipso2
  cfg%Lclcalipso = lclcalipso
  cfg%Lclhcalipso = lclhcalipso
  cfg%Lclisccp2 = lclisccp2
  cfg%Lcllcalipso = lcllcalipso
  cfg%Lclmcalipso = lclmcalipso
  cfg%Lcltcalipso = lcltcalipso
  cfg%Lcltlidarradar = lcltlidarradar
  cfg%Lctpisccp = lctpisccp
  cfg%Ldbze94 = ldbze94
  cfg%Ltauisccp = ltauisccp
  cfg%Ltclisccp = ltclisccp
  cfg%Llongitude = llongitude
  cfg%Llatitude = llatitude
  cfg%Lparasol_refl = lparasol_refl
  cfg%LclMISR = lclmisr
  cfg%Lmeantbisccp = lmeantbisccp
  cfg%Lmeantbclrisccp = lmeantbclrisccp
  cfg%Lfrac_out = lfrac_out
  cfg%Lbeta_mol532 = lbeta_mol532
  cfg%Ltbrttov = ltbrttov

 END SUBROUTINE read_cosp_output_nl