Changeset 5093 for LMDZ6/branches/Amaury_dev/libf/phylmdiso
- Timestamp:
- Jul 21, 2024, 1:07:18 PM (4 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmdiso
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/add_wake_tend.F90
r5082 r5093 5 5 ) 6 6 !=================================================================== 7 ! Ajoute les tendances li �es aux diverses parametrisations physiques aux7 ! Ajoute les tendances liées aux diverses parametrisations physiques aux 8 8 ! variables d'etat des poches froides. 9 9 !=================================================================== -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90
r5086 r5093 127 127 ! wdens_ref: initial number of wakes per unit area (3D) or per 128 128 ! unit length (2D), at the beginning of each time step 129 ! Tgw : 1 sur la p �riode de onde de gravit�130 ! Cgw : vitesse de propagation de onde de gravit �129 ! Tgw : 1 sur la période de onde de gravité 130 ! Cgw : vitesse de propagation de onde de gravité 131 131 ! LL : distance entre 2 poches 132 132 133 133 ! ------------------------------------------------------------------------- 134 ! D �claration de variables134 ! Déclaration de variables 135 135 ! ------------------------------------------------------------------------- 136 136 … … 196 196 ! ------------------- 197 197 198 ! Variables �fixer198 ! Variables à fixer 199 199 INTEGER, SAVE :: igout 200 200 !$OMP THREADPRIVATE(igout) … … 383 383 ! Configuration de coefgw,stark,wdens (22/02/06 by YU Jingmei) 384 384 385 ! coefgw : Coefficient pour les ondes de gravit �385 ! coefgw : Coefficient pour les ondes de gravité 386 386 ! stark : Coefficient k dans Cstar=k*sqrt(2*WAPE) 387 ! wdens : Densit �surfacique de poche froide387 ! wdens : Densité surfacique de poche froide 388 388 ! ------------------------------------------------------------------------- 389 389 … … 1088 1088 1089 1089 ! cc nrlmd Ajout d'un recalcul de wdens dans le cas d'un entrainement 1090 ! n �gatif de ktop �kupper --------1091 ! cc On calcule pour cela une densit �wdens0 pour laquelle on1090 ! négatif de ktop à kupper -------- 1091 ! cc On calcule pour cela une densité wdens0 pour laquelle on 1092 1092 ! aurait un entrainement nul --- 1093 1093 !jyg< … … 1096 1096 ! des descentes unsaturees. Nous faisons alors l'hypothese que la 1097 1097 ! convection profonde cree directement de nouvelles poches, sans passer 1098 ! par les thermiques. La nouvelle valeur de wdens est alors impos �e.1098 ! par les thermiques. La nouvelle valeur de wdens est alors imposée. 1099 1099 1100 1100 DO i = 1, klon … … 1195 1195 DO i = 1, klon 1196 1196 IF (wk_adv(i)) THEN 1197 ! cc nrlmd Introduction du taux de mortalit �des poches et1197 ! cc nrlmd Introduction du taux de mortalité des poches et 1198 1198 ! test sur sigmaw_max=0.4 1199 1199 ! cc d_sigmaw(i) = gfl(i)*Cstar(i)*dtimesub … … 1308 1308 1309 1309 ! c DO i=1,klon 1310 ! c print*,'Pente entre 0 et kupper (r �f�rence)'1310 ! c print*,'Pente entre 0 et kupper (référence)' 1311 1311 ! c $ ,omg(i,kupper(i)+1)/(pupper(i)-ph(i,1)) 1312 1312 ! c print*,'Pente entre ktop et kupper' … … 1602 1602 1603 1603 1604 ! Coefficient de r �partition1604 ! Coefficient de répartition 1605 1605 1606 1606 crep(i, k) = crep_sol*(ph(i,kupper(i))-ph(i,k))/ & … … 1646 1646 ! 1647 1647 1648 ! cc nrlmd Prise en compte du taux de mortalit �1649 ! cc D �finitions de entr, detr1648 ! cc nrlmd Prise en compte du taux de mortalité 1649 ! cc Définitions de entr, detr 1650 1650 !jyg< 1651 1651 !! detr(i, k) = 0. … … 1664 1664 1665 1665 1666 ! ajout d'un effet onde de gravit �-Tgw(k)*deltatw(k) 03/02/06 YU1666 ! ajout d'un effet onde de gravité -Tgw(k)*deltatw(k) 03/02/06 YU 1667 1667 ! Jingmei 1668 1668 … … 2601 2601 discrim = b*b - 4.*a*c 2602 2602 ! print*, 'x, a, b, c, discrim', x, a, b, c, discrim 2603 IF (a+b>=0.) THEN !! Condition suffisante pour la positivit �de ovap2603 IF (a+b>=0.) THEN !! Condition suffisante pour la positivité de ovap 2604 2604 alpha1(i) = 1. 2605 2605 ELSE -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyaqua_mod.F90
r5090 r5093 963 963 if (type_profil==20) then 964 964 print*,'Profile SST 20' 965 ! Méthode 13 "Qmax2K" plateau réel �| l'Equateur augmenté +2K965 ! Méthode 13 "Qmax2K" plateau réel é| l'Equateur augmenté +2K 966 966 967 967 do j=1,klon … … 972 972 if (type_profil==21) then 973 973 print*,'Profile SST 21' 974 ! Méthode 13 "Qmax2K" plateau réel �| l'Equateur augmenté +2K974 ! Méthode 13 "Qmax2K" plateau réel é| l'Equateur augmenté +2K 975 975 do j=1,klon 976 976 phy_sst(j,i)=252.+55.*(1-sin(rlatd(j))**2)
Note: See TracChangeset
for help on using the changeset viewer.