!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% !--------------- 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 !! AI ! 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 !! logical, save :: Lradar_sim,Llidar_sim,Lisccp_sim,Lmodis_sim,Lmisr_sim,Lrttov_sim, & Lalbisccp,Latb532,Lboxptopisccp,Lboxtauisccp,LcfadDbze94, & LcfadLidarsr532,Lclcalipso2,Lclcalipso,Lclhcalipso,Lclisccp,Lcllcalipso, & Lclmcalipso,Lcltcalipso,Lcltlidarradar,Lpctisccp,Ldbze94,Ltauisccp,Lcltisccp, & Lclcalipsoliq,Lclcalipsoice,Lclcalipsoun, & Lclcalipsotmp,Lclcalipsotmpliq,Lclcalipsotmpice,Lclcalipsotmpun, & Lcltcalipsoliq,Lcltcalipsoice,Lcltcalipsoun, & Lclhcalipsoliq,Lclhcalipsoice,Lclhcalipsoun, & Lclmcalipsoliq,Lclmcalipsoice,Lclmcalipsoun, & Lcllcalipsoliq,Lcllcalipsoice,Lcllcalipsoun, & Ltoffset,LparasolRefl,LclMISR,Lmeantbisccp,Lmeantbclrisccp, & Lfracout,LlidarBetaMol532,Ltbrttov, & Lcltmodis,Lclwmodis,Lclimodis,Lclhmodis,Lclmmodis,Lcllmodis,Ltautmodis,Ltauwmodis,Ltauimodis,Ltautlogmodis, & Ltauwlogmodis,Ltauilogmodis,Lreffclwmodis,Lreffclimodis,Lpctmodis,Llwpmodis, & Liwpmodis,Lclmodis,Lcrimodis,Lcrlmodis,Lclopaquecalipso,Lclthincalipso, & !OPAQ (2) Lclzopaquecalipso,Lclcalipsoopaque,Lclcalipsothin,Lclcalipsozopaque,Lclcalipsoopacity, & !OPAQ (5) LprofSR,Lproftemp !TIBO (2) namelist/COSP_OUTPUT/Lradar_sim,Llidar_sim,Lisccp_sim,Lmodis_sim,Lmisr_sim,Lrttov_sim, & Lalbisccp,Latb532,Lboxptopisccp,Lboxtauisccp,LcfadDbze94, & LcfadLidarsr532,Lclcalipso2,Lclcalipso,Lclhcalipso,Lclisccp, & Lcllcalipso,Lclmcalipso,Lcltcalipso,Lcltlidarradar,Lpctisccp,Ldbze94,Ltauisccp, & Lclcalipsoliq,Lclcalipsoice,Lclcalipsoun, & Lclcalipsotmp,Lclcalipsotmpliq,Lclcalipsotmpice,Lclcalipsotmpun, & Lcltcalipsoliq,Lcltcalipsoice,Lcltcalipsoun, & Lclhcalipsoliq,Lclhcalipsoice,Lclhcalipsoun, & Lclmcalipsoliq,Lclmcalipsoice,Lclmcalipsoun, & Lcllcalipsoliq,Lcllcalipsoice,Lcllcalipsoun, & Lcltisccp,Ltoffset,LparasolRefl,LclMISR,Lmeantbisccp,Lmeantbclrisccp, & Lfracout,LlidarBetaMol532,Ltbrttov, & Lcltmodis,Lclwmodis,Lclimodis,Lclhmodis,Lclmmodis,Lcllmodis,Ltautmodis,Ltauwmodis,Ltauimodis,Ltautlogmodis, & Ltauwlogmodis,Ltauilogmodis,Lreffclwmodis,Lreffclimodis,Lpctmodis,Llwpmodis, & Liwpmodis,Lclmodis,Lcrimodis,Lcrlmodis,Lclopaquecalipso,Lclthincalipso, & !OPAQ (2) Lclzopaquecalipso,Lclcalipsoopaque,Lclcalipsothin,Lclcalipsozopaque,Lclcalipsoopacity, & !OPAQ (5) LprofSR,Lproftemp !TIBO (2) 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(Lmodis_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(Lclcalipsoliq) CALL bcast(Lclcalipsoice) CALL bcast(Lclcalipsoun) CALL bcast(Lclcalipsotmp) CALL bcast(Lclcalipsotmpliq) CALL bcast(Lclcalipsotmpice) CALL bcast(Lclcalipsotmpun) CALL bcast(Lcltcalipsoliq) CALL bcast(Lcltcalipsoice) CALL bcast(Lcltcalipsoun) CALL bcast(Lclhcalipsoliq) CALL bcast(Lclhcalipsoice) CALL bcast(Lclhcalipsoun) CALL bcast(Lclmcalipsoliq) CALL bcast(Lclmcalipsoice) CALL bcast(Lclmcalipsoun) CALL bcast(Lcllcalipsoliq) CALL bcast(Lcllcalipsoice) CALL bcast(Lcllcalipsoun) 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(Lcltisccp) CALL bcast(Lparasol_refl) CALL bcast(LclMISR) CALL bcast(Lmeantbisccp) CALL bcast(Lmeantbclrisccp) CALL bcast(Lfrac_out) CALL bcast(Lbeta_mol532) CALL bcast(Lcltmodis) CALL bcast(Lclwmodis) CALL bcast(Lclimodis) CALL bcast(Lclhmodis) CALL bcast(Lclmmodis) CALL bcast(Lcllmodis) CALL bcast(Ltautmodis) CALL bcast(Ltauwmodis) CALL bcast(Ltauimodis) CALL bcast(Ltautlogmodis) CALL bcast(Ltauwlogmodis) CALL bcast(Ltauilogmodis) CALL bcast(Lreffclwmodis) CALL bcast(Lreffclimodis) CALL bcast(Lpctmodis) CALL bcast(Llwpmodis) CALL bcast(Liwpmodis) CALL bcast(Lclmodis) CALL bcast(Ltbrttov) CALL bcast(Lcrimodis) CALL bcast(Lcrlmodis) CALL bcast(Lclopaquecalipso) !OPAQ CALL bcast(Lclthincalipso) !OPAQ CALL bcast(Lclzopaquecalipso) !OPAQ CALL bcast(Lclcalipsoopaque) !OPAQ CALL bcast(Lclcalipsothin) !OPAQ CALL bcast(Lclcalipsozopaque) !OPAQ CALL bcast(Lclcalipsoopacity) !OPAQ CALL bcast(LprofSR) !TIBO CALL bcast(Lproftemp) !TIBO !$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. !! AI Lclcalipsoliq = .false. Lclcalipsoice = .false. Lclcalipsoun = .false. Lclcalipsotmp = .false. Lclcalipsotmpun = .false. Lclcalipsotmpliq = .false. Lclcalipsotmpice = .false. Lclhcalipsoliq = .false. Lcllcalipsoliq = .false. Lclmcalipsoliq = .false. Lcltcalipsoliq = .false. Lclhcalipsoice = .false. Lcllcalipsoice = .false. Lclmcalipsoice = .false. Lcltcalipsoice = .false. Lclhcalipsoun = .false. Lcllcalipsoun = .false. Lclmcalipsoun = .false. Lcltcalipsoun = .false. Lclopaquecalipso = .false. !OPAQ Lclthincalipso = .false. !OPAQ Lclzopaquecalipso = .false. !OPAQ Lclcalipsoopaque = .false. !OPAQ Lclcalipsothin = .false. !OPAQ Lclcalipsozopaque = .false. !OPAQ Lclcalipsoopacity = .false. !OPAQ LprofSR = .false. !TIBO Lproftemp = .false. !TIBO 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 if (.not.Lmodis_sim) then Lcltmodis=.false. Lclwmodis=.false. Lclimodis=.false. Lclhmodis=.false. Lclmmodis=.false. Lcllmodis=.false. Ltautmodis=.false. Ltauwmodis=.false. Ltauimodis=.false. Ltautlogmodis=.false. Ltauwlogmodis=.false. Ltauilogmodis=.false. Lreffclwmodis=.false. Lreffclimodis=.false. Lpctmodis=.false. Llwpmodis=.false. Liwpmodis=.false. Lclmodis=.false. Lcrimodis=.false. Lcrlmodis=.false. endif if (Lmodis_sim) Lisccp_sim = .true. ! 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%Lmodis_sim = Lmodis_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 (LcfadDbze94) cfg%out_list(i) = 'cfadDbze94' i = i+1 if (LcfadLidarsr532) cfg%out_list(i) = 'cfadLidarsr532' 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 (Lclisccp) cfg%out_list(i) = 'clisccp' 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 (Lcllcalipsoice) cfg%out_list(i) = 'cllcalipsoice' i = i+1 if (Lclmcalipsoice) cfg%out_list(i) = 'clmcalipsoice' i = i+1 if (Lclhcalipsoice) cfg%out_list(i) = 'clhcalipsoice' i = i+1 if (Lcltcalipsoice) cfg%out_list(i) = 'cltcalipsoice' i = i+1 if (Lcllcalipsoliq) cfg%out_list(i) = 'cllcalipsoliq' i = i+1 if (Lclmcalipsoliq) cfg%out_list(i) = 'clmcalipsoliq' i = i+1 if (Lclhcalipsoliq) cfg%out_list(i) = 'clhcalipsoliq' i = i+1 if (Lcltcalipsoliq) cfg%out_list(i) = 'cltcalipsoliq' i = i+1 if (Lcllcalipsoun) cfg%out_list(i) = 'cllcalipsoun' i = i+1 if (Lclmcalipsoun) cfg%out_list(i) = 'clmcalipsoun' i = i+1 if (Lclhcalipsoun) cfg%out_list(i) = 'clhcalipsoun' i = i+1 if (Lcltcalipsoun) cfg%out_list(i) = 'cltcalipsoun' i = i+1 if (Lclcalipsoice) cfg%out_list(i) = 'clcalipsoice' i = i+1 if (Lclcalipsoliq) cfg%out_list(i) = 'clcalipsoliq' i = i+1 if (Lclcalipsoun) cfg%out_list(i) = 'clcalipsoun' i = i+1 if (Lclcalipsotmp) cfg%out_list(i) = 'clcalipsotmp' i = i+1 if (Lclcalipsotmpice) cfg%out_list(i) = 'clcalipsotmpice' i = i+1 if (Lclcalipsotmpliq) cfg%out_list(i) = 'clcalipsotmpliq' i = i+1 if (Lclcalipsotmpun) cfg%out_list(i) = 'clcalipsotmpun' i = i+1 if (Lcltlidarradar) cfg%out_list(i) = 'cltlidarradar' i = i+1 if (Lpctisccp) cfg%out_list(i) = 'pctisccp' i = i+1 if (Ldbze94) cfg%out_list(i) = 'dbze94' i = i+1 if (Ltauisccp) cfg%out_list(i) = 'tauisccp' i = i+1 if (Lcltisccp) cfg%out_list(i) = 'cltisccp' i = i+1 if (Ltoffset) cfg%out_list(i) = 'toffset' i = i+1 if (LparasolRefl) cfg%out_list(i) = 'parasolRefl' 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 (Lfracout) cfg%out_list(i) = 'fracout' i = i+1 if (LlidarBetaMol532) cfg%out_list(i) = 'lidarBetaMol532' i = i+1 if (Ltbrttov) cfg%out_list(i) = 'tbrttov' i = i+1 if (Lcltmodis) cfg%out_list(i) = 'cltmodis' i = i+1 if (Lclwmodis) cfg%out_list(i) = 'clwmodis' i = i+1 if (Lclimodis) cfg%out_list(i) = 'climodis' i = i+1 if (Lclhmodis) cfg%out_list(i) = 'clhmodis' i = i+1 if (Lclmmodis) cfg%out_list(i) = 'clmmodis' i = i+1 if (Lcllmodis) cfg%out_list(i) = 'cllmodis' i = i+1 if (Ltautmodis) cfg%out_list(i) = 'tautmodis' i = i+1 if (Ltauwmodis) cfg%out_list(i) = 'tauwmodis' i = i+1 if (Ltauimodis) cfg%out_list(i) = 'tauimodis' i = i+1 if (Ltautlogmodis) cfg%out_list(i) = 'tautlogmodis' i = i+1 if (Ltauwlogmodis) cfg%out_list(i) = 'tauwlogmodis' i = i+1 if (Ltauilogmodis) cfg%out_list(i) = 'tauilogmodis' i = i+1 if (Lreffclwmodis) cfg%out_list(i) = 'reffclwmodis' i = i+1 if (Lreffclimodis) cfg%out_list(i) = 'reffclimodis' i = i+1 if (Lpctmodis) cfg%out_list(i) = 'pctmodis' i = i+1 if (Llwpmodis) cfg%out_list(i) = 'lwpmodis' i = i+1 if (Liwpmodis) cfg%out_list(i) = 'iwpmodis' i = i+1 if (Lclmodis) cfg%out_list(i) = 'clmodis' i = i+1 if (Lcrimodis) cfg%out_list(i) = 'crimodis' i = i+1 if (Lcrlmodis) cfg%out_list(i) = 'crlmodis' i = i+1 !OPAQ if (Lclopaquecalipso) cfg%out_list(i) = 'clopaquecalipso' !OPAQ i = i+1 !OPAQ if (Lclthincalipso) cfg%out_list(i) = 'clthincalipso' !OPAQ i = i+1 !OPAQ if (Lclzopaquecalipso) cfg%out_list(i) = 'clzopaquecalipso' !OPAQ i = i+1 !OPAQ if (Lclcalipsoopaque) cfg%out_list(i) = 'clcalipsoopaque' !OPAQ i = i+1 !OPAQ if (Lclcalipsothin) cfg%out_list(i) = 'clcalipsothin' !OPAQ i = i+1 !OPAQ if (Lclcalipsozopaque) cfg%out_list(i) = 'clcalipsozopaque' !OPAQ i = i+1 !OPAQ if (Lclcalipsoopacity) cfg%out_list(i) = 'clcalipsoopacity' !OPAQ i = i+1 !TIBO if (LprofSR) cfg%out_list(i) = 'profSR' !TIBO i = i+1 !TIBO if (Lproftemp) cfg%out_list(i) = 'proftemp' !TIBO if (i /= N_OUT_LIST) then print *, 'COSP_IO: wrong number of output diagnostics' print *, i,N_OUT_LIST stop endif ! Copy diagnostic flags to cfg structure ! ISCCP simulator cfg%Lalbisccp = Lalbisccp cfg%Latb532 = Latb532 cfg%Lboxptopisccp = Lboxptopisccp cfg%Lboxtauisccp = Lboxtauisccp cfg%Lmeantbisccp = Lmeantbisccp cfg%Lmeantbclrisccp = Lmeantbclrisccp cfg%Lclisccp = Lclisccp cfg%Lpctisccp = Lpctisccp cfg%Ltauisccp = Ltauisccp cfg%Lcltisccp = Lcltisccp ! CloudSat simulator cfg%Ldbze94 = Ldbze94 cfg%LcfadDbze94 = LcfadDbze94 ! CALIPSO/PARASOL simulator cfg%LcfadLidarsr532 = LcfadLidarsr532 cfg%Lclcalipso2 = Lclcalipso2 cfg%Lclcalipso = Lclcalipso cfg%Lclhcalipso = Lclhcalipso cfg%Lcllcalipso = Lcllcalipso cfg%Lclmcalipso = Lclmcalipso cfg%Lcltcalipso = Lcltcalipso cfg%Lclhcalipsoice = Lclhcalipsoice cfg%Lcllcalipsoice = Lcllcalipsoice cfg%Lclmcalipsoice = Lclmcalipsoice cfg%Lcltcalipsoice = Lcltcalipsoice cfg%Lclhcalipsoliq = Lclhcalipsoliq cfg%Lcllcalipsoliq = Lcllcalipsoliq cfg%Lclmcalipsoliq = Lclmcalipsoliq cfg%Lcltcalipsoliq = Lcltcalipsoliq cfg%Lclhcalipsoun = Lclhcalipsoun cfg%Lcllcalipsoun = Lcllcalipsoun cfg%Lclmcalipsoun = Lclmcalipsoun cfg%Lcltcalipsoun = Lcltcalipsoun cfg%Lclcalipsoice = Lclcalipsoice cfg%Lclcalipsoliq = Lclcalipsoliq cfg%Lclcalipsoun = Lclcalipsoun cfg%Lclcalipsotmp = Lclcalipsotmp cfg%Lclcalipsotmpice = Lclcalipsotmpice cfg%Lclcalipsotmpliq = Lclcalipsotmpliq cfg%Lclcalipsotmpun = Lclcalipsotmpun cfg%Lcltlidarradar = Lcltlidarradar cfg%LparasolRefl = LparasolRefl cfg%Lclopaquecalipso = Lclopaquecalipso !OPAQ cfg%Lclthincalipso = Lclthincalipso !OPAQ cfg%Lclzopaquecalipso = Lclzopaquecalipso !OPAQ cfg%Lclcalipsoopaque = Lclcalipsoopaque !OPAQ cfg%Lclcalipsothin = Lclcalipsothin !OPAQ cfg%Lclcalipsozopaque = Lclcalipsozopaque !OPAQ cfg%Lclcalipsoopacity = Lclcalipsoopacity !OPAQ cfg%LprofSR = LprofSR !TIBO cfg%Lproftemp = Lproftemp !TIBO ! MISR simulator cfg%LclMISR = LclMISR ! Other cfg%Ltoffset = Ltoffset cfg%Lfracout = Lfracout cfg%LlidarBetaMol532 = LlidarBetaMol532 ! RTTOV cfg%Ltbrttov = Ltbrttov ! MODIS simulator cfg%Lcltmodis=Lcltmodis cfg%Lclwmodis=Lclwmodis cfg%Lclimodis=Lclimodis cfg%Lclhmodis=Lclhmodis cfg%Lclmmodis=Lclmmodis cfg%Lcllmodis=Lcllmodis cfg%Ltautmodis=Ltautmodis cfg%Ltauwmodis=Ltauwmodis cfg%Ltauimodis=Ltauimodis cfg%Ltautlogmodis=Ltautlogmodis cfg%Ltauwlogmodis=Ltauwlogmodis cfg%Ltauilogmodis=Ltauilogmodis cfg%Lreffclwmodis=Lreffclwmodis cfg%Lreffclimodis=Lreffclimodis cfg%Lpctmodis=Lpctmodis cfg%Llwpmodis=Llwpmodis cfg%Liwpmodis=Liwpmodis cfg%Lclmodis=Lclmodis cfg%Lcrimodis=Lcrimodis cfg%Lcrlmodis=Lcrlmodis END SUBROUTINE READ_COSP_OUTPUT_NL