Changeset 2844
- Timestamp:
- 2011-09-20T15:40:17+02:00 (13 years ago)
- Location:
- branches/2011/dev_r2802_MERCATOR9_floats/NEMOGCM/NEMO/OPA_SRC/FLO
- Files:
-
- 3 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 ! -
branches/2011/dev_r2802_MERCATOR9_floats/NEMOGCM/NEMO/OPA_SRC/FLO/florst.F90
r2843 r2844 17 17 !! * Modules used 18 18 USE flo_oce ! ocean drifting floats 19 !USE oce ! ocean dynamics and tracers20 19 USE dom_oce ! ocean space and time domain 21 20 USE lib_mpp ! distribued memory computing library 22 21 USE in_out_manager ! I/O manager 23 !USE phycst ! physic constants24 !USE dianam ! build name of file (routine)25 !USE ioipsl26 22 27 23 IMPLICIT NONE … … 29 25 30 26 PUBLIC flo_rst ! routine called by floats.F90 31 32 !INTEGER :: jfl ! number of floats33 !INTEGER :: numfl ! logical unit for floats netcdf output34 27 35 28 !! * Substitutions … … 45 38 SUBROUTINE flo_rst( kt ) 46 39 !!--------------------------------------------------------------------- 47 !! *** ROUTINE flo_ wri***40 !! *** ROUTINE flo_rst *** 48 41 !! 49 42 !! ** Purpose : … … 85 78 END DO 86 79 87 CALL ctl_opn( inum, clname, 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 88 REWIND inum 80 inum=0 81 IF( lwp )THEN 82 CALL ctl_opn( inum, clname, 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 83 REWIND inum 84 ENDIF 89 85 ! 90 86 DO jpn = 1, jpnij -
branches/2011/dev_r2802_MERCATOR9_floats/NEMOGCM/NEMO/OPA_SRC/FLO/flowri.F90
r2841 r2844 52 52 CONTAINS 53 53 54 INTEGER FUNCTION flo_wri_alloc 54 INTEGER FUNCTION flo_wri_alloc() 55 55 !!------------------------------------------------------------------- 56 56 !! *** FUNCTION flo_wri_alloc *** … … 220 220 221 221 #if defined key_iomput 222 IF(lwp)WRITE(numout,*)"zlon ",zlon ; call FLUSH(numout)223 222 CALL iom_put( "traj_lon" , zlon ) 224 223 CALL iom_put( "traj_lat" , zlat )
Note: See TracChangeset
for help on using the changeset viewer.