Changeset 2764


Ignore:
Timestamp:
2011-05-09T13:20:38+02:00 (10 years ago)
Author:
cetlod
Message:

update the dtadyn.F90 routine, see ticket #823

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90

    r2762 r2764  
    7575   REAL(wp), ALLOCATABLE, SAVE,   DIMENSION(:,:,:) :: empdta     ! E-P at two consecutive times 
    7676   REAL(wp), ALLOCATABLE, SAVE,   DIMENSION(:,:,:) :: qsrdta     ! short wave heat flux at two consecutive times 
    77    REAL(wp), ALLOCATABLE, SAVE,   DIMENSION(:,:,:) :: bblxdta    ! frequency of bbl in the x direction at 2 consecutive times  
    78    REAL(wp), ALLOCATABLE, SAVE,   DIMENSION(:,:,:) :: bblydta    ! frequency of bbl in the y direction at 2 consecutive times  
     77   REAL(wp), ALLOCATABLE, SAVE,   DIMENSION(:,:,:) :: bblxdta    ! bbl diffusive coef. in the x direction at 2 consecutive times  
     78   REAL(wp), ALLOCATABLE, SAVE,   DIMENSION(:,:,:) :: bblydta    ! bbl diffusive coef. in the y direction at 2 consecutive times  
    7979   LOGICAL :: l_offbbl 
    8080#if defined key_ldfslp && ! defined key_c1d 
     
    360360#if defined key_trabbl 
    361361      IF( .NOT. lk_c1d .AND. nn_bbl_ldf == 1 ) THEN 
    362          IF( iom_varid( numfl_u, 'sobblcox', ldstop = .FALSE. ) > 0  .AND. & 
    363          &   iom_varid( numfl_v, 'sobblcoy', ldstop = .FALSE. ) > 0 ) THEN 
    364              CALL iom_get( numfl_u, jpdom_data, 'sobblcox', zbblx(:,:), jkenr ) 
    365              CALL iom_get( numfl_v, jpdom_data, 'sobblcoy', zbbly(:,:), jkenr ) 
     362         IF( iom_varid( numfl_u, 'ahu_bbl', ldstop = .FALSE. ) > 0  .AND. & 
     363         &   iom_varid( numfl_v, 'ahv_bbl', ldstop = .FALSE. ) > 0 ) THEN 
     364             CALL iom_get( numfl_u, jpdom_data, 'ahu_bbl', zbblx(:,:), jkenr ) 
     365             CALL iom_get( numfl_v, jpdom_data, 'ahv_bbl', zbbly(:,:), jkenr ) 
    366366             l_offbbl = .TRUE. 
    367367         ENDIF 
     
    440440      ! 
    441441      wspddta(:,:,2)  = zwspd(:,:) * tmask(:,:,1) 
    442       frlddta(:,:,2)  = MIN( 1., zice(:,:) ) * tmask(:,:,1) 
    443       empdta (:,:,2)  = zemp(:,:) * tmask(:,:,1) 
    444       qsrdta (:,:,2)  = zqsr(:,:) * tmask(:,:,1) 
    445       hmlddta(:,:,2)  = zmld(:,:) * tmask(:,:,1) 
     442      frlddta(:,:,2)  = zice (:,:) * tmask(:,:,1) 
     443      empdta (:,:,2)  = zemp (:,:) * tmask(:,:,1) 
     444      qsrdta (:,:,2)  = zqsr (:,:) * tmask(:,:,1) 
     445      hmlddta(:,:,2)  = zmld (:,:) * tmask(:,:,1) 
    446446 
    447447#if defined key_trabbl 
    448448      IF( l_offbbl ) THEN  
    449          bblxdta(:,:,2) = MAX( 0., zbblx(:,:) ) 
    450          bblydta(:,:,2) = MAX( 0., zbbly(:,:) ) 
    451          WHERE( bblxdta(:,:,2) > 2. ) bblxdta(:,:,2) = 0. 
    452          WHERE( bblydta(:,:,2) > 2. ) bblydta(:,:,2) = 0. 
     449         bblxdta(:,:,2) = zbblx(:,:)  * umask(:,:,1) 
     450         bblydta(:,:,2) = zbbly(:,:)  * vmask(:,:,1) 
    453451      ENDIF 
    454452#endif 
Note: See TracChangeset for help on using the changeset viewer.