Index: LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h =================================================================== --- LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h (revision 2681) +++ LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h (revision 2683) @@ -55,7 +55,7 @@ !Config Key = prt_level -!Config Desc = niveau d'impressions de débogage +!Config Desc = niveau d'impressions de d?bogage !Config Def = 0 -!Config Help = Niveau d'impression pour le débogage +!Config Help = Niveau d'impression pour le d?bogage !Config (0 = minimum d'impression) ! prt_level = 0 @@ -134,5 +134,5 @@ ENDIF -!Paramètres de forçage +!Param?tres de for?age !Config Key = tend_t !Config Desc = forcage ou non par advection de T @@ -1123,6 +1123,6 @@ !---------------------------------------------------------------------- ! Calcul de l'advection verticale (ascendance et subsidence) de -! température et d'humidité. Hypothèse : ce qui rentre de l'extérieur -! a les mêmes caractéristiques que l'air de la colonne 1D (WTG) ou +! temp?rature et d'humidit?. Hypoth?se : ce qui rentre de l'ext?rieur +! a les m?mes caract?ristiques que l'air de la colonne 1D (WTG) ou ! sans WTG rajouter une advection horizontale !---------------------------------------------------------------------- @@ -1197,6 +1197,6 @@ !---------------------------------------------------------------------- ! Calcul de l'advection verticale (ascendance et subsidence) de -! température et d'humidité. Hypothèse : ce qui rentre de l'extérieur -! a les mêmes caractéristiques que l'air de la colonne 1D (WTG) ou +! temp?rature et d'humidit?. Hypoth?se : ce qui rentre de l'ext?rieur +! a les m?mes caract?ristiques que l'air de la colonne 1D (WTG) ou ! sans WTG rajouter une advection horizontale !---------------------------------------------------------------------- @@ -2951,5 +2951,5 @@ endif if (annee_ref.eq.1992 .and. day1.lt.day_ini_toga) then - print*,'TOGA-COARE a débuté le 1er Nov 1992 (jour julien=306)' + print*,'TOGA-COARE a d?but? le 1er Nov 1992 (jour julien=306)' print*,'Changer dayref dans run.def' stop @@ -3783,5 +3783,5 @@ !===================================================================== subroutine read_dice(fich_dice,nlevel,ntime & - & ,zz,pres,th,qv,u,v,o3 & + & ,zz,pres,t,qv,u,v,o3 & & ,shf,lhf,lwup,swup,tg,ustar,psurf,ug,vg & & ,hadvt,hadvq,hadvu,hadvv,w,omega) @@ -3793,4 +3793,5 @@ #include "netcdf.inc" +#include "YOMCST.h" integer ntime,nlevel @@ -3800,5 +3801,5 @@ real*8 zz(nlevel) - real*8 th(nlevel),pres(nlevel) + real*8 th(nlevel),pres(nlevel),t(nlevel) real*8 qv(nlevel),u(nlevel),v(nlevel),o3(nlevel) real*8 shf(ntime),lhf(ntime),lwup(ntime),swup(ntime),tg(ntime) @@ -3806,4 +3807,5 @@ real*8 hadvt(nlevel,ntime),hadvq(nlevel,ntime),hadvu(nlevel,ntime) real*8 hadvv(nlevel,ntime),w(nlevel,ntime),omega(nlevel,ntime) + real*8 pzero integer nid, ierr @@ -3812,4 +3814,5 @@ integer var3didin(nbvar3d) + pzero=100000. ierr = NF_OPEN(fich_dice,NF_NOWRITE,nid) if (ierr.NE.NF_NOERR) then @@ -3986,4 +3989,7 @@ endif ! write(*,*)'lecture th ok',th + do k=1,nlevel + t(k)=th(k)*(pres(k)/pzero)**rkappa + enddo #ifdef NC_DOUBLE @@ -4798,5 +4804,5 @@ ! ! Cette formule remplace d_q = (1/tau) [rh_targ - rh] qsat(T_new) -! qui n'était pas correcte. +! qui n'?tait pas correcte. ! IF (tnew.LT.RTT) THEN @@ -4873,2 +4879,3 @@ END + Index: LMDZ5/trunk/libf/phylmd/dyn1d/1D_decl_cases.h =================================================================== --- LMDZ5/trunk/libf/phylmd/dyn1d/1D_decl_cases.h (revision 2681) +++ LMDZ5/trunk/libf/phylmd/dyn1d/1D_decl_cases.h (revision 2683) @@ -146,5 +146,5 @@ real zz_dice(nlev_dice) - real th_dice(nlev_dice),qv_dice(nlev_dice) + real t_dice(nlev_dice),qv_dice(nlev_dice) real u_dice(nlev_dice), v_dice(nlev_dice),o3_dice(nlev_dice) real ht_dice(nlev_dice,nt_dice) @@ -153,5 +153,5 @@ real w_dice(nlev_dice,nt_dice),omega_dice(nlev_dice,nt_dice) real o3_mod(llm),hu_mod(llm),hv_mod(llm) - real th_dicei(nlev_dice),qv_dicei(nlev_dice) + real t_dicei(nlev_dice),qv_dicei(nlev_dice) real u_dicei(nlev_dice), v_dicei(nlev_dice),o3_dicei(nlev_dice) real ht_dicei(nlev_dice) @@ -287,2 +287,3 @@ !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! + Index: LMDZ5/trunk/libf/phylmd/dyn1d/1D_interp_cases.h =================================================================== --- LMDZ5/trunk/libf/phylmd/dyn1d/1D_interp_cases.h (revision 2681) +++ LMDZ5/trunk/libf/phylmd/dyn1d/1D_interp_cases.h (revision 2683) @@ -118,7 +118,7 @@ ! vertical interpolation: CALL interp_dice_vertical(play,nlev_dice,nt_dice,plev_dice & - & ,th_dice,qv_dice,u_dice,v_dice,o3_dice & + & ,t_dice,qv_dice,u_dice,v_dice,o3_dice & & ,ht_profd,hq_profd,hu_profd,hv_profd,w_profd,omega_profd & - & ,th_mod,qv_mod,u_mod,v_mod,o3_mod & + & ,t_mod,qv_mod,u_mod,v_mod,o3_mod & & ,ht_mod,hq_mod,hu_mod,hv_mod,w_mod,omega_mod,mxcalc) ! do l = 1, llm @@ -810,2 +810,3 @@ !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! + Index: LMDZ5/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h =================================================================== --- LMDZ5/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h (revision 2681) +++ LMDZ5/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h (revision 2683) @@ -367,5 +367,5 @@ fich_dice='dice_driver.nc' call read_dice(fich_dice,nlev_dice,nt_dice & - & ,zz_dice,plev_dice,th_dice,qv_dice,u_dice,v_dice,o3_dice & + & ,zz_dice,plev_dice,t_dice,qv_dice,u_dice,v_dice,o3_dice & & ,shf_dice,lhf_dice,lwup_dice,swup_dice,tg_dice,ustar_dice& & ,psurf_dice,ug_dice,vg_dice,ht_dice,hq_dice & @@ -376,5 +376,5 @@ !champs initiaux: do k=1,nlev_dice - th_dicei(k)=th_dice(k) + t_dicei(k)=t_dice(k) qv_dicei(k)=qv_dice(k) u_dicei(k)=u_dice(k) @@ -405,7 +405,7 @@ CALL interp_dice_vertical(play,nlev_dice,nt_dice,plev_dice & - & ,th_dicei,qv_dicei,u_dicei,v_dicei,o3_dicei & + & ,t_dicei,qv_dicei,u_dicei,v_dicei,o3_dicei & & ,ht_dicei,hq_dicei,hu_dicei,hv_dicei,w_dicei,omega_dicei& - & ,th_mod,qv_mod,u_mod,v_mod,o3_mod & + & ,t_mod,qv_mod,u_mod,v_mod,o3_mod & & ,ht_mod,hq_mod,hu_mod,hv_mod,w_mod,omega_mod,mxcalc) @@ -425,6 +425,6 @@ do l = 1, llm ! Ligne du dessous ?? decommenter si on lit theta au lieu de temp - temp(l) = th_mod(l)*(play(l)/pzero)**rkappa -! temp(l) = t_mod(l) +! temp(l) = th_mod(l)*(play(l)/pzero)**rkappa + temp(l) = t_mod(l) q(l,1) = qv_mod(l) q(l,2) = 0.0 @@ -476,5 +476,5 @@ !--------------------------------------------------------------------- -!!!! Si la temperature de surface n'est pas imposée: +!!!! Si la temperature de surface n'est pas impos??e: if (forcing_gabls4) then @@ -531,5 +531,5 @@ write(*,*) 'SST initiale mxcalc: ',tsurf,mxcalc do l = 1, llm -! Ligne du dessous à decommenter si on lit theta au lieu de temp +! Ligne du dessous ?? decommenter si on lit theta au lieu de temp ! temp(l) = th_mod(l)*(play(l)/pzero)**rkappa temp(l) = t_mod(l) @@ -909,2 +909,3 @@ endif !forcing_case !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! +