- Timestamp:
- Jul 24, 2024, 4:23:34 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_plume_6A.f90
r5116 r5117 25 25 IMPLICIT NONE 26 26 27 integer, intent(in) :: itap, lev_out, lunout1, igout, ngrid, nlay28 real, intent(in) :: ptimestep29 real, intent(in), dimension(ngrid, nlay) :: ztv30 real, intent(in), dimension(ngrid, nlay) :: zthl31 real, intent(in), dimension(ngrid, nlay) :: po32 real, intent(in), dimension(ngrid, nlay) :: zl33 real, intent(in), dimension(ngrid, nlay) :: rhobarz34 real, intent(in), dimension(ngrid, nlay + 1) :: zlev35 real, intent(in), dimension(ngrid, nlay + 1) :: pplev36 real, intent(in), dimension(ngrid, nlay) :: pphi37 real, intent(in), dimension(ngrid, nlay) :: zpspsk38 real, intent(in), dimension(ngrid) :: f039 40 integer, intent(out) :: lalim(ngrid)41 real, intent(out), dimension(ngrid, nlay) :: alim_star42 real, intent(out), dimension(ngrid) :: alim_star_tot43 real, intent(out), dimension(ngrid, nlay) :: detr_star44 real, intent(out), dimension(ngrid, nlay) :: entr_star45 real, intent(out), dimension(ngrid, nlay + 1) :: f_star46 real, intent(out), dimension(ngrid, nlay) :: csc47 real, intent(out), dimension(ngrid, nlay) :: ztva48 real, intent(out), dimension(ngrid, nlay) :: ztla49 real, intent(out), dimension(ngrid, nlay) :: zqla50 real, intent(out), dimension(ngrid, nlay) :: zqta51 real, intent(out), dimension(ngrid, nlay) :: zha52 real, intent(out), dimension(ngrid, nlay + 1) :: zw253 real, intent(out), dimension(ngrid, nlay + 1) :: w_est54 real, intent(out), dimension(ngrid, nlay) :: ztva_est55 real, intent(out), dimension(ngrid, nlay) :: zqsatth56 integer, intent(out), dimension(ngrid) :: lmix57 integer, intent(out), dimension(ngrid) :: lmix_bis58 real, intent(out), dimension(ngrid) :: linter27 INTEGER, INTENT(IN) :: itap, lev_out, lunout1, igout, ngrid, nlay 28 REAL, INTENT(IN) :: ptimestep 29 REAL, INTENT(IN), DIMENSION(ngrid, nlay) :: ztv 30 REAL, INTENT(IN), DIMENSION(ngrid, nlay) :: zthl 31 REAL, INTENT(IN), DIMENSION(ngrid, nlay) :: po 32 REAL, INTENT(IN), DIMENSION(ngrid, nlay) :: zl 33 REAL, INTENT(IN), DIMENSION(ngrid, nlay) :: rhobarz 34 REAL, INTENT(IN), DIMENSION(ngrid, nlay + 1) :: zlev 35 REAL, INTENT(IN), DIMENSION(ngrid, nlay + 1) :: pplev 36 REAL, INTENT(IN), DIMENSION(ngrid, nlay) :: pphi 37 REAL, INTENT(IN), DIMENSION(ngrid, nlay) :: zpspsk 38 REAL, INTENT(IN), DIMENSION(ngrid) :: f0 39 40 INTEGER, INTENT(OUT) :: lalim(ngrid) 41 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: alim_star 42 REAL, INTENT(OUT), DIMENSION(ngrid) :: alim_star_tot 43 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: detr_star 44 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: entr_star 45 REAL, INTENT(OUT), DIMENSION(ngrid, nlay + 1) :: f_star 46 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: csc 47 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: ztva 48 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: ztla 49 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: zqla 50 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: zqta 51 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: zha 52 REAL, INTENT(OUT), DIMENSION(ngrid, nlay + 1) :: zw2 53 REAL, INTENT(OUT), DIMENSION(ngrid, nlay + 1) :: w_est 54 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: ztva_est 55 REAL, INTENT(OUT), DIMENSION(ngrid, nlay) :: zqsatth 56 INTEGER, INTENT(OUT), DIMENSION(ngrid) :: lmix 57 INTEGER, INTENT(OUT), DIMENSION(ngrid) :: lmix_bis 58 REAL, INTENT(OUT), DIMENSION(ngrid) :: linter 59 59 60 60 REAL zdw2, zdw2bis 61 61 REAL zw2modif 62 62 REAL zw2fact, zw2factbis 63 REAL, dimension(ngrid, nlay) :: zeps64 65 REAL, dimension(ngrid) :: wmaxa63 REAL, DIMENSION(ngrid, nlay) :: zeps 64 65 REAL, DIMENSION(ngrid) :: wmaxa 66 66 67 67 INTEGER ig, l, k, lt, it, lm 68 integernbpb69 70 real, dimension(ngrid, nlay) :: detr71 real, dimension(ngrid, nlay) :: entr72 real, dimension(ngrid, nlay + 1) :: wa_moy73 real, dimension(ngrid, nlay) :: ztv_est74 real, dimension(ngrid) :: ztemp, zqsat75 real, dimension(ngrid, nlay) :: zqla_est76 real, dimension(ngrid, nlay) :: zta_est77 78 real, dimension(ngrid, nlay) :: zbuoy, gamma, zdqt79 realzdz, zalpha, zw2m80 real, dimension(ngrid, nlay) :: zbuoyjam, zdqtjam81 realzbuoybis, zdz2, zdz3, lmel, entrbis, zdzbis82 real, dimension(ngrid) :: d_temp83 realztv1, ztv2, factinv, zinv, zlmel84 realzlmelup, zlmeldwn, zlt, zltdwn, zltup85 realatv1, atv2, btv1, btv286 realztv_est1, ztv_est287 realzcor, zdelta, zcvm5, qlbef88 realzbetalpha, coefzlmel89 realeps90 logicalZsat91 LOGICAL, dimension(ngrid) :: active, activetmp68 INTEGER nbpb 69 70 REAL, DIMENSION(ngrid, nlay) :: detr 71 REAL, DIMENSION(ngrid, nlay) :: entr 72 REAL, DIMENSION(ngrid, nlay + 1) :: wa_moy 73 REAL, DIMENSION(ngrid, nlay) :: ztv_est 74 REAL, DIMENSION(ngrid) :: ztemp, zqsat 75 REAL, DIMENSION(ngrid, nlay) :: zqla_est 76 REAL, DIMENSION(ngrid, nlay) :: zta_est 77 78 REAL, DIMENSION(ngrid, nlay) :: zbuoy, gamma, zdqt 79 REAL zdz, zalpha, zw2m 80 REAL, DIMENSION(ngrid, nlay) :: zbuoyjam, zdqtjam 81 REAL zbuoybis, zdz2, zdz3, lmel, entrbis, zdzbis 82 REAL, DIMENSION(ngrid) :: d_temp 83 REAL ztv1, ztv2, factinv, zinv, zlmel 84 REAL zlmelup, zlmeldwn, zlt, zltdwn, zltup 85 REAL atv1, atv2, btv1, btv2 86 REAL ztv_est1, ztv_est2 87 REAL zcor, zdelta, zcvm5, qlbef 88 REAL zbetalpha, coefzlmel 89 REAL eps 90 LOGICAL Zsat 91 LOGICAL, DIMENSION(ngrid) :: active, activetmp 92 92 REAL fact_gamma, fact_gamma2, fact_epsilon2 93 93 REAL coefc 94 REAL, dimension(ngrid, nlay) :: c295 96 if(ngrid==1) PRINT*, 'THERMCELL PLUME MODIFIE 2014/07/11'94 REAL, DIMENSION(ngrid, nlay) :: c2 95 96 IF (ngrid==1) PRINT*, 'THERMCELL PLUME MODIFIE 2014/07/11' 97 97 Zsat = .FALSE. 98 98 ! Initialisation … … 102 102 103 103 ! Initialisations des variables r?elles 104 if(1==1) THEN104 IF (1==1) THEN 105 105 ztva(:, :) = ztv(:, :) 106 106 ztva_est(:, :) = ztva(:, :) … … 173 173 ! Le panache va prendre au debut les caracteristiques de l'air contenu 174 174 ! dans cette couche. 175 if(active(ig)) THEN175 IF (active(ig)) THEN 176 176 ztla(ig, 1) = zthl(ig, 1) 177 177 zqta(ig, 1) = po(ig, 1) … … 197 197 do ig = 1, ngrid 198 198 active(ig) = active(ig) & 199 & . and. zw2(ig, l)>1.e-10 &200 & . and. f_star(ig, l) + alim_star(ig, l)>1.e-10199 & .AND. zw2(ig, l)>1.e-10 & 200 & .AND. f_star(ig, l) + alim_star(ig, l)>1.e-10 201 201 enddo 202 202 … … 247 247 248 248 !-------------------------------------------------- 249 if(iflag_thermals_ed<8) THEN249 IF (iflag_thermals_ed<8) THEN 250 250 !-------------------------------------------------- 251 251 !AJ052014: J'ai remplac?? la boucle do par un do while … … 284 284 ztv2 = atv2 * zlt + btv2 285 285 286 if (ztv(ig, lt)>ztv1.and.ztv(ig, lt)<ztv2) THEN286 IF (ztv(ig, lt)>ztv1.AND.ztv(ig, lt)<ztv2) THEN 287 287 !-------------------------------------------------- 288 288 !AJ052014: D??calage de zinv qui est entre le haut … … 292 292 zinv = zltdwn + zdz3 * factinv 293 293 294 if(zlmeldwn>=zinv) THEN294 IF (zlmeldwn>=zinv) THEN 295 295 ztv_est(ig, l) = atv2 * zlmel + btv2 296 296 zbuoyjam(ig, l) = fact_shell * RG * (ztva_est(ig, l) - ztv_est(ig, l)) / ztv_est(ig, l) & … … 311 311 endif 312 312 313 else ! if (ztv(ig,lt).gt.ztv1. and.ztv(ig,lt).lt.ztv2) THEN314 if(zlmeldwn>zltdwn) THEN313 else ! if (ztv(ig,lt).gt.ztv1.AND.ztv(ig,lt).lt.ztv2) THEN 314 IF (zlmeldwn>zltdwn) THEN 315 315 zbuoyjam(ig, l) = fact_shell * RG * ((ztva_est(ig, l) - & 316 316 ztv(ig, lt)) / ztv(ig, lt)) + (1. - fact_shell) * zbuoy(ig, l) … … 328 328 ! & po(ig,lt))/po(ig,lt)+((zdz2-lmel)/zdz3)*(zqta(ig,l-1)- & 329 329 ! & po(ig,lt-1))/po(ig,lt-1)) 330 endif ! if (ztv(ig,lt).gt.ztv1. and.ztv(ig,lt).lt.ztv2) THEN330 endif ! if (ztv(ig,lt).gt.ztv1.AND.ztv(ig,lt).lt.ztv2) THEN 331 331 else ! if (iflag_thermals_ed.lt.8) THEN 332 332 lt = l + 1 … … 373 373 !AJ052014: J'ai remplac? w_est(ig,l) par zw2(ig,l) 374 374 !-------------------------------------------------- 375 if(iflag_thermals_ed==8) THEN375 IF (iflag_thermals_ed==8) THEN 376 376 ! Ancienne version 377 377 ! w_est(ig,l+1)=Max(0.0001,(zdz/zdzbis)*(exp(-zw2fact)* & … … 404 404 endif 405 405 406 if(iflag_thermals_ed<6) THEN406 IF (iflag_thermals_ed<6) THEN 407 407 zalpha = f0(ig) * f_star(ig, l) / sqrt(w_est(ig, l + 1)) / rhobarz(ig, l) 408 408 ! fact_epsilon=0.0005/(zalpha+0.025)**0.5 … … 441 441 442 442 do ig = 1, ngrid 443 if(active(ig)) THEN443 IF (active(ig)) THEN 444 444 ! zw2m=max(0.5*(w_est(ig,l)+w_est(ig,l+1)),0.1) 445 445 zw2m = w_est(ig, l + 1) … … 471 471 ! entrbis=entr_star(ig,l) 472 472 473 if(iflag_thermals_ed<6) THEN473 IF (iflag_thermals_ed<6) THEN 474 474 fact_epsilon = 0.0002 / (zalpha + 0.1) 475 475 endif … … 508 508 ! En dessous de lalim, on prend le max de alim_star et entr_star pour 509 509 ! alim_star et 0 sinon 510 if(l<lalim(ig)) THEN510 IF (l<lalim(ig)) THEN 511 511 alim_star(ig, l) = max(alim_star(ig, l), entr_star(ig, l)) 512 512 entr_star(ig, l) = 0. 513 513 endif 514 ! if (l.lt.lalim(ig). and.alim_star(ig,l)>alim_star(ig,l-1)) THEN514 ! if (l.lt.lalim(ig).AND.alim_star(ig,l)>alim_star(ig,l-1)) THEN 515 515 ! alim_star(ig,l)=entrbis 516 516 ! endif … … 529 529 !=========================================================================== 530 530 531 activetmp(:) = active(:) . and. f_star(:, l + 1)>1.e-10532 do ig = 1, ngrid 533 if(activetmp(ig)) THEN531 activetmp(:) = active(:) .AND. f_star(:, l + 1)>1.e-10 532 do ig = 1, ngrid 533 IF (activetmp(ig)) THEN 534 534 Zsat = .FALSE. 535 535 ztla(ig, l) = (f_star(ig, l) * ztla(ig, l - 1) + & … … 546 546 CALL thermcell_qsat(ngrid, activetmp, pplev(:, l), ztemp, zqta(:, l), zqsatth(:, l)) 547 547 do ig = 1, ngrid 548 if(activetmp(ig)) THEN548 IF (activetmp(ig)) THEN 549 549 ! on ecrit de maniere conservative (sat ou non) 550 550 ! T = Tl +Lv/Cp ql … … 577 577 ! & (zw2(ig,l)-zdw2)+zdw2)+(zdzbis-zdz)/zdzbis* & 578 578 ! & (exp(-zw2factbis)*(zw2(ig,l-1)-zdw2bis)+zdw2)) 579 if(iflag_thermals_ed==8) THEN579 IF (iflag_thermals_ed==8) THEN 580 580 zw2(ig, l + 1) = Max(0.0001, exp(-zw2fact) * (zw2(ig, l) - zdw2) + zdw2) 581 581 else … … 586 586 ! & (exp(-zw2factbis)*(zw2(ig,l-1)-zdw2bis)+zdw2bis)) 587 587 588 if(iflag_thermals_ed<6) THEN588 IF (iflag_thermals_ed<6) THEN 589 589 zalpha = f0(ig) * f_star(ig, l) / sqrt(zw2(ig, l + 1)) / rhobarz(ig, l) 590 590 ! fact_epsilon=0.0005/(zalpha+0.025)**0.5 … … 607 607 enddo 608 608 609 if(prt_level>=20) PRINT*, 'coucou calcul detr 460: ig, l', ig, l609 IF (prt_level>=20) PRINT*, 'coucou calcul detr 460: ig, l', ig, l 610 610 611 611 !=========================================================================== … … 615 615 nbpb = 0 616 616 do ig = 1, ngrid 617 if (zw2(ig, l + 1)>0. .and. zw2(ig, l + 1)<1.e-10) THEN617 IF (zw2(ig, l + 1)>0. .AND. zw2(ig, l + 1)<1.e-10) THEN 618 618 ! stop 'On tombe sur le cas particulier de thermcell_dry' 619 619 ! PRINT*,'On tombe sur le cas particulier de thermcell_plume' … … 623 623 endif 624 624 625 if(zw2(ig, l + 1)<0.) THEN625 IF (zw2(ig, l + 1)<0.) THEN 626 626 linter(ig) = (l * (zw2(ig, l + 1) - zw2(ig, l)) & 627 627 - zw2(ig, l)) / (zw2(ig, l + 1) - zw2(ig, l)) … … 637 637 wa_moy(ig, l + 1) = sqrt(zw2(ig, l + 1)) 638 638 639 if(wa_moy(ig, l + 1)>wmaxa(ig)) THEN639 IF (wa_moy(ig, l + 1)>wmaxa(ig)) THEN 640 640 ! lmix est le niveau de la couche ou w (wa_moy) est maximum 641 641 !on rajoute le calcul de lmix_bis 642 if(zqla(ig, l)<1.e-10) THEN642 IF (zqla(ig, l)<1.e-10) THEN 643 643 lmix_bis(ig) = l + 1 644 644 endif … … 648 648 enddo 649 649 650 if(nbpb>0) THEN650 IF (nbpb>0) THEN 651 651 PRINT*, 'WARNING on tombe ', nbpb, ' x sur un pb pour l=', l, ' dans thermcell_plume' 652 652 endif … … 667 667 enddo 668 668 669 if(prt_level>=20) PRINT*, 'coucou calcul detr 470: ig, l', ig, l669 IF (prt_level>=20) PRINT*, 'coucou calcul detr 470: ig, l', ig, l 670 670 RETURN 671 671 end … … 705 705 REAL f0(ngrid) 706 706 INTEGER lalim(ngrid) 707 integerlev_out ! niveau pour les print708 integernbpb709 710 realalim_star_tot(ngrid)707 INTEGER lev_out ! niveau pour les print 708 INTEGER nbpb 709 710 REAL alim_star_tot(ngrid) 711 711 712 712 REAL ztva(ngrid, nlay) … … 747 747 INTEGER ig, l, k 748 748 749 realzdz, zbuoy(ngrid, nlay), zalpha, gamma(ngrid, nlay), zdqt(ngrid, nlay), zw2m750 realzbuoybis751 realzcor, zdelta, zcvm5, qlbef, zdz2752 realbetalpha, zbetalpha753 realeps, afact754 logicalZsat749 REAL zdz, zbuoy(ngrid, nlay), zalpha, gamma(ngrid, nlay), zdqt(ngrid, nlay), zw2m 750 REAL zbuoybis 751 REAL zcor, zdelta, zcvm5, qlbef, zdz2 752 REAL betalpha, zbetalpha 753 REAL eps, afact 754 LOGICAL Zsat 755 755 LOGICAL active(ngrid), activetmp(ngrid) 756 756 REAL fact_gamma, fact_epsilon, fact_gamma2, fact_epsilon2 … … 767 767 768 768 ! Initialisations des variables reeles 769 if(1==1) THEN769 IF (1==1) THEN 770 770 ztva(:, :) = ztv(:, :) 771 771 ztva_est(:, :) = ztva(:, :) … … 819 819 do l = 1, nlay - 1 820 820 do ig = 1, ngrid 821 if (ztv(ig, l)> ztv(ig, l + 1) .and. ztv(ig, 1)>=ztv(ig, l)) THEN821 IF (ztv(ig, l)> ztv(ig, l + 1) .AND. ztv(ig, 1)>=ztv(ig, l)) THEN 822 822 alim_star(ig, l) = MAX((ztv(ig, l) - ztv(ig, l + 1)), 0.) & 823 823 * sqrt(zlev(ig, l + 1)) … … 829 829 do l = 1, nlay 830 830 do ig = 1, ngrid 831 if(alim_star_tot(ig) > 1.e-10) THEN831 IF (alim_star_tot(ig) > 1.e-10) THEN 832 832 alim_star(ig, l) = alim_star(ig, l) / alim_star_tot(ig) 833 833 endif … … 848 848 ! Le panache va prendre au debut les caracteristiques de l'air contenu 849 849 ! dans cette couche. 850 if(active(ig)) THEN850 IF (active(ig)) THEN 851 851 ztla(ig, 1) = zthl(ig, 1) 852 852 zqta(ig, 1) = po(ig, 1) … … 872 872 do ig = 1, ngrid 873 873 active(ig) = active(ig) & 874 & . and. zw2(ig, l)>1.e-10 &875 & . and. f_star(ig, l) + alim_star(ig, l)>1.e-10874 & .AND. zw2(ig, l)>1.e-10 & 875 & .AND. f_star(ig, l) + alim_star(ig, l)>1.e-10 876 876 enddo 877 877 … … 911 911 w_est(ig, l + 1) = Max(0.0001, exp(-zw2fact) * (w_est(ig, l) - zdw2) + zdw2) 912 912 913 if(w_est(ig, l + 1)<0.) THEN913 IF (w_est(ig, l + 1)<0.) THEN 914 914 w_est(ig, l + 1) = zw2(ig, l) 915 915 endif … … 923 923 924 924 do ig = 1, ngrid 925 if(active(ig)) THEN925 IF (active(ig)) THEN 926 926 zw2m = max(0.5 * (w_est(ig, l) + w_est(ig, l + 1)), 0.1) 927 927 zw2m = w_est(ig, l + 1) … … 942 942 ! En dessous de lalim, on prend le max de alim_star et entr_star pour 943 943 ! alim_star et 0 sinon 944 if(l<lalim(ig)) THEN944 IF (l<lalim(ig)) THEN 945 945 alim_star(ig, l) = max(alim_star(ig, l), entr_star(ig, l)) 946 946 entr_star(ig, l) = 0. … … 958 958 !calcul de la vitesse verticale en melangeant Tl et qt du thermique 959 959 !--------------------------------------------------------------------------- 960 activetmp(:) = active(:) . and. f_star(:, l + 1)>1.e-10961 do ig = 1, ngrid 962 if(activetmp(ig)) THEN960 activetmp(:) = active(:) .AND. f_star(:, l + 1)>1.e-10 961 do ig = 1, ngrid 962 IF (activetmp(ig)) THEN 963 963 Zsat = .FALSE. 964 964 ztla(ig, l) = (f_star(ig, l) * ztla(ig, l - 1) + & … … 976 976 977 977 do ig = 1, ngrid 978 if(activetmp(ig)) THEN978 IF (activetmp(ig)) THEN 979 979 ! on ecrit de maniere conservative (sat ou non) 980 980 ! T = Tl +Lv/Cp ql … … 996 996 enddo 997 997 998 if(prt_level>=20) PRINT*, 'coucou calcul detr 460: ig, l', ig, l998 IF (prt_level>=20) PRINT*, 'coucou calcul detr 460: ig, l', ig, l 999 999 1000 1000 !--------------------------------------------------------------------------- … … 1004 1004 nbpb = 0 1005 1005 do ig = 1, ngrid 1006 if (zw2(ig, l + 1)>0. .and. zw2(ig, l + 1)<1.e-10) THEN1006 IF (zw2(ig, l + 1)>0. .AND. zw2(ig, l + 1)<1.e-10) THEN 1007 1007 ! stop 'On tombe sur le cas particulier de thermcell_dry' 1008 1008 ! PRINT*,'On tombe sur le cas particulier de thermcell_plume' … … 1012 1012 endif 1013 1013 1014 if(zw2(ig, l + 1)<0.) THEN1014 IF (zw2(ig, l + 1)<0.) THEN 1015 1015 linter(ig) = (l * (zw2(ig, l + 1) - zw2(ig, l)) & 1016 1016 - zw2(ig, l)) / (zw2(ig, l + 1) - zw2(ig, l)) … … 1025 1025 wa_moy(ig, l + 1) = sqrt(zw2(ig, l + 1)) 1026 1026 1027 if(wa_moy(ig, l + 1)>wmaxa(ig)) THEN1027 IF (wa_moy(ig, l + 1)>wmaxa(ig)) THEN 1028 1028 ! lmix est le niveau de la couche ou w (wa_moy) est maximum 1029 1029 !on rajoute le calcul de lmix_bis 1030 if(zqla(ig, l)<1.e-10) THEN1030 IF (zqla(ig, l)<1.e-10) THEN 1031 1031 lmix_bis(ig) = l + 1 1032 1032 endif … … 1036 1036 enddo 1037 1037 1038 if(nbpb>0) THEN1038 IF (nbpb>0) THEN 1039 1039 PRINT*, 'WARNING on tombe ', nbpb, ' x sur un pb pour l=', l, ' dans thermcell_plume' 1040 1040 endif … … 1055 1055 enddo 1056 1056 1057 if(prt_level>=20) PRINT*, 'coucou calcul detr 470: ig, l', ig, l1057 IF (prt_level>=20) PRINT*, 'coucou calcul detr 470: ig, l', ig, l 1058 1058 end 1059 1059 END MODULE lmdz_thermcell_plume_6A
Note: See TracChangeset
for help on using the changeset viewer.
