Ignore:
Timestamp:
2015-03-06T12:30:35+01:00 (7 years ago)
Author:
cetlod
Message:

bugfix:define wmask and rn2b for offline configuration, see ticket#1480

File:
1 edited

Legend:

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

    r5120 r5131  
    245245      tsn(:,:,:,jp_sal) = sf_dyn(jf_sal)%fnow(:,:,:)  * tmask(:,:,:)    ! salinity 
    246246      ! 
    247       CALL eos    ( tsn, rhd, rhop, gdept_0(:,:,:) )                                       ! In any case, we need rhop 
     247      ! 
     248      CALL eos    ( tsn, rhd, rhop, gdept_0(:,:,:) ) ! In any case, we need rhop 
     249      CALL eos_rab( tsn, rab_n )       ! now    local thermal/haline expension ratio at T-points 
     250      CALL bn2    ( tsn, rab_n, rn2 ) ! before Brunt-Vaisala frequency need for zdfmxl 
     251 
     252      rn2b(:,:,:) = rn2(:,:,:)         ! need for zdfmxl 
    248253      CALL zdf_mxl( kt )                                                   ! In any case, we need mxl  
    249254      ! 
     
    535540      !!--------------------------------------------------------------------- 
    536541#if defined key_ldfslp && ! defined key_c1d 
     542      CALL eos    ( pts, rhd, rhop, gdept_0(:,:,:) ) 
    537543      CALL eos_rab( pts, rab_n )       ! now local thermal/haline expension ratio at T-points 
    538544      CALL bn2    ( pts, rab_n, rn2  ) ! now    Brunt-Vaisala 
    539         ! Partial steps: before Horizontal DErivative 
    540         ! only gtsu, gtsv, rhd, gru , grv are used  
     545 
     546      ! Partial steps: before Horizontal DErivative 
    541547      IF( ln_zps  .AND. .NOT. ln_isfcav)                            & 
    542548         &            CALL zps_hde    ( kt, jpts, pts, gtsu, gtsv,  &  ! Partial steps: before horizontal gradient 
     
    547553         &                                 gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the first ocean level 
    548554 
    549  
    550  
    551          !                                                            ! of t, s, rd at the bottom ocean level 
     555      rn2b(:,:,:) = rn2(:,:,:)         ! need for zdfmxl 
    552556      CALL zdf_mxl( kt )            ! mixed layer depth 
    553557      CALL ldf_slp( kt, rhd, rn2 )  ! slopes 
Note: See TracChangeset for help on using the changeset viewer.