Changeset 894
- Timestamp:
- 2008-04-21T14:36:04+02:00 (16 years ago)
- Location:
- trunk/NEMO/LIM_SRC_3
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_3/ice.F90
r888 r894 455 455 !!-------------------------------------------------------------------------- 456 456 ! REMOVE 457 INTEGER (wp), PUBLIC, DIMENSION(jpl) :: & !:457 INTEGER, PUBLIC, DIMENSION(jpl) :: & !: 458 458 ice_types !: Vector making the connection between types and categories 459 459 460 INTEGER (wp), PUBLIC, DIMENSION(jpm,2) :: & !:460 INTEGER, PUBLIC, DIMENSION(jpm,2) :: & !: 461 461 ice_cat_bounds !: Matrix containing the integer upper and 462 462 !: lower boundaries of ice thickness categories 463 463 464 464 ! REMOVE 465 INTEGER (wp), PUBLIC, DIMENSION(jpm) :: & !:465 INTEGER, PUBLIC, DIMENSION(jpm) :: & !: 466 466 ice_ncat_types !: Vector containing the number of thickness categories in each ice type 467 467 -
trunk/NEMO/LIM_SRC_3/limdia.F90
r888 r894 421 421 ! oututs on file ice_evolu 422 422 !MV IF( MOD( numit , ninfo ) == 0 ) THEN 423 WRITE( 90,fmtw) ( titvar(jv), vinfom(jv)/naveg, jv = 1, nvinfo )423 WRITE(numevo_ice,fmtw) ( titvar(jv), vinfom(jv)/naveg, jv = 1, nvinfo ) 424 424 naveg = 0 425 425 DO jv = 1, nvinfo … … 453 453 454 454 CHARACTER(len=jpchinf) :: titinf 455 CHARACTER(len=50) :: clname 455 456 !!------------------------------------------------------------------- 456 457 … … 576 577 577 578 ! opening "ice_evolu" file 579 clname = 'ice.evolu' 578 580 irecl = ( jpchinf + 1 ) * nvinfo 579 OPEN( 90, file='ice.evolu', status='unknown', RECL = irecl)580 OPEN( 90, file='ice.evolu', status='unknown')581 CALL ctlopn( numevo_ice, clname, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', 582 & irecl, numout, lwp, 1 ) 581 583 582 584 !- ecriture de 2 lignes d''entete : 583 WRITE( 90,1000) fmtr, fmtw, fmtitr, nvinfo, ntot, 0, nfrinf585 WRITE(numevo_ice,1000) fmtr, fmtw, fmtitr, nvinfo, ntot, 0, nfrinf 584 586 zxx0 = 0.001 * REAL(ninfo) 585 587 zxx1 = 0.001 * REAL(ndeb) 586 WRITE( 90,1111) REAL(jpchinf), 0., zxx1, zxx0, 0., 0., 0588 WRITE(numevo_ice,1111) REAL(jpchinf), 0., zxx1, zxx0, 0., 0., 0 587 589 588 590 !- ecriture de 2 lignes de titre : 589 WRITE( 90,'(A,I8,A,I8,A,I5)') &591 WRITE(numevo_ice,'(A,I8,A,I8,A,I5)') & 590 592 'Evolution chronologique - Experience '//cexper & 591 593 //' de', ndeb, ' a', nferme, ' pas', ninfo 592 WRITE( 90,fmtitr) ( titvar(jv), jv = 1, nvinfo )594 WRITE(numevo_ice,fmtitr) ( titvar(jv), jv = 1, nvinfo ) 593 595 594 596 -
trunk/NEMO/LIM_SRC_3/limitd_me.F90
r888 r894 625 625 DO ji = 1, jpi 626 626 627 IF & 628 ( ( a_i(ji,jj,jl) .GT. epsi11 ) & 627 IF( ( a_i(ji,jj,jl) .GT. epsi11 ) & 629 628 .AND. ( athorn(ji,jj,jl) .GT. 0.0 ) ) THEN 630 629 hi = v_i(ji,jj,jl) / a_i(ji,jj,jl) -
trunk/NEMO/LIM_SRC_3/limupdate.F90
r888 r894 942 942 !------------ 943 943 944 !RB had to split the loop for mpp reproducibility, why ??? 945 DO jj = 1, jpj 946 DO ji = 1, jpim1 944 DO jj = 2, jpjm1 945 DO ji = fs_2, fs_jpim1 947 946 IF ( at_i(ji,jj) .EQ. 0.0 ) THEN ! what to do if there is no ice 948 ! mask u949 947 IF ( at_i(ji+1,jj) .EQ. 0.0 ) u_ice(ji,jj) = 0.0 ! right side 950 ENDIF951 END DO952 END DO953 DO jj = 1, jpj954 DO ji = 2, jpi955 IF ( at_i(ji,jj) .EQ. 0.0 ) THEN ! what to do if there is no ice956 ! mask u957 948 IF ( at_i(ji-1,jj) .EQ. 0.0 ) u_ice(ji-1,jj) = 0.0 ! left side 958 ENDIF959 END DO960 END DO961 DO jj = 1, jpj - 1962 DO ji = 1, jpi963 IF ( at_i(ji,jj) .EQ. 0.0 ) THEN ! what to do if there is no ice964 ! mask u965 949 IF ( at_i(ji,jj+1) .EQ. 0.0 ) v_ice(ji,jj) = 0.0 ! upper side 966 ENDIF967 END DO968 END DO969 DO jj = 2, jpj970 DO ji = 1, jpi971 IF ( at_i(ji,jj) .EQ. 0.0 ) THEN ! what to do if there is no ice972 ! mask u973 950 IF ( at_i(ji,jj-1) .EQ. 0.0 ) v_ice(ji-1,jj) = 0.0 ! bottom side 974 951 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.