Changeset 5160 for LMDZ6/branches/Amaury_dev/libf/phylmdiso
- Timestamp:
- Aug 3, 2024, 2:56:58 PM (3 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmdiso
- Files:
-
- 2 added
- 2 deleted
- 12 edited
- 4 moved
-
add_phys_tend_mod.F90 (modified) (4 diffs)
-
calcul_REGDYN.h (deleted)
-
calwake.F90 (modified) (1 diff)
-
cv3_enthalpmix.F90 (modified) (5 diffs)
-
cv3_estatmix.F90 (modified) (5 diffs)
-
cv3_routines.F90 (modified) (17 diffs)
-
cv3p_mixing.F90 (modified) (4 diffs)
-
cv_driver.F90 (modified) (2 diffs)
-
cva_driver.F90 (modified) (5 diffs)
-
lmdz_chem.f90 (added)
-
lmdz_chem.h (moved) (moved from LMDZ6/branches/Amaury_dev/libf/phylmdiso/chem.h)
-
lmdz_concvl.F90 (modified) (2 diffs)
-
lmdz_ini_undefstd.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/phylmdiso/ini_undefSTD.F90) (1 diff)
-
lmdz_moy_undefstd.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/phylmdiso/moy_undefSTD.F90) (1 diff)
-
lmdz_plevel_new.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/phylmdiso/plevel_new.F90) (1 diff)
-
lmdz_undefstd.f90 (added)
-
lmdz_wake.F90 (modified) (4 diffs)
-
physiq_mod.F90 (modified) (22 diffs)
-
reevap.F90 (modified) (1 diff)
-
undefSTD.F90 (deleted)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/add_phys_tend_mod.F90
r5158 r5160 269 269 debug_level=10 270 270 IF (first) THEN 271 print*,"TestJLD rcpv, rcw, rcs",rcpv, rcw, rcs271 PRINT *,"TestJLD rcpv, rcw, rcs",rcpv, rcw, rcs 272 272 first=.FALSE. 273 273 endif 274 274 275 ! print*,'add_phys_tend: paprs ',paprs275 ! PRINT *,'add_phys_tend: paprs ',paprs 276 276 ! When in diagnostic mode, save initial values of out variables 277 277 IF (diag_mode == 1) THEN … … 746 746 debug_level=10 747 747 IF (first) THEN 748 print*,"TestJLD rcpv, rcw, rcs",rcpv, rcw, rcs748 PRINT *,"TestJLD rcpv, rcw, rcs",rcpv, rcw, rcs 749 749 first=.FALSE. 750 750 endif 751 751 752 ! print*,'add_phys_tend: paprs ',paprs752 ! PRINT *,'add_phys_tend: paprs ',paprs 753 753 !====================================================================== 754 754 ! Diagnostics for energy conservation tests … … 813 813 d_ek_col(:) = (zek_col(:,2)-zek_col(:,1))/phys_tstep 814 814 815 print*,'zdu ', zdu816 print*,'zdv ', zdv817 print*,'d_ek_col, zek_col(2), zek_col(1) ',d_ek_col(1), zek_col(1,2), zek_col(1,1)815 PRINT *,'zdu ', zdu 816 PRINT *,'zdv ', zdv 817 PRINT *,'d_ek_col, zek_col(2), zek_col(1) ',d_ek_col(1), zek_col(1,2), zek_col(1,1) 818 818 819 819 d_h_dair_col(:) = (zh_dair_col(:,2)-zh_dair_col(:,1))/phys_tstep … … 937 937 bilh_ok=0 938 938 939 !! print*,'prt_level:',prt_level,' fl_ebil:',fl_ebil,' fl_cor_ebil:',fl_cor_ebil939 !!PRINT *,'prt_level:',prt_level,' fl_ebil:',fl_ebil,' fl_cor_ebil:',fl_cor_ebil 940 940 IF ( (fl_ebil > 0) .AND. (klon == 1)) THEN 941 941 bilq_bnd = 0. -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/calwake.F90
r5158 r5160 143 143 144 144 IF (prt_level >= 10) THEN 145 print*, '-> calwake, wake_s, wgen input ', wake_s(1), wgen(1)145 PRINT *, '-> calwake, wake_s, wgen input ', wake_s(1), wgen(1) 146 146 ENDIF 147 147 -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_enthalpmix.F90
r5158 r5160 78 78 79 79 80 !! print*,' ->cv3_vertmix, plim1,plim2 ', plim1,plim2 !jyg80 !! PRINT *,' ->cv3_vertmix, plim1,plim2 ', plim1,plim2 !jyg 81 81 plim2p(:) = min(plim2(:),plim1(:)-dpmin) 82 82 j1(:)=nd … … 122 122 END DO 123 123 124 !! print*,'cv3_vertmix, j1,j2,coef ', j1,j2,coef !jyg124 !! PRINT *,'cv3_vertmix, j1,j2,coef ', j1,j2,coef !jyg 125 125 126 126 DO j = 1, nd … … 135 135 dpw(i) = dpw(i) + wi(i, j) 136 136 137 !! print*,'cv3_vertmix, j, wi(1,j),dpw ', j, wi(1,j),dpw !jyg137 !! PRINT *,'cv3_vertmix, j, wi(1,j),dpw ', j, wi(1,j),dpw !jyg 138 138 139 139 END IF … … 167 167 168 168 169 !! print*,'cv3_vertmix, rdcp ', rdcp !jyg169 !! PRINT *,'cv3_vertmix, rdcp ', rdcp !jyg 170 170 171 171 … … 233 233 END DO 234 234 235 !! print*,' cv3_vertmix->' !jyg235 !! PRINT *,' cv3_vertmix->' !jyg 236 236 237 237 -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_estatmix.F90
r5158 r5160 81 81 82 82 83 !! print*,' ->cv3_vertmix, plim1,plim2 ', plim1,plim2 !jyg83 !! PRINT *,' ->cv3_vertmix, plim1,plim2 ', plim1,plim2 !jyg 84 84 plim2p(:) = min(plim2(:),plim1(:)-dpmin) 85 85 j1(:)=nd … … 122 122 END DO 123 123 124 !! print*,'cv3_vertmix, j1,j2,coef ', j1,j2,coef !jyg124 !! PRINT *,'cv3_vertmix, j1,j2,coef ', j1,j2,coef !jyg 125 125 126 126 DO j = 1, nd … … 135 135 dpw(i) = dpw(i) + wi(i, j) 136 136 137 !! print*,'cv3_vertmix, j, wi(1,j),dpw ', j, wi(1,j),dpw !jyg137 !! PRINT *,'cv3_vertmix, j, wi(1,j),dpw ', j, wi(1,j),dpw !jyg 138 138 139 139 END IF … … 167 167 168 168 169 !! print*,'cv3_vertmix, rdcp ', rdcp !jyg169 !! PRINT *,'cv3_vertmix, rdcp ', rdcp !jyg 170 170 171 171 DO i = 1, len … … 216 216 END DO 217 217 218 !! print*,' cv3_vertmix->' !jyg218 !! PRINT *,' cv3_vertmix->' !jyg 219 219 220 220 -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_routines.F90
r5158 r5160 202 202 INTEGER il 203 203 204 ! print*,'cv3_incrcount : noconv_stop ',noconv_stop205 ! print*,'cv3_incrcount in, sig(1,nd) ',sig(1,nd)204 ! PRINT *,'cv3_incrcount : noconv_stop ',noconv_stop 205 ! PRINT *,'cv3_incrcount in, sig(1,nd) ',sig(1,nd) 206 206 IF(ok_convstop) THEN 207 207 DO il = 1, len … … 215 215 END DO 216 216 ENDIF ! (ok_convstop) 217 ! print*,'cv3_incrcount out, sig(1,nd) ',sig(1,nd)217 ! PRINT *,'cv3_incrcount out, sig(1,nd) ',sig(1,nd) 218 218 219 219 … … 280 280 (p(i,k-1)-p(i,k))/ph(i, k) !convect3 281 281 282 ! c print*,' gz(',k,')',gz(i,k),' tvx',tvx,' tvy ',tvy282 ! c PRINT *,' gz(',k,')',gz(i,k),' tvx',tvx,' tvy ',tvy 283 283 284 284 ! ori gz(i,k)=gz(i,k-1)+hrd*(tv(i,k-1)+tv(i,k)) … … 611 611 ENDIF 612 612 IF (prt_level >= 10) THEN 613 print*,'cv3_feed : iflag(1), pfeed(1), plcl(1), wghti(1,k) ', &613 PRINT *,'cv3_feed : iflag(1), pfeed(1), plcl(1), wghti(1,k) ', & 614 614 iflag(1), pfeed(1), plcl(1), (wghti(1,k),k=1,10) 615 615 ENDIF … … 1582 1582 #endif 1583 1583 IF (prt_level >= 10) THEN 1584 print*,'cv3_undilute2.0. icvflag_Tpa, t(1,k), q(1,k), qs(1,k) ', &1584 PRINT *,'cv3_undilute2.0. icvflag_Tpa, t(1,k), q(1,k), qs(1,k) ', & 1585 1585 icvflag_Tpa, (k, t(1,k), q(1,k), qs(1,k), k = 1,nl) 1586 1586 ENDIF … … 1775 1775 ENDIF ! (tg .gt. Tx) 1776 1776 1777 !! print*,' j, k, Um, U, Ux, aa, bb, discr, dd, ddelta ', j, k, Um, U, Ux, aa, bb, discr, dd, ddelta1778 !! print*,' j, k, ah0(i), ahg, tg, qg, tp(i,k), ff ', j, k, ah0(i), ahg, tg, qg, tp(i,k), ff1777 !! PRINT *,' j, k, Um, U, Ux, aa, bb, discr, dd, ddelta ', j, k, Um, U, Ux, aa, bb, discr, dd, ddelta 1778 !! PRINT *,' j, k, ah0(i), ahg, tg, qg, tp(i,k), ff ', j, k, ah0(i), ahg, tg, qg, tp(i,k), ff 1779 1779 END IF ! (k>=(icbs(i)+1)) 1780 1780 END DO ! i = 1, ncum … … 1839 1839 Tp(i,k) = tg + (ah0(i) - ahg)/ & 1840 1840 (cpd + (cl-cpd)*qta(i,k-1) + (lv(i,k)+frac(i,k)*lf(i,k))*dqgdT) 1841 !! print*,'undilute2 iterations k, Tp(i,k), ah0(i), ahg ', &1841 !! PRINT *,'undilute2 iterations k, Tp(i,k), ah0(i), ahg ', & 1842 1842 !! k, Tp(i,k), ah0(i), ahg 1843 1843 END IF ! (k>=(icbs(i)+1)) … … 2216 2216 ! ========================================================================= 2217 2217 IF (prt_level >= 10) THEN 2218 print*,'cv3_undilute2.1. tp(1,k), tvp(1,k) ', &2218 PRINT *,'cv3_undilute2.1. tp(1,k), tvp(1,k) ', & 2219 2219 (k, tp(1,k), tvp(1,k), k = 1,nl) 2220 2220 ENDIF … … 3696 3696 3697 3697 ! ------------------------------------------------------ 3698 IF (prt_level >= 10) print*,' ->cv3_unsat, iflag(1) ', iflag(1)3698 IF (prt_level >= 10) PRINT *,' ->cv3_unsat, iflag(1) ', iflag(1) 3699 3699 3700 3700 smallestreal=tiny(smallestreal) … … 4114 4114 ! water(il,i)=revap*revap 4115 4115 ! prec(il,i)=revap*revap 4116 !! print*,' i,b6,c6,revap,evap(il,i),water(il,i),wdtrain(il) ', &4116 !! PRINT *,' i,b6,c6,revap,evap(il,i),water(il,i),wdtrain(il) ', & 4117 4117 !! i,b6,c6,revap,evap(il,i),water(il,i),wdtrain(il) 4118 4118 !!---end jyg--- … … 4139 4139 b6 = bfac*50.*sigd(il)*(ph(il,i)-ph(il,i+1))*sigt*afac1 4140 4140 c6 = prec(il, i+1) + 0.5*bfac*wdtrain(il) 4141 ! print*,'bfac,sigd(il),sigt,afac1 ',bfac,sigd(il),sigt,afac14142 ! print*,'prec(il,i+1),wdtrain(il) ',prec(il,i+1),wdtrain(il)4143 ! print*,'b6,c6,b6*b6+4.*c6 ',b6,c6,b6*b6+4.*c64141 ! PRINT *,'bfac,sigd(il),sigt,afac1 ',bfac,sigd(il),sigt,afac1 4142 ! PRINT *,'prec(il,i+1),wdtrain(il) ',prec(il,i+1),wdtrain(il) 4143 ! PRINT *,'b6,c6,b6*b6+4.*c6 ',b6,c6,b6*b6+4.*c6 4144 4144 IF (c6>b6*b6+1.E-20) THEN 4145 4145 revap = 2.*c6/(b6+sqrt(b6*b6+4.*c6)) … … 4163 4163 !JYG--- Dans tous les cas, evaporation = [tt ce qui entre dans la couche i] 4164 4164 ! moins [tt ce qui sort de la couche i] 4165 ! print*, 'evap avec ice'4165 ! PRINT *, 'evap avec ice' 4166 4166 evap(il, i) = (wdtrain(il)+sigd(il)*wt(il,i)*(prec(il,i+1)-prec(il,i))) / & 4167 4167 (sigd(il)*(ph(il,i)-ph(il,i+1))*100.) … … 4227 4227 water(il, i) = 0. 4228 4228 END IF 4229 ! print*, 'evap sans ice'4229 ! PRINT *, 'evap sans ice' 4230 4230 evap(il, i) = (wdtrain(il)+sigd(il)*wt(il,i)*(water(il,i+1)-water(il,i)))/ & 4231 4231 (sigd(il)*(ph(il,i)-ph(il,i+1))*100.) … … 5645 5645 !----------------------------------------------------------- 5646 5646 5647 !! print*,'yield, i, amp1, ad', i, amp1(1), ad(1)5647 !! PRINT *,'yield, i, amp1, ad', i, amp1(1), ad(1) 5648 5648 5649 5649 DO il = 1, ncum … … 6798 6798 6799 6799 IF (prt_level >= 5) THEN 6800 print*,' CV3_YIELD : alpha_qpos ',alpha_qpos(1)6800 PRINT *,' CV3_YIELD : alpha_qpos ',alpha_qpos(1) 6801 6801 ENDIF 6802 6802 … … 7022 7022 END IF 7023 7023 END IF 7024 ! c print*,'cbmf',il,i,k,cbmf(il),wghti(il,k)7024 ! c PRINT *,'cbmf',il,i,k,cbmf(il),wghti(il,k) 7025 7025 END DO 7026 7026 END DO … … 7035 7035 dnwd(il, i) = dnwd(il, i) + dn1(il, k, i) 7036 7036 END IF 7037 ! c print*,'upwd',il,i,k,inb(il),upwd(il,i),m(il,k),up1(il,k,i)7037 ! c PRINT *,'upwd',il,i,k,inb(il),upwd(il,i),m(il,k),up1(il,k,i) 7038 7038 END DO 7039 7039 END DO -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3p_mixing.F90
r5158 r5160 278 278 ELSE IF (j > i) THEN 279 279 IF (prt_level >= 10) THEN 280 print*,'cv3p_mixing i, j, Sij given by the no-precip eq. ', i, j, Sij(il,i,j)280 PRINT *,'cv3p_mixing i, j, Sij given by the no-precip eq. ', i, j, Sij(il,i,j) 281 281 ENDIF 282 282 END IF ! new ! IF ((i>=icb(il)) .AND. (i<=inb(il)) .AND. (j>=(icb(il)-1)) & … … 476 476 !jygdebug< 477 477 IF (prt_level >= 10) THEN 478 print*,'cv3p_mixing i, nent(i), icb, inb ',i, nent(igout,i), icb(igout), inb(igout)478 PRINT *,'cv3p_mixing i, nent(i), icb, inb ',i, nent(igout,i), icb(igout), inb(igout) 479 479 IF (nent(igout,i) > 0) THEN 480 print*,'i,(j,Sij(i,j),j=icb-1,inb) ',i,(j,Sij(igout,i,j),j=icb(igout)-1,inb(igout))480 PRINT *,'i,(j,Sij(i,j),j=icb-1,inb) ',i,(j,Sij(igout,i,j),j=icb(igout)-1,inb(igout)) 481 481 ENDIF 482 482 ENDIF … … 848 848 END IF 849 849 850 ! print*,h(il,j)-hent(il,i,j),LV(il,j)*rs(il,j)/(cpd*rrv*t(il,j)*850 ! PRINT *,h(il,j)-hent(il,i,j),LV(il,j)*rs(il,j)/(cpd*rrv*t(il,j)* 851 851 ! : t(il,j)) 852 852 … … 865 865 !IM 301008 end 866 866 867 ! print*,'mix : i,j,hent(il,i,j),Sigij(il,i,j) ',867 ! PRINT *,'mix : i,j,hent(il,i,j),Sigij(il,i,j) ', 868 868 ! : i,j,hent(il,i,j),Sigij(il,i,j) 869 869 -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv_driver.F90
r5158 r5160 491 491 ! --- SET CONSTANTS AND PARAMETERS 492 492 ! ------------------------------------------------------------------- 493 ! print*, '-> cv_driver' !jyg493 ! PRINT *, '-> cv_driver' !jyg 494 494 ! -- set simulation flags: 495 495 ! (common cvflag) … … 1227 1227 END IF ! ncum>0 1228 1228 1229 ! print*, 'fin cv_driver ->' !jyg1229 ! PRINT *, 'fin cv_driver ->' !jyg 1230 1230 1231 1231 END SUBROUTINE cv_driver -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cva_driver.F90
r5158 r5160 666 666 667 667 668 ! print*, 't1, t1_wake ',(k,t1(1,k),t1_wake(1,k),k=1,nd)669 ! print*, 'q1, q1_wake ',(k,q1(1,k),q1_wake(1,k),k=1,nd)668 ! PRINT *, 't1, t1_wake ',(k,t1(1,k),t1_wake(1,k),k=1,nd) 669 ! PRINT *, 'q1, q1_wake ',(k,q1(1,k),q1_wake(1,k),k=1,nd) 670 670 671 671 ! ------------------------------------------------------------------- … … 912 912 END IF 913 913 914 ! print*, 'cv3_feed-> iflag1, plcl1 ',iflag1(1),plcl1(1)914 ! PRINT *, 'cv3_feed-> iflag1, plcl1 ',iflag1(1),plcl1(1) 915 915 916 916 ! -------------------------------------------------------------------- … … 1013 1013 ! ------------------------------------------------------------------- 1014 1014 1015 ! print*,' avant triggering, iflag_con ',iflag_con1015 ! PRINT *,' avant triggering, iflag_con ',iflag_con 1016 1016 1017 1017 IF (iflag_con==3) THEN … … 1344 1344 END IF 1345 1345 1346 ! print*,'cv_closure-> cape ',cape(1)1346 ! PRINT *,'cv_closure-> cape ',cape(1) 1347 1347 1348 1348 ! ------------------------------------------------------------------- … … 1559 1559 #endif 1560 1560 1561 ! print*,'cv_unsat-> mp ',mp1562 ! print*,'cv_unsat-> water ',water1561 ! PRINT *,'cv_unsat-> mp ',mp 1562 ! PRINT *,'cv_unsat-> water ',water 1563 1563 ! ------------------------------------------------------------------- 1564 1564 ! --- YIELD -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_concvl.F90
r5159 r5160 500 500 #endif 501 501 502 ! print*, '-> cv_driver' !jyg502 ! PRINT *, '-> cv_driver' !jyg 503 503 CALL cv_driver(klon, klev, klevp1, ntra, iflag_con, & 504 504 t, q, qs, u, v, tra, & … … 524 524 #endif 525 525 ) 526 ! print*, 'cv_driver ->' !jyg526 ! PRINT *, 'cv_driver ->' !jyg 527 527 528 528 #ifdef ISO -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_ini_undefstd.f90
r5159 r5160 1 link ../phylmd/ ini_undefSTD.F901 link ../phylmd/lmdz_ini_undefstd.f90 -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_moy_undefstd.f90
r5159 r5160 1 link ../phylmd/ moy_undefSTD.F901 link ../phylmd/lmdz_moy_undefstd.f90 -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_plevel_new.f90
r5159 r5160 1 link ../phylmd/ plevel_new.F901 link ../phylmd/lmdz_plevel_new.f90 -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90
r5158 r5160 1096 1096 1097 1097 DO i = 1, klon 1098 ! c print*,' isubstep,wk_adv(i),cstar(i),wape(i) ',1098 ! c PRINT *,' isubstep,wk_adv(i),cstar(i),wape(i) ', 1099 1099 ! c $ isubstep,wk_adv(i),cstar(i),wape(i) 1100 1100 IF (wk_adv(i) .AND. cstar(i)>0.01) THEN … … 1176 1176 1177 1177 IF (prt_level >= 10) THEN 1178 print*,'wake, cstar(1), cstar(1)/cstart, rad_wk(1), tau_wk_inv(1), drdt(1) ', &1178 PRINT *,'wake, cstar(1), cstar(1)/cstart, rad_wk(1), tau_wk_inv(1), drdt(1) ', & 1179 1179 cstar(1), cstar(1)/cstart, rad_wk(1), tau_wk_inv(1), drdt(1) 1180 print*,'wake, wdens(1), awdens(1), act(1), d_awdens(1) ', &1180 PRINT *,'wake, wdens(1), awdens(1), act(1), d_awdens(1) ', & 1181 1181 wdens(1), awdens(1), act(1), d_awdens(1) 1182 print*,'wake, wgen, -(wdens-awdens)*tau_wk_inv, -2.*wdens*gfl*drdt_pos, d_wdens ', &1182 PRINT *,'wake, wgen, -(wdens-awdens)*tau_wk_inv, -2.*wdens*gfl*drdt_pos, d_wdens ', & 1183 1183 wgen(1), -(wdens(1)-awdens(1))*tau_wk_inv(1), -2.*wdens(1)*gfl(1)*drdt_pos, d_wdens(1) 1184 print*,'wake, d_sig_gen(1), d_sig_death(1), d_sig_col(1), d_sigmaw(1) ', &1184 PRINT *,'wake, d_sig_gen(1), d_sig_death(1), d_sig_col(1), d_sigmaw(1) ', & 1185 1185 d_sig_gen(1), d_sig_death(1), d_sig_col(1), d_sigmaw(1) 1186 1186 ENDIF … … 1319 1319 END DO 1320 1320 END DO 1321 !! print*,'omg(igout,k) ', (k,omg(igout,k),k=1,klev)1321 !! PRINT *,'omg(igout,k) ', (k,omg(igout,k),k=1,klev) 1322 1322 ! cc nrlmd 1323 1323 ! c DO i=1,klon … … 2327 2327 DO i = 1, klon 2328 2328 IF (ok_qx_qw(i)) THEN 2329 !! print*,'wake, hw0(i), dthmin(i) ', hw0(i), dthmin(i)2330 !! print*,'wake, 2.*sum_dth(i)/(hw0(i)*dthmin(i)) ', &2329 !! PRINT *,'wake, hw0(i), dthmin(i) ', hw0(i), dthmin(i) 2330 !! PRINT *,'wake, 2.*sum_dth(i)/(hw0(i)*dthmin(i)) ', & 2331 2331 !! 2.*sum_dth(i)/(hw0(i)*dthmin(i)) 2332 !! print*,'wake, sum_half_dth(i), sum_dth(i) ', &2332 !! PRINT *,'wake, sum_half_dth(i), sum_dth(i) ', & 2333 2333 !! sum_half_dth(i), sum_dth(i) 2334 2334 IF ((hw0(i) < 1.) .OR. (dthmin(i) >= -delta_t_min) ) THEN 2335 2335 wape2(i) = -1. 2336 !! print*,'wake, rej 1'2336 !! PRINT *,'wake, rej 1' 2337 2337 ELSE IF (iflag_wk_check_trgl==1.AND.abs(2.*sum_dth(i)/(hw0(i)*dthmin(i)) - 1.) > 0.5) THEN 2338 2338 wape2(i) = -1. 2339 !! print*,'wake, rej 2'2339 !! PRINT *,'wake, rej 2' 2340 2340 ELSE IF (abs(sum_half_dth(i)) < 0.5*abs(sum_dth(i)) ) THEN 2341 2341 wape2(i) = -1. 2342 !! print*,'wake, rej 3'2342 !! PRINT *,'wake, rej 3' 2343 2343 END IF 2344 2344 END IF -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/physiq_mod.F90
r5158 r5160 93 93 USE phys_output_var_mod, ONLY: cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv 94 94 USE phys_output_var_mod, ONLY: cloud_cover_sw, cloud_cover_sw_s2 95 USE lmdz_plevel_new, ONLY: plevel_new 96 USE lmdz_moy_undefstd, ONLY: moy_undefstd 97 USE lmdz_undefstd, ONLY: undefstd 98 USE lmdz_ini_undefstd, ONLY: ini_undefstd 95 99 96 100 USE lmdz_geometry, ONLY: longitude, latitude, boundslon, boundslat, ind_cell_glo … … 106 110 USE paramLMDZ_phy_mod 107 111 108 USE phytracr_spl_mod, ONLY: phytracr_spl, phytracr_spl_out_init112 USE lmdz_phytracr_spl, ONLY: phytracr_spl, phytracr_spl_out_init 109 113 USE phys_output_write_spl_mod 110 114 USE phytrac_mod, ONLY: phytrac_init, phytrac … … 428 432 429 433 USE lmdz_yomcst 434 USE lmdz_phys_constants, ONLY: dobson_u 430 435 431 436 IMPLICIT NONE … … 902 907 REAL zmasse(klon, nbp_lev),exner(klon, nbp_lev) 903 908 ! (column-density of mass of air in a cell, in kg m-2) 904 REAL, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2905 909 906 910 !IM sorties … … 925 929 EXTERNAL suphel ! initialiser certaines constantes 926 930 EXTERNAL transp ! transport total de l'eau et de l'energie 927 !IM 931 928 932 EXTERNAL haut2bas !variables de haut en bas 929 EXTERNAL ini_undefSTD !initialise a 0 une variable a 1 niveau de pression930 EXTERNAL undefSTD !somme les valeurs definies d'1 var a 1 niveau de pression931 ! EXTERNAL moy_undefSTD !moyenne d'1 var a 1 niveau de pression932 933 ! EXTERNAL moyglo_aire 933 934 ! moyenne globale d'1 var ponderee par l'aire de la maille (moyglo_pondaire) … … 1655 1656 ENDIF 1656 1657 1657 IF (prt_level>=1) print*,'CONVERGENCE PHYSIQUE THERM 1 '1658 IF (prt_level>=1) PRINT *,'CONVERGENCE PHYSIQUE THERM 1 ' 1658 1659 1659 1660 !====================================================================== … … 1804 1805 IF (nbapp_cv == 0) nbapp_cv=86400./phys_tstep 1805 1806 IF (nbapp_wk == 0) nbapp_wk=86400./phys_tstep 1806 print*,'physiq, nbapp_cv, nbapp_wk ',nbapp_cv,nbapp_wk1807 PRINT *,'physiq, nbapp_cv, nbapp_wk ',nbapp_cv,nbapp_wk 1807 1808 IF (MOD(NINT(86400./phys_tstep),nbapp_cv)==0) THEN 1808 1809 cvpas_0 = NINT( 86400./phys_tstep)/nbapp_cv 1809 1810 cvpas = cvpas_0 1810 print*,'physiq, cvpas ',cvpas1811 PRINT *,'physiq, cvpas ',cvpas 1811 1812 ELSE 1812 1813 WRITE(lunout,*) 'le nombre de pas de temps physique doit etre un ', & … … 1820 1821 IF (MOD(NINT(86400./phys_tstep),nbapp_wk)==0) THEN 1821 1822 wkpas = NINT( 86400./phys_tstep)/nbapp_wk 1822 ! print*,'physiq, wkpas ',wkpas1823 ! PRINT *,'physiq, wkpas ',wkpas 1823 1824 ELSE 1824 1825 WRITE(lunout,*) 'le nombre de pas de temps physique doit etre un ', & … … 3214 3215 !!jyg IF (prt_level .ge. 2 .AND. mod(iflag_pbl_split,2) .EQ. 1) THEN 3215 3216 IF (prt_level >= 2 .AND. mod(iflag_pbl_split,10) >= 1) THEN 3216 print*,'debut du splitting de la PBL, wake_s = ', wake_s(:)3217 print*,'debut du splitting de la PBL, wake_deltat = ', wake_deltat(:,1)3218 print*,'debut du splitting de la PBL, wake_deltaq = ', wake_deltaq(:,1)3217 PRINT *,'debut du splitting de la PBL, wake_s = ', wake_s(:) 3218 PRINT *,'debut du splitting de la PBL, wake_deltat = ', wake_deltat(:,1) 3219 PRINT *,'debut du splitting de la PBL, wake_deltaq = ', wake_deltaq(:,1) 3219 3220 ENDIF 3220 3221 ! !! … … 3574 3575 3575 3576 !!jyg IF (MOD(itapcv,cvpas).EQ.0) THEN 3576 !! print*,' physiq : itapcv, cvpas, itap-1, cvpas_0 ', &3577 !! PRINT *,' physiq : itapcv, cvpas, itap-1, cvpas_0 ', & 3577 3578 !! itapcv, cvpas, itap-1, cvpas_0 3578 3579 IF (MOD(itapcv,cvpas)==0 .OR. MOD(itap-1,cvpas_0)==0) THEN … … 3805 3806 !>jyg 3806 3807 3807 !! print*,'physiq. q_w(1,k), q_x(1,k) ', &3808 !! PRINT *,'physiq. q_w(1,k), q_x(1,k) ', & 3808 3809 !! (k, q_w(1,k), q_x(1,k),k=1,25) 3809 3810 … … 4025 4026 abs(d_q_con(i,k)) > dqcon_multistep_max) THEN 4026 4027 cvpas = 1 4027 !! print*,'physiq1, i,k,d_t_con(i,k),d_q_con(i,k) ', &4028 !! PRINT *,'physiq1, i,k,d_t_con(i,k),d_q_con(i,k) ', & 4028 4029 !! i,k,d_t_con(i,k),d_q_con(i,k) 4029 4030 ENDIF … … 4177 4178 ENDIF 4178 4179 IF (mod(itap,1000) == 0) THEN 4179 print*,' physiq, nombre de steps ou cvpas = 1 : ', Ncvpaseq14180 PRINT *,' physiq, nombre de steps ou cvpas = 1 : ', Ncvpaseq1 4180 4181 ENDIF 4181 4182 … … 4451 4452 4452 4453 IF (prt_level >= 10) THEN 4453 print*,' physiq, after calwake, wake_s: ',wake_s(:)4454 print*,' physiq, after calwake, wake_deltat: ',wake_deltat(:,1)4455 print*,' physiq, after calwake, wake_deltaq: ',wake_deltaq(:,1)4454 PRINT *,' physiq, after calwake, wake_s: ',wake_s(:) 4455 PRINT *,' physiq, after calwake, wake_deltat: ',wake_deltat(:,1) 4456 PRINT *,' physiq, after calwake, wake_deltaq: ',wake_deltaq(:,1) 4456 4457 ENDIF 4457 4458 … … 4933 4934 !------------------------------------------------------------------------- 4934 4935 IF (prt_level >=10) THEN 4935 print*,'itap, ->fisrtilp ',itap4936 PRINT *,'itap, ->fisrtilp ',itap 4936 4937 ENDIF 4937 4938 … … 5932 5933 5933 5934 IF (prt_level >=10) THEN 5934 print*,' ->radlwsw, number 1 '5935 PRINT *,' ->radlwsw, number 1 ' 5935 5936 ENDIF 5936 5937 … … 6011 6012 6012 6013 IF (prt_level >=10) THEN 6013 print*,' ->radlwsw, number 2 '6014 PRINT *,' ->radlwsw, number 2 ' 6014 6015 ENDIF 6015 6016 … … 6195 6196 6196 6197 IF (prt_level >=10) THEN 6197 print*,' CALL orography ? ', ok_orodr6198 PRINT *,' CALL orography ? ', ok_orodr 6198 6199 ENDIF 6199 6200 … … 6447 6448 flh2o=1 6448 6449 IF(flag_verbose_strataer) THEN 6449 print*,'IN physiq_mod: ok_qemiss =yes (',ok_qemiss,'), flh2o=',flh2o6450 print*,'IN physiq_mod: flag_emit=',flag_emit,', nErupt=',nErupt6451 print*,'IN physiq_mod: nAerErupt=',nAerErupt6450 PRINT *,'IN physiq_mod: ok_qemiss =yes (',ok_qemiss,'), flh2o=',flh2o 6451 PRINT *,'IN physiq_mod: flag_emit=',flag_emit,', nErupt=',nErupt 6452 PRINT *,'IN physiq_mod: nAerErupt=',nAerErupt 6452 6453 ENDIF 6453 6454 … … 6460 6461 day_cur<(day_emit_vol(ieru)+injdur)) THEN 6461 6462 6462 IF(flag_verbose_strataer) print*,'IN physiq_mod: date=',year_cur,mth_cur,day_cur6463 IF(flag_verbose_strataer) PRINT *,'IN physiq_mod: date=',year_cur,mth_cur,day_cur 6463 6464 ! initialisation of q tendency emission 6464 6465 d_q_emiss(:,:)=0. … … 6475 6476 nAerErupt+1,0) 6476 6477 6477 IF(flag_verbose_strataer) print*,'IN physiq_mod: min max d_q_emiss=',&6478 IF(flag_verbose_strataer) PRINT *,'IN physiq_mod: min max d_q_emiss=',& 6478 6479 minval(d_q_emiss),maxval(d_q_emiss) 6479 6480 … … 7057 7058 7058 7059 IF (prt_level >=10) THEN 7059 print*,'Convertir les incrementations en tendances '7060 PRINT *,'Convertir les incrementations en tendances ' 7060 7061 ENDIF 7061 7062 -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/reevap.F90
r5153 r5160 41 41 ivapcur= iqWIsoPha(ixt,ivap) 42 42 isolcur= iqWIsoPha(ixt,isol) 43 ! print*,'rrevap ; fl_cor_ebil:',fl_cor_ebil,' iflag_ice_thermo:',iflag_ice_thermo,' RVTMP2',RVTMP243 !PRINT *,'rrevap ; fl_cor_ebil:',fl_cor_ebil,' iflag_ice_thermo:',iflag_ice_thermo,' RVTMP2',RVTMP2 44 44 DO k = 1, klev ! re-evaporation de l'eau liquide nuageuse 45 45 DO i = 1, klon
Note: See TracChangeset
for help on using the changeset viewer.
