Ignore:
Timestamp:
2020-04-09T19:30:08+02:00 (8 months ago)
Author:
techene
Message:

replace h. and gde. in case key_qco is activated - quick and dirty

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DOM/domain.F90

    r12724 r12731  
    149149      CALL dom_msk( ik_top, ik_bot )    ! Masks 
    150150      ! 
    151       ht_0(:,:) = 0._wp  ! Reference ocean thickness 
     151!!st      ht_0(:,:) = 0._wp  ! Reference ocean thickness 
    152152      hu_0(:,:) = 0._wp 
    153153      hv_0(:,:) = 0._wp 
     
    166166 
    167167      ! 
     168#if defined key_qco 
     169      !           !==  initialisation of time varying coordinate  ==!   Quasi-Euerian coordinate case 
     170      ! 
     171      IF( .NOT.l_offline )   CALL dom_qe_init( Kbb, Kmm, Kaa ) 
     172      ! 
     173      IF( ln_linssh )        CALL ctl_stop('STOP','domain: key_qco and ln_linssh = T are incompatible') 
     174      ! 
     175#else 
    168176      !           !==  time varying part of coordinate system  ==! 
    169177      ! 
     
    176184            gde3w(:,:,:)    = gde3w_0(:,:,:)    ! = gdept as the sum of e3t 
    177185         ! 
    178 #if defined key_qco 
    179          !  Quasi-Euerian coordinate : no initialisation of e3. scale factors 
    180 #else 
    181186         DO jt = 1, jpt                         ! vertical scale factors 
    182187            e3t(:,:,:,jt) =  e3t_0(:,:,:) 
     
    189194            e3f(:,:,:)    =  e3f_0(:,:,:) 
    190195         ! 
    191 #endif 
    192          ! 
    193196         DO jt = 1, jpt                         ! water column thickness and its inverse 
    194197            hu(:,:,jt)    =    hu_0(:,:) 
     
    201204      ELSE                       != time varying : initialize before/now/after variables 
    202205         ! 
    203 #if defined key_qco 
    204          IF( .NOT.l_offline )   CALL dom_qe_init( Kbb, Kmm, Kaa ) 
    205 #else 
    206206         IF( .NOT.l_offline )   CALL dom_vvl_init( Kbb, Kmm, Kaa ) 
     207         ! 
     208      ENDIF 
    207209#endif 
    208          ! 
    209       ENDIF 
     210 
    210211      ! 
    211212      IF( lk_c1d         )   CALL cor_c1d       ! 1D configuration: Coriolis set at T-point 
Note: See TracChangeset for help on using the changeset viewer.