Changeset 2844 for branches/2011/dev_r2802_MERCATOR9_floats/NEMOGCM/NEMO/OPA_SRC/FLO/flodom.F90

Timestamp:

2011-09-20T15:40:17+02:00 (13 years ago)

Author:

cbricaud

Message:

minor corrections

File:

• branches/2011/dev_r2802_MERCATOR9_floats/NEMOGCM/NEMO/OPA_SRC/FLO/flodom.F90 (modified) (2 diffs)

branches/2011/dev_r2802_MERCATOR9_floats/NEMOGCM/NEMO/OPA_SRC/FLO/flodom.F90

@@ -319 +319 @@
 
          jfl1=0
+         ierr=0
          CALL ctl_opn( inum, clname2, 'OLD', 'FORMATTED', 'SEQUENTIAL',  1, numout, .TRUE., 1 )
-         DO WHILE (ierr .GE. 0)
+         DO WHILE (ierr .EQ. 0)
             jfl1=jfl1+1
-            READ (inum,*, iostat=ierr)
+            READ(inum,*, iostat=ierr)
          END DO
          CLOSE(inum)
          IF( (jfl1-1) .NE. ifl )THEN 
-            WRITE(cltmp,'(A20,A20,A3,i4.4,A10,i4.4)')"the number of floats in",TRIM(clname2), &
+            WRITE(cltmp,'(A25,A20,A3,i4.4,A10,i4.4)')"the number of floats in ",TRIM(clname2), &
                                                      " = ",jfl1," is not equal to jfl= ",ifl
             CALL ctl_stop('STOP',TRIM(cltmp) )
@@ -437 +438 @@
       IF( ABS(dlx) > 1.0 ) dlx = 1.0
       !
-      dld = ATAN(DSQRT( ( 1-dlx )/( 1+dlx ) )) * 222.24 / dls
+      dld = ATAN(DSQRT( 1.d0 * ( 1.-dlx )/( 1.+dlx ) )) * 222.24 / dls
       dstnce = dld * 1000.
       !