Changeset 1617


Ignore:
Timestamp:
2009-08-12T16:59:14+02:00 (11 years ago)
Author:
ctlod
Message:

remove the TKE profile of penetration case (2) for the parameter nn_htau, see ticket: #521

Location:
trunk/NEMO/OPA_SRC/ZDF
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/ZDF/zdftke.F90

    r1601 r1617  
    7676   REAL(wp) ::   rn_bshear= 1.e-20          ! background shear (>0) 
    7777   INTEGER  ::   nn_etau  = 0               ! type of depth penetration of surface tke (=0/1/2) 
    78    INTEGER  ::   nn_htau  = 0               ! type of tke profile of penetration (=0/1/2) 
     78   INTEGER  ::   nn_htau  = 0               ! type of tke profile of penetration (=0/1) 
    7979   REAL(wp) ::   rn_efr   = 1.0_wp          ! fraction of TKE surface value which penetrates in the ocean 
    8080   LOGICAL  ::   ln_lc    = .FALSE.         ! Langmuir cells (LC) as a source term of TKE or not 
     
    664664      IF( nn_mxl  < 0    .OR. nn_mxl  > 3   )   CALL ctl_stop( 'bad flag: nn_mxl is  0, 1 or 2 ' ) 
    665665      IF( nn_pdl  < 0    .OR. nn_pdl  > 1   )   CALL ctl_stop( 'bad flag: nn_pdl is  0 or 1    ' ) 
    666       IF( nn_htau < 0    .OR. nn_htau > 2   )   CALL ctl_stop( 'bad flag: nn_htau is 0, 1 or 2 ' ) 
     666      IF( nn_htau < 0    .OR. nn_htau > 1   )   CALL ctl_stop( 'bad flag: nn_htau is 0 or 1    ' ) 
    667667      IF( rn_lc   < 0.15 .OR. rn_lc   > 0.2 )   CALL ctl_stop( 'bad value: rn_lc must be between 0.15 and 0.2 ' ) 
    668668 
     
    674674         CASE( 0 )                                    ! constant depth penetration (here 10 meters) 
    675675            htau(:,:) = 10.e0 
    676          CASE( 1 )                                    ! F(latitude) : 5m to 40m at high lat. 
     676         CASE( 1 )                                    ! F(latitude) : 0.5m to 30m at high lat. 
    677677            DO jj = 1, jpj 
    678678               DO ji = 1, jpi   
    679                   htau(ji,jj) = MAX( 5., MIN( 40., 60.*ABS( SIN( rpi/180. * gphit(ji,jj) ) ) )   ) 
    680                END DO 
    681             END DO 
    682          CASE( 2 )                                    ! F(latitude) : 5m to 60m at high lat. 
    683             DO jj = 1, jpj 
    684                DO ji = 1, jpi 
    685                   htau(ji,jj) =  MAX( 5.,6./4.* MIN( 40., 60.*ABS( SIN( rpi/180. * gphit(ji,jj) ) ) )   ) 
     679                  htau(ji,jj) = MAX( 0.5, 3./4. * MIN( 40., 60.*ABS( SIN( rpi/180. * gphit(ji,jj) ) ) )   ) 
    686680               END DO 
    687681            END DO 
  • trunk/NEMO/OPA_SRC/ZDF/zdftke_old.F90

    r1601 r1617  
    8181   REAL(wp) ::   rn_cri   = 2._wp / 9._wp   ! critic Richardson number 
    8282   INTEGER  ::   nn_etau  = 0               ! type of depth penetration of surface tke (=0/1/2) 
    83    INTEGER  ::   nn_htau  = 0               ! type of tke profile of penetration (=0/1/2) 
     83   INTEGER  ::   nn_htau  = 0               ! type of tke profile of penetration (=0/1) 
    8484   REAL(wp) ::   rn_lmin0 = 0.4_wp          ! surface  min value of mixing length 
    8585   REAL(wp) ::   rn_lmin  = 0.1_wp          ! interior min value of mixing length 
     
    555555            END DO 
    556556         CASE( 1 )                                    ! meridional profile  1 
    557             DO jj = 2, jpjm1                          ! ( 5m in the tropics to a maximum of 40 m at high lat.) 
    558                DO ji = fs_2, fs_jpim1   ! vector opt. 
    559                   zhtau(ji,jj) = MAX( 5., MIN( 40., 60.*ABS( SIN( rpi/180. * gphit(ji,jj) ) ) )   ) 
    560                END DO 
    561             END DO 
    562          CASE( 2 )                                    ! meridional profile 2 
    563             DO jj = 2, jpjm1                          ! ( 5m in the tropics to a maximum of 60 m at high lat.) 
    564                DO ji = fs_2, fs_jpim1   ! vector opt. 
    565                   zhtau(ji,jj) =  MAX( 5.,6./4.* MIN( 40., 60.*ABS( SIN( rpi/180. * gphit(ji,jj) ) ) )   ) 
     557            DO jj = 2, jpjm1                          ! ( 0.5m in the tropics to a maximum of 30 m at high lat.) 
     558               DO ji = fs_2, fs_jpim1   ! vector opt. 
     559                  zhtau(ji,jj) = MAX( 0.5, 3./4. * MIN( 40., 60.*ABS( SIN( rpi/180. * gphit(ji,jj) ) ) )   ) 
    566560               END DO 
    567561            END DO 
     
    781775      IF( nn_pdl  < 0    .OR. nn_pdl  > 1   )   CALL ctl_stop( 'bad flag: nn_pdl is  0 or 1    ' ) 
    782776      IF( nn_ave  < 0    .OR. nn_ave  > 1   )   CALL ctl_stop( 'bad flag: nn_ave is  0 or 1    ' ) 
    783       IF( nn_htau < 0    .OR. nn_htau > 2   )   CALL ctl_stop( 'bad flag: nn_htau is 0, 1 or 2 ' ) 
     777      IF( nn_htau < 0    .OR. nn_htau > 1   )   CALL ctl_stop( 'bad flag: nn_htau is 0 or 1    ' ) 
    784778      IF( rn_lc   < 0.15 .OR. rn_lc   > 0.2 )   CALL ctl_stop( 'bad value: rn_lc must be between 0.15 and 0.2 ' ) 
    785779 
Note: See TracChangeset for help on using the changeset viewer.