Changeset 3852 for LMDZ6/branches/LMDZ-tracers/libf/dyn3d/leapfrog.F
- Timestamp:
- Feb 22, 2021, 5:28:31 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ-tracers/libf/dyn3d/leapfrog.F
r3416 r3852 11 11 use IOIPSL 12 12 #endif 13 USE infotrac, ONLY: nqtot ,ok_iso_verif13 USE infotrac, ONLY: nqtot 14 14 USE guide_mod, ONLY : guide_main 15 15 USE write_field, ONLY: writefield … … 237 237 jH_cur = jH_cur - int(jH_cur) 238 238 239 if (ok_iso_verif) then 240 call check_isotopes_seq(q,ip1jmp1,'leapfrog 321') 241 endif !if (ok_iso_verif) then 239 call check_isotopes_seq(q,ip1jmp1,'leapfrog 321') 242 240 243 241 #ifdef CPP_IOIPSL … … 271 269 ! CALL filtreg ( finvmaold ,jjp1, llm, -2,2, .TRUE., 1 ) 272 270 273 if (ok_iso_verif) then 274 call check_isotopes_seq(q,ip1jmp1,'leapfrog 400') 275 endif !if (ok_iso_verif) then 271 call check_isotopes_seq(q,ip1jmp1,'leapfrog 400') 276 272 277 273 2 CONTINUE ! Matsuno backward or leapfrog step begins here … … 323 319 endif 324 320 325 326 if (ok_iso_verif) then 327 call check_isotopes_seq(q,ip1jmp1,'leapfrog 589') 328 endif !if (ok_iso_verif) then 321 call check_isotopes_seq(q,ip1jmp1,'leapfrog 589') 329 322 330 323 c----------------------------------------------------------------------- … … 345 338 c ------------------------------------------------------------- 346 339 347 if (ok_iso_verif) then 348 call check_isotopes_seq(q,ip1jmp1, 340 call check_isotopes_seq(q,ip1jmp1, 349 341 & 'leapfrog 686: avant caladvtrac') 350 endif !if (ok_iso_verif) then351 342 352 343 IF( forward. OR . leapf ) THEN … … 376 367 c ---------------------------------- 377 368 378 if (ok_iso_verif) then 379 write(*,*) 'leapfrog 720' 380 call check_isotopes_seq(q,ip1jmp1,'leapfrog 756') 381 endif !if (ok_iso_verif) then 369 call check_isotopes_seq(q,ip1jmp1,'leapfrog 756') 382 370 383 371 CALL integrd ( nqtot,vcovm1,ucovm1,tetam1,psm1,massem1 , … … 385 373 ! $ finvmaold ) 386 374 387 if (ok_iso_verif) then 388 write(*,*) 'leapfrog 724' 389 call check_isotopes_seq(q,ip1jmp1,'leapfrog 762') 390 endif !if (ok_iso_verif) then 375 call check_isotopes_seq(q,ip1jmp1,'leapfrog 762') 391 376 392 377 c .P.Le Van (26/04/94 ajout de finvpold dans l'appel d'integrd) … … 552 537 CALL massdair(p,masse) 553 538 554 if (ok_iso_verif) then 555 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1196') 556 endif !if (ok_iso_verif) then 539 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1196') 557 540 558 541 c----------------------------------------------------------------------- … … 639 622 c preparation du pas d'integration suivant ...... 640 623 641 if (ok_iso_verif) then 642 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1509') 643 endif !if (ok_iso_verif) then 624 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1509') 644 625 645 626 IF ( .NOT.purmats ) THEN … … 703 684 ENDIF ! of IF((MOD(itau,iperiod).EQ.0).OR.(itau.EQ.itaufin)) 704 685 705 if (ok_iso_verif) then 706 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1584') 707 endif !if (ok_iso_verif) then 686 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1584') 708 687 709 688 c----------------------------------------------------------------------- … … 785 764 ELSE ! of IF (.not.purmats) 786 765 787 if (ok_iso_verif) then 788 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1664') 789 endif !if (ok_iso_verif) then 766 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1664') 790 767 791 768 c ........................................................ … … 812 789 ELSE ! of IF(forward) i.e. backward step 813 790 814 if (ok_iso_verif) then 815 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1698') 816 endif !if (ok_iso_verif) then 791 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1698') 817 792 818 793 IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin) THEN
Note: See TracChangeset
for help on using the changeset viewer.