Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/alpale_th.F90
- Timestamp:
- Aug 2, 2024, 2:12:03 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/alpale_th.F90
r5142 r5158 135 135 136 136 !----Initialisations 137 doi=1,klon137 DO i=1,klon 138 138 !!jyg proba_notrig(i)=1. 139 139 random_notrig(i)=1e6*ale_bl_stat(i)-int(1e6*ale_bl_stat(i)) … … 156 156 IF (iflag_trig_bl==1) THEN 157 157 !----Tirage al\'eatoire et calcul de ale_bl_trig 158 doi=1,klon158 DO i=1,klon 159 159 IF ( (ale_bl_stat(i) > abs(cin(i))+1.e-10) ) THEN 160 160 proba_notrig(i)=proba_notrig(i)* & … … 177 177 178 178 ELSE IF (iflag_trig_bl>=2) THEN 179 doi=1,klon179 DO i=1,klon 180 180 IF ( (Ale_bl(i) > abs(cin(i))+1.e-10) ) THEN 181 181 proba_notrig(i)=proba_notrig(i)* & … … 208 208 !-----------Statistical closure----------- 209 209 IF (iflag_clos_bl==1) THEN 210 doi=1,klon210 DO i=1,klon 211 211 !CR: alp probabiliste 212 212 IF (ale_bl_trig(i)>0.) THEN … … 217 217 ELSE IF (iflag_clos_bl==2) THEN 218 218 !CR: alp calculee dans thermcell_main 219 doi=1,klon219 DO i=1,klon 220 220 alp_bl(i)=alp_bl_stat(i) 221 221 enddo … … 244 244 ! iflag_trig_bl value. 245 245 IF (iflag_trig_bl==1) then ! use ale_bl_stat (Rochetin computation) 246 doi=1,klon246 DO i=1,klon 247 247 ale_bl_ref(i)=ale_bl_stat(i) 248 248 enddo 249 249 ELSE IF (iflag_trig_bl>=2) then ! use ale_bl (old computation) 250 doi=1,klon250 DO i=1,klon 251 251 ale_bl_ref(i)=Ale_bl(i) 252 252 enddo … … 255 255 256 256 !----Initializations and random number generation 257 doi=1,klon257 DO i=1,klon 258 258 !!jyg proba_notrig(i)=1. 259 259 random_notrig(i)=1e6*ale_bl_stat(i)-int(1e6*ale_bl_stat(i)) … … 273 273 274 274 !----alp_bl computation 275 doi=1,klon275 DO i=1,klon 276 276 IF ( (ale_bl_ref(i) > abs(cin(i))+1.e-10) ) THEN 277 277 birth_number = n2(i)*exp(-strig(i)/s2(i)) … … 292 292 293 293 !----ale_bl_trig computation 294 doi=1,klon294 DO i=1,klon 295 295 IF (random_notrig(i) >= proba_notrig(i)) THEN 296 296 ale_bl_trig(i)=ale_bl_ref(i) … … 324 324 WRITE(lunout,*)'Couplage Thermiques/Emanuel seulement si T<0' 325 325 ENDIF 326 doi=1,klon326 DO i=1,klon 327 327 IF (t_seri(i,lmax_th(i))>273.) THEN 328 328 Ale_bl(i)=0.
Note: See TracChangeset
for help on using the changeset viewer.