Changeset 5270 for LMDZ6/trunk/libf/phylmdiso/phyaqua_mod.F90
- Timestamp:
- Oct 24, 2024, 1:55:38 PM (2 days ago)
- File:
-
- 1 edited
-
LMDZ6/trunk/libf/phylmdiso/phyaqua_mod.F90 (modified) (15 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmdiso/phyaqua_mod.F90
r5249 r5270 539 539 IMPLICIT NONE 540 540 541 include "netcdf.inc"542 543 541 INTEGER, INTENT (IN) :: klon 544 542 REAL, INTENT (IN) :: phy_nat(klon, 360) … … 593 591 USE mod_phys_lmdz_transfert_para, ONLY: gather 594 592 USE phys_cal_mod, ONLY: year_len 595 USE netcdf, ONLY: nf90_put_var, nf90_def_var 593 USE netcdf, ONLY: nf90_clobber, nf90_close, nf90_noerr, nf90_strerror, nf90_put_att, nf90_def_var, & 594 nf90_def_dim, nf90_create, nf90_put_var, nf90_unlimited, nf90_global, nf90_64bit_offset, & 595 nf90_enddef 596 596 USE lmdz_cppkeys_wrapper, ONLY: nf90_format 597 597 IMPLICIT NONE 598 include "netcdf.inc"599 598 600 599 INTEGER, INTENT (IN) :: klon … … 627 626 PRINT *, 'writelim: Ecriture du fichier limit' 628 627 629 ierr = nf_create('limit.nc', IOR(NF_CLOBBER,NF_64BIT_OFFSET), nid) 630 631 ierr = nf_put_att_text(nid, nf_global, 'title', 30, & 632 'Fichier conditions aux limites') 633 ! ! ierr = NF_DEF_DIM (nid, "points_physiques", klon, ndim) 634 ierr = nf_def_dim(nid, 'points_physiques', klon_glo, ndim) 635 ierr = nf_def_dim(nid, 'time', nf_unlimited, ntim) 628 ierr = nf90_create('limit.nc', IOR(nf90_clobber,nf90_64bit_offset), nid) 629 630 ierr = nf90_put_att(nid, nf90_global, 'title', 'Fichier conditions aux limites') 631 ! ! ierr = nf90_def_dim (nid, "points_physiques", klon, ndim) 632 ierr = nf90_def_dim(nid, 'points_physiques', klon_glo, ndim) 633 ierr = nf90_def_dim(nid, 'time', nf90_unlimited, ntim) 636 634 637 635 dims(1) = ndim … … 639 637 640 638 ierr = nf90_def_var(nid, 'TEMPS', nf90_format, ntim, id_tim) 641 ierr = nf _put_att_text(nid, id_tim, 'title', 17, 'Jour dans l annee')639 ierr = nf90_put_att(nid, id_tim, 'title', 'Jour dans l annee') 642 640 643 641 ierr = nf90_def_var(nid, 'NAT', nf90_format, dims, id_nat) 644 ierr = nf_put_att_text(nid, id_nat, 'title', 23, & 645 'Nature du sol (0,1,2,3)') 642 ierr = nf90_put_att(nid, id_nat, 'title', 'Nature du sol (0,1,2,3)') 646 643 647 644 ierr = nf90_def_var(nid, 'SST', nf90_format, dims, id_sst) 648 ierr = nf_put_att_text(nid, id_sst, 'title', 35, & 649 'Temperature superficielle de la mer') 645 ierr = nf90_put_att(nid, id_sst, 'title', 'Temperature superficielle de la mer') 650 646 651 647 ierr = nf90_def_var(nid, 'BILS', nf90_format, dims, id_bils) 652 ierr = nf_put_att_text(nid, id_bils, 'title', 32, & 653 'Reference flux de chaleur au sol') 648 ierr = nf90_put_att(nid, id_bils, 'title', 'Reference flux de chaleur au sol') 654 649 655 650 ierr = nf90_def_var(nid, 'ALB', nf90_format, dims, id_alb) 656 ierr = nf _put_att_text(nid, id_alb, 'title', 19, 'Albedo a la surface')651 ierr = nf90_put_att(nid, id_alb, 'title', 'Albedo a la surface') 657 652 658 653 ierr = nf90_def_var(nid, 'RUG', nf90_format, dims, id_rug) 659 ierr = nf _put_att_text(nid, id_rug, 'title', 8, 'Rugosite')654 ierr = nf90_put_att(nid, id_rug, 'title', 'Rugosite') 660 655 661 656 ierr = nf90_def_var(nid, 'FTER', nf90_format, dims, id_fter) 662 ierr = nf _put_att_text(nid, id_fter, 'title',10,'Frac. Land')657 ierr = nf90_put_att(nid, id_fter, 'title','Frac. Land') 663 658 ierr = nf90_def_var(nid, 'FOCE', nf90_format, dims, id_foce) 664 ierr = nf _put_att_text(nid, id_foce, 'title',11,'Frac. Ocean')659 ierr = nf90_put_att(nid, id_foce, 'title','Frac. Ocean') 665 660 ierr = nf90_def_var(nid, 'FSIC', nf90_format, dims, id_fsic) 666 ierr = nf _put_att_text(nid, id_fsic, 'title',13,'Frac. Sea Ice')661 ierr = nf90_put_att(nid, id_fsic, 'title','Frac. Sea Ice') 667 662 ierr = nf90_def_var(nid, 'FLIC', nf90_format, dims, id_flic) 668 ierr = nf _put_att_text(nid, id_flic, 'title',14,'Frac. Land Ice')669 670 ierr = nf _enddef(nid)671 IF (ierr/=nf _noerr) THEN663 ierr = nf90_put_att(nid, id_flic, 'title','Frac. Land Ice') 664 665 ierr = nf90_enddef(nid) 666 IF (ierr/=nf90_noerr) THEN 672 667 WRITE (*, *) 'writelim error: failed to end define mode' 673 WRITE (*, *) nf _strerror(ierr)668 WRITE (*, *) nf90_strerror(ierr) 674 669 END IF 675 670 … … 678 673 DO k = 1, year_len 679 674 ierr = nf90_put_var(nid, id_tim, k, [k]) 680 IF (ierr/=nf _noerr) THEN675 IF (ierr/=nf90_noerr) THEN 681 676 WRITE (*, *) 'writelim error with temps(k),k=', k 682 WRITE (*, *) nf _strerror(ierr)677 WRITE (*, *) nf90_strerror(ierr) 683 678 END IF 684 679 END DO … … 691 686 IF (is_master) THEN 692 687 ierr = nf90_put_var(nid, id_nat, phy_glo) 693 IF (ierr/=nf _noerr) THEN688 IF (ierr/=nf90_noerr) THEN 694 689 WRITE (*, *) 'writelim error with phy_nat' 695 WRITE (*, *) nf _strerror(ierr)690 WRITE (*, *) nf90_strerror(ierr) 696 691 END IF 697 692 END IF … … 700 695 IF (is_master) THEN 701 696 ierr = nf90_put_var(nid, id_sst, phy_glo) 702 IF (ierr/=nf _noerr) THEN697 IF (ierr/=nf90_noerr) THEN 703 698 WRITE (*, *) 'writelim error with phy_sst' 704 WRITE (*, *) nf _strerror(ierr)699 WRITE (*, *) nf90_strerror(ierr) 705 700 END IF 706 701 END IF … … 709 704 IF (is_master) THEN 710 705 ierr = nf90_put_var(nid, id_bils, phy_glo) 711 IF (ierr/=nf _noerr) THEN706 IF (ierr/=nf90_noerr) THEN 712 707 WRITE (*, *) 'writelim error with phy_bil' 713 WRITE (*, *) nf _strerror(ierr)708 WRITE (*, *) nf90_strerror(ierr) 714 709 END IF 715 710 END IF … … 718 713 IF (is_master) THEN 719 714 ierr = nf90_put_var(nid, id_alb, phy_glo) 720 IF (ierr/=nf _noerr) THEN715 IF (ierr/=nf90_noerr) THEN 721 716 WRITE (*, *) 'writelim error with phy_alb' 722 WRITE (*, *) nf _strerror(ierr)717 WRITE (*, *) nf90_strerror(ierr) 723 718 END IF 724 719 END IF … … 727 722 IF (is_master) THEN 728 723 ierr = nf90_put_var(nid, id_rug, phy_glo) 729 IF (ierr/=nf _noerr) THEN724 IF (ierr/=nf90_noerr) THEN 730 725 WRITE (*, *) 'writelim error with phy_rug' 731 WRITE (*, *) nf _strerror(ierr)726 WRITE (*, *) nf90_strerror(ierr) 732 727 END IF 733 728 END IF … … 736 731 IF (is_master) THEN 737 732 ierr = nf90_put_var(nid, id_fter, phy_glo) 738 IF (ierr/=nf _noerr) THEN733 IF (ierr/=nf90_noerr) THEN 739 734 WRITE (*, *) 'writelim error with phy_fter' 740 WRITE (*, *) nf _strerror(ierr)735 WRITE (*, *) nf90_strerror(ierr) 741 736 END IF 742 737 END IF … … 745 740 IF (is_master) THEN 746 741 ierr = nf90_put_var(nid, id_foce, phy_glo) 747 IF (ierr/=nf _noerr) THEN742 IF (ierr/=nf90_noerr) THEN 748 743 WRITE (*, *) 'writelim error with phy_foce' 749 WRITE (*, *) nf _strerror(ierr)744 WRITE (*, *) nf90_strerror(ierr) 750 745 END IF 751 746 END IF … … 754 749 IF (is_master) THEN 755 750 ierr = nf90_put_var(nid, id_fsic, phy_glo) 756 IF (ierr/=nf _noerr) THEN751 IF (ierr/=nf90_noerr) THEN 757 752 WRITE (*, *) 'writelim error with phy_fsic' 758 WRITE (*, *) nf _strerror(ierr)753 WRITE (*, *) nf90_strerror(ierr) 759 754 END IF 760 755 END IF … … 763 758 IF (is_master) THEN 764 759 ierr = nf90_put_var(nid, id_flic, phy_glo) 765 IF (ierr/=nf _noerr) THEN760 IF (ierr/=nf90_noerr) THEN 766 761 WRITE (*, *) 'writelim error with phy_flic' 767 WRITE (*, *) nf _strerror(ierr)762 WRITE (*, *) nf90_strerror(ierr) 768 763 END IF 769 764 END IF … … 771 766 ! close file: 772 767 IF (is_master) THEN 773 ierr = nf _close(nid)768 ierr = nf90_close(nid) 774 769 END IF 775 770
Note: See TracChangeset
for help on using the changeset viewer.
