Changeset 15659
- Timestamp:
- 2022-01-19T19:36:56+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.4_GO8_paquage_branch/src/OCE/ZDF/zdftke.F90
r14075 r15659 802 802 IF( nn_mxl < 0 .OR. nn_mxl > 3 ) CALL ctl_stop( 'bad flag: nn_mxl is 0, 1 or 2 ' ) 803 803 IF( nn_pdl < 0 .OR. nn_pdl > 1 ) CALL ctl_stop( 'bad flag: nn_pdl is 0 or 1 ' ) 804 IF( nn_htau < 0 .OR. nn_htau > 1 ) CALL ctl_stop( 'bad flag: nn_htau is 0, 1 or 2' )804 IF( ( nn_htau < 0 .OR. nn_htau > 1 ) .AND. nn_htau .NE. 4 ) CALL ctl_stop( 'bad flag: nn_htau is 0, 1 or 4 ' ) 805 805 IF( nn_etau == 3 .AND. .NOT. ln_cpl ) CALL ctl_stop( 'nn_etau == 3 : HF taum only known in coupled mode' ) 806 806 ! … … 820 820 CASE( 1 ) ! F(latitude) : 0.5m to 30m poleward of 40 degrees 821 821 htau(:,:) = MAX( 0.5_wp, MIN( 30._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(:,:) ) ) ) ) 822 CASE( 4 ) ! F(latitude) : 0.5m to 10m/30m poleward of 13/40 degrees north/south 823 DO jj = 2, jpjm1 824 DO ji = fs_2, fs_jpim1 ! vector opt. 825 IF( gphit(ji,jj) <= 0._wp ) THEN 826 htau(ji,jj) = MAX( 0.5_wp, MIN( 30._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) ) ) 827 ELSE 828 htau(ji,jj) = MAX( 0.5_wp, MIN( 10._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) ) ) 829 ENDIF 830 END DO 831 END DO 822 832 END SELECT 823 833 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.