Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat
- Timestamp:
- Aug 2, 2024, 2:12:03 PM (7 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.f90
r5136 r5158 364 364 ! -------------------------------------------------------------- 365 365 DO l = 1, llm 366 doi = 1, iim367 doj = 1, jjm366 DO i = 1, iim 367 DO j = 1, jjm 368 368 zrot(i, j, l) = (pvcov(i + 1, j, l) - pvcov(i, j, l) & 369 369 + pucov(i, j + 1, l) - pucov(i, j, l)) & … … 471 471 IF (CPPKEY_PHYS) THEN 472 472 473 doisplit = 1, nsplit_phys473 DO isplit = 1, nsplit_phys 474 474 475 475 jH_cur_split = jH_cur + (isplit - 1) * dtvr / (daysec * nsplit_phys) -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_calfis_loc.F90
r5136 r5158 285 285 !CDIR ON_ADB(index_i) 286 286 !CDIR ON_ADB(index_j) 287 doig0=1,klon287 DO ig0=1,klon 288 288 i=index_i(ig0) 289 289 j=index_j(ig0) … … 309 309 !CDIR ON_ADB(index_i) 310 310 !CDIR ON_ADB(index_j) 311 doig0=1,klon311 DO ig0=1,klon 312 312 i=index_i(ig0) 313 313 j=index_j(ig0) … … 319 319 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 320 320 DO l=1,llm 321 doig0=1,klon321 DO ig0=1,klon 322 322 i=index_i(ig0) 323 323 j=index_j(ig0) … … 334 334 !CDIR ON_ADB(index_i) 335 335 !CDIR ON_ADB(index_j) 336 doig0=1,klon336 DO ig0=1,klon 337 337 i=index_i(ig0) 338 338 j=index_j(ig0) … … 356 356 !CDIR ON_ADB(index_i) 357 357 !CDIR ON_ADB(index_j) 358 doig0=1,klon358 DO ig0=1,klon 359 359 i=index_i(ig0) 360 360 j=index_j(ig0) … … 373 373 !CDIR ON_ADB(index_i) 374 374 !CDIR ON_ADB(index_j) 375 doig0=1,klon375 DO ig0=1,klon 376 376 i=index_i(ig0) 377 377 j=index_j(ig0) … … 386 386 !CDIR ON_ADB(index_i) 387 387 !CDIR ON_ADB(index_j) 388 doig0=1,klon388 DO ig0=1,klon 389 389 i=index_i(ig0) 390 390 j=index_j(ig0) … … 420 420 !CDIR ON_ADB(index_j) 421 421 !CDIR SPARSE 422 doig0=kstart,kend422 DO ig0=kstart,kend 423 423 i=index_i(ig0) 424 424 j=index_j(ig0) … … 446 446 447 447 DO l=1,llm 448 doi=1,iim449 doj=jjb,jje448 DO i=1,iim 449 DO j=jjb,jje 450 450 zrot(i,j,l) = (pvcov(i+1,j,l) - pvcov(i,j,l) & 451 451 + pucov(i,j+1,l) - pucov(i,j,l)) & … … 545 545 !CDIR ON_ADB(index_i) 546 546 !CDIR ON_ADB(index_j) 547 doig0=1,klon547 DO ig0=1,klon 548 548 i=index_i(ig0) 549 549 j=index_j(ig0) … … 592 592 offset=klon_omp_begin-1 593 593 594 dol=1,llm+1595 doi=1,klon594 DO l=1,llm+1 595 DO i=1,klon 596 596 zplev_omp(i,l)=zplev(offset+i,l) 597 597 enddo 598 598 enddo 599 599 600 dol=1,llm601 doi=1,klon600 DO l=1,llm 601 DO i=1,klon 602 602 zplay_omp(i,l)=zplay(offset+i,l) 603 603 enddo 604 604 enddo 605 605 606 dol=1,llm607 doi=1,klon606 DO l=1,llm 607 DO i=1,klon 608 608 zpk_omp(i,l)=zpk(offset+i,l) 609 609 enddo 610 610 enddo 611 611 612 dol=1,llm613 doi=1,klon612 DO l=1,llm 613 DO i=1,klon 614 614 zphi_omp(i,l)=zphi(offset+i,l) 615 615 enddo 616 616 enddo 617 617 618 doi=1,klon618 DO i=1,klon 619 619 zphis_omp(i)=zphis(offset+i) 620 620 enddo 621 621 622 622 623 dol=1,llm623 DO l=1,llm 624 624 presnivs_omp(l)=presnivs(l) 625 625 enddo 626 626 627 dol=1,llm628 doi=1,klon627 DO l=1,llm 628 DO i=1,klon 629 629 zufi_omp(i,l)=zufi(offset+i,l) 630 630 enddo 631 631 enddo 632 632 633 dol=1,llm634 doi=1,klon633 DO l=1,llm 634 DO i=1,klon 635 635 zvfi_omp(i,l)=zvfi(offset+i,l) 636 636 enddo 637 637 enddo 638 638 639 dol=1,llm640 doi=1,klon639 DO l=1,llm 640 DO i=1,klon 641 641 zrfi_omp(i,l)=zrfi(offset+i,l) 642 642 enddo 643 643 enddo 644 644 645 dol=1,llm646 doi=1,klon645 DO l=1,llm 646 DO i=1,klon 647 647 ztfi_omp(i,l)=ztfi(offset+i,l) 648 648 enddo 649 649 enddo 650 650 651 doiq=1,nqtot652 dol=1,llm653 doi=1,klon651 DO iq=1,nqtot 652 DO l=1,llm 653 DO i=1,klon 654 654 zqfi_omp(i,l,iq)=zqfi(offset+i,l,iq) 655 655 enddo … … 657 657 enddo 658 658 659 dol=1,llm660 doi=1,klon659 DO l=1,llm 660 DO i=1,klon 661 661 zdufi_omp(i,l)=zdufi(offset+i,l) 662 662 enddo 663 663 enddo 664 664 665 dol=1,llm666 doi=1,klon665 DO l=1,llm 666 DO i=1,klon 667 667 zdvfi_omp(i,l)=zdvfi(offset+i,l) 668 668 enddo 669 669 enddo 670 670 671 dol=1,llm672 doi=1,klon671 DO l=1,llm 672 DO i=1,klon 673 673 zdtfi_omp(i,l)=zdtfi(offset+i,l) 674 674 enddo 675 675 enddo 676 676 677 doiq=1,nqtot678 dol=1,llm679 doi=1,klon677 DO iq=1,nqtot 678 DO l=1,llm 679 DO i=1,klon 680 680 zdqfi_omp(i,l,iq)=zdqfi(offset+i,l,iq) 681 681 enddo … … 683 683 enddo 684 684 685 doi=1,klon685 DO i=1,klon 686 686 zdpsrf_omp(i)=zdpsrf(offset+i) 687 687 enddo 688 688 689 dol=1,llm690 doi=1,klon689 DO l=1,llm 690 DO i=1,klon 691 691 flxwfi_omp(i,l)=flxwfi(offset+i,l) 692 692 enddo … … 705 705 zdqfic_omp(:,:,:)=0. 706 706 707 doisplit=1,nsplit_phys707 DO isplit=1,nsplit_phys 708 708 709 709 jH_cur_split=jH_cur+(isplit-1) * dtvr / (daysec *nsplit_phys) … … 743 743 !$OMP BARRIER 744 744 745 dol=1,llm+1746 doi=1,klon745 DO l=1,llm+1 746 DO i=1,klon 747 747 zplev(offset+i,l)=zplev_omp(i,l) 748 748 enddo 749 749 enddo 750 750 751 dol=1,llm752 doi=1,klon751 DO l=1,llm 752 DO i=1,klon 753 753 zplay(offset+i,l)=zplay_omp(i,l) 754 754 enddo 755 755 enddo 756 756 757 dol=1,llm758 doi=1,klon757 DO l=1,llm 758 DO i=1,klon 759 759 zphi(offset+i,l)=zphi_omp(i,l) 760 760 enddo … … 762 762 763 763 764 doi=1,klon764 DO i=1,klon 765 765 zphis(offset+i)=zphis_omp(i) 766 766 enddo 767 767 768 768 769 dol=1,llm769 DO l=1,llm 770 770 presnivs(l)=presnivs_omp(l) 771 771 enddo 772 772 773 dol=1,llm774 doi=1,klon773 DO l=1,llm 774 DO i=1,klon 775 775 zufi(offset+i,l)=zufi_omp(i,l) 776 776 enddo 777 777 enddo 778 778 779 dol=1,llm780 doi=1,klon779 DO l=1,llm 780 DO i=1,klon 781 781 zvfi(offset+i,l)=zvfi_omp(i,l) 782 782 enddo 783 783 enddo 784 784 785 dol=1,llm786 doi=1,klon785 DO l=1,llm 786 DO i=1,klon 787 787 ztfi(offset+i,l)=ztfi_omp(i,l) 788 788 enddo 789 789 enddo 790 790 791 doiq=1,nqtot792 dol=1,llm793 doi=1,klon791 DO iq=1,nqtot 792 DO l=1,llm 793 DO i=1,klon 794 794 zqfi(offset+i,l,iq)=zqfi_omp(i,l,iq) 795 795 enddo … … 797 797 enddo 798 798 799 dol=1,llm800 doi=1,klon799 DO l=1,llm 800 DO i=1,klon 801 801 zdufi(offset+i,l)=zdufi_omp(i,l) 802 802 enddo 803 803 enddo 804 804 805 dol=1,llm806 doi=1,klon805 DO l=1,llm 806 DO i=1,klon 807 807 zdvfi(offset+i,l)=zdvfi_omp(i,l) 808 808 enddo 809 809 enddo 810 810 811 dol=1,llm812 doi=1,klon811 DO l=1,llm 812 DO i=1,klon 813 813 zdtfi(offset+i,l)=zdtfi_omp(i,l) 814 814 enddo 815 815 enddo 816 816 817 doiq=1,nqtot818 dol=1,llm819 doi=1,klon817 DO iq=1,nqtot 818 DO l=1,llm 819 DO i=1,klon 820 820 zdqfi(offset+i,l,iq)=zdqfi_omp(i,l,iq) 821 821 enddo … … 823 823 enddo 824 824 825 doi=1,klon825 DO i=1,klon 826 826 zdpsrf(offset+i)=zdpsrf_omp(i) 827 827 enddo … … 943 943 !CDIR ON_ADB(index_j) 944 944 !cdir NODEP 945 doig0=kstart,kend945 DO ig0=kstart,kend 946 946 i=index_i(ig0) 947 947 j=index_j(ig0) … … 979 979 !CDIR ON_ADB(index_j) 980 980 !cdir NODEP 981 doig0=kstart,kend981 DO ig0=kstart,kend 982 982 i=index_i(ig0) 983 983 j=index_j(ig0) … … 1079 1079 !CDIR ON_ADB(index_j) 1080 1080 !cdir NODEP 1081 doig0=kstart,kend1081 DO ig0=kstart,kend 1082 1082 i=index_i(ig0) 1083 1083 j=index_j(ig0) … … 1124 1124 !CDIR ON_ADB(index_j) 1125 1125 !cdir NODEP 1126 doig0=kstart,kend1126 DO ig0=kstart,kend 1127 1127 i=index_i(ig0) 1128 1128 j=index_j(ig0) -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90
r5139 r5158 538 538 coslat0 = 0.5 539 539 phiso = phis 540 doj = 2, jmp1 - 1540 DO j = 2, jmp1 - 1 541 541 PRINT*, 'avant if ', cos(rlatu(j)), coslat0 542 542 IF (cos(rlatu(j))<coslat0) THEN … … 544 544 ifiltre = (coslat0 / cos(rlatu(j)) - 1.) / 2. 545 545 wwf = 0. 546 doi = 1, ifiltre546 DO i = 1, ifiltre 547 547 wwf(i) = 1. 548 548 enddo 549 549 wwf(ifiltre + 1) = (coslat0 / cos(rlatu(j)) - 1.) / 2. - ifiltre 550 doi = 1, imp1 - 1550 DO i = 1, imp1 - 1 551 551 IF (masque(i, j)>0.9) THEN 552 552 ssz = phis(i, j) 553 doifi = 1, ifiltre + 1553 DO ifi = 1, ifiltre + 1 554 554 ii = i + ifi 555 555 IF (ii>imp1 - 1) ii = ii - imp1 + 1 -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90
r5137 r5158 102 102 103 103 !!$OMP PARALLEL DEFAULT(SHARED) COPYIN(/temps/) 104 ! 104 ! Copy all threadprivate variables in temps_mod 105 105 !$OMP PARALLEL DEFAULT(SHARED) COPYIN(annee_ref,day_ini,day_ref,start_time) 106 106 -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/test_disvert_m.F90
r5134 r5158 44 44 IF (any(p(:, :llm) <= p_lay .OR. p_lay <= p(:, 2:))) THEN 45 45 ! List details and stop: 46 dol = 1, llm47 doi = 1, ngrid46 DO l = 1, llm 47 DO i = 1, ngrid 48 48 IF (p(i, l) <= p_lay(i, l)) THEN 49 49 print 1000, "ps = ", ps(i) / 100., "hPa, p(level ", l, &
Note: See TracChangeset
for help on using the changeset viewer.