Changeset 1415 for trunk/NEMO/OPA_SRC/SOL/solisl.F90
- Timestamp:
- 2009-04-29T16:45:11+02:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SOL/solisl.F90
r1152 r1415 558 558 !! * Local declarations 559 559 INTEGER :: ji, jj, jni, jnj, jl ! dummy loop indices 560 INTEGER :: ios ! temporary integers 561 INTEGER :: & 562 inum ! temporary logical unit 560 INTEGER :: ios ! temporary integers 561 INTEGER :: inum ! temporary logical unit 563 562 REAL(wp), DIMENSION(jpi,jpj,2) :: zwx 564 563 REAL(wp), DIMENSION(jpisl*jpisl) :: ztab … … 599 598 ENDIF 600 599 601 CALL restclo(numisl)602 603 600 ELSE 604 601 605 CALL iom_close( inum )606 607 602 ! II. Island matrix computation 608 603 ! ============================= … … 685 680 ! -------------------------------------- 686 681 687 CALL restput( numisl, 'aisl' , jpisl, jpisl, 1, 0, aisl )688 CALL restput( numisl, 'aislm1', jpisl, jpisl, 1, 0, aislm1 )689 CALL restclo( numisl)682 CALL iom_rstput( 0, 0, inum, 'aisl' , aisl ) 683 CALL iom_rstput( 0, 0, inum, 'aislm1', aislm1 ) 684 CALL iom_close ( inum ) 690 685 691 686 ENDIF … … 741 736 INTEGER :: iimlu, ijmlu, inmlu, iju 742 737 INTEGER :: ii, ij, icile, icut, inmax, indic 743 INTEGER :: itime744 738 REAL(wp) :: zepsr, zeplu, zgwgt 745 REAL(wp) :: zep(jpisl), z dept(1), zprec(4)746 REAL(wp) :: zd ate0, zdt739 REAL(wp) :: zep(jpisl), zprec(4) 740 REAL(wp) :: zdt 747 741 REAL(wp) :: t2p1(jpi,1,1) 748 742 INTEGER :: iloc … … 796 790 END DO 797 791 ENDIF 792 CALL iom_close( inum ) 793 ! 798 794 ELSE 799 795 ! islands FILE does not EXIST : icut=999 800 796 icut = 999 801 797 ENDIF 802 803 CALL iom_close( inum )804 798 805 799 ! the read precision is not the required one : icut=888 … … 1071 1065 ! =============================== 1072 1066 1073 CALL ymds2ju( 0, 1, 1, 0.e0, zdate0 )1074 1067 zprec(1) = FLOAT(jpi) 1075 1068 zprec(2) = FLOAT(jpj) 1076 1069 zprec(3) = FLOAT(jpisl) 1077 1070 IF(lwp) WRITE(numout,*) clname 1078 zdept(1) = 0. 1079 itime = 0 1080 CALL restini( 'NONE', jpi, jpj, glamt, gphit, 1, zdept, & 1081 & clname, itime, zdate0, rdt, numisl, domain_id=nidom ) 1071 CALL iom_open ( clname, numisl, ldwrt = .TRUE., kiolib = jprstlib ) 1082 1072 IF( icile == 0 .AND. icut /= 0 ) THEN 1083 1073 IF(lwp) THEN … … 1087 1077 ENDIF 1088 1078 zprec(4) = epsisl 1089 CALL restput(numisl,'PRECISION',1,1,4,0,zprec)1079 CALL iom_rstput( 0, 0, numisl, 'PRECISION' , zprec ) 1090 1080 DO jni = 1, jpisl 1091 1081 IF(jni < 10) THEN … … 1096 1086 WRITE(clisl,'("island",I3)') jni 1097 1087 ENDIF 1098 CALL restput( numisl, clisl, jpi, jpj, 1, 0, bsfisl(:,:,jni) ) 1099 END DO 1088 CALL iom_rstput( 0, 0, numisl, clisl , bsfisl(:,:,jni) ) 1089 END DO 1090 CALL iom_close( numisl ) 1100 1091 ENDIF 1101 1092 … … 1120 1111 1121 1112 zprec(4) = zepsr 1122 CALL restput( numisl, 'PRECISION', 1, 1, 1, 0, zprec )1113 CALL iom_rstput( 0, 0, numisl, 'PRECISION' , zprec ) 1123 1114 DO jni = 1, jpisl 1124 1115 IF( jni < 10 ) THEN … … 1129 1120 WRITE(clisl,'("island",I3)') jni 1130 1121 ENDIF 1131 CALL restput( numisl, clisl, jpi, jpj, 1, 0, bsfisl(:,:,jni) )1132 END DO 1133 CALL restclo(numisl)1122 CALL iom_rstput( 0, 0, numisl, clisl , bsfisl(:,:,jni) ) 1123 END DO 1124 CALL iom_close( numisl ) 1134 1125 CALL ctl_stop( ' ' ) 1135 1126 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.