Changeset 12804


Ignore:
Timestamp:
2020-04-23T11:06:57+02:00 (6 months ago)
Author:
dancopsey
Message:

Merged in NEMO_4.0.1_GO6_mixing branch

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/r4.0-HEAD_r12713_dan_test_clems_branch/src/OCE/ZDF/zdftke.F90

    r12801 r12804  
    710710      IF( nn_mxl  < 0   .OR.  nn_mxl  > 3 )   CALL ctl_stop( 'bad flag: nn_mxl is  0, 1 or 2 ' ) 
    711711      IF( nn_pdl  < 0   .OR.  nn_pdl  > 1 )   CALL ctl_stop( 'bad flag: nn_pdl is  0 or 1    ' ) 
    712       IF( nn_htau < 0   .OR.  nn_htau > 1 )   CALL ctl_stop( 'bad flag: nn_htau is 0, 1 or 2 ' ) 
     712      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 ' ) 
    713713      IF( nn_etau == 3 .AND. .NOT. ln_cpl )   CALL ctl_stop( 'nn_etau == 3 : HF taum only known in coupled mode' ) 
    714714      ! 
     
    728728         CASE( 1 )                                 ! F(latitude) : 0.5m to 30m poleward of 40 degrees 
    729729            htau(:,:) = MAX(  0.5_wp, MIN( 30._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(:,:) ) ) )   )             
     730         CASE( 4 )                                 ! F(latitude) : 0.5m to 10m/30m poleward of 13/40 degrees north/south 
     731            DO jj = 2, jpjm1 
     732               DO ji = fs_2, fs_jpim1   ! vector opt. 
     733                  IF( gphit(ji,jj) <= 0._wp ) THEN 
     734                     htau(ji,jj) = MAX(  0.5_wp, MIN( 30._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) )   ) 
     735                  ELSE 
     736                     htau(ji,jj) = MAX(  0.5_wp, MIN( 10._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) )   ) 
     737                  ENDIF 
     738               END DO 
     739            END DO 
    730740         END SELECT 
    731741      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.