** Warning **
Issuing rollback() due to DESTROY without explicit disconnect() of DBD::mysql::db handle dbname=MITgcm at /usr/local/share/lxr/lib/LXR/Common.pm line 1224.
Last-Modified: Fri, 16 May 2024 05:11:24 GMT
Content-Type: text/html; charset=utf-8
MITgcm/MITgcm/pkg/zonal_filt/zonal_filt_readparms.F
File indexing completed on 2018-03-02 18:44:41 UTC
view on github raw 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