Context Navigation

cosp_output_mod.F90 @ 4352

Last change on this file since 4352 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

File size: 25.6 KB

Rev	Line
[3358]	1	! A.Idelkadi sept 2013 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	2	!! Module pour declarer et initialiser les parametres de controle des fichiers de sorties et des champs a sortir
	3	!! La routine cosp_output_open (appelee 1 seule fois dans phy_cosp.F90) permet :
	4	!! de creer les fichiers avec leurs grilles horizontales et verticales
	5	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	6	MODULE cosp_output_mod
	7
	8	! USE MOD_COSP_CONSTANTS
	9	! USE MOD_COSP_TYPES
	10
	11	! use MOD_COSP_INTERFACE_v1p4
	12	! use MOD_COSP_CONFIG
	13	! use MOD_COSP_Modis_Simulator, only : cosp_modis
	14	! use mod_modis_sim, only : numMODISReffIceBins, reffICE_binCenters, &
	15	! numMODISReffLiqBins, reffLIQ_binCenters
	16	USE COSP_KINDS, ONLY: wp,dp
	17	IMPLICIT NONE
	18	! cosp_output_mod
	19	INTEGER :: i
	20	!!!!!!! Controle des fichier de sorties Cosp !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	21	LOGICAL, DIMENSION(3), SAVE :: cosp_outfilekeys
	22	INTEGER, DIMENSION(3), SAVE :: cosp_nidfiles
	23	!$OMP THREADPRIVATE(cosp_outfilekeys, cosp_nidfiles)
	24	INTEGER, DIMENSION(3), SAVE :: nhoricosp,nvert,nvertmcosp,nvertcol,nvertbze, &
	25	nvertsratio,nvertisccp,nvertp,nverttemp,nvertmisr, &
	26	nvertReffIce,nvertReffLiq,nverttau
	27	REAL, DIMENSION(3), SAVE :: zoutm_cosp
	28	!$OMP THREADPRIVATE(nhoricosp, nvert,nvertmcosp,nvertcol,nvertsratio,nvertbze,nvertisccp,nvertp,zoutm_cosp,nverttemp,nvertmisr)
	29	!$OMP THREADPRIVATE(nvertReffIce,nvertReffLiq,nverttau)
	30	REAL, SAVE :: zdtimemoy_cosp
	31	!$OMP THREADPRIVATE(zdtimemoy_cosp)
	32	CHARACTER(LEN=20), DIMENSION(3), SAVE :: cosp_outfiletypes
	33	CHARACTER(LEN=20), DIMENSION(3), SAVE :: cosp_outfilenames
	34	REAL, DIMENSION(3), SAVE :: cosp_ecritfiles
	35	!$OMP THREADPRIVATE(cosp_outfiletypes, cosp_outfilenames, cosp_ecritfiles)
	36
	37	!!!! Controle des variables a sortir dans les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	38	TYPE ctrl_outcosp
	39	LOGICAL,DIMENSION(3) :: cles !!! Sortir ou non le champs
	40	CHARACTER(len=20) :: name
	41	CHARACTER(len=150) :: description !!! Nom
	42	CHARACTER(len=20) :: unit !!! Unite
	43	CHARACTER(len=20),DIMENSION(3) :: cosp_typeecrit !!! Operation (ave, inst, ...)
	44	END TYPE ctrl_outcosp
	45
	46	! CALIPSO vars
	47	TYPE(ctrl_outcosp), SAVE :: o_cllcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	48	"cllcalipso", "Lidar Low-level Cloud Fraction", "1", (/ ('', i=1, 3) /))
	49	TYPE(ctrl_outcosp), SAVE :: o_clmcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	50	"clmcalipso", "Lidar Mid-level Cloud Fraction", "1", (/ ('', i=1, 3) /))
	51	TYPE(ctrl_outcosp), SAVE :: o_clhcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	52	"clhcalipso", "Lidar High-level Cloud Fraction", "1", (/ ('', i=1, 3) /))
	53	TYPE(ctrl_outcosp), SAVE :: o_cltcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	54	"cltcalipso", "Lidar Total Cloud Fraction", "1", (/ ('', i=1, 3) /))
	55	TYPE(ctrl_outcosp), SAVE :: o_clcalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	56	"clcalipso", "Lidar Cloud Fraction (532 nm)", "1", (/ ('', i=1, 3) /))
	57	TYPE(ctrl_outcosp), SAVE :: o_cfad_lidarsr532 = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
	58	"cfad_lidarsr532", "Lidar Scattering Ratio CFAD (532 nm)", "1", (/ ('', i=1, 3) /))
	59	TYPE(ctrl_outcosp), SAVE :: o_parasol_refl = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	60	"parasol_refl", "PARASOL-like mono-directional reflectance","1", (/ ('', i=1, 3) /))
	61	TYPE(ctrl_outcosp), SAVE :: o_parasol_crefl = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	62	"parasol_crefl", "PARASOL-like mono-directional reflectance (integral)","1", (/ ('', i=1, 3) /))
	63	TYPE(ctrl_outcosp), SAVE :: o_Ncrefl = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	64	"Ncrefl", "Nb PARASOL-like mono-directional reflectance (integral)","1", (/ ('', i=1, 3) /))
	65	TYPE(ctrl_outcosp), SAVE :: o_atb532 = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
	66	"atb532", "Lidar Attenuated Total Backscatter (532 nm)","1", (/ ('', i=1, 3) /))
	67	TYPE(ctrl_outcosp), SAVE :: o_beta_mol532 = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
	68	"beta_mol532", "Lidar Molecular Backscatter (532 nm)","m-1 sr-1", (/ ('', i=1, 3) /))
	69	!! AI 11 2015
	70	TYPE(ctrl_outcosp), SAVE :: o_cllcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	71	"cllcalipsoice", "CALIPSO Ice-Phase Low Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	72	TYPE(ctrl_outcosp), SAVE :: o_cllcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	73	"cllcalipsoliq", "CALIPSO Liq-Phase Low Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	74	TYPE(ctrl_outcosp), SAVE :: o_clmcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	75	"clmcalipsoice", "CALIPSO Ice-Phase Mid Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	76	TYPE(ctrl_outcosp), SAVE :: o_clmcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	77	"clmcalipsoliq", "CALIPSO Liq-Phase Mid Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	78	TYPE(ctrl_outcosp), SAVE :: o_clhcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	79	"clhcalipsoice", "CALIPSO Ice-Phase High Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	80	TYPE(ctrl_outcosp), SAVE :: o_clhcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	81	"clhcalipsoliq", "CALIPSO Liq-Phase High Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	82	TYPE(ctrl_outcosp), SAVE :: o_cltcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	83	"cltcalipsoice", "CALIPSO Ice-Phase Tot Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	84	TYPE(ctrl_outcosp), SAVE :: o_cltcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	85	"cltcalipsoliq", "CALIPSO Liq-Phase Tot Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	86	TYPE(ctrl_outcosp), SAVE :: o_cllcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	87	"cllcalipsoun", "CALIPSO Undefined-Phase Low Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	88	TYPE(ctrl_outcosp), SAVE :: o_clmcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	89	"clmcalipsoun", "CALIPSO Undefined-Phase Mid Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	90	TYPE(ctrl_outcosp), SAVE :: o_clhcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	91	"clhcalipsoun", "CALIPSO Undefined-Phase High Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	92	TYPE(ctrl_outcosp), SAVE :: o_cltcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	93	"cltcalipsoun", "CALIPSO Undefined-Phase Tot Level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	94	TYPE(ctrl_outcosp), SAVE :: o_clcalipsoice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	95	"clcalipsoice", "Lidar Ice-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
	96	TYPE(ctrl_outcosp), SAVE :: o_clcalipsoliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	97	"clcalipsoliq", "Lidar Liq-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
	98	TYPE(ctrl_outcosp), SAVE :: o_clcalipsoun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	99	"clcalipsoun", "Lidar Undef-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
	100	TYPE(ctrl_outcosp), SAVE :: o_clcalipsotmpice = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	101	"clcalipsotmpice", "Lidar Ice-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
	102	TYPE(ctrl_outcosp), SAVE :: o_clcalipsotmpliq = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	103	"clcalipsotmpliq", "Lidar Liq-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
	104	TYPE(ctrl_outcosp), SAVE :: o_clcalipsotmpun = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	105	"clcalipsotmpun", "Lidar Undef-Phase Cloud Fraction", "%", (/ ('', i=1, 3) /))
	106	TYPE(ctrl_outcosp), SAVE :: o_clcalipsotmp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	107	"clcalipsotmp", "Lidar Cloud Fraction", "%", (/ ('', i=1, 3) /))
	108
	109	TYPE(ctrl_outcosp), SAVE :: o_clopaquecalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
	110	"clopaquecalipso", "Lidar Opaque Cloud Fraction", "%", (/ ('', i=1, 3) /)) !OPAQ
	111	TYPE(ctrl_outcosp), SAVE :: o_clthincalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
	112	"clthincalipso", "Lidar Thin Cloud Fraction", "%", (/ ('', i=1, 3) /)) !OPAQ
	113	TYPE(ctrl_outcosp), SAVE :: o_clzopaquecalipso = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
	114	"clzopaquecalipso", "Lidar mean opacity altitude", "m", (/ ('', i=1, 3) /)) !OPAQ
	115	TYPE(ctrl_outcosp), SAVE :: o_clcalipsoopaque = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
	116	"clcalipsoopaque", "Lidar Opaque profile Cloud Fraction", "%", (/ ('', i=1, 3) /)) !OPAQ
	117	TYPE(ctrl_outcosp), SAVE :: o_clcalipsothin = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
	118	"clcalipsothin", "Lidar Thin profile Cloud Fraction", "%", (/ ('', i=1, 3) /)) !OPAQ
	119	TYPE(ctrl_outcosp), SAVE :: o_clcalipsozopaque = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
	120	"clcalipsozopaque", "Lidar z_opaque Fraction", "%", (/ ('', i=1, 3) /)) !OPAQ
	121	TYPE(ctrl_outcosp), SAVE :: o_clcalipsoopacity = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !OPAQ
	122	"clcalipsoopacity", "Lidar opacity Fraction", "%", (/ ('', i=1, 3) /)) !OPAQ
	123
	124	TYPE(ctrl_outcosp), SAVE :: o_proftemp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !TIBO
	125	"proftemp", "Temperature profiles (40 lev)", "K", (/ ('', i=1, 3) /)) !TIBO
	126	TYPE(ctrl_outcosp), SAVE :: o_profSR = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), & !TIBO
	127	"profSR", "Lidar Scattering Ratio profiles (532 nm)", "1", (/ ('', i=1, 3) /)) !TIBO
	128
	129	! Radar Cloudsat
	130	TYPE(ctrl_outcosp), SAVE :: o_cfadDbze94 = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	131	"cfadDbze94", "CloudSat Radar Reflectivity CFAD", "%", (/ ('', i=1, 3) /))
	132	TYPE(ctrl_outcosp), SAVE :: o_dbze94 = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	133	"dbze94", "CloudSat Radar Reflectivity", "%", (/ ('', i=1, 3) /))
	134
	135	! Calipso + Cloudsat
	136	TYPE(ctrl_outcosp), SAVE :: o_clcalipso2 = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	137	"clcalipso2", "CALIPSO Cloud Fraction Undetected by CloudSat", "1", (/ ('', i=1, 3) /))
	138	TYPE(ctrl_outcosp), SAVE :: o_cltlidarradar = ctrl_outcosp((/ .TRUE., .TRUE.,.TRUE. /), &
	139	"cltlidarradar", "Lidar and Radar Total Cloud Fraction", "%", (/ ('', i=1, 3) /))
	140
	141	! ISCCP vars
	142	TYPE(ctrl_outcosp), SAVE :: o_sunlit = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	143	"sunlit", "1 for day points, 0 for nightime","1",(/ ('', i=1, 3) /))
	144	TYPE(ctrl_outcosp), SAVE :: o_clisccp2 = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	145	"clisccp2", "Cloud Fraction as Calculated by the ISCCP Simulator","%", (/ ('', i=1, 3) /))
	146	TYPE(ctrl_outcosp), SAVE :: o_boxtauisccp = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
	147	"boxtauisccp", "Optical Depth in Each Column as Calculated by the ISCCP Simulator","1", (/ ('', i=1, 3) /))
	148	TYPE(ctrl_outcosp), SAVE :: o_boxptopisccp = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
	149	"boxptopisccp", "Cloud Top Pressure in Each Column as Calculated by the ISCCP Simulator","Pa", (/ ('', i=1, 3) /))
	150	TYPE(ctrl_outcosp), SAVE :: o_tclisccp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	151	"tclisccp", "Total Cloud Fraction as Calculated by the ISCCP Simulator", "%", (/ ('', i=1, 3) /))
	152	TYPE(ctrl_outcosp), SAVE :: o_ctpisccp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	153	"ctpisccp", "Mean Cloud Top Pressure as Calculated by the ISCCP Simulator", "Pa", (/ ('', i=1, 3) /))
	154	TYPE(ctrl_outcosp), SAVE :: o_tauisccp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	155	"tauisccp", "Optical Depth as Calculated by the ISCCP Simulator", "1", (/ ('', i=1, 3) /))
	156	TYPE(ctrl_outcosp), SAVE :: o_albisccp = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	157	"albisccp", "Mean Cloud Albedo as Calculated by the ISCCP Simulator", "1", (/ ('', i=1, 3) /))
	158	TYPE(ctrl_outcosp), SAVE :: o_meantbisccp = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
	159	"meantbisccp", " Mean all-sky 10.5 micron brightness temperature as calculated &
	160	by the ISCCP Simulator","K", (/ ('', i=1, 3) /))
	161	TYPE(ctrl_outcosp), SAVE :: o_meantbclrisccp = ctrl_outcosp((/ .FALSE., .FALSE., .FALSE. /), &
	162	"meantbclrisccp", "Mean clear-sky 10.5 micron brightness temperature as calculated &
	163	by the ISCCP Simulator","K", (/ ('', i=1, 3) /))
	164
	165	! MISR simulator
	166	TYPE(ctrl_outcosp), SAVE :: o_clMISR = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	167	"clMISR", "Cloud Fraction as Calculated by the MISR Simulator","%", (/ ('', i=1, 3) /))
	168
	169	! MODIS simulator
	170	TYPE(ctrl_outcosp), SAVE :: o_cllmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	171	"cllmodis", "MODIS Low-level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	172	TYPE(ctrl_outcosp), SAVE :: o_clmmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	173	"clmmodis", "MODIS Mid-level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	174	TYPE(ctrl_outcosp), SAVE :: o_clhmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	175	"clhmodis", "MODIS High-level Cloud Fraction", "%", (/ ('', i=1, 3) /))
	176	TYPE(ctrl_outcosp), SAVE :: o_cltmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	177	"cltmodis", "MODIS Total Cloud Fraction", "%", (/ ('', i=1, 3) /))
	178	TYPE(ctrl_outcosp), SAVE :: o_clwmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	179	"clwmodis", "MODIS Cloud Fraction water mean", "%", (/ ('', i=1, 3) /))
	180	TYPE(ctrl_outcosp), SAVE :: o_climodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	181	"climodis", "MODIS Cloud Fraction ice mean", "%", (/ ('', i=1, 3) /))
	182	TYPE(ctrl_outcosp), SAVE :: o_tautmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	183	"tautmodis", "MODIS Optical_Thickness_Total_Mean", "1", (/ ('', i=1, 3) /))
	184	TYPE(ctrl_outcosp), SAVE :: o_tauwmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	185	"tauwmodis", "MODIS Optical_Thickness_Water_Mean", "1", (/ ('', i=1, 3) /))
	186	TYPE(ctrl_outcosp), SAVE :: o_tauimodis= ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	187	"tauimodis", "MODIS Optical_Thickness_Ice_Mean", "1", (/ ('', i=1, 3) /))
	188	TYPE(ctrl_outcosp), SAVE :: o_tautlogmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	189	"tautlogmodis", "MODIS Optical_Thickness_Total_logMean", "1", (/ ('', i=1, 3) /))
	190	TYPE(ctrl_outcosp), SAVE :: o_tauwlogmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	191	"tauwlogmodis", "MODIS Optical_Thickness_Water_logMean", "1", (/ ('', i=1, 3) /))
	192	TYPE(ctrl_outcosp), SAVE :: o_tauilogmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	193	"tauilogmodis", "MODIS Optical_Thickness_Ice_logMean", "1", (/ ('', i=1, 3) /))
	194	TYPE(ctrl_outcosp), SAVE :: o_reffclwmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	195	"reffclwmodis", "Modis Cloud_Particle_Size_Water_Mean", "m", (/ ('', i=1, 3) /))
	196	TYPE(ctrl_outcosp), SAVE :: o_reffclimodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	197	"reffclimodis", "Modis Cloud_Particle_Size_Ice_Mean", "m", (/ ('', i=1, 3) /))
	198	TYPE(ctrl_outcosp), SAVE :: o_pctmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	199	"pctmodis", "Modis Cloud_Top_Pressure_Total_Mean", "Pa", (/ ('', i=1, 3) /))
	200	TYPE(ctrl_outcosp), SAVE :: o_lwpmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	201	"lwpmodis", "Modis Liquid_Water_Path_Mean", "kg m-2", (/ ('', i=1, 3) /))
	202	TYPE(ctrl_outcosp), SAVE :: o_iwpmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	203	"iwpmodis", "Modis Ice_Water_Path_Mean", "kg m-2", (/ ('', i=1, 3) /))
	204	TYPE(ctrl_outcosp), SAVE :: o_clmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	205	"clmodis", "MODIS Cloud Area Fraction", "%", (/ ('', i=1, 3) /))
	206	TYPE(ctrl_outcosp), SAVE :: o_crimodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	207	"crimodis", "Optical_Thickness_vs_ReffIce from Modis", "%", (/ ('',i=1, 3) /))
	208	TYPE(ctrl_outcosp), SAVE :: o_crlmodis = ctrl_outcosp((/ .TRUE., .TRUE.,.TRUE. /), &
	209	"crlmodis", "Optical_Thickness_vs_ReffLiq from Modis", "%", (/ ('',i=1, 3) /))
	210
	211	! Rttovs simulator
	212	TYPE(ctrl_outcosp), SAVE :: o_tbrttov = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	213	"tbrttov", "Rttovs Cloud Area Fraction", "%", (/ ('', i=1, 3) /))
	214
	215	! Scops and others
	216	TYPE(ctrl_outcosp), SAVE :: o_fracout = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
	217	"fracout", "Subcolumn output from SCOPS", "%", (/ ('', i=1, 3) /))
	218
	219	LOGICAL, SAVE :: cosp_varsdefined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL
	220	REAL, SAVE :: Cosp_fill_value
	221	!$OMP THREADPRIVATE(Cosp_fill_value)
	222
	223
	224	CONTAINS
	225
	226	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	227	!!!!!!!!! Ouverture des fichier et definition des axes!!!!!!!!
	228	!! histbeg, histvert
	229	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	230
	231	SUBROUTINE cosp_output_open(Nlevlmdz, Ncolumns, presnivs, dtime, freq_cosp, &
	232	ok_mensuelCOSP, ok_journeCOSP, ok_hfCOSP, ok_all_xml, &
[3369]	233	ecrit_mth, ecrit_day, ecrit_hf, use_vgrid, vgrid)
[3358]	234	use MOD_COSP_INTERFACE_v1p4, only : cosp_vgrid
	235	! use MOD_COSP
	236	use mod_cosp_config, only : DBZE_BINS, SR_BINS, CFAD_ZE_MIN, PARASOL_NREFL, &
	237	CFAD_ZE_WIDTH,vgrid_zl,vgrid_zu,vgrid_z,PARASOL_SZA, &
	238	isccp_histPresCenters,tau_binCenters, LIDAR_NTEMP, &
	239	LIDAR_PHASE_TEMP,misr_histHgtCenters,numMISRHgtBins, &
	240	numMODISReffIceBins,reffICE_binCenters, &
	241	numMODISReffLiqBins, reffLIQ_binCenters, pres_binCenters
	242
	243	USE iophy
	244	USE ioipsl
	245	USE phys_cal_mod
	246	USE time_phylmdz_mod, ONLY: day_ref, annee_ref, day_ini, start_time, itau_phy
	247	USE print_control_mod, ONLY: lunout
	248
	249	#ifdef CPP_XIOS
	250	! ug Pour les sorties XIOS
	251	USE wxios
	252	#endif
	253
	254	IMPLICIT NONE
	255
	256	!!! Variables d'entree
	257	integer :: Nlevlmdz, Ncolumns ! Number of levels
	258	real,dimension(Nlevlmdz) :: presnivs
	259	real :: dtime, freq_cosp, ecrit_day, ecrit_hf, ecrit_mth
[3369]	260	logical :: use_vgrid
[3358]	261	logical :: ok_mensuelCOSP, ok_journeCOSP, ok_hfCOSP, ok_all_xml
	262	type(cosp_vgrid) :: vgrid ! Information on vertical grid of stats
	263	! type(cosp_lidarstats) :: stlidar ! Summary statistics from lidar simulator
	264
	265	!!! Variables locales
	266	integer :: idayref, iff, ii
	267	real :: zjulian,zjulian_start
	268	real(wp),dimension(Ncolumns) :: column_ax
	269	real(wp),dimension(DBZE_BINS) :: dbze_ax
	270	CHARACTER(LEN=20), DIMENSION(3) :: chfreq = (/ '1day', '1d ', '3h ' /)
	271	real(wp),parameter,dimension(SR_BINS) :: sratio_ax = (/0.005, &
	272	0.605,2.09,4.,6., &
	273	8.5,12.5,17.5,22.5,27.5,35.,45.,55.,70.,50040./)
	274
	275	!!! Variables d'entree
	276
	277	#ifdef CPP_XIOS
	278	! ug Variables utilisées pour récupérer le calendrier pour xios
	279	INTEGER :: x_an, x_mois, x_jour
	280	REAL :: x_heure
	281	INTEGER :: ini_an, ini_mois, ini_jour
	282	REAL :: ini_heure
	283	#endif
	284
	285	WRITE(lunout,*) 'Debut cosp_output_mod.F90'
	286	print*,'cosp_varsdefined',cosp_varsdefined
	287	! Initialisations (Valeurs par defaut)
	288
	289	!! Definition valeurs axes
	290	do ii=1,Ncolumns
	291	column_ax(ii) = real(ii)
	292	enddo
	293
	294	do i=1,DBZE_BINS
	295	dbze_ax(i) = CFAD_ZE_MIN + CFAD_ZE_WIDTH*(i - 0.5)
	296	enddo
	297
	298	cosp_outfilenames(1) = 'histmthCOSP'
	299	cosp_outfilenames(2) = 'histdayCOSP'
	300	cosp_outfilenames(3) = 'histhfCOSP'
	301
	302	cosp_outfiletypes(1) = 'ave(X)'
	303	cosp_outfiletypes(2) = 'ave(X)'
	304	cosp_outfiletypes(3) = 'ave(X)'
	305
	306	cosp_outfilekeys(1) = ok_mensuelCOSP
	307	cosp_outfilekeys(2) = ok_journeCOSP
	308	cosp_outfilekeys(3) = ok_hfCOSP
	309
	310	cosp_ecritfiles(1) = mth_len*86400.
	311	cosp_ecritfiles(2) = 1.*86400.
	312	cosp_ecritfiles(3) = 0.125*86400.
	313
	314	! Lecture des parametres dans output.def ou config.def
	315
	316	CALL getin('cosp_outfilenames',cosp_outfilenames)
	317	CALL getin('cosp_outfilekeys',cosp_outfilekeys)
	318	CALL getin('cosp_ecritfiles',cosp_ecritfiles)
	319	CALL getin('cosp_outfiletypes',cosp_outfiletypes)
	320
	321	WRITE(lunout,*)'cosp_outfilenames=',cosp_outfilenames
	322	WRITE(lunout,*)'cosp_outfilekeys=',cosp_outfilekeys
	323	WRITE(lunout,*)'cosp_ecritfiles=',cosp_ecritfiles
	324	WRITE(lunout,*)'cosp_outfiletypes=',cosp_outfiletypes
	325
	326	idayref = day_ref
	327	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	328	CALL ymds2ju(annee_ref, 1, day_ini, start_time, zjulian_start)
	329
	330	#ifdef CPP_XIOS
	331
	332	! recuperer la valeur indefine Xios
	333	! CALL xios_get_field_attr("clcalipso",default_value=Cosp_fill_value)
	334	! Cosp_fill_value=missing_val
	335	Cosp_fill_value=0.
	336	print*,'Cosp_fill_value=',Cosp_fill_value
	337
	338	CALL wxios_add_vaxis("height", vgrid%Nlvgrid, vgrid%z)
	339	print*,'wxios_add_vaxis vgrid%Nlvgrid, vgrid%z',vgrid%Nlvgrid,vgrid%z
	340
	341	WRITE(lunout,*) 'wxios_add_vaxis height_mlev, Nlevlmdz vgrid%mz ', &
	342	Nlevlmdz,vgrid%mz
	343	CALL wxios_add_vaxis("height_mlev", Nlevlmdz, vgrid%mz)
	344
	345	WRITE(lunout,*) 'wxios_add_vaxis sza, PARASOL_NREFL ', &
	346	PARASOL_NREFL, PARASOL_SZA
	347	CALL wxios_add_vaxis("sza", PARASOL_NREFL, PARASOL_SZA)
	348
	349	WRITE(lunout,*) 'wxios_add_vaxis pressure2 ',7,pres_binCenters
	350	CALL wxios_add_vaxis("pressure2", 7, pres_binCenters)
	351
	352	WRITE(lunout,*) 'wxios_add_vaxis column ',Ncolumns,column_ax
	353	CALL wxios_add_vaxis("column", Ncolumns, column_ax)
	354
	355	WRITE(lunout,*) 'wxios_add_vaxis temp LIDAR_NTEMP, LIDAR_PHASE_TEMP ', &
	356	LIDAR_NTEMP, LIDAR_PHASE_TEMP
	357	CALL wxios_add_vaxis("temp", LIDAR_NTEMP, LIDAR_PHASE_TEMP)
	358
	359	WRITE(lunout,*) 'wxios_add_vaxis cth16 numMISRHgtBins, misr_histHgtCenters ', &
	360	numMISRHgtBins, misr_histHgtCenters
[3435]	361	CALL wxios_add_vaxis("cth16", numMISRHgtBins, misr_histHgtCenters)
[3358]	362
	363	WRITE(lunout,*) 'wxios_add_vaxis dbze DBZE_BINS, dbze_ax ', &
	364	DBZE_BINS, dbze_ax
	365	CALL wxios_add_vaxis("dbze", DBZE_BINS, dbze_ax)
	366
	367	WRITE(lunout,*) 'wxios_add_vaxis scatratio SR_BINS, sratio_ax', &
	368	SR_BINS, sratio_ax
	369	CALL wxios_add_vaxis("scatratio", SR_BINS, sratio_ax)
	370
	371	WRITE(lunout,*) 'wxios_add_vaxis ReffIce numMODISReffIceBins, &
	372	reffICE_binCenters',numMODISReffIceBins, reffICE_binCenters
	373	CALL wxios_add_vaxis("ReffIce", numMODISReffIceBins, reffICE_binCenters)
	374
	375	WRITE(lunout,*) 'wxios_add_vaxis ReffLiq numMODISReffLiqBins, &
	376	reffLIQ_binCenters', numMODISReffLiqBins, reffLIQ_binCenters
	377	CALL wxios_add_vaxis("ReffLiq", numMODISReffLiqBins, reffLIQ_binCenters)
	378
	379	WRITE(lunout,*) 'wxios_add_vaxis 7, tau_binCenters', &
	380	7, tau_binCenters
	381	CALL wxios_add_vaxis("tau", 7, tau_binCenters)
	382
	383	#endif
	384
	385	zdtimemoy_cosp = freq_COSP ! Frequence ou l on moyenne
	386
	387	DO iff=1,3
	388	zoutm_cosp(iff) = cosp_ecritfiles(iff) ! Frequence ou l on ecrit en seconde
	389
	390	IF (cosp_outfilekeys(iff)) THEN
	391	CALL histbeg_phy_all(cosp_outfilenames(iff),itau_phy,zjulian,&
	392	dtime,nhoricosp(iff),cosp_nidfiles(iff))
	393	! print*,'histbeg_phy nhoricosp(iff),cosp_nidfiles(iff)', &
	394	! nhoricosp(iff),cosp_nidfiles(iff)
	395
	396	#ifdef CPP_XIOS
	397	IF (.not. ok_all_xml) then
	398	WRITE(lunout,*) 'wxios_add_file ',cosp_outfilenames(iff)
	399	CALL wxios_add_file(cosp_outfilenames(iff),chfreq(iff),10)
	400	ENDIF
	401	#endif
	402
	403	#ifndef CPP_IOIPSL_NO_OUTPUT
	404	! Definition de l'axe vertical
	405	if (use_vgrid) then
	406	! Axe vertical Cosp 40 niveaux (en m)
	407	CALL histvert(cosp_nidfiles(iff),"height","height","m",vgrid%Nlvgrid,vgrid%z,nvert(iff))
	408	else
	409	! Axe vertical modele LMDZ presnivs
	410	CALL histvert(cosp_nidfiles(iff),"presnivs","Vertical levels","Pa",vgrid%Nlvgrid,presnivs,nvert(iff),"down")
	411	endif
	412	! Axe vertical niveaux modele (en m)
	413	CALL histvert(cosp_nidfiles(iff),"height_mlev","height_mlev","m",Nlevlmdz,vgrid%mz,nvertmcosp(iff))
	414
	415	CALL histvert(cosp_nidfiles(iff),"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp(iff))
	416
[3369]	417	CALL histvert(cosp_nidfiles(iff),"pressure2","pressure","mb",7,pres_binCenters,nvertisccp(iff),"down")
[3358]	418
	419	CALL histvert(cosp_nidfiles(iff),"column","column","count",Ncolumns,column_ax,nvertcol(iff)) !DBUG
	420
	421	CALL histvert(cosp_nidfiles(iff),"temp","temperature","C",LIDAR_NTEMP,LIDAR_PHASE_TEMP,nverttemp(iff))
	422
[3369]	423	CALL histvert(cosp_nidfiles(iff),"cth","altitude","m",numMISRHgtBins,misr_histHgtCenters,nvertmisr(iff))
[3358]	424
	425	CALL histvert(cosp_nidfiles(iff),"ReffIce","Effective_particle_size_Ice","microns",numMODISReffIceBins, reffICE_binCenters, &
	426	nvertReffIce(iff))
	427
	428	CALL histvert(cosp_nidfiles(iff),"ReffLiq","Effective_particle_size_Liq","microns",numMODISReffLiqBins, reffLIQ_binCenters, &
	429	nvertReffLiq(iff))
	430
	431	CALL histvert(cosp_nidfiles(iff),"dbze","equivalent_reflectivity_factor","dBZ",DBZE_BINS,dbze_ax,nvertbze(iff))
	432
	433	CALL histvert(cosp_nidfiles(iff),"scatratio","backscattering_ratio","1",SR_BINS,sratio_ax,nvertsratio(iff))
	434
[3369]	435	CALL histvert(cosp_nidfiles(iff),"tau","cloud optical depth","1",7,tau_binCenters,nverttau(iff))
[3358]	436
	437	!!! Valeur indefinie en cas IOIPSL
	438	Cosp_fill_value=0.
	439
	440	#endif
	441
	442	ENDIF
	443	ENDDO
	444
	445	end SUBROUTINE cosp_output_open
	446
	447	END MODULE cosp_output_mod

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

Download in other formats:

Original Format