Context Navigation

phyetat0_mod.F90 @ 5117

Last change on this file since 5117 was 5117, checked in by abarral, 4 months ago
rename modules properly lmdz_* move some unused files to obsolete/ (lint) uppercase fortran keywords
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 24.6 KB

Rev	Line
[1403]	1	! $Id: phyetat0_mod.F90 5117 2024-07-24 14:23:34Z abarral $
[782]	2
[4358]	3	MODULE phyetat0_mod
[5111]	4	USE lmdz_abort_physic, ONLY: abort_physic
[4358]	5	PRIVATE
[4367]	6	PUBLIC :: phyetat0
[4358]	7
	8	CONTAINS
	9
[5111]	10	SUBROUTINE phyetat0(fichnom, clesphy0, tabcntr0)
[1279]	11
[5116]	12	USE dimphy, ONLY: klon, zmasq, klev
[5111]	13	USE iophy, ONLY: init_iophy_new
	14	USE ocean_cpl_mod, ONLY: ocean_cpl_init
	15	USE fonte_neige_mod, ONLY: fonte_neige_init
	16	USE pbl_surface_mod, ONLY: pbl_surface_init
	17	USE surface_data, ONLY: type_ocean, version_ocean
	18	USE phyetat0_get_mod, ONLY: phyetat0_get, phyetat0_srf
	19	USE phys_state_var_mod, ONLY: ancien_ok, clwcon, detr_therm, phys_tstep, &
	20	qsol, fevap, z0m, z0h, agesno, &
	21	du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
	22	falb_dir, falb_dif, prw_ancien, prlw_ancien, prsw_ancien, prbsw_ancien, &
	23	ftsol, pbl_tke, pctsrf, q_ancien, ql_ancien, qs_ancien, qbs_ancien, rneb_ancien, radpas, radsol, rain_fall, ratqs, &
	24	rnebcon, rugoro, sig1, snow_fall, bs_fall, solaire_etat0, sollw, sollwdown, &
	25	solsw, solswfdiff, t_ancien, u_ancien, v_ancien, w01, wake_cstar, wake_deltaq, &
	26	wake_deltat, wake_delta_pbl_TKE, delta_tsurf, beta_aridity, wake_fip, wake_pe, &
	27	wake_s, awake_s, wake_dens, awake_dens, cv_gen, zgam, zmax0, zmea, zpic, zsig, &
	28	zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, u10m, v10m, treedrg, &
	29	ale_wake, ale_bl_stat, ds_ns, dt_ns, delta_sst, delta_sal, dter, dser, &
	30	dt_ds, ratqs_inter_
	31	!FC
[5112]	32	USE lmdz_geometry, ONLY: longitude_deg, latitude_deg
[5111]	33	USE iostart, ONLY: close_startphy, get_field, get_var, open_startphy
	34	USE infotrac_phy, ONLY: nqtot, nbtr, type_trac, tracers
[5117]	35	USE lmdz_readTracFiles, ONLY: maxlen, new2oldH2O
[5111]	36	USE traclmdz_mod, ONLY: traclmdz_from_restart
	37	USE carbon_cycle_mod, ONLY: carbon_cycle_init, carbon_cycle_cpl, carbon_cycle_tr, carbon_cycle_rad, co2_send, RCO2_glo
	38	USE indice_sol_mod, ONLY: nbsrf, is_ter, epsfra, is_lic, is_oce, is_sic
	39	USE ocean_slab_mod, ONLY: nslay, tslab, seaice, tice, ocean_slab_init
	40	USE time_phylmdz_mod, ONLY: init_iteration, pdtphys, itau_phy
[5117]	41	USE lmdz_wxios, ONLY: missing_val_xios => missing_val, using_xios
	42	USE netcdf, ONLY: missing_val_netcdf => nf90_fill_real
	43	USE config_ocean_skin_m, ONLY: activate_ocean_skin
[967]	44
[5111]	45	IMPLICIT none
	46	!======================================================================
	47	! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
	48	! Objet: Lecture de l'etat initial pour la physique
	49	!======================================================================
	50	include "dimsoil.h"
	51	include "clesphys.h"
	52	include "alpale.h"
	53	include "compbl.h"
	54	include "YOMCST.h"
	55	!======================================================================
	56	CHARACTER() fichnom
[524]	57
[5111]	58	! les variables globales lues dans le fichier restart
[1001]	59
[5111]	60	REAL tsoil(klon, nsoilmx, nbsrf)
	61	REAL qsurf(klon, nbsrf)
	62	REAL snow(klon, nbsrf)
[5117]	63	REAL fder(klon)
[5111]	64	REAL run_off_lic_0(klon)
	65	REAL fractint(klon)
	66	REAL trs(klon, nbtr)
	67	REAL zts(klon)
	68	! pour drag arbres FC
	69	REAL drg_ter(klon, klev)
[651]	70
[5111]	71	CHARACTER*6 ocean_in
	72	LOGICAL ok_veget_in
[879]	73
[5111]	74	INTEGER longcles
	75	PARAMETER (longcles = 20)
	76	REAL clesphy0(longcles)
[766]	77
[5111]	78	REAL xmin, xmax
[766]	79
[5111]	80	INTEGER nid, nvarid
	81	INTEGER ierr, i, nsrf, isoil, k
	82	INTEGER length
	83	PARAMETER (length = 100)
	84	INTEGER it, iq, isw
	85	REAL tab_cntrl(length), tabcntr0(length)
	86	CHARACTER*7 str7
	87	CHARACTER*2 str2
	88	LOGICAL :: found
	89	REAL :: lon_startphy(klon), lat_startphy(klon)
	90	CHARACTER(LEN = maxlen) :: tname, t(2)
	91	REAL :: missing_val
[1827]	92
[5111]	93	IF (using_xios) THEN
	94	missing_val = missing_val_xios
	95	ELSE
	96	missing_val = missing_val_netcdf
	97	ENDIF
[1827]	98
[5111]	99	! FH1D
	100	! real iolat(jjm+1)
	101	!real iolat(jjm+1-1/(iim*jjm))
[1827]	102
[5111]	103	! Ouvrir le fichier contenant l'etat initial:
[1827]	104
[5111]	105	CALL open_startphy(fichnom)
[1827]	106
[5111]	107	! Lecture des parametres de controle:
[1827]	108
[5111]	109	CALL get_var("controle", tab_cntrl)
[956]	110
[5111]	111	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	112	! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
	113	! Les constantes de la physiques sont lues dans la physique seulement.
	114	! Les egalites du type
	115	! tab_cntrl( 5 )=clesphy0(1)
	116	! sont remplacees par
	117	! clesphy0(1)=tab_cntrl( 5 )
	118	! On inverse aussi la logique.
	119	! On remplit les tab_cntrl avec les parametres lus dans les .def
	120	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1279]	121
[5111]	122	DO i = 1, length
	123	tabcntr0(i) = tab_cntrl(i)
	124	ENDDO
[1279]	125
[5111]	126	tab_cntrl(1) = pdtphys
	127	tab_cntrl(2) = radpas
[4298]	128
[5111]	129	! co2_ppm : value from the previous time step
[4298]	130
[5111]	131	! co2_ppm0 : initial value of atmospheric CO2 (from create_etat0_limit.e .def)
	132	co2_ppm0 = 284.32
	133	! when no initial value is available e.g., from a restart
	134	! this variable must be set in a .def file which will then be
	135	! used by the conf_phys_m.F90 routine.
	136	! co2_ppm0 = 284.32 (illustrative example on how to set the variable in .def
	137	! file, for a pre-industrial CO2 concentration value)
[1279]	138
[5111]	139	IF (carbon_cycle_tr .OR. carbon_cycle_cpl) THEN
	140	co2_ppm = tab_cntrl(3)
	141	RCO2 = co2_ppm * 1.0e-06 * RMCO2 / RMD
	142	IF (tab_cntrl(17) > 0. .AND. carbon_cycle_rad) THEN
	143	RCO2_glo = tab_cntrl(17)
	144	ELSE
	145	RCO2_glo = co2_ppm0 * 1.0e-06 * RMCO2 / RMD
	146	ENDIF
	147	! ELSE : keep value from .def
	148	ENDIF
[524]	149
[5111]	150	solaire_etat0 = tab_cntrl(4)
	151	tab_cntrl(5) = iflag_con
	152	tab_cntrl(6) = nbapp_rad
[956]	153
[5111]	154	IF (iflag_cycle_diurne>=1) tab_cntrl(7) = iflag_cycle_diurne
	155	IF (soil_model) tab_cntrl(8) = 1.
	156	IF (new_oliq) tab_cntrl(9) = 1.
	157	IF (ok_orodr) tab_cntrl(10) = 1.
	158	IF (ok_orolf) tab_cntrl(11) = 1.
	159	IF (ok_limitvrai) tab_cntrl(12) = 1.
[956]	160
[5111]	161	itau_phy = tab_cntrl(15)
[956]	162
[5111]	163	clesphy0(1) = tab_cntrl(5)
	164	clesphy0(2) = tab_cntrl(6)
	165	clesphy0(3) = tab_cntrl(7)
	166	clesphy0(4) = tab_cntrl(8)
	167	clesphy0(5) = tab_cntrl(9)
	168	clesphy0(6) = tab_cntrl(10)
	169	clesphy0(7) = tab_cntrl(11)
	170	clesphy0(8) = tab_cntrl(12)
	171	clesphy0(9) = tab_cntrl(17)
[2344]	172
[5111]	173	! set time iteration
	174	CALL init_iteration(itau_phy)
[766]	175
[5111]	176	! read latitudes and make a sanity check (because already known from dyn)
	177	CALL get_field("latitude", lat_startphy)
	178	DO i = 1, klon
	179	IF (ABS(lat_startphy(i) - latitude_deg(i))>=1) THEN
	180	WRITE(, ) "phyetat0: Error! Latitude discrepancy wrt startphy file:", &
	181	" i=", i, " lat_startphy(i)=", lat_startphy(i), &
	182	" latitude_deg(i)=", latitude_deg(i)
[3505]	183	! This is presumably serious enough to abort run
[5111]	184	CALL abort_physic("phyetat0", "discrepancy in latitudes!", 1)
[3505]	185	ENDIF
[5111]	186	IF (ABS(lat_startphy(i) - latitude_deg(i))>=0.0001) THEN
	187	WRITE(, ) "phyetat0: Warning! Latitude discrepancy wrt startphy file:", &
	188	" i=", i, " lat_startphy(i)=", lat_startphy(i), &
	189	" latitude_deg(i)=", latitude_deg(i)
	190	ENDIF
	191	ENDDO
[1001]	192
[5111]	193	! read longitudes and make a sanity check (because already known from dyn)
	194	CALL get_field("longitude", lon_startphy)
	195	DO i = 1, klon
	196	IF (ABS(lon_startphy(i) - longitude_deg(i))>=1) THEN
	197	IF (ABS(360 - ABS(lon_startphy(i) - longitude_deg(i)))>=1) THEN
	198	WRITE(, ) "phyetat0: Error! Longitude discrepancy wrt startphy file:", &
	199	" i=", i, " lon_startphy(i)=", lon_startphy(i), &
	200	" longitude_deg(i)=", longitude_deg(i)
	201	! This is presumably serious enough to abort run
	202	CALL abort_physic("phyetat0", "discrepancy in longitudes!", 1)
	203	ENDIF
	204	ENDIF
	205	IF (ABS(lon_startphy(i) - longitude_deg(i))>=1) THEN
	206	IF (ABS(360 - ABS(lon_startphy(i) - longitude_deg(i))) > 0.0001) THEN
	207	WRITE(, ) "phyetat0: Warning! Longitude discrepancy wrt startphy file:", &
	208	" i=", i, " lon_startphy(i)=", lon_startphy(i), &
	209	" longitude_deg(i)=", longitude_deg(i)
	210	ENDIF
	211	ENDIF
	212	ENDDO
[766]	213
[5111]	214	! Lecture du masque terre mer
[1001]	215
[5111]	216	CALL get_field("masque", zmasq, found)
	217	IF (.NOT. found) THEN
	218	PRINT*, 'phyetat0: Le champ <masque> est absent'
	219	PRINT *, 'fichier startphy non compatible avec phyetat0'
	220	ENDIF
[766]	221
[5111]	222	! Lecture des fractions pour chaque sous-surface
[766]	223
[5111]	224	! initialisation des sous-surfaces
[766]	225
[5111]	226	pctsrf = 0.
[766]	227
[5111]	228	! fraction de terre
[766]	229
[5111]	230	CALL get_field("FTER", pctsrf(:, is_ter), found)
	231	IF (.NOT. found) PRINT*, 'phyetat0: Le champ <FTER> est absent'
[766]	232
[5111]	233	! fraction de glace de terre
[1001]	234
[5111]	235	CALL get_field("FLIC", pctsrf(:, is_lic), found)
	236	IF (.NOT. found) PRINT*, 'phyetat0: Le champ <FLIC> est absent'
[1001]	237
[5111]	238	! fraction d'ocean
[1001]	239
[5111]	240	CALL get_field("FOCE", pctsrf(:, is_oce), found)
	241	IF (.NOT. found) PRINT*, 'phyetat0: Le champ <FOCE> est absent'
[1001]	242
[5111]	243	! fraction glace de mer
[1001]	244
[5111]	245	CALL get_field("FSIC", pctsrf(:, is_sic), found)
	246	IF (.NOT. found) PRINT*, 'phyetat0: Le champ <FSIC> est absent'
[1827]	247
[5111]	248	! Verification de l'adequation entre le masque et les sous-surfaces
	249
	250	fractint(1:klon) = pctsrf(1:klon, is_ter) &
	251	+ pctsrf(1:klon, is_lic)
	252	DO i = 1, klon
	253	IF (abs(fractint(i) - zmasq(i)) > EPSFRA) THEN
	254	WRITE(, ) 'phyetat0: attention fraction terre pas ', &
	255	'coherente ', i, zmasq(i), pctsrf(i, is_ter) &
	256	, pctsrf(i, is_lic)
[1827]	257	WRITE(, ) 'Je force la coherence zmasq=fractint'
	258	zmasq(i) = fractint(i)
[5111]	259	ENDIF
	260	ENDDO
	261	fractint (1:klon) = pctsrf(1:klon, is_oce) &
	262	+ pctsrf(1:klon, is_sic)
	263	DO i = 1, klon
	264	IF (abs(fractint(i) - (1. - zmasq(i))) > EPSFRA) THEN
	265	WRITE(, ) 'phyetat0 attention fraction ocean pas ', &
	266	'coherente ', i, zmasq(i), pctsrf(i, is_oce) &
	267	, pctsrf(i, is_sic)
[2053]	268	WRITE(, ) 'Je force la coherence zmasq=1.-fractint'
[2052]	269	zmasq(i) = 1. - fractint(i)
[5111]	270	ENDIF
	271	ENDDO
[1827]	272
[5111]	273	!===================================================================
	274	! Lecture des temperatures du sol:
	275	!===================================================================
[1827]	276
[5111]	277	found = phyetat0_get(ftsol(:, 1), "TS", "Surface temperature", 283.)
	278	IF (found) THEN
	279	DO nsrf = 2, nbsrf
	280	ftsol(:, nsrf) = ftsol(:, 1)
	281	ENDDO
	282	ELSE
	283	found = phyetat0_srf(ftsol, "TS", "Surface temperature", 283.)
	284	ENDIF
[524]	285
[5111]	286	!===================================================================
	287	! Lecture des albedo difus et direct
	288	!===================================================================
[2237]	289
[5111]	290	DO nsrf = 1, nbsrf
	291	DO isw = 1, nsw
[5075]	292	IF (isw>99) THEN
[5111]	293	PRINT*, "Trop de bandes SW"
	294	CALL abort_physic("phyetat0", "", 1)
[2237]	295	ENDIF
[2252]	296	WRITE(str2, '(i2.2)') isw
[5111]	297	found = phyetat0_srf(falb_dir(:, isw, :), "A_dir_SW" // str2 // "srf", "Direct Albedo", 0.2)
	298	found = phyetat0_srf(falb_dif(:, isw, :), "A_dif_SW" // str2 // "srf", "Direct Albedo", 0.2)
	299	ENDDO
	300	ENDDO
[2237]	301
[5111]	302	found = phyetat0_srf(u10m, "U10M", "u a 10m", 0.)
	303	found = phyetat0_srf(v10m, "V10M", "v a 10m", 0.)
[2569]	304
[5111]	305	!===================================================================
	306	! Lecture dans le cas iflag_pbl_surface =1
	307	!===================================================================
[4537]	308
[5117]	309	IF (iflag_physiq <= 1) THEN
[5111]	310	!===================================================================
	311	! Lecture des temperatures du sol profond:
	312	!===================================================================
[524]	313
[5111]	314	DO isoil = 1, nsoilmx
[5075]	315	IF (isoil>99) THEN
[5111]	316	PRINT*, "Trop de couches "
	317	CALL abort_physic("phyetat0", "", 1)
[1827]	318	ENDIF
[5111]	319	WRITE(str2, '(i2.2)') isoil
	320	found = phyetat0_srf(tsoil(:, isoil, :), "Tsoil" // str2 // "srf", "Temp soil", 0.)
[1827]	321	IF (.NOT. found) THEN
[5111]	322	PRINT*, "phyetat0: Le champ <Tsoil" // str7 // "> est absent"
	323	PRINT*, " Il prEND DOnc la valeur de surface"
	324	tsoil(:, isoil, :) = ftsol(:, :)
[1827]	325	ENDIF
[5111]	326	ENDDO
[524]	327
[5111]	328	!=======================================================================
	329	! Lecture precipitation/evaporation
	330	!=======================================================================
[1001]	331
[5111]	332	found = phyetat0_srf(qsurf, "QS", "Near surface hmidity", 0.)
	333	found = phyetat0_get(qsol, "QSOL", "Surface hmidity / bucket", 0.)
	334	found = phyetat0_srf(snow, "SNOW", "Surface snow", 0.)
	335	found = phyetat0_srf(fevap, "EVAP", "evaporation", 0.)
	336	found = phyetat0_get(snow_fall, "snow_f", "snow fall", 0.)
	337	found = phyetat0_get(rain_fall, "rain_f", "rain fall", 0.)
	338	IF (ok_bs) THEN
	339	found = phyetat0_get(bs_fall, "bs_f", "blowing snow fall", 0.)
	340	ELSE
	341	bs_fall(:) = 0.
	342	ENDIF
	343	!=======================================================================
	344	! Radiation
	345	!=======================================================================
[1001]	346
[5111]	347	found = phyetat0_get(solsw, "solsw", "net SW radiation surf", 0.)
	348	found = phyetat0_get(solswfdiff, "solswfdiff", "fraction of SW radiation surf that is diffuse", 1.)
	349	found = phyetat0_get(sollw, "sollw", "net LW radiation surf", 0.)
	350	found = phyetat0_get(sollwdown, "sollwdown", "down LW radiation surf", 0.)
	351	IF (.NOT. found) THEN
	352	sollwdown(:) = 0. ; zts(:) = 0.
	353	DO nsrf = 1, nbsrf
	354	zts(:) = zts(:) + ftsol(:, nsrf) * pctsrf(:, nsrf)
	355	ENDDO
	356	sollwdown(:) = sollw(:) + RSIGMA * zts(:)**4
	357	ENDIF
[2188]	358
[5111]	359	found = phyetat0_get(radsol, "RADS", "Solar radiation", 0.)
	360	found = phyetat0_get(fder, "fder", "Flux derivative", 0.)
[2188]	361
[1827]	362
[5111]	363	! Lecture de la longueur de rugosite
	364	found = phyetat0_srf(z0m, "RUG", "Z0m ancien", 0.001)
	365	IF (found) THEN
	366	z0h(:, 1:nbsrf) = z0m(:, 1:nbsrf)
	367	ELSE
	368	found = phyetat0_srf(z0m, "Z0m", "Roughness length, momentum ", 0.001)
	369	found = phyetat0_srf(z0h, "Z0h", "Roughness length, enthalpy ", 0.001)
	370	ENDIF
	371	!FC
	372	IF (ifl_pbltree>0) THEN
	373	!CALL get_field("FTER", pctsrf(:, is_ter), found)
	374	treedrg(:, 1:klev, 1:nbsrf) = 0.0
	375	CALL get_field("treedrg_ter", drg_ter(:, :), found)
	376	! found=phyetat0_srf(treedrg,"treedrg","drag from vegetation" , 0.)
	377	!lecture du profile de freinage des arbres
[5117]	378	IF (.NOT. found) THEN
[5111]	379	treedrg(:, 1:klev, 1:nbsrf) = 0.0
	380	ELSE
	381	treedrg(:, 1:klev, is_ter) = drg_ter(:, :)
	382	! found=phyetat0_get(treedrg,"treedrg","freinage arbres",0.)
	383	ENDIF
	384	ELSE
	385	! initialize treedrg(), because it will be written in restartphy.nc
	386	treedrg(:, :, :) = 0.0
	387	ENDIF
[1827]	388
[5111]	389	endif ! iflag_physiq <= 1
[4537]	390
[5111]	391	! Lecture de l'age de la neige:
	392	found = phyetat0_srf(agesno, "AGESNO", "SNOW AGE", 0.001)
[1827]	393
[5111]	394	ancien_ok = .TRUE.
	395	ancien_ok = ancien_ok.AND.phyetat0_get(t_ancien, "TANCIEN", "TANCIEN", 0.)
	396	ancien_ok = ancien_ok.AND.phyetat0_get(q_ancien, "QANCIEN", "QANCIEN", 0.)
	397	ancien_ok = ancien_ok.AND.phyetat0_get(ql_ancien, "QLANCIEN", "QLANCIEN", 0.)
	398	ancien_ok = ancien_ok.AND.phyetat0_get(qs_ancien, "QSANCIEN", "QSANCIEN", 0.)
	399	ancien_ok = ancien_ok.AND.phyetat0_get(rneb_ancien, "RNEBANCIEN", "RNEBANCIEN", 0.)
	400	ancien_ok = ancien_ok.AND.phyetat0_get(u_ancien, "UANCIEN", "UANCIEN", 0.)
	401	ancien_ok = ancien_ok.AND.phyetat0_get(v_ancien, "VANCIEN", "VANCIEN", 0.)
	402	ancien_ok = ancien_ok.AND.phyetat0_get(prw_ancien, "PRWANCIEN", "PRWANCIEN", 0.)
	403	ancien_ok = ancien_ok.AND.phyetat0_get(prlw_ancien, "PRLWANCIEN", "PRLWANCIEN", 0.)
	404	ancien_ok = ancien_ok.AND.phyetat0_get(prsw_ancien, "PRSWANCIEN", "PRSWANCIEN", 0.)
[1827]	405
[5111]	406	! cas specifique des variables de la neige soufflee
	407	IF (ok_bs) THEN
	408	ancien_ok = ancien_ok.AND.phyetat0_get(qbs_ancien, "QBSANCIEN", "QBSANCIEN", 0.)
	409	ancien_ok = ancien_ok.AND.phyetat0_get(prbsw_ancien, "PRBSWANCIEN", "PRBSWANCIEN", 0.)
	410	ELSE
	411	qbs_ancien(:, :) = 0.
	412	prbsw_ancien(:) = 0.
	413	ENDIF
[4578]	414
[5111]	415	! Ehouarn: addtional tests to check if t_ancien, q_ancien contain
	416	! dummy values (as is the case when generated by ce0l,
	417	! or by iniaqua)
	418	IF ((maxval(q_ancien)==minval(q_ancien)) .OR. &
	419	(maxval(ql_ancien)==minval(ql_ancien)) .OR. &
	420	(maxval(qs_ancien)==minval(qs_ancien)) .OR. &
	421	(maxval(rneb_ancien)==minval(rneb_ancien)) .OR. &
	422	(maxval(prw_ancien)==minval(prw_ancien)) .OR. &
	423	(maxval(prlw_ancien)==minval(prlw_ancien)) .OR. &
	424	(maxval(prsw_ancien)==minval(prsw_ancien)) .OR. &
	425	(maxval(t_ancien)==minval(t_ancien))) THEN
	426	ancien_ok = .FALSE.
	427	ENDIF
[2494]	428
[5111]	429	IF (ok_bs) THEN
	430	IF ((maxval(qbs_ancien)==minval(qbs_ancien)) .OR. &
	431	(maxval(prbsw_ancien)==minval(prbsw_ancien))) THEN
	432	ancien_ok = .FALSE.
	433	ENDIF
[4579]	434	ENDIF
	435
[5111]	436	found = phyetat0_get(clwcon, "CLWCON", "CLWCON", 0.)
	437	found = phyetat0_get(rnebcon, "RNEBCON", "RNEBCON", 0.)
	438	found = phyetat0_get(ratqs, "RATQS", "RATQS", 0.)
[1827]	439
[5111]	440	found = phyetat0_get(run_off_lic_0, "RUNOFFLIC0", "RUNOFFLIC0", 0.)
[1827]	441
[5111]	442	!==================================
	443	! TKE
	444	!==================================
[5099]	445
[5116]	446	IF (iflag_pbl>1) THEN
[5111]	447	found = phyetat0_srf(pbl_tke, "TKE", "Turb. Kinetic. Energ. ", 1.e-8)
	448	ENDIF
[1403]	449
[5116]	450	IF (iflag_pbl>1 .AND. iflag_wake>=1 .AND. iflag_pbl_split >=1) THEN
[5111]	451	found = phyetat0_srf(wake_delta_pbl_tke, "DELTATKE", "Del TKE wk/env", 0.)
	452	!! found=phyetat0_srf(delta_tsurf,"DELTA_TSURF","Delta Ts wk/env ",0.)
	453	found = phyetat0_srf(delta_tsurf, "DELTATS", "Delta Ts wk/env ", 0.)
	454	!! found=phyetat0_srf(beta_aridity,"BETA_S","Aridity factor ",1.)
	455	found = phyetat0_srf(beta_aridity, "BETAS", "Aridity factor ", 1.)
	456	ENDIF !(iflag_pbl>1 .AND. iflag_wake>=1 .AND. iflag_pbl_split >=1 )
[2159]	457
[5111]	458	!==================================
	459	! thermiques, poches, convection
	460	!==================================
[1403]	461
[5111]	462	! Emanuel
	463	found = phyetat0_get(sig1, "sig1", "sig1", 0.)
	464	found = phyetat0_get(w01, "w01", "w01", 0.)
[1403]	465
[5111]	466	! Wake
	467	found = phyetat0_get(wake_deltat, "WAKE_DELTAT", "Delta T wake/env", 0.)
	468	found = phyetat0_get(wake_deltaq, "WAKE_DELTAQ", "Delta hum. wake/env", 0.)
	469	found = phyetat0_get(wake_s, "WAKE_S", "Wake frac. area", 0.)
	470	found = phyetat0_get(awake_s, "AWAKE_S", "Active Wake frac. area", 0.)
	471	!jyg<
	472	! Set wake_dens to -1000. when there is no restart so that the actual
	473	! initialization is made in calwake.
	474	!! found=phyetat0_get(1,wake_dens,"WAKE_DENS","Wake num. /unit area",0.)
	475	found = phyetat0_get(wake_dens, "WAKE_DENS", "Wake num. /unit area", -1000.)
	476	found = phyetat0_get(awake_dens, "AWAKE_DENS", "Active Wake num. /unit area", 0.)
	477	found = phyetat0_get(cv_gen, "CV_GEN", "CB birth rate", 0.)
	478	!>jyg
	479	found = phyetat0_get(wake_cstar, "WAKE_CSTAR", "WAKE_CSTAR", 0.)
	480	found = phyetat0_get(wake_pe, "WAKE_PE", "WAKE_PE", 0.)
	481	found = phyetat0_get(wake_fip, "WAKE_FIP", "WAKE_FIP", 0.)
[879]	482
[5111]	483	! Thermiques
	484	found = phyetat0_get(zmax0, "ZMAX0", "ZMAX0", 40.)
	485	found = phyetat0_get(f0, "F0", "F0", 1.e-5)
	486	found = phyetat0_get(fm_therm, "FM_THERM", "Thermals mass flux", 0.)
	487	found = phyetat0_get(entr_therm, "ENTR_THERM", "Thermals Entrain.", 0.)
	488	found = phyetat0_get(detr_therm, "DETR_THERM", "Thermals Detrain.", 0.)
[782]	489
[5111]	490	! ALE/ALP
	491	found = phyetat0_get(ale_bl, "ALE_BL", "ALE BL", 0.)
	492	found = phyetat0_get(ale_bl_trig, "ALE_BL_TRIG", "ALE BL_TRIG", 0.)
	493	found = phyetat0_get(alp_bl, "ALP_BL", "ALP BL", 0.)
	494	found = phyetat0_get(ale_wake, "ALE_WAKE", "ALE_WAKE", 0.)
	495	found = phyetat0_get(ale_bl_stat, "ALE_BL_STAT", "ALE_BL_STAT", 0.)
[1279]	496
[5111]	497	! fisrtilp/Clouds 0.002 could be ratqsbas. But can stay like this as well
	498	found = phyetat0_get(ratqs_inter_, "RATQS_INTER", "Relative width of the lsc sugrid scale water", 0.002)
[3856]	499
[5111]	500	!===========================================
	501	! Read and send field trs to traclmdz
	502	!===========================================
[1827]	503
[5111]	504	!--OB now this is for co2i - ThL: and therefore also for inco
	505	IF (ANY(type_trac == ['co2i', 'inco'])) THEN
	506	IF (carbon_cycle_cpl) THEN
	507	ALLOCATE(co2_send(klon), stat = ierr)
[4170]	508	IF (ierr /= 0) CALL abort_physic('phyetat0', 'pb allocation co2_send', 1)
[5111]	509	found = phyetat0_get(co2_send, "co2_send", "co2 send", co2_ppm0)
	510	ENDIF
	511	ELSE IF (type_trac == 'lmdz') THEN
	512	it = 0
	513	DO iq = 1, nqtot
[4071]	514	IF(.NOT.(tracers(iq)%isAdvected .AND. tracers(iq)%isInPhysics)) CYCLE
[5111]	515	it = it + 1
[4358]	516	tname = tracers(iq)%name
[5111]	517	t(1) = 'trs_' // TRIM(tname); t(2) = 'trs_' // TRIM(new2oldH2O(tname))
	518	found = phyetat0_get(trs(:, it), t(:), "Surf trac" // TRIM(tname), 0.)
	519	END DO
	520	CALL traclmdz_from_restart(trs)
	521	ENDIF
[1827]	522
	523
[5111]	524	!===========================================
	525	! ondes de gravite / relief
	526	!===========================================
[2251]	527
[5111]	528	! ondes de gravite non orographiques
	529	IF (ok_gwd_rando) found = &
	530	phyetat0_get(du_gwd_rando, "du_gwd_rando", "du_gwd_rando", 0.)
	531	IF (.NOT. ok_hines .AND. ok_gwd_rando) found &
	532	= phyetat0_get(du_gwd_front, "du_gwd_front", "du_gwd_front", 0.)
[1938]	533
[5111]	534	! prise en compte du relief sous-maille
	535	found = phyetat0_get(zmea, "ZMEA", "sub grid orography", 0.)
	536	found = phyetat0_get(zstd, "ZSTD", "sub grid orography", 0.)
	537	found = phyetat0_get(zsig, "ZSIG", "sub grid orography", 0.)
	538	found = phyetat0_get(zgam, "ZGAM", "sub grid orography", 0.)
	539	found = phyetat0_get(zthe, "ZTHE", "sub grid orography", 0.)
	540	found = phyetat0_get(zpic, "ZPIC", "sub grid orography", 0.)
	541	found = phyetat0_get(zval, "ZVAL", "sub grid orography", 0.)
	542	found = phyetat0_get(zmea, "ZMEA", "sub grid orography", 0.)
	543	found = phyetat0_get(rugoro, "RUGSREL", "sub grid orography", 0.)
[2252]	544
[5111]	545	!===========================================
	546	! Initialize ocean
	547	!===========================================
[2251]	548
[5111]	549	IF (type_ocean == 'slab') THEN
[3435]	550	CALL ocean_slab_init(phys_tstep, pctsrf)
[5075]	551	IF (nslay==1) THEN
[5111]	552	found = phyetat0_get(tslab, ["tslab01", "tslab "], "tslab", 0.)
[2656]	553	ELSE
[5111]	554	DO i = 1, nslay
	555	WRITE(str2, '(i2.2)') i
	556	found = phyetat0_get(tslab(:, i), "tslab" // str2, "tslab", 0.)
	557	ENDDO
[3449]	558	ENDIF
[5111]	559	IF (.NOT. found) THEN
	560	PRINT*, "phyetat0: Le champ <tslab> est absent"
	561	PRINT*, "Initialisation a tsol_oce"
	562	DO i = 1, nslay
	563	tslab(:, i) = MAX(ftsol(:, is_oce), 271.35)
	564	ENDDO
	565	ENDIF
[2251]	566
[2209]	567	! Sea ice variables
	568	IF (version_ocean == 'sicINT') THEN
[5111]	569	found = phyetat0_get(tice, "slab_tice", "slab_tice", 0.)
	570	IF (.NOT. found) THEN
	571	PRINT*, "phyetat0: Le champ <tice> est absent"
	572	PRINT*, "Initialisation a tsol_sic"
	573	tice(:) = ftsol(:, is_sic)
	574	ENDIF
	575	found = phyetat0_get(seaice, "seaice", "seaice", 0.)
	576	IF (.NOT. found) THEN
	577	PRINT*, "phyetat0: Le champ <seaice> est absent"
	578	PRINT*, "Initialisation a 0/1m suivant fraction glace"
	579	seaice(:) = 0.
	580	WHERE (pctsrf(:, is_sic)>EPSFRA)
	581	seaice = 917.
	582	ENDWHERE
	583	ENDIF
[3449]	584	ENDIF !sea ice INT
[5111]	585	ENDIF ! Slab
[2057]	586
[5117]	587	IF (activate_ocean_skin >= 1) THEN
	588	IF (activate_ocean_skin == 2 .AND. type_ocean == 'couple') THEN
[4358]	589	found = phyetat0_get(delta_sal, "delta_sal", &
[5111]	590	"ocean-air interface salinity minus bulk salinity", 0.)
[4358]	591	found = phyetat0_get(delta_sst, "delta_SST", &
[5111]	592	"ocean-air interface temperature minus bulk SST", 0.)
[4370]	593	found = phyetat0_get(dter, "dter", &
[5111]	594	"ocean-air interface temperature minus subskin temperature", 0.)
[4370]	595	found = phyetat0_get(dser, "dser", &
[5111]	596	"ocean-air interface salinity minus subskin salinity", 0.)
[4370]	597	found = phyetat0_get(dt_ds, "dt_ds", "(tks / tkt) * dTer", 0.)
	598
	599	where (pctsrf(:, is_oce) == 0.)
[5111]	600	delta_sst = missing_val
	601	delta_sal = missing_val
	602	dter = missing_val
	603	dser = missing_val
	604	dt_ds = missing_val
[4370]	605	end where
[5111]	606	end if
[3815]	607
[5111]	608	found = phyetat0_get(ds_ns, "dS_ns", "delta salinity near surface", 0.)
	609	found = phyetat0_get(dt_ns, "dT_ns", "delta temperature near surface", &
	610	0.)
	611
	612	where (pctsrf(:, is_oce) == 0.)
[3815]	613	ds_ns = missing_val
	614	dt_ns = missing_val
	615	delta_sst = missing_val
	616	delta_sal = missing_val
[5111]	617	end where
	618	end if
[3815]	619
[5111]	620	! on ferme le fichier
	621	CALL close_startphy
[1827]	622
[5111]	623	! Initialize module pbl_surface_mod
[1827]	624
[5117]	625	IF (iflag_physiq <= 1) THEN
[5111]	626	CALL pbl_surface_init(fder, snow, qsurf, tsoil)
	627	endif
[1827]	628
[5111]	629	! Initialize module ocean_cpl_mod for the case of coupled ocean
	630	IF (type_ocean == 'couple') THEN
	631	CALL ocean_cpl_init(phys_tstep, longitude_deg, latitude_deg)
	632	ENDIF
[1827]	633
[5111]	634	! CALL init_iophy_new(latitude_deg, longitude_deg)
[2054]	635
[5111]	636	! Initilialize module fonte_neige_mod
	637	CALL fonte_neige_init(run_off_lic_0)
[1827]	638
[5111]	639	END SUBROUTINE phyetat0
[2243]	640
[4358]	641	END MODULE phyetat0_mod
[2243]	642

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

Download in other formats:

Original Format