Changeset 735 for LMDZ4/branches/V3_test/libf/phylmd/thermcell.F
- Timestamp:
- Nov 10, 2006, 3:50:16 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/V3_test/libf/phylmd/thermcell.F
r704 r735 14 14 c de "thermiques" explicitement representes 15 15 c 16 c R �criture �partir d'un listing papier �Habas, le 14/02/0017 c 18 c le thermique est suppos �homog�e et dissip�par m�ange avec19 c son environnement. la longueur l_mix contr �e l'efficacit�du20 c m �ange21 c 22 c Le calcul du transport des diff �entes esp�es se fait en prenant16 c Reecriture a partir d'un listing papier à Habas, le 14/02/00 17 c 18 c le thermique est suppose homogene et dissipe par melange avec 19 c son environnement. la longueur l_mix controle l'efficacite du 20 c melange 21 c 22 c Le calcul du transport des differentes especes se fait en prenant 23 23 c en compte: 24 24 c 1. un flux de masse montant … … 512 512 c calcul de la largeur de chaque ascendance dans le cas conservatif. 513 513 c dans ce cas simple, on suppose que la largeur de l'ascendance provenant 514 c d'une couche est �ale �la hauteur de la couche alimentante.514 c d'une couche est egale a la hauteur de la couche alimentante. 515 515 c La vitesse maximale dans l'ascendance est aussi prise comme estimation 516 516 c de la vitesse d'entrainement horizontal dans la couche alimentante. … … 551 551 c print*,'10 OK convect8' 552 552 c print*,'WA2 ',wa_moy 553 c calcul de la fraction de la maille concern �par l'ascendance en tenant553 c calcul de la fraction de la maille concerne par l'ascendance en tenant 554 554 c compte de l'epluchage du thermique. 555 555 c … … 1119 1119 if ((fm(ig,k+1)+detr(ig,k))*ptimestep.gt. 1120 1120 s 1.e-5*masse(ig,k)) then 1121 c On it �e sur la valeur du coeff de freinage.1121 c On itere sur la valeur du coeff de freinage. 1122 1122 c gamma0=rho(ig,k)*(zlev(ig,k+1)-zlev(ig,k)) 1123 1123 gamma0=masse(ig,k) … … 1125 1125 s *0.5/larga(ig) 1126 1126 c gamma0=0. 1127 c la premi �e fois on multiplie le coefficient de freinage1127 c la premiere fois on multiplie le coefficient de freinage 1128 1128 c par le module du vent dans la couche en dessous. 1129 1129 dua=ua(ig,k-1)-u(ig,k-1) … … 1317 1317 if ((fm(ig,k+1)+detr(ig,k))*ptimestep.gt. 1318 1318 s 1.e-5*masse(ig,k)) then 1319 c On it �e sur la valeur du coeff de freinage.1319 c On itere sur la valeur du coeff de freinage. 1320 1320 c gamma0=rho(ig,k)*(zlev(ig,k+1)-zlev(ig,k)) 1321 1321 gamma0=masse(ig,k) … … 1328 1328 zf=0. 1329 1329 zf2=1./(1.-zf) 1330 c la premi �e fois on multiplie le coefficient de freinage1330 c la premiere fois on multiplie le coefficient de freinage 1331 1331 c par le module du vent dans la couche en dessous. 1332 1332 dua=ua(ig,k-1)-u(ig,k-1)
Note: See TracChangeset
for help on using the changeset viewer.