- Timestamp:
- 2011-09-20T15:40:17+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2802_MERCATOR9_floats/NEMOGCM/NEMO/OPA_SRC/FLO/flodom.F90
r2843 r2844 319 319 320 320 jfl1=0 321 ierr=0 321 322 CALL ctl_opn( inum, clname2, 'OLD', 'FORMATTED', 'SEQUENTIAL', 1, numout, .TRUE., 1 ) 322 DO WHILE (ierr . GE. 0)323 DO WHILE (ierr .EQ. 0) 323 324 jfl1=jfl1+1 324 READ 325 READ(inum,*, iostat=ierr) 325 326 END DO 326 327 CLOSE(inum) 327 328 IF( (jfl1-1) .NE. ifl )THEN 328 WRITE(cltmp,'(A2 0,A20,A3,i4.4,A10,i4.4)')"the number of floats in",TRIM(clname2), &329 WRITE(cltmp,'(A25,A20,A3,i4.4,A10,i4.4)')"the number of floats in ",TRIM(clname2), & 329 330 " = ",jfl1," is not equal to jfl= ",ifl 330 331 CALL ctl_stop('STOP',TRIM(cltmp) ) … … 437 438 IF( ABS(dlx) > 1.0 ) dlx = 1.0 438 439 ! 439 dld = ATAN(DSQRT( ( 1-dlx )/( 1+dlx ) )) * 222.24 / dls440 dld = ATAN(DSQRT( 1.d0 * ( 1.-dlx )/( 1.+dlx ) )) * 222.24 / dls 440 441 dstnce = dld * 1000. 441 442 !
Note: See TracChangeset
for help on using the changeset viewer.