Changeset 895


Ignore:
Timestamp:
2008-04-21T17:32:02+02:00 (13 years ago)
Author:
rblod
Message:

Fix small bugs in sea-ice (previous commit) and old flyspray: 189, 185, 161

Location:
trunk/NEMO
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/limdia.F90

    r894 r895  
    579579       clname = 'ice.evolu' 
    580580       irecl = ( jpchinf + 1 ) * nvinfo  
    581        CALL ctlopn( numevo_ice, clname, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', 
     581       CALL ctlopn( numevo_ice, clname, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL',    & 
    582582          &         irecl, numout, lwp, 1 ) 
    583583 
  • trunk/NEMO/LIM_SRC_3/limupdate.F90

    r894 r895  
    5252   !! * Accessibility 
    5353   PUBLIC lim_update ! routine called by ice_step 
     54 
     55   !! * Substitutions 
     56#  include "vectopt_loop_substitute.h90" 
    5457 
    5558   !! * Module variables 
  • trunk/NEMO/OPA_SRC/DIA/diawri.F90

    r888 r895  
    169169 
    170170         ! WRITE root name in date.file for use by postpro 
    171          CALL dia_nam( clhstnam, nwrite,' ' ) 
    172          clname = 'date.file' 
    173          CALL ctlopn( inum, clname,  'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', 1, numout, lwp, 1 ) 
    174          WRITE(inum,*) clhstnam 
    175          CLOSE(inum) 
     171         IF( lwp) THEN 
     172            CALL dia_nam( clhstnam, nwrite,' ' ) 
     173            clname = 'date.file' 
     174            CALL ctlopn( inum, clname,  'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', 1, numout, lwp, 1 ) 
     175            WRITE(inum,*) clhstnam 
     176            CLOSE(inum) 
     177         ENDIF 
    176178 
    177179         ! Define the T grid FILE ( nid_T ) 
  • trunk/NEMO/OPA_SRC/LDF/ldfeiv.F90

    r888 r895  
    137137            ! Compute aeiw by multiplying Ro^2 and T^-1 
    138138            aeiw(ji,jj) = zross(ji,jj) * zross(ji,jj) * SQRT( zah(ji,jj) / zhw(ji,jj) ) * tmask(ji,jj,1) 
    139             IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN   ! ORCA R02 
     139         END DO 
     140      END DO 
     141 
     142      IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN   ! ORCA R02 
     143         DO jj = 2, jpjm1 
     144!CDIR NOVERRCHK  
     145            DO ji = fs_2, fs_jpim1   ! vector opt. 
    140146               ! Take the minimum between aeiw and aeiv0 for depth levels 
    141147               ! lower than 20 (21 in w- point) 
    142148               IF( mbathy(ji,jj) <= 21. ) aeiw(ji,jj) = MIN( aeiw(ji,jj), 1000. ) 
    143             ENDIF 
    144          END DO 
    145       END DO 
     149            END DO 
     150         END DO 
     151      ENDIF 
    146152 
    147153      ! Decrease the coefficient in the tropics (20N-20S)  
  • trunk/NEMO/OPA_SRC/SBC/sbcfwb.F90

    r888 r895  
    6060      INTEGER, INTENT( in ) ::   kn_fwb   ! ocean time-step index 
    6161      !! 
    62       INTEGER  ::   inum = 11             ! temporary logical unit 
     62      INTEGER  ::   inum                  ! temporary logical unit 
    6363      INTEGER  ::   ikty, iyear           !  
    6464      LOGICAL  ::   llbon 
    65       CHARACTER (len=32) ::   clname = 'EMPave_old.dat' 
     65      CHARACTER (len=32) ::   clname  
    6666      REAL(wp) ::   z_emp                 ! temporary scalars 
    6767      !!---------------------------------------------------------------------- 
     
    9898         IF( kt == nit000 ) THEN 
    9999            ! Read the corrective factor on precipitations (empold) 
    100             INQUIRE( FILE=clname, EXIST=llbon ) 
    101             IF( llbon ) THEN 
    102                OPEN ( inum, FILE=clname) 
    103                READ ( inum, "(24X,I8,2ES24.16)" ) iyear, a_fwb_b, a_fwb 
    104                CLOSE( inum ) 
    105                empold = a_fwb                  ! current year freshwater budget correction 
    106                !                               ! estimate from the previous year budget 
    107                IF(lwp)WRITE(numout,*) 
    108                IF(lwp)WRITE(numout,*)'sbc_fwb : year = ',iyear  , ' freshwater budget correction = ', empold 
    109                IF(lwp)WRITE(numout,*)'          year = ',iyear-1, ' freshwater budget read       = ', a_fwb 
    110                IF(lwp)WRITE(numout,*)'          year = ',iyear-2, ' freshwater budget read       = ', a_fwb_b 
    111             ELSE 
    112                WRITE(ctmp1,*)'sbc_fwb : unable to read the file', clname 
    113                CALL ctl_stop( ctmp1 ) 
    114             ENDIF 
    115          ENDIF 
     100            clname = 'EMPave_old.dat' 
     101            CALL ctlopn( inum, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL',   & 
     102               &           1, numout, .FALSE., 1 ) 
     103 
     104            READ ( inum, "(24X,I8,2ES24.16)" ) iyear, a_fwb_b, a_fwb 
     105            CLOSE( inum ) 
     106            empold = a_fwb                  ! current year freshwater budget correction 
     107            !                               ! estimate from the previous year budget 
     108            IF(lwp)WRITE(numout,*) 
     109            IF(lwp)WRITE(numout,*)'sbc_fwb : year = ',iyear  , ' freshwater budget correction = ', empold 
     110            IF(lwp)WRITE(numout,*)'          year = ',iyear-1, ' freshwater budget read       = ', a_fwb 
     111            IF(lwp)WRITE(numout,*)'          year = ',iyear-2, ' freshwater budget read       = ', a_fwb_b 
     112         ENDIF    
    116113         !  
    117114         ! Update empold if new year start 
     
    134131         ! 
    135132         ! save empold value in a file 
    136          IF( kt == nitend ) THEN 
    137             OPEN( inum, FILE='EMPave.dat' ) 
     133         IF( kt == nitend .AND. lwp ) THEN 
     134            clname = 'EMPav.dat'     
     135            CALL ctlopn( inum, clname, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL',   & 
     136               &           1, numout, .FALSE., 0 ) 
    138137            WRITE(inum, "(24X,I8,2ES24.16)" ) nyear, a_fwb_b, a_fwb 
    139138         ENDIF 
Note: See TracChangeset for help on using the changeset viewer.