Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/phylmd/cv30_routines.F90
- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cv30_routines.F90
r4593 r5082 3082 3082 3083 3083 ! fraction deau condensee dans les melanges convertie en precip : epm 3084 ! et eau condens ée précipitée dans masse d'air saturé: l_m*dM_m/dzdz.dzdz3084 ! et eau condens�e pr�cipit�e dans masse d'air satur� : l_m*dM_m/dzdz.dzdz 3085 3085 DO j = 1, nam1 3086 3086 DO k = 1, j - 1 … … 3277 3277 3278 3278 ! On fait varier epmax en fn de la cape 3279 ! Il faut donc recalculer ep, et hp qui a d éjà été calculéet3280 ! qui en d épend3281 ! Toutes les autres variables fn de ep sont calcul ées plus bas.3279 ! Il faut donc recalculer ep, et hp qui a d�j� �t� calcul� et 3280 ! qui en d�pend 3281 ! Toutes les autres variables fn de ep sont calcul�es plus bas. 3282 3282 3283 3283 INCLUDE "cvthermo.h" … … 3304 3304 ! on recalcule ep et hp 3305 3305 3306 if (coef_epmax_cape .gt.1e-12) then3306 if (coef_epmax_cape>1e-12) then 3307 3307 do i=1,ncum 3308 3308 epmax_diag(i)=epmax-coef_epmax_cape*sqrt(cape(i)) … … 3327 3327 do k=minorig+1,nl 3328 3328 do i=1,ncum 3329 if((k .ge.icb(i)).and.(k.le.inb(i)))then3329 if((k>=icb(i)).and.(k<=inb(i)))then 3330 3330 hp(i,k)=h(i,nk(i))+(lv(i,k)+(cpd-cpv)*t(i,k))*ep(i,k)*clw(i,k) 3331 3331 endif … … 3335 3335 do i=1,ncum 3336 3336 do k=1,nl 3337 if (abs(hp_bak(i,k)-hp(i,k)) .gt.0.01) then3337 if (abs(hp_bak(i,k)-hp(i,k))>0.01) then 3338 3338 write(*,*) 'i,k=',i,k 3339 3339 write(*,*) 'coef_epmax_cape=',coef_epmax_cape
Note: See TracChangeset
for help on using the changeset viewer.