Changeset 8769


Ignore:
Timestamp:
2017-11-21T16:31:59+01:00 (3 years ago)
Author:
flavoni
Message:

add control ISNAN function, and fix in domain.F90 output domain file name

Location:
branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/OPA_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

    r8568 r8769  
    592592      !                       ! ============================= ! 
    593593      !          
    594       clnam = 'domcfg_out'  ! filename (configuration information) 
     594      clnam = cn_domcfg_out  ! filename (configuration information) 
    595595      CALL iom_open( TRIM(clnam), inum, ldwrt = .TRUE., kiolib = jprstlib ) 
    596596       
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/OPA_SRC/stpctl.F90

    r8568 r8769  
    9191      ENDIF 
    9292      ! 
    93       IF ( zmax(1) > 10._wp .OR.   &                     ! too large sea surface height ( > 10 m) 
    94          & zmax(2) > 10._wp .OR.   &                     ! too large velocity ( > 10 m/s) 
    95          & zmax(3) >  0._wp ) THEN                       ! negative sea surface salinity 
     93      IF (  zmax(1) >  10._wp .OR.   &                     ! too large sea surface height ( > 10 m) 
     94         &  zmax(2) >  10._wp .OR.   &                     ! too large velocity ( > 10 m/s) 
     95         &  zmax(3) >=  0._wp .OR.   &                     ! negative or zero sea surface salinity 
     96         &  ISNAN( zmax(1) + zmax(2) + zmax(3) )  ) THEN   ! NaN encounter in the tests 
    9697         IF( lk_mpp ) THEN 
    9798            CALL mpp_maxloc( ABS(sshn)        , tmask(:,:,1), zzz, iih, ijh ) 
     
    108109         IF(lwp) THEN 
    109110            WRITE(numout,cform_err) 
    110             WRITE(numout,*) ' stpctl: |ssh| > 10 m   or   |U| > 10 m/s   or   SSS < 0' 
     111            WRITE(numout,*) ' stpctl: |ssh| > 10 m   or   |U| > 10 m/s   or   SSS < 0   or   NaN encounter in the tests' 
    111112            WRITE(numout,*) ' ====== ' 
    112113            WRITE(numout,9100) kt,   zmax(1), iih, ijh 
Note: See TracChangeset for help on using the changeset viewer.