Changeset 938
- Timestamp:
- Mar 31, 2008, 6:14:45 PM (17 years ago)
- Location:
- LMDZ4/trunk/libf
- Files:
-
- 17 edited
-
dyn3d/guide.F (modified) (2 diffs)
-
dyn3d/wrgrads.F (modified) (4 diffs)
-
phylmd/calltherm.F90 (modified) (5 diffs)
-
phylmd/concvl.F (modified) (1 diff)
-
phylmd/cva_driver.F (modified) (6 diffs)
-
phylmd/thermcell.F (modified) (6 diffs)
-
phylmd/thermcell_closure.F90 (modified) (2 diffs)
-
phylmd/thermcell_dq.F90 (modified) (2 diffs)
-
phylmd/thermcell_dry.F90 (modified) (4 diffs)
-
phylmd/thermcell_dv2.F90 (modified) (1 diff)
-
phylmd/thermcell_env.F90 (modified) (2 diffs)
-
phylmd/thermcell_flux.F90 (modified) (17 diffs)
-
phylmd/thermcell_height.F90 (modified) (1 diff)
-
phylmd/thermcell_init.F90 (modified) (1 diff)
-
phylmd/thermcell_main.F90 (modified) (25 diffs)
-
phylmd/thermcell_plume.F90 (modified) (5 diffs)
-
phylmd/yamada4.F (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/dyn3d/guide.F
r762 r938 122 122 c calcul de l'humidite saturante 123 123 C----------------------------------------------------------------------- 124 124 c print*,'OK0' 125 125 CALL pression( ip1jmp1, ap, bp, ps, p ) 126 126 call massdair(p,masse) 127 127 c print*,'OK1' 128 128 CALL exner_hyb(ip1jmp1,ps,p,alpha,beta,pks,pk,pkf) 129 129 c print*,'OK2' 130 130 tnat(:,:)=pk(:,:)*teta(:,:)/cpp 131 131 c print*,'OK3' 132 132 unskap = 1./ kappa 133 133 pres(:,:)=preff*(pk(:,:)/cpp)**unskap 134 134 c print*,'OK4' 135 135 call q_sat(iip1*jjp1*llm,tnat,pres,qsat) 136 136 … … 144 144 c----------------------------------------------------------------------- 145 145 146 146 c print*,'ONLINE=',online 147 147 if(online.eq.-1) then 148 148 return -
LMDZ4/trunk/libf/dyn3d/wrgrads.F
r541 r938 29 29 writectl=.false. 30 30 31 31 c print*,if,iid(if),jid(if),ifd(if),jfd(if) 32 32 iii=iid(if) 33 33 iji=jid(if) … … 38 38 lm=lmd(if) 39 39 40 41 40 c print*,'im,jm,lm,name,firsttime(if)' 41 c print*,im,jm,lm,name,firsttime(if) 42 42 43 43 if(firsttime(if)) then … … 57 57 tvar(ivar(if),if)=titlevar(1:lnblnk(titlevar)) 58 58 nld(ivar(if),if)=nl 59 60 59 c print*,'initialisation ecriture de ',var(ivar(if),if) 60 c print*,'if ivar(if) nld ',if,ivar(if),nld(ivar(if),if) 61 61 endif 62 62 writectl=.true. … … 81 81 endif 82 82 83 84 83 c print*,'ivar(if),nvar(if),var(ivar(if),if),writectl' 84 c print*,ivar(if),nvar(if),var(ivar(if),if),writectl 85 85 do l=1,nl 86 86 irec(if)=irec(if)+1 -
LMDZ4/trunk/libf/phylmd/calltherm.F90
r927 r938 14 14 #include "dimphy.h" 15 15 #include "thermcell.h" 16 #include "iniprint.h" 16 17 17 18 ! A inclure eventuellement dans les fichiers de configuration … … 73 74 74 75 integer i,k 75 76 76 !******************************************************** 77 77 … … 85 85 Ale_bl(:)=0. 86 86 Alp_bl(:)=0. 87 print*,'thermV4 nsplit: ',nsplit_thermals,' weak_inversion' 88 87 if (prt_level.ge.10) then 88 print*,'thermV4 nsplit: ',nsplit_thermals,' weak_inversion' 89 endif 89 90 90 91 ! tests sur les valeurs negatives de l'eau … … 92 93 do i=1,klon 93 94 if (.not.q_seri(i,k).ge.0.) then 95 if (prt_level.ge.10) then 94 96 print*,'WARN eau<0 avant therm i=',i,' k=',k & 95 97 & ,' dq,q',d_q_the(i,k),q_seri(i,k) 96 q_seri(i,k)=1.e-15 98 endif 99 q_seri(i,k)=1.e-15 97 100 endif 98 101 enddo … … 224 227 DO i = 1, klon 225 228 if (.not.q_seri(i,k).ge.0.) then 229 if (prt_level.ge.10) then 226 230 print*,'WARN eau<0 apres therm i=',i,' k=',k & 227 231 & ,' dq,q',d_q_the(i,k),q_seri(i,k), & 228 232 & 'fm=',zfm_therm(i,k),'entr=',entr_therm(i,k) 229 q_seri(i,k)=1.e-15 233 endif 234 q_seri(i,k)=1.e-15 230 235 ! stop 231 236 endif -
LMDZ4/trunk/libf/phylmd/concvl.F
r879 r938 382 382 ENDDO 383 383 endif 384 384 c print*, 'concvl->: dd_t,dd_q ',dd_t(1,1),dd_q(1,1) 385 385 386 386 RETURN -
LMDZ4/trunk/libf/phylmd/cva_driver.F
r879 r938 305 305 306 306 logical ok_inhib ! True => possible inhibition of convection by dryness 307 logical, save :: debut=.true. 307 308 308 309 real plcl1(klon) … … 507 508 508 509 if (iflag_con.eq.3) then 509 510 print*,'Emanuel version 3 nouvelle' 510 511 if (debut) THEN 512 print*,'Emanuel version 3 nouvelle' 513 endif 514 511 515 CALL cv3_prelim(len,nd,ndp1,t1,q1,p1,ph1 ! nd->na 512 516 o ,lv1,cpn1,tv1,gz1,h1,hm1,th1) … … 774 778 ELSE 775 779 CALL cv3_mixscale(nloc,ncum,nd,ment,m) 776 print *,' cv3_mixscale-> ' 780 if (debut) THEN 781 print *,' cv3_mixscale-> ' 782 endif !(debut) THEN 777 783 ENDIF 778 784 endif … … 785 791 endif 786 792 c 787 print *,' cv_mixing ->' 793 if (debut) THEN 794 print *,' cv_mixing ->' 795 endif !(debut) THEN 788 796 c do i = 1,klev 789 797 c print*,'cv_mixing-> i,ment ',i,(ment(1,i,j),j=1,klev) … … 794 802 !------------------------------------------------------------------- 795 803 if (iflag_con.eq.3) then 796 print *,' cva_driver -> cv3_unsat ' 804 if (debut) THEN 805 print *,' cva_driver -> cv3_unsat ' 806 endif !(debut) THEN 797 807 798 808 CALL cv3_unsat(nloc,ncum,nd,nd,ntra,icb,inb,iflag ! na->nd … … 810 820 endif 811 821 c 812 print *,'cv_unsat-> ' 822 if (debut) THEN 823 print *,'cv_unsat-> ' 824 debut=.FALSE. 825 endif !(debut) THEN 813 826 ! 814 827 c print *,'cv_unsat-> mp ',mp -
LMDZ4/trunk/libf/phylmd/thermcell.F
r878 r938 991 991 #ifdef troisD 992 992 c if (sorties) then 993 print*,'Debut des wrgradsfi'993 c print*,'Debut des wrgradsfi' 994 994 995 995 c print*,'16 OK convect8' … … 1068 1068 c print*,'18 OK convect8' 1069 1069 c endif 1070 print*,'Fin des wrgradsfi'1070 c print*,'Fin des wrgradsfi' 1071 1071 #endif 1072 1072 … … 1167 1167 c detr_star(ig,l)=0. 1168 1168 endif 1169 1169 c print*,'ok detr_star' 1170 1170 cprise en compte du detrainement dans le calcul du flux 1171 1171 f_star(ig,l+1)=f_star(ig,l)+alim_star(ig,l)+entr_star(ig,l) … … 1188 1188 s -zw2(ig,l))/(zw2(ig,l+1)-zw2(ig,l)) 1189 1189 zw2(ig,l+1)=0. 1190 1190 c print*,'linter=',linter(ig) 1191 1191 else 1192 1192 wa_moy(ig,l+1)=sqrt(zw2(ig,l+1)) … … 1199 1199 enddo 1200 1200 enddo 1201 1201 c print*,'fin calcul zw2' 1202 1202 c 1203 1203 c Calcul de la couche correspondant a la hauteur du thermique … … 1230 1230 if (l.le.lmax(ig)) then 1231 1231 if (zw2(ig,l).lt.0.)then 1232 1232 c print*,'pb2 zw2<0' 1233 1233 endif 1234 1234 zw2(ig,l)=sqrt(zw2(ig,l)) -
LMDZ4/trunk/libf/phylmd/thermcell_closure.F90
r878 r938 7 7 IMPLICIT NONE 8 8 9 #include "iniprint.h" 9 10 INTEGER ngrid,nlay 10 11 INTEGER ig,k … … 61 62 f0(ig)=f(ig) 62 63 enddo 63 if ( lev_out.ge.1) print*,'apres fermeture'64 if (prt_level.ge.1) print*,'apres fermeture' 64 65 65 66 ! -
LMDZ4/trunk/libf/phylmd/thermcell_dq.F90
r878 r938 3 3 implicit none 4 4 5 #include "iniprint.h" 5 6 !======================================================================= 6 7 ! … … 26 27 ! calcul du detrainement 27 28 28 if ( lev_out.ge.1) print*,'Q2 THERMCEL_DQ 0'29 if (prt_level.ge.1) print*,'Q2 THERMCEL_DQ 0' 29 30 30 31 do k=1,nlay -
LMDZ4/trunk/libf/phylmd/thermcell_dry.F90
r878 r938 7 7 IMPLICIT NONE 8 8 #include "YOMCST.h" 9 #include "iniprint.h" 9 10 INTEGER l,ig 10 11 … … 120 121 if (zw2(ig,l+1)>0. .and. zw2(ig,l+1).lt.1.e-10) then 121 122 ! stop'On tombe sur le cas particulier de thermcell_dry' 122 123 ! print*,'On tombe sur le cas particulier de thermcell_dry' 123 124 zw2(ig,l+1)=0. 124 125 linter(ig)=l+1 … … 142 143 enddo 143 144 enddo 144 if ( lev_out.ge.1) print*,'fin calcul zw2'145 if (prt_level.ge.1) print*,'fin calcul zw2' 145 146 ! 146 147 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 201 202 do ig=1,ngrid 202 203 if(lalim(ig)>lmax(ig)) then 204 if ( prt_level > 1 ) THEN 203 205 print*,'WARNING thermcell_dry ig=',ig,' lalim=',lalim(ig),' lmax(ig)=',lmax(ig) 204 lmax(ig)=lalim(ig) 206 endif 207 lmax(ig)=lalim(ig) 205 208 endif 206 209 enddo -
LMDZ4/trunk/libf/phylmd/thermcell_dv2.F90
r878 r938 4 4 implicit none 5 5 6 #include "iniprint.h" 6 7 !======================================================================= 7 8 ! -
LMDZ4/trunk/libf/phylmd/thermcell_env.F90
r878 r938 12 12 #include "YOETHF.h" 13 13 #include "FCTTRE.h" 14 #include "iniprint.h" 14 15 15 16 INTEGER ngrid,nlay … … 113 114 ! --------------------------------------------------- 114 115 115 if ( lev_out.ge.1) print*,'0 OK convect8'116 if (prt_level.ge.1) print*,'0 OK convect8' 116 117 117 118 DO 1010 l=1,nlay -
LMDZ4/trunk/libf/phylmd/thermcell_flux.F90
r883 r938 7 7 & lalim,lmax,alim_star, & 8 8 & entr_star,detr_star,f,rhobarz,zlev,zw2,fm,entr, & 9 & detr,zqla,zmax,lev_out,lunout ,igout)9 & detr,zqla,zmax,lev_out,lunout1,igout) 10 10 11 11 … … 15 15 16 16 IMPLICIT NONE 17 #include "iniprint.h" 17 18 18 19 INTEGER ig,l … … 42 43 integer igout 43 44 integer lev_out 44 integer lunout 45 integer lunout1 45 46 46 47 REAL f_old,ddd0,eee0,ddd,eee,zzz … … 65 66 fm(:,:)=0. 66 67 67 if ( lev_out.ge.10) then68 if (prt_level.ge.10) then 68 69 write(lunout,*) 'Dans thermcell_flux 0' 69 70 write(lunout,*) 'flux base ',f(igout) … … 81 82 ! de lmax(ig) 82 83 !------------------------------------------------------------------------- 83 84 if ( prt_level > 1 ) THEN 84 85 do l=1,klev 85 86 do ig=1,ngrid … … 103 104 enddo 104 105 enddo 105 106 endif !( prt_level > 1 ) THEN 106 107 !------------------------------------------------------------------------- 107 108 ! Multiplication par le flux de masse issu de la femreture … … 113 114 enddo 114 115 115 if ( lev_out.ge.10) then116 if (prt_level.ge.10) then 116 117 write(lunout,*) 'Dans thermcell_flux 1' 117 118 write(lunout,*) 'flux base ',f(igout) … … 153 154 enddo 154 155 155 if ( lev_out.ge.10) &156 if (prt_level.ge.10) & 156 157 & call printflux(ngrid,klev,lunout,igout,f,lmax,lalim, & 157 158 & ptimestep,masse,entr,detr,fm,'2 ') … … 197 198 ! enddo 198 199 199 if ( lev_out.ge.10) &200 if (prt_level.ge.10) & 200 201 & write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, & 201 202 & entr(igout,l),detr(igout,l),fm(igout,l+1) … … 223 224 ! enddo 224 225 225 if ( lev_out.ge.10) &226 if (prt_level.ge.10) & 226 227 & write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, & 227 228 & entr(igout,l),detr(igout,l),fm(igout,l+1) 228 229 else 230 if (l.eq.1) then 229 231 print*,'Test sur les fractions croissantes inhibe dans thermcell_flux2' 232 endif 230 233 endif 231 234 … … 246 249 ! enddo 247 250 248 if ( lev_out.ge.10) &251 if (prt_level.ge.10) & 249 252 & write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, & 250 253 & entr(igout,l),detr(igout,l),fm(igout,l+1) … … 295 298 296 299 297 if ( lev_out.ge.10) &300 if (prt_level.ge.10) & 298 301 & write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, & 299 302 & entr(igout,l),detr(igout,l),fm(igout,l+1) … … 320 323 ! enddo 321 324 322 if ( lev_out.ge.10) &325 if (prt_level.ge.10) & 323 326 & write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, & 324 327 & entr(igout,l),detr(igout,l),fm(igout,l+1) … … 364 367 365 368 366 if ( lev_out.ge.10) &369 if (prt_level.ge.10) & 367 370 & write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, & 368 371 & entr(igout,l),detr(igout,l),fm(igout,l+1) … … 371 374 enddo 372 375 373 if ( lev_out.ge.10) &376 if (prt_level.ge.10) & 374 377 & call printflux(ngrid,klev,lunout,igout,f,lmax,lalim, & 375 378 & ptimestep,masse,entr,detr,fm,'8 ') … … 442 445 443 446 if (ncorecfm1+ncorecfm2+ncorecfm3+ncorecfm4+ncorecfm5+ncorecalpha > 0 ) then 444 447 if (prt_level.ge.10) then 445 448 print*,'PB thermcell : on a du coriger ',ncorecfm1,'x fm1',& 446 449 & ncorecfm2,'x fm2',ncorecfm3,'x fm3 et', & … … 450 453 & ncorecfm8,'x fm8', & 451 454 & ncorecalpha,'x alpha' 455 endif 452 456 endif 453 457 454 if ( lev_out.ge.10) &458 if (prt_level.ge.10) & 455 459 & call printflux(ngrid,klev,lunout,igout,f,lmax,lalim, & 456 460 & ptimestep,masse,entr,detr,fm,'fin') -
LMDZ4/trunk/libf/phylmd/thermcell_height.F90
r878 r938 6 6 !----------------------------------------------------------------------------- 7 7 IMPLICIT NONE 8 #include "iniprint.h" 8 9 9 10 INTEGER ig,l -
LMDZ4/trunk/libf/phylmd/thermcell_init.F90
r878 r938 6 6 !---------------------------------------------------------------------- 7 7 IMPLICIT NONE 8 #include "iniprint.h" 9 8 10 INTEGER l,ig 9 11 !arguments d entree -
LMDZ4/trunk/libf/phylmd/thermcell_main.F90
r927 r938 44 44 #include "YOETHF.h" 45 45 #include "FCTTRE.h" 46 #include "iniprint.h" 46 47 47 48 ! arguments: … … 61 62 62 63 integer,save :: igout=1 63 integer,save :: lunout =664 integer,save :: lunout1=6 64 65 integer,save :: lev_out=10 65 66 … … 162 163 seuil=0.25 163 164 164 if ( lev_out.ge.1) print*,'thermcell_main V4'165 if (prt_level.ge.1) print*,'thermcell_main V4' 165 166 166 167 sorties=.true. … … 191 192 & pplev,zo,zh,zl,ztv,zthl,zu,zv,zpspsk,zqsat,lev_out) 192 193 193 if ( lev_out.ge.1) print*,'thermcell_main apres thermcell_env'194 if (prt_level.ge.1) print*,'thermcell_main apres thermcell_env' 194 195 195 196 !------------------------------------------------------------------------ … … 246 247 enddo 247 248 248 if ( lev_out.ge.1) print*,'thermcell_main apres initialisation'249 if (prt_level.ge.1) print*,'thermcell_main apres initialisation' 249 250 250 251 !------------------------------------------------------------------ … … 314 315 315 316 316 if ( lev_out.ge.1) print*,'thermcell_main apres thermcell_init'317 if ( lev_out.ge.10) then317 if (prt_level.ge.1) print*,'thermcell_main apres thermcell_init' 318 if (prt_level.ge.10) then 318 319 write(lunout,*) 'Dans thermcell_main 1' 319 320 write(lunout,*) 'lmin ',lmin(igout) … … 324 325 endif 325 326 326 do ig=1,klon327 !v1d do ig=1,klon 327 328 !v1d if (alim_star(ig,1).gt.1.e-10) then 328 329 !v1d therm=.true. 329 330 !v1d endif 330 enddo331 !v1d enddo 331 332 !----------------------------------------------------------------------------- 332 333 ! 3. wmax_sec et zmax_sec sont les vitesses et altitudes maximum d'un … … 342 343 call test_ltherm(ngrid,nlay,pplev,pplay,lalim,seuil,ztv,po,ztva,zqla,f_star,zw2,'thermcell_dry lalim ') 343 344 344 if ( lev_out.ge.1) print*,'thermcell_main apres thermcell_dry'345 if ( lev_out.ge.10) then345 if (prt_level.ge.1) print*,'thermcell_main apres thermcell_dry' 346 if (prt_level.ge.10) then 346 347 write(lunout,*) 'Dans thermcell_main 1b' 347 348 write(lunout,*) 'lmin ',lmin(igout) … … 365 366 call test_ltherm(ngrid,nlay,pplev,pplay,lmix ,seuil,ztv,po,ztva,zqla,f_star,zw2,'thermcell_plum lmix ') 366 367 367 if ( lev_out.ge.1) print*,'thermcell_main apres thermcell_plume'368 if ( lev_out.ge.10) then368 if (prt_level.ge.1) print*,'thermcell_main apres thermcell_plume' 369 if (prt_level.ge.10) then 369 370 write(lunout,*) 'Dans thermcell_main 2' 370 371 write(lunout,*) 'lmin ',lmin(igout) … … 388 389 call test_ltherm(ngrid,nlay,pplev,pplay,lmax ,seuil,ztv,po,ztva,zqla,f_star,zw2,'thermcell_heig lmax ') 389 390 390 if ( lev_out.ge.1) print*,'thermcell_main apres thermcell_height'391 if (prt_level.ge.1) print*,'thermcell_main apres thermcell_height' 391 392 392 393 !------------------------------------------------------------------------------- … … 397 398 & zlev,lalim,alim_star,zmax_sec,wmax_sec,zmax,wmax,f,f0,lev_out) 398 399 399 if( lev_out.ge.1)print*,'thermcell_closure apres thermcell_closure'400 if(prt_level.ge.1)print*,'thermcell_closure apres thermcell_closure' 400 401 401 402 !------------------------------------------------------------------------------- … … 408 409 & detr,zqla,zmax,lev_out,lunout,igout) 409 410 410 if ( lev_out.ge.1) print*,'thermcell_main apres thermcell_flux'411 if (prt_level.ge.1) print*,'thermcell_main apres thermcell_flux' 411 412 call test_ltherm(ngrid,nlay,pplev,pplay,lalim,seuil,ztv,po,ztva,zqla,f_star,zw2,'thermcell_flux lalim ') 412 413 call test_ltherm(ngrid,nlay,pplev,pplay,lmax ,seuil,ztv,po,ztva,zqla,f_star,zw2,'thermcell_flux lmax ') … … 476 477 enddo 477 478 478 479 ! print*,'14 OK convect8' 479 480 !------------------------------------------------------------------ 480 481 ! Calculs de diagnostiques pour les sorties … … 483 484 484 485 if (sorties) then 485 print*,'14a OK convect8'486 ! print*,'14a OK convect8' 486 487 ! calcul du niveau de condensation 487 488 ! initialisation … … 503 504 enddo 504 505 enddo 505 506 ! print*,'14b OK convect8' 506 507 do k=nlay,1,-1 507 508 do ig=1,ngrid … … 512 513 enddo 513 514 enddo 514 515 ! print*,'14c OK convect8' 515 516 !calcul des moments 516 517 !initialisation … … 524 525 enddo 525 526 enddo 526 527 ! print*,'14d OK convect8' 527 528 do l=1,nlay 528 529 do ig=1,ngrid … … 636 637 637 638 !calcul du ratqscdiff 638 639 ! print*,'14e OK convect8' 639 640 var=0. 640 641 vardiff=0. … … 645 646 enddo 646 647 enddo 647 648 ! print*,'14f OK convect8' 648 649 do ig=1,ngrid 649 650 do l=1,lalim(ig) … … 656 657 enddo 657 658 enddo 658 659 ! print*,'14g OK convect8' 659 660 do l=1,nlay 660 661 do ig=1,ngrid … … 672 673 673 674 #ifdef und 674 if ( lev_out.ge.1) print*,'thermcell_main sorties 1D'675 if (prt_level.ge.1) print*,'thermcell_main sorties 1D' 675 676 #include "thermcell_out1d.h" 676 677 #endif … … 678 679 679 680 ! #define troisD 680 if ( lev_out.ge.1) print*,'thermcell_main sorties 3D'681 if (prt_level.ge.1) print*,'thermcell_main sorties 3D' 681 682 #ifdef troisD 682 683 #include "thermcell_out3d.h" … … 685 686 endif 686 687 687 if ( lev_out.ge.1) print*,'thermcell_main FIN OK'688 if (prt_level.ge.1) print*,'thermcell_main FIN OK' 688 689 689 690 return … … 693 694 694 695 subroutine test_ltherm(klon,klev,pplev,pplay,long,seuil,ztv,po,ztva,zqla,f_star,zw2,comment) 695 696 integer klon,klev 696 IMPLICIT NONE 697 #include "iniprint.h" 698 699 integer i, k, klon,klev 697 700 real pplev(klon,klev+1),pplay(klon,klev) 698 701 real ztv(klon,klev) … … 706 709 character*21 comment 707 710 708 print*,'WARNING !!! TEST ',comment 711 if (prt_level.ge.1) THEN 712 print*,'WARNING !!! TEST ',comment 713 endif 709 714 return 710 715 -
LMDZ4/trunk/libf/phylmd/thermcell_plume.F90
r878 r938 13 13 #include "YOETHF.h" 14 14 #include "FCTTRE.h" 15 #include "iniprint.h" 15 16 16 17 INTEGER ngrid,klev … … 94 95 enddo 95 96 enddo 96 if ( lev_out.ge.1) print*,'7 OK convect8'97 if (prt_level.ge.1) print*,'7 OK convect8' 97 98 do k=1,klev+1 98 99 do ig=1,ngrid … … 104 105 enddo 105 106 106 if ( lev_out.ge.1) print*,'8 OK convect8'107 if (prt_level.ge.1) print*,'8 OK convect8' 107 108 do ig=1,ngrid 108 109 linter(ig)=1. … … 210 211 & -rhobarz(ig,l)*sqrt(w_est(ig,l))*sqrt(l_mix*zlev(ig,l))) & 211 212 & /(r_aspect*zmax_sec(ig))) 212 if ( lev_out.ge.20) print*,'coucou calcul detr 1'213 if (prt_level.ge.20) print*,'coucou calcul detr 1' 213 214 else if ((zlev(ig,l+1).lt.zmax_sec(ig)).and. & 214 215 & (zqla(ig,l-1).lt.1.e-10)) then … … 221 222 & *((zmax_sec(ig)-zlev(ig,l))/ & 222 223 & ((zmax_sec(ig)-zlev(ig,lmix(ig)))))**2.) 223 if ( lev_out.ge.20) print*,'coucou calcul detr 2'224 if (prt_level.ge.20) print*,'coucou calcul detr 2' 224 225 else 225 226 detr_star(ig,l)=0.002*f0(ig)*f_star(ig,l) & 226 227 & *(zlev(ig,l+1)-zlev(ig,l)) 227 if ( lev_out.ge.20) print*,'coucou calcul detr 3'228 if (prt_level.ge.20) print*,'coucou calcul detr 3' 228 229 229 230 endif -
LMDZ4/trunk/libf/phylmd/yamada4.F
r926 r938 7 7 use dimphy 8 8 IMPLICIT NONE 9 #include "iniprint.h" 9 10 c....................................................................... 10 11 cym#include "dimensions.h" … … 423 424 c D'apres Holtslag Boville. 424 425 425 print*,'YAMADA4 0' 426 426 if (prt_level>1) THEN 427 print*,'YAMADA4 0' 428 endif !(prt_level>1) THEN 427 429 do ig=1,ngrid 428 430 coriol(ig)=1.e-4 … … 456 458 endif 457 459 458 print*,'YAMADA4 1' 460 if (prt_level>1) THEN 461 print*,'YAMADA4 1' 462 endif !(prt_level>1) THEN 459 463 c Diagnostique pour stokage 460 464
Note: See TracChangeset
for help on using the changeset viewer.
