Changeset 5086 for LMDZ6/branches/Amaury_dev/libf/phylmdiso
- Timestamp:
- Jul 19, 2024, 7:54:50 PM (5 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmdiso
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/concvl.F90
r4613 r5086 415 415 ! em_wght(k)=wght_th(i,k) 416 416 ! print*,'em_wght=',em_wght(k),wght_th(i,k) 417 ! end do417 ! END DO 418 418 ! END DO 419 419 -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv30_routines.F90
r4491 r5086 2032 2032 ! !!! traent(il,i,j,k)=sij(il,i,j)*tra(il,i,k) 2033 2033 ! !!! : +(1.-sij(il,i,j))*tra(il,nk(il),k) 2034 ! !!! end do2034 ! !!! END DO 2035 2035 elij(il, i, j) = altem 2036 2036 elij(il, i, j) = amax1(0.0, elij(il,i,j)) … … 3155 3155 ! : +tra(il,i,j)*(mp(il,i)-mp(il,i+1)) 3156 3156 ! trap(il,i,j)=trap(il,i,j)/mp(il,i) 3157 ! end do3157 ! END DO 3158 3158 3159 3159 ELSE … … 3172 3172 ! do j=1,ntra 3173 3173 ! trap(il,i,j)=trap(il,i+1,j) 3174 ! end do3174 ! END DO 3175 3175 3176 3176 END IF -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_routines.F90
r5082 r5086 2965 2965 !!!! traent(il,i,j,k)=sij(il,i,j)*tra(il,i,k) 2966 2966 !!!! : +(1.-sij(il,i,j))*tra(il,nk(il),k) 2967 !!!! end do2967 !!!! END DO 2968 2968 elij(il, i, j) = altem 2969 2969 elij(il, i, j) = max(0.0, elij(il,i,j)) … … 4485 4485 !AC! endif ! (i.lt.inb(il) .and. lwork(il)) 4486 4486 !AC! enddo 4487 !AC! end do4487 !AC! END DO 4488 4488 4489 4489 #ifdef ISO -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90
r4594 r5086 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 … … 1748 1748 ! c do i=1,klon 1749 1749 ! c print*,alpha(i) 1750 ! c end do1750 ! c END DO 1751 1751 ! cc 1752 1752 DO k = 1, klev … … 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/phyetat0_mod.F90
r5082 r5086 342 342 IF (.NOT. found) THEN 343 343 PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent" 344 PRINT*, " Il pr end donc la valeur de surface"344 PRINT*, " Il prEND DOnc la valeur de surface" 345 345 tsoil(:, isoil, :)=ftsol(:, :) 346 346 ENDIF -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/physiq_mod.F90
r5082 r5086 1869 1869 998 CONTINUE 1870 1870 CLOSE(98) 1871 CONTINUE1872 1871 IF(nCFMIP>npCFMIP) THEN 1873 1872 print*,'nCFMIP > npCFMIP : augmenter npCFMIP et recompiler'
Note: See TracChangeset
for help on using the changeset viewer.