Changeset 2286 for LMDZ5/trunk/libf/dyn3dmem
- Timestamp:
- May 20, 2015, 3:27:07 PM (9 years ago)
- Location:
- LMDZ5/trunk/libf/dyn3dmem
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3dmem/advtrac_loc.F
r2270 r2286 155 155 c$OMP BARRIER 156 156 157 write(*,*) 'advtrac 157: appel de vlspltgen_loc'157 !write(*,*) 'advtrac 157: appel de vlspltgen_loc' 158 158 call vlspltgen_loc( q,iadv, 2., massem, wg , 159 159 * pbarug,pbarvg,dtvr,p, 160 160 * pk,teta ) 161 161 162 write(*,*) 'advtrac 162: apres appel vlspltgen_loc'162 !write(*,*) 'advtrac 162: apres appel vlspltgen_loc' 163 163 if (ok_iso_verif) then 164 164 call check_isotopes(q,ijb_u,ije_u,'advtrac 162') -
LMDZ5/trunk/libf/dyn3dmem/caladvtrac_loc.F
r2270 r2286 56 56 !$OMP THREADPRIVATE(Request_vanleer) 57 57 58 write(*,*) 'caladvtrac 58: entree'58 !write(*,*) 'caladvtrac 58: entree' 59 59 ijbu=ij_begin 60 60 ijeu=ij_end … … 109 109 110 110 IF ( iadvtr.EQ.iapp_tracvl ) THEN 111 write(*,*) 'caladvtrac 133'111 !write(*,*) 'caladvtrac 133' 112 112 c$OMP MASTER 113 113 call suspend_timer(timer_caldyn) … … 184 184 CALL WriteField_u('wg1',wg_adv) 185 185 #endif 186 write(*,*) 'caladvtrac 185'186 !write(*,*) 'caladvtrac 185' 187 187 CALL advtrac_loc( pbarug_adv,pbarvg_adv,wg_adv, 188 188 * p_adv, massem_adv,q_adv, teta_adv, 189 189 . pk_adv) 190 write(*,*) 'caladvtrac 189'190 !write(*,*) 'caladvtrac 189' 191 191 192 192 -
LMDZ5/trunk/libf/dyn3dmem/integrd_loc.F
r2281 r2286 164 164 c$OMP END MASTER 165 165 c$OMP BARRIER 166 write(*,*) 'integrd 170'166 !write(*,*) 'integrd 170' 167 167 IF (.NOT. Checksum_all) THEN 168 168 call WriteField_v('int_vcov',vcov) … … 194 194 195 195 c 196 write(*,*) 'integrd 200'196 !write(*,*) 'integrd 200' 197 197 C$OMP MASTER 198 198 if (pole_nord) THEN … … 285 285 c 286 286 c 287 write(*,*) 'integrd 291'287 !write(*,*) 'integrd 291' 288 288 IF (pole_nord) THEN 289 289 … … 352 352 endif !if (ok_iso_verif) then 353 353 354 write(*,*) 'integrd 341'354 !write(*,*) 'integrd 341' 355 355 CALL qminimum_loc( q, nq, deltap ) 356 write(*,*) 'integrd 343'356 !write(*,*) 'integrd 343' 357 357 358 358 if (ok_iso_verif) then -
LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F
r2270 r2286 308 308 309 309 1 CONTINUE ! Matsuno Forward step begins here 310 write(*,*) 'leapfrog 298: itau=',itau310 !write(*,*) 'leapfrog 298: itau=',itau 311 311 jD_cur = jD_ref + day_ini - day_ref + & 312 312 & itau/day_step … … 713 713 IF( forward. OR . leapf ) THEN 714 714 ! Ehouarn: NB: fields sent to advtrac are those at the beginning of the time step 715 write(*,*) 'leapfrog 679: avant CALL caladvtrac_loc'715 !write(*,*) 'leapfrog 679: avant CALL caladvtrac_loc' 716 716 CALL caladvtrac_loc(q,pbaru,pbarv, 717 717 * p, masse, dq, teta, 718 718 . flxw,pk, iapptrac) 719 719 720 write(*,*) 'leapfrog 719'720 !write(*,*) 'leapfrog 719' 721 721 if (ok_iso_verif) then 722 722 call check_isotopes(q,ijb_u,ije_u, … … 755 755 ! CALL FTRACE_REGION_BEGIN("integrd") 756 756 757 write(*,*) 'leapfrog 720'757 !write(*,*) 'leapfrog 720' 758 758 if (ok_iso_verif) then 759 759 call check_isotopes(q,ijb_u,ije_u,'leapfrog 756') … … 765 765 ! $ finvmaold ) 766 766 767 write(*,*) 'leapfrog 724'767 !write(*,*) 'leapfrog 724' 768 768 if (ok_iso_verif) then 769 769 call check_isotopes(q,ijb_u,ije_u,'leapfrog 762') … … 1148 1148 call check_isotopes(q,ijb_u,ije_u,'leapfrog 1132') 1149 1149 endif !if (ok_iso_verif) then 1150 write(*,*) 'leapfrog 1134: iflag_phys=',iflag_phys1150 !write(*,*) 'leapfrog 1134: iflag_phys=',iflag_phys 1151 1151 1152 1152 IF(iflag_phys.EQ.2) THEN ! "Newtonian" case … … 1221 1221 c dissipation horizontale et verticale des petites echelles: 1222 1222 c ---------------------------------------------------------- 1223 write(*,*) 'leapfrog 1163: apdiss=',apdiss1223 !write(*,*) 'leapfrog 1163: apdiss=',apdiss 1224 1224 IF(apdiss) THEN 1225 1225 -
LMDZ5/trunk/libf/dyn3dmem/qminimum_loc.F
r2285 r2286 50 50 c 51 51 52 write(*,*) 'qminimum 52: entree'52 !write(*,*) 'qminimum 52: entree' 53 53 if (ok_iso_verif) then 54 54 call check_isotopes(q,ij_begin,ij_end,'qminimum 52') -
LMDZ5/trunk/libf/dyn3dmem/vlsplt_loc.F
r2281 r2286 56 56 INTEGER ijb,ije,ijb_x,ije_x 57 57 58 write(*,*) 'vlsplt 58: entree dans vlx_loc, iq,ijb_x=',59 & iq,ijb_x58 !write(*,*) 'vlsplt 58: entree dans vlx_loc, iq,ijb_x=', 59 ! & iq,ijb_x 60 60 c calcul de la pente a droite et a gauche de la maille 61 61 … … 154 154 ENDIF ! (pente_max.lt.-1.e-5) 155 155 156 write(*,*) 'vlx 156: iq,ijb_x=',iq,ijb_x156 !write(*,*) 'vlx 156: iq,ijb_x=',iq,ijb_x 157 157 158 158 c bouclage de la pente en iip1: … … 273 273 ENDDO 274 274 niju=iju 275 PRINT*,'vlx 278, niju,nl',niju,nl(l)275 !PRINT*,'vlx 278, niju,nl',niju,nl(l) 276 276 277 277 c traitement des mailles … … 331 331 ! Il faut faire ça avant d'avoir mis à jour q et masse 332 332 333 write(*,*) 'vlsplt 326: iq,ijb_x,nqfils(iq)=',iq,ijb_x,nqfils(iq)333 !write(*,*) 'vlsplt 326: iq,ijb_x,nqfils(iq)=',iq,ijb_x,nqfils(iq) 334 334 335 335 if (nqfils(iq).gt.0) then … … 354 354 ! end CRisi 355 355 356 write(*,*) 'vlsplt 360: iq,ijb_x=',iq,ijb_x356 !write(*,*) 'vlsplt 360: iq,ijb_x=',iq,ijb_x 357 357 358 358 c calcul des tENDances … … 373 373 ENDDO 374 374 c$OMP END DO NOWAIT 375 write(*,*) 'vlsplt 380: iq,ijb_x=',iq,ijb_x375 !write(*,*) 'vlsplt 380: iq,ijb_x=',iq,ijb_x 376 376 377 377 ! retablir les fils en rapport de melange par rapport a l'air: … … 394 394 endif !if (nqfils(iq).gt.0) then 395 395 396 write(*,*) 'vlsplt 399: iq,ijb_x=',iq,ijb_x396 !write(*,*) 'vlsplt 399: iq,ijb_x=',iq,ijb_x 397 397 c CALL SCOPY((jjm-1)*llm,q(iip1+iip1,1),iip1,q(iip2,1),iip1) 398 398 c CALL SCOPY((jjm-1)*llm,masse(iip1+iip1,1),iip1,masse(iip2,1),iip1) … … 726 726 ! CRisi: appel récursif de l'advection sur les fils. 727 727 ! Il faut faire ça avant d'avoir mis à jour q et masse 728 write(*,*) 'vly 689: iq,nqfils(iq)=',iq,nqfils(iq)728 !write(*,*) 'vly 689: iq,nqfils(iq)=',iq,nqfils(iq) 729 729 730 730 ijb=ij_begin-2*iip1 … … 926 926 c sens de W 927 927 928 write(*,*) 'vlsplt 926: entree dans vlz_loc, iq=',iq928 !write(*,*) 'vlsplt 926: entree dans vlz_loc, iq=',iq 929 929 #ifdef BIDON 930 930 IF(testcpu) THEN … … 982 982 c calcul de - d( q * w )/ d(sigma) qu'on ajoute a dq pour calculer dq 983 983 984 write(*,*) 'vlz 982,ijb,ije=',ijb,ije984 !write(*,*) 'vlz 982,ijb,ije=',ijb,ije 985 985 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 986 986 DO l = 1,llm-1 … … 998 998 ENDDO 999 999 c$OMP END DO NOWAIT 1000 write(*,*) 'vlz 1001'1000 !write(*,*) 'vlz 1001' 1001 1001 1002 1002 c$OMP MASTER … … 1010 1010 ! CRisi: appel récursif de l'advection sur les fils. 1011 1011 ! Il faut faire ça avant d'avoir mis à jour q et masse 1012 write(*,*) 'vlsplt 942: iq,nqfils(iq)=',iq,nqfils(iq)1012 !write(*,*) 'vlsplt 942: iq,nqfils(iq)=',iq,nqfils(iq) 1013 1013 if (nqfils(iq).gt.0) then 1014 1014 do ifils=1,nqdesc(iq) -
LMDZ5/trunk/libf/dyn3dmem/vlspltgen_loc.F
r2281 r2286 202 202 ! DO iq=1,nqtot 203 203 DO iq=1,nqperes ! CRisi: on ne boucle que sur les pères= ceux qui sont transportés directement par l'air 204 write(*,*) 'vlspltgen 192: iq,iadv=',iq,iadv(iq)204 !write(*,*) 'vlspltgen 192: iq,iadv=',iq,iadv(iq) 205 205 #ifdef DEBUG_IO 206 206 CALL WriteField_u('zq',zq(:,:,iq)) … … 298 298 299 299 do iq=1,nqperes 300 write(*,*) 'vlspltgen 279: iq=',iq300 !write(*,*) 'vlspltgen 279: iq=',iq 301 301 302 302 if(iadv(iq) == 0) then … … 351 351 352 352 do iq=1,nqperes 353 write(*,*) 'vlspltgen 321: iq=',iq353 !write(*,*) 'vlspltgen 321: iq=',iq 354 354 #ifdef DEBUG_IO 355 355 CALL WriteField_u('zq',zq(:,:,iq)) … … 383 383 384 384 do iq=1,nqperes 385 write(*,*) 'vlspltgen 349: iq=',iq385 !write(*,*) 'vlspltgen 349: iq=',iq 386 386 #ifdef DEBUG_IO 387 387 CALL WriteField_u('zq',zq(:,:,iq)) … … 448 448 c$OMP BARRIER 449 449 do iq=1,nqperes 450 write(*,*) 'vlspltgen 409: iq=',iq450 !write(*,*) 'vlspltgen 409: iq=',iq 451 451 452 452 if(iadv(iq) == 0) then … … 470 470 471 471 enddo 472 write(*,*) 'vlspltgen_loc 476'473 474 c$OMP BARRIER 475 write(*,*) 'vlspltgen_loc 477'472 !write(*,*) 'vlspltgen_loc 476' 473 474 c$OMP BARRIER 475 !write(*,*) 'vlspltgen_loc 477' 476 476 c$OMP MASTER 477 477 call VTb(VTHallo) … … 489 489 490 490 491 write(*,*) 'vlspltgen_loc 494'491 !write(*,*) 'vlspltgen_loc 494' 492 492 if (ok_iso_verif) then 493 493 call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 461') … … 495 495 496 496 do iq=1,nqperes 497 write(*,*) 'vlspltgen 449: iq=',iq497 !write(*,*) 'vlspltgen 449: iq=',iq 498 498 #ifdef DEBUG_IO 499 499 CALL WriteField_u('zq',zq(:,:,iq)) … … 526 526 527 527 do iq=1,nqperes 528 write(*,*) 'vlspltgen 477: iq=',iq528 !write(*,*) 'vlspltgen 477: iq=',iq 529 529 #ifdef DEBUG_IO 530 530 CALL WriteField_u('zq',zq(:,:,iq)) … … 553 553 enddo !do iq=1,nqperes 554 554 555 write(*,*) 'vlspltgen 550: apres derniere serie de call vlx'555 !write(*,*) 'vlspltgen 550: apres derniere serie de call vlx' 556 556 if (ok_iso_verif) then 557 557 call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 521') … … 560 560 ijb=ij_begin 561 561 ije=ij_end 562 write(*,*) 'vlspltgen_loc 557'562 !write(*,*) 'vlspltgen_loc 557' 563 563 c$OMP BARRIER 564 564 565 write(*,*) 'vlspltgen_loc 559'565 !write(*,*) 'vlspltgen_loc 559' 566 566 DO iq=1,nqtot 567 write(*,*) 'vlspltgen_loc 561, iq=',iq567 !write(*,*) 'vlspltgen_loc 561, iq=',iq 568 568 #ifdef DEBUG_IO 569 569 CALL WriteField_u('zq',zq(:,:,iq)) … … 579 579 ENDDO 580 580 c$OMP END DO NOWAIT 581 write(*,*) 'vlspltgen_loc 575'581 !write(*,*) 'vlspltgen_loc 575' 582 582 583 583 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) … … 588 588 ENDDO 589 589 c$OMP END DO NOWAIT 590 write(*,*) 'vlspltgen_loc 583'590 !write(*,*) 'vlspltgen_loc 583' 591 591 ENDDO !DO iq=1,nqtot 592 592 … … 603 603 cc$OMP BARRIER 604 604 605 write(*,*) 'vlspltgen 597: sortie'605 !write(*,*) 'vlspltgen 597: sortie' 606 606 RETURN 607 607 END -
LMDZ5/trunk/libf/dyn3dmem/vlspltqs_loc.F
r2281 r2286 47 47 INTEGER ijb,ije,ijb_x,ije_x 48 48 49 write(*,*) 'vlspltqs 58: entree vlxqs_loc, iq,ijb_x=',50 & iq,ijb_x49 !write(*,*) 'vlspltqs 58: entree vlxqs_loc, iq,ijb_x=', 50 ! & iq,ijb_x 51 51 52 52 c calcul de la pente a droite et a gauche de la maille … … 280 280 ENDDO 281 281 niju=iju 282 PRINT*,'vlxqs 280: niju,nl',niju,nl(l)282 !PRINT*,'vlxqs 280: niju,nl',niju,nl(l) 283 283 284 284 c traitement des mailles … … 336 336 ! CRisi: appel récursif de l'advection sur les fils. 337 337 ! Il faut faire ça avant d'avoir mis à jour q et masse 338 write(*,*) 'vlspltqs 336: iq,ijb_x,nqfils(iq)=',339 & iq,ijb_x,nqfils(iq)338 !write(*,*) 'vlspltqs 336: iq,ijb_x,nqfils(iq)=', 339 ! & iq,ijb_x,nqfils(iq) 340 340 341 341 if (nqfils(iq).gt.0) then … … 353 353 do ifils=1,nqfils(iq) 354 354 iq2=iqfils(ifils,iq) 355 write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2355 !write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2 356 356 call vlx_loc(Ratio,pente_max,masse,u_mq,ijb_x,ije_x,iq2) 357 357 enddo !do ifils=1,nqfils(iq) … … 359 359 ! end CRisi 360 360 361 write(*,*) 'vlspltqs 360: iq,ijb_x=',iq,ijb_x361 !write(*,*) 'vlspltqs 360: iq,ijb_x=',iq,ijb_x 362 362 363 363 c calcul des tendances … … 379 379 c$OMP END DO NOWAIT 380 380 381 write(*,*) 'vlspltqs 380: iq,ijb_x=',iq,ijb_x381 !write(*,*) 'vlspltqs 380: iq,ijb_x=',iq,ijb_x 382 382 383 383 ! retablir les fils en rapport de melange par rapport a l'air: … … 398 398 endif !if (nqfils(iq).gt.0) then 399 399 400 write(*,*) 'vlspltqs 399: iq,ijb_x=',iq,ijb_x400 !write(*,*) 'vlspltqs 399: iq,ijb_x=',iq,ijb_x 401 401 402 402 c CALL SCOPY((jjm-1)*llm,q(iip1+iip1,1),iip1,q(iip2,1),iip1) … … 477 477 l=3 478 478 if ((ij.ge.ijb).and.(ij.le.ije)) then 479 write(*,*) 'vlyqs 480: ij,l,iq,ijb,q(ij,l,:)=',480 & ij,l,iq,ijb,q(ij,l,:)479 !write(*,*) 'vlyqs 480: ij,l,iq,ijb,q(ij,l,:)=', 480 ! & ij,l,iq,ijb,q(ij,l,:) 481 481 endif 482 482 … … 725 725 ! CRisi: appel récursif de l'advection sur les fils. 726 726 ! Il faut faire ça avant d'avoir mis à jour q et masse 727 write(*,*) 'vlyqs 689: iq,nqfils(iq)=',iq,nqfils(iq)727 !write(*,*) 'vlyqs 689: iq,nqfils(iq)=',iq,nqfils(iq) 728 728 729 729 ijb=ij_begin-2*iip1
Note: See TracChangeset
for help on using the changeset viewer.