Changeset 69 for trunk/NEMO/OPA_SRC
- Timestamp:
- 2004-04-22T14:17:26+02:00 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/mppini_2.h90
r15 r69 40 40 !!---------------------------------------------------------------------- 41 41 !! Local variables 42 CHARACTER (len=25) :: clexp ! temporary name 42 CHARACTER (len=25) :: clexp, & ! temporary name 43 & clname ! " " 44 LOGICAL :: llbon ! check the existence of bathy files 43 45 INTEGER :: ji, jj, jn, jproc, jarea ! dummy loop indices 44 46 INTEGER :: inum = 11 ! temporary logical unit … … 98 100 99 101 ! open the file 102 clname = 'bathy_level' 100 103 IF(lwp) WRITE(numout,*) ' bathymetry file in levels' 101 CALL ctlopn( inum, 'bathymetry', 'OLD', 'FORMATTED', 'SEQUENTIAL', &102 1 , numout , lwp , 1 )103 104 ! read bathymetry file105 REWIND inum106 READ(inum,9001) clexp, iim, ijm107 READ(inum,'(/)') 108 ifreq=40109 il1=1110 DO jn=1,jpidta/ifreq+1104 INQUIRE( FILE=clname, EXIST=llbon ) 105 IF( llbon ) THEN 106 IF(lwp) WRITE(numout,*) 107 IF(lwp) WRITE(numout,*) ' read level bathymetry in ', clname 108 IF(lwp) WRITE(numout,*) 109 OPEN( UNIT=inum, FILE=clname, FORM='FORMATTED', ACCESS='SEQUENTIAL', STATUS='OLD' ) 110 111 ! read bathymetry file 112 REWIND inum 113 READ(inum,9001) clexp, iim, ijm 111 114 READ(inum,'(/)') 112 il2 = MIN( jpidta, il1+ifreq-1 ) 113 READ(inum,9002) ( ii, ji = il1, il2, 5 ) 114 READ(inum,'(/)') 115 DO jj = jpjdta, 1, -1 116 READ(inum,9003) ij, ( idata(ji,jj), ji = il1, il2 ) 117 END DO 118 il1 = il1 + ifreq 119 END DO 120 CLOSE(inum) 121 122 9001 FORMAT(1x,a15,2i8) 123 9002 FORMAT(3x,13(i3,12x)) 124 9003 FORMAT(i3,41i3) 115 ifreq=40 116 il1=1 117 DO jn=1,jpidta/ifreq+1 118 READ(inum,'(/)') 119 il2 = MIN( jpidta, il1+ifreq-1 ) 120 READ(inum,9002) ( ii, ji = il1, il2, 5 ) 121 READ(inum,'(/)') 122 DO jj = jpjdta, 1, -1 123 READ(inum,9003) ij, ( idata(ji,jj), ji = il1, il2 ) 124 END DO 125 il1 = il1 + ifreq 126 END DO 127 CLOSE(inum) 128 129 9001 FORMAT(1x,a15,2i8) 130 9002 FORMAT(3x,13(i3,12x)) 131 9003 FORMAT(i3,41i3) 132 133 ELSE 134 IF(lwp) WRITE(numout,cform_err) 135 IF(lwp) WRITE(numout,*)' mppini_2 : unable to read the file', clname 136 nstop = nstop + 1 137 ENDIF 125 138 126 139 ! land/sea mask over the global/zoom domain
Note: See TracChangeset
for help on using the changeset viewer.