File indexing completed on 2018-03-02 18:44:41 UTC
view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
d9332b57eb Jean*0001 #include "ZONAL_FILT_OPTIONS.h"
0002
0003 SUBROUTINE ZONAL_FILT_READPARMS( myThid )
ae4c29e0db Jean*0004
d9332b57eb Jean*0005
0006
ae4c29e0db Jean*0007
0008
d9332b57eb Jean*0009 IMPLICIT NONE
0010
0011
0012 #include "SIZE.h"
0013 #include "EEPARAMS.h"
0014 #include "PARAMS.h"
0015 #include "ZONAL_FILT.h"
0016
0017
0018 INTEGER myThid
0019
0020 #ifdef ALLOW_ZONAL_FILT
ae4c29e0db Jean*0021
0022
0023
0024 CHARACTER*(MAX_LEN_MBUF) msgBuf
0025 INTEGER iUnit
d9332b57eb Jean*0026
ae4c29e0db Jean*0027 NAMELIST /ZONFILT_PARM01/
d9332b57eb Jean*0028 & zonal_filt_uvStar, zonal_filt_TrStagg,
ae4c29e0db Jean*0029 & zonal_filt_lat,
d9332b57eb Jean*0030 & zonal_filt_cospow, zonal_filt_sinpow,
0031 & zonal_filt_mode2dx
0032
ae4c29e0db Jean*0033 IF ( .NOT.useZONAL_FILT ) THEN
0034
0035 _BEGIN_MASTER(myThid)
0036
0037
0038
0039 CALL PACKAGES_UNUSED_MSG( 'useZONAL_FILT', ' ', 'zonfilt' )
0040 _END_MASTER(myThid)
0041 RETURN
0042 ENDIF
d9332b57eb Jean*0043
0044
0045
ae4c29e0db Jean*0046
d9332b57eb Jean*0047
0048 _BEGIN_MASTER(myThid)
0049
0050 WRITE(msgBuf,'(A)')' ZONAL_FILT_READPARMS: opening data.zonfilt'
0051 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
0052 & SQUEEZE_RIGHT , 1)
0053
0054 CALL OPEN_COPY_DATA_FILE(
0055 I 'data.zonfilt', 'ZONAL_FILT_READPARMS',
0056 O iUnit,
0057 I myThid )
0058
0059
ae4c29e0db Jean*0060 zonal_filt_uvStar = .TRUE.
0061 zonal_filt_TrStagg = .TRUE.
d9332b57eb Jean*0062 zonal_filt_lat = 90.
0063 zonal_filt_sinpow = 2
ae4c29e0db Jean*0064 zonal_filt_cospow = 2
d9332b57eb Jean*0065 zonal_filt_mode2dx = 0
0066
0067
0068 READ(UNIT=iUnit,NML=ZONFILT_PARM01)
0069
ae4c29e0db Jean*0070 WRITE(msgBuf,'(A)')
d9332b57eb Jean*0071 & ' ZONAL_FILT_READPARMS: finished reading data.zonfilt'
0072 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
0073 & SQUEEZE_RIGHT , 1)
0074
0075
0076
0077
0078
0079
7a77863887 Mart*0080 #ifdef SINGLE_DISK_IO
d9332b57eb Jean*0081 CLOSE(iUnit)
7a77863887 Mart*0082 #else
0083 CLOSE(iUnit,STATUS='DELETE')
0084 #endif /* SINGLE_DISK_IO */
d9332b57eb Jean*0085 _END_MASTER(myThid)
0086
0087
0088 _BARRIER
0089
0090
0091
0092 IF ( .NOT.zonal_filt_uvStar ) THEN
0093
0094
ae4c29e0db Jean*0095
d9332b57eb Jean*0096
0097
0098 IF ( rigidLid ) THEN
0099 WRITE(msgBuf,'(2A)') 'ZONAL_FILT with rigidLid ',
0100 & 'needs zonal_filt_uvStar=.true.'
0101 CALL PRINT_ERROR( msgBuf , 1)
0102 STOP 'ABNORMAL END: S/R ZONAL_FILT_READPARMS'
0103 ELSEIF ( .NOT.exactConserv ) THEN
0104 WRITE(msgBuf,'(2A)') 'S/R ZONAL_FILT_READPARMS: WARNING <<< ',
0105 & 'applying Filter after SOLVE_FOR_P (zonal_filt_uvStar=FALSE)'
1a18fc3009 Jean*0106 CALL PRINT_MESSAGE(msgBuf, errorMessageUnit, SQUEEZE_RIGHT,1)
d9332b57eb Jean*0107 WRITE(msgBuf,'(2A)') 'S/R ZONAL_FILT_READPARMS: WARNING <<< ',
0108 & 'requires to recompute Eta after ==> turn on exactConserv !'
1a18fc3009 Jean*0109 CALL PRINT_MESSAGE(msgBuf, errorMessageUnit, SQUEEZE_RIGHT,1)
ae4c29e0db Jean*0110 ENDIF
d9332b57eb Jean*0111
ae4c29e0db Jean*0112 ENDIF
d9332b57eb Jean*0113
0114 #endif /* ALLOW_ZONAL_FILT */
ae4c29e0db Jean*0115
d9332b57eb Jean*0116 RETURN
0117 END