Changeset 4143 for LMDZ6/trunk/libf/dyn3d/leapfrog.F
- Timestamp:
- May 9, 2022, 12:35:40 PM (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d/leapfrog.F
r4120 r4143 11 11 use IOIPSL 12 12 #endif 13 USE infotrac, ONLY: nqtot, ok_iso_verif13 USE infotrac, ONLY: nqtot, isoCheck 14 14 USE guide_mod, ONLY : guide_main 15 15 USE write_field, ONLY: writefield … … 26 26 USE temps_mod, ONLY: jD_ref,jH_ref,itaufin,day_ini,day_ref, 27 27 & start_time,dt 28 USE strings_mod, ONLY: msg 28 29 29 30 IMPLICIT NONE … … 237 238 jH_cur = jH_cur - int(jH_cur) 238 239 239 if (ok_iso_verif) then 240 call check_isotopes_seq(q,ip1jmp1,'leapfrog 321') 241 endif !if (ok_iso_verif) then 240 call check_isotopes_seq(q,ip1jmp1,'leapfrog 321') 242 241 243 242 #ifdef CPP_IOIPSL … … 271 270 ! CALL filtreg ( finvmaold ,jjp1, llm, -2,2, .TRUE., 1 ) 272 271 273 if (ok_iso_verif) then 274 call check_isotopes_seq(q,ip1jmp1,'leapfrog 400') 275 endif !if (ok_iso_verif) then 272 call check_isotopes_seq(q,ip1jmp1,'leapfrog 400') 276 273 277 274 2 CONTINUE ! Matsuno backward or leapfrog step begins here … … 324 321 325 322 326 if (ok_iso_verif) then 327 call check_isotopes_seq(q,ip1jmp1,'leapfrog 589') 328 endif !if (ok_iso_verif) then 323 call check_isotopes_seq(q,ip1jmp1,'leapfrog 589') 329 324 330 325 c----------------------------------------------------------------------- … … 345 340 c ------------------------------------------------------------- 346 341 347 if (ok_iso_verif) then 348 call check_isotopes_seq(q,ip1jmp1, 342 call check_isotopes_seq(q,ip1jmp1, 349 343 & 'leapfrog 686: avant caladvtrac') 350 endif !if (ok_iso_verif) then351 344 352 345 IF( forward. OR . leapf ) THEN … … 376 369 c ---------------------------------- 377 370 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 371 CALL msg('720', modname, isoCheck) 372 call check_isotopes_seq(q,ip1jmp1,'leapfrog 756') 382 373 383 374 CALL integrd ( nqtot,vcovm1,ucovm1,tetam1,psm1,massem1 , … … 385 376 ! $ finvmaold ) 386 377 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 378 CALL msg('724', modname, isoCheck) 379 call check_isotopes_seq(q,ip1jmp1,'leapfrog 762') 391 380 392 381 c .P.Le Van (26/04/94 ajout de finvpold dans l'appel d'integrd) … … 552 541 CALL massdair(p,masse) 553 542 554 if (ok_iso_verif) then 555 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1196') 556 endif !if (ok_iso_verif) then 543 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1196') 557 544 558 545 c----------------------------------------------------------------------- … … 639 626 c preparation du pas d'integration suivant ...... 640 627 641 if (ok_iso_verif) then 642 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1509') 643 endif !if (ok_iso_verif) then 628 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1509') 644 629 645 630 IF ( .NOT.purmats ) THEN … … 703 688 ENDIF ! of IF((MOD(itau,iperiod).EQ.0).OR.(itau.EQ.itaufin)) 704 689 705 if (ok_iso_verif) then 706 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1584') 707 endif !if (ok_iso_verif) then 690 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1584') 708 691 709 692 c----------------------------------------------------------------------- … … 790 773 ELSE ! of IF (.not.purmats) 791 774 792 if (ok_iso_verif) then 793 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1664') 794 endif !if (ok_iso_verif) then 775 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1664') 795 776 796 777 c ........................................................ … … 817 798 ELSE ! of IF(forward) i.e. backward step 818 799 819 if (ok_iso_verif) then 820 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1698') 821 endif !if (ok_iso_verif) then 800 call check_isotopes_seq(q,ip1jmp1,'leapfrog 1698') 822 801 823 802 IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin) THEN
Note: See TracChangeset
for help on using the changeset viewer.