Changeset 5862


Ignore:
Timestamp:
2015-11-05T15:03:28+01:00 (5 years ago)
Author:
gm
Message:

#1613: vvl by default: non-vvl: initialize _b,n,a scale factors with _0 arrays

Location:
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90

    r5845 r5862  
    469469         sshn(:,:) = sshn(:,:) - snwice_mass(:,:) * r1_rau0 
    470470         sshb(:,:) = sshb(:,:) - snwice_mass(:,:) * r1_rau0 
     471!!gm  like for lim3 : 
     472#if defined key_vvl             
     473!!gm 
    471474         do jk = 1,jpkm1                     ! adjust initial vertical scale factors 
    472475          e3t_n(:,:,jk) = e3t_0(:,:,jk)*( 1._wp + sshn(:,:)*tmask(:,:,1)/(ht_0(:,:) + 1.0 - tmask(:,:,1)) ) 
     
    500503            gde3w_n(:,:,jk) = gdept_n(:,:,jk  ) - sshn   (:,:) 
    501504         END DO 
     505!!gm 
     506#endif 
     507!!gm 
    502508      ENDIF 
    503509      ! 
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90

    r5845 r5862  
    182182#if defined key_vvl 
    183183   LOGICAL, PUBLIC, PARAMETER ::   lk_vvl = .TRUE.    !: variable grid flag 
     184   LOGICAL, PUBLIC ::   lk_linssh = .FALSE.    !: variable grid flag 
    184185#else 
    185186   LOGICAL, PUBLIC, PARAMETER ::   lk_vvl = .FALSE.   !: variable grid flag 
     187   LOGICAL, PUBLIC ::   lk_linssh = .TRUE.    !: variable grid flag 
    186188#endif 
    187189!!gm 
     
    202204    
    203205   !                                                      !  ref. ! before  !   now   !  after  ! 
    204    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   ht_0            ,    ht_n ,    ht_a   !: t-depth              [m] 
     206   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   ht_0            ,    ht_n             !: t-depth              [m] 
    205207   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   hu_0  ,    hu_b ,    hu_n ,    hu_a   !: u-depth              [m] 
    206208   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   hv_0  ,    hv_b ,    hv_n ,    hv_a   !: u-depth              [m] 
     
    355357         &                      hu_b(jpi,jpj) , hv_b(jpi,jpj) , r1_hu_b(jpi,jpj) , r1_hv_b(jpi,jpj) ,     & 
    356358         &      ht_n(jpi,jpj) , hu_n(jpi,jpj) , hv_n(jpi,jpj) , r1_hu_n(jpi,jpj) , r1_hv_n(jpi,jpj) ,     & 
    357          &      ht_a(jpi,jpj) , hu_a(jpi,jpj) , hv_a(jpi,jpj) , r1_hu_a(jpi,jpj) , r1_hv_a(jpi,jpj) , STAT=ierr(6)  ) 
     359         &                      hu_a(jpi,jpj) , hv_a(jpi,jpj) , r1_hu_a(jpi,jpj) , r1_hv_a(jpi,jpj) , STAT=ierr(6)  ) 
    358360         ! 
    359361         ! 
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

    r5845 r5862  
    8383      ENDIF 
    8484      ! 
    85       !              !==  Reference coordinate system  ==! 
     85      !                       !==  Reference coordinate system  ==! 
    8686      ! 
    8787                     CALL dom_nam               ! read namelist ( namrun, namdom ) 
     
    101101      END DO 
    102102      ! 
    103       !              !==  time varying part of coordinate system  ==! 
    104       ! 
    105       IF( lk_vvl ) THEN                         ! time varying : initialize before/now/after variables 
    106          CALL dom_vvl_init  
    107          ! 
    108       ELSE                                      ! Fix in time : set to the reference one for all 
    109          !    before         !          now          !       after         ! 
    110          gdept_b = gdept_0   ;   gdept_n = gdept_0   !        ---          ! depth of grid-points 
    111          gdepw_b = gdepw_0   ;   gdepw_n = gdepw_0   !        ---          ! 
    112                                  gde3w_n = gde3w_0   !        ---          ! 
     103      IF( lk_vvl ) THEN       !==  time varying coordinate system  ==! 
     104         ! 
     105         CALL dom_vvl_init                      ! set before/now/after variables 
     106         ! 
     107      ELSE                    !==  Fix in time  ==!   set everything to the reference one for all 
     108         ! 
     109         !       before        !          now          !       after         ! 
     110         ;  gdept_b = gdept_0  ;   gdept_n = gdept_0   !        ---          ! depth of grid-points 
     111         ;  gdepw_b = gdepw_0  ;   gdepw_n = gdepw_0   !        ---          ! 
     112         ;                     ;   gde3w_n = gde3w_0   !        ---          ! 
    113113         !                                                                   
    114            e3t_b =   e3t_0   ;     e3t_n =   e3t_0   ;   e3t_a =  e3t_0    ! scale factors 
    115            e3u_b =   e3u_0   ;     e3u_n =   e3u_0   ;   e3u_a =  e3u_0    ! 
    116            e3v_b =   e3v_0   ;     e3v_n =   e3u_0   ;   e3v_a =  e3v_0    ! 
    117                              ;     e3f_n =   e3f_0   !        ---          ! 
    118            e3w_b =   e3w_0   ;     e3w_n =   e3w_0   !        ---          ! 
    119           e3uw_b =  e3uw_0   ;    e3uw_n =  e3uw_0   !        ---          ! 
    120           e3vw_b =  e3vw_0   ;    e3vw_n =  e3vw_0   !        ---          ! 
    121          ! 
    122          !                                            !  
     114         ;    e3t_b =   e3t_0  ;     e3t_n =   e3t_0   ;   e3t_a =  e3t_0    ! scale factors 
     115         ;    e3u_b =   e3u_0  ;     e3u_n =   e3u_0   ;   e3u_a =  e3u_0    ! 
     116         ;    e3v_b =   e3v_0  ;     e3v_n =   e3v_0   ;   e3v_a =  e3v_0    ! 
     117         ;                     ;     e3f_n =   e3f_0   !        ---          ! 
     118         ;    e3w_b =   e3w_0  ;     e3w_n =   e3w_0   !        ---          ! 
     119         ;   e3uw_b =  e3uw_0  ;    e3uw_n =  e3uw_0   !        ---          ! 
     120         ;   e3vw_b =  e3vw_0  ;    e3vw_n =  e3vw_0   !        ---          ! 
     121         ! 
    123122         CALL wrk_alloc( jpi,jpj,   z1_hu_0, z1_hv_0 ) 
    124123         ! 
    125          z1_hu_0(:,:) = 1._wp / ( hu_0(:,:) + 1._wp - umask_i(:,:) ) * umask_i(:,:)    ! _i mask due to ISF 
    126          z1_hv_0(:,:) = 1._wp / ( hv_0(:,:) + 1._wp - vmask_i(:,:) ) * vmask_i(:,:) 
    127          ! 
    128          !        before       !         now         !       after         ! 
    129          ;                     ;     ht_n =    hu_0  ;    ht_a =    hu_0   ! water column thickness 
    130          ;     hu_b =    hu_0  ;     hu_n =    hu_0  ;    hu_a =    hu_0   !  
    131          ;     hv_b =    hv_0  ;     hv_n =    hv_0  ;    hv_a =    hv_0   ! 
    132          ;  r1_hu_b = z1_hu_0  ;  r1_hu_n = z1_hu_0  ; r1_hu_a = z1_hu_0   ! inverse of water column thickness 
    133          ;  r1_hv_b = z1_hv_0  ;  r1_hv_n = z1_hv_0  ; r1_hv_a = z1_hv_0   ! 
     124         z1_hu_0(:,:) = umask_i(:,:) / ( hu_0(:,:) + 1._wp - umask_i(:,:) )     ! _i mask due to ISF 
     125         z1_hv_0(:,:) = vmask_i(:,:) / ( hv_0(:,:) + 1._wp - vmask_i(:,:) ) 
     126         ! 
     127         !        before       !          now          !       after         ! 
     128         ;                     ;      ht_n =    ht_0   !                     ! water column thickness 
     129         ;     hu_b =    hu_0  ;      hu_n =    hu_0   ;    hu_a =    hu_0   !  
     130         ;     hv_b =    hv_0  ;      hv_n =    hv_0   ;    hv_a =    hv_0   ! 
     131         ;  r1_hu_b = z1_hu_0  ;   r1_hu_n = z1_hu_0   ; r1_hu_a = z1_hu_0   ! inverse of water column thickness 
     132         ;  r1_hv_b = z1_hv_0  ;   r1_hv_n = z1_hv_0   ; r1_hv_a = z1_hv_0   ! 
    134133         ! 
    135134         CALL wrk_dealloc( jpi,jpj,   z1_hu_0, z1_hv_0 ) 
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r5845 r5862  
    7474      !!                ***  FUNCTION dom_vvl_alloc  *** 
    7575      !!---------------------------------------------------------------------- 
    76       IF( ln_vvl_zstar ) dom_vvl_alloc = 0 
     76      IF( ln_vvl_zstar )   dom_vvl_alloc = 0 
    7777      IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 
    7878         ALLOCATE( tilde_e3t_b(jpi,jpj,jpk)  , tilde_e3t_n(jpi,jpj,jpk) , tilde_e3t_a(jpi,jpj,jpk) ,   & 
     
    539539      !                                        ! Inverse of the local depth 
    540540!!gm BUG ?  don't understand the use of umask_i here ..... 
    541       r1_hu_a(:,:) = 1._wp / ( hu_a(:,:) + 1._wp - umask_i(:,:) ) * umask_i(:,:) 
    542       r1_hv_a(:,:) = 1._wp / ( hv_a(:,:) + 1._wp - vmask_i(:,:) ) * vmask_i(:,:) 
     541      r1_hu_a(:,:) = umask_i(:,:) / ( hu_a(:,:) + 1._wp - umask_i(:,:) ) 
     542      r1_hv_a(:,:) = vmask_i(:,:) / ( hv_a(:,:) + 1._wp - vmask_i(:,:) ) 
    543543      ! 
    544544      CALL wrk_dealloc( jpi,jpj,       zht, z_scale, zwu, zwv, zhdiv ) 
     
    630630      gdepw_n(:,:,1) = 0.0_wp 
    631631      gde3w_n(:,:,1) = gdept_n(:,:,1) - sshn(:,:) 
    632  
    633632      DO jk = 2, jpk 
    634633         DO jj = 1,jpj 
     
    645644      END DO 
    646645 
    647       ! ---------------------------------------------------------------------------------- 
     646      ! Local depth and Inverse of the local depth of the water 
     647      ! ------------------------------------------------------- 
    648648      hu_n(:,:) = hu_a(:,:)   ;   r1_hu_n(:,:) = r1_hu_a(:,:) 
    649649      hv_n(:,:) = hv_a(:,:)   ;   r1_hv_n(:,:) = r1_hv_a(:,:) 
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r5845 r5862  
    875875      !!                  ***  ROUTINE zgr_zco  *** 
    876876      !! 
    877       !! ** Purpose :   define the z-coordinate system 
     877      !! ** Purpose :   define the reference z-coordinate system 
    878878      !! 
    879879      !! ** Method  :   set 3D coord. arrays to reference 1D array  
     
    907907      !!                      
    908908      !! ** Purpose :   the depth and vertical scale factor in partial step 
    909       !!      z-coordinate case 
     909      !!      reference z-coordinate case 
    910910      !! 
    911911      !! ** Method  :   Partial steps : computes the 3D vertical scale factors 
     
    959959      IF( nn_timing == 1 )  CALL timing_start('zgr_zps') 
    960960      ! 
    961       CALL wrk_alloc( jpi, jpj, jpk, zprt ) 
     961      CALL wrk_alloc( jpi,jpj,jpk,  zprt ) 
    962962      ! 
    963963      IF(lwp) WRITE(numout,*) 
     
    11971197      END IF 
    11981198      ! 
    1199       CALL wrk_dealloc( jpi, jpj, jpk, zprt ) 
     1199      CALL wrk_dealloc( jpi,jpj,jpk,  zprt ) 
    12001200      ! 
    12011201      IF( nn_timing == 1 )  CALL timing_stop('zgr_zps') 
     
    12341234      IF( nn_timing == 1 )  CALL timing_start('zgr_isf') 
    12351235      ! 
    1236       CALL wrk_alloc( jpi, jpj, zbathy, zmask, zrisfdep) 
    1237       CALL wrk_alloc( jpi, jpj, zmisfdep, zmbathy ) 
     1236      CALL wrk_alloc( jpi,jpj,  zbathy, zmask, zrisfdep) 
     1237      CALL wrk_alloc( jpi,jpj,  zmisfdep, zmbathy ) 
    12381238 
    12391239 
     
    12511251         WHERE( 0._wp < risfdep(:,:) .AND. risfdep(:,:) >= zdepth )   misfdep(:,:) = jk+1  
    12521252      END DO  
    1253       WHERE (risfdep(:,:) <= e3t_1d(1) .AND. risfdep(:,:) .GT. 0._wp) 
    1254          risfdep(:,:) = 0. ; misfdep(:,:) = 1 
     1253      WHERE (risfdep(:,:) <= e3t_1d(1) .AND. risfdep(:,:) > 0._wp) 
     1254         risfdep(:,:) = 0.   ;  misfdep(:,:) = 1 
    12551255      END WHERE 
    12561256  
     
    12591259! run the bathy check 10 times to be sure all the modif in the bathy or iceshelf draft are compatible together 
    12601260      DO jl = 1, 10      
    1261          WHERE (bathy(:,:) .EQ. risfdep(:,:) ) 
     1261         WHERE (bathy(:,:) == risfdep(:,:) ) 
    12621262            misfdep(:,:) = 0 ; risfdep(:,:) = 0._wp 
    12631263            mbathy (:,:) = 0 ; bathy  (:,:) = 0._wp 
     
    12661266            misfdep(:,:) = 0; risfdep(:,:) = 0._wp  
    12671267            mbathy (:,:) = 0; bathy  (:,:) = 0._wp 
    1268          ENDWHERE 
     1268         END WHERE 
    12691269         IF( lk_mpp ) THEN 
    12701270            zbathy(:,:) = FLOAT( misfdep(:,:) ) 
     
    13111311               ! find the minimum change option: 
    13121312               ! test bathy 
    1313                IF (risfdep(ji,jj) .GT. 1) THEN 
     1313               IF (risfdep(ji,jj) > 1) THEN 
    13141314               zbathydiff =ABS(bathy(ji,jj)   - (gdepw_1d(mbathy (ji,jj)+1) & 
    13151315                 &   + MIN( e3zps_min, e3t_1d(mbathy (ji,jj)+1)*e3zps_rat ))) 
     
    13171317                 &   - MIN( e3zps_min, e3t_1d(misfdep(ji,jj)-1)*e3zps_rat ))) 
    13181318  
    1319                   IF (bathy(ji,jj) .GT. risfdep(ji,jj) .AND. mbathy(ji,jj) .LT. misfdep(ji,jj)) THEN 
     1319                  IF (bathy(ji,jj) > risfdep(ji,jj) .AND. mbathy(ji,jj) < misfdep(ji,jj)) THEN 
    13201320                     IF (zbathydiff .LE. zrisfdepdiff) THEN 
    13211321                        bathy(ji,jj) = gdepw_1d(mbathy(ji,jj)) + MIN( e3zps_min, e3t_1d(mbathy(ji,jj)+1)*e3zps_rat ) 
     
    17681768      IF( nn_timing == 1 )  CALL timing_start('zgr_sco') 
    17691769      ! 
    1770       CALL wrk_alloc( jpi, jpj, zenv, ztmp, zmsk, zri, zrj, zhbat , ztmpi1, ztmpi2, ztmpj1, ztmpj2 ) 
     1770      CALL wrk_alloc( jpi,jpj,  zenv, ztmp, zmsk, zri, zrj, zhbat , ztmpi1, ztmpi2, ztmpj1, ztmpj2 ) 
    17711771      ! 
    17721772      REWIND( numnam_ref )              ! Namelist namzgr_sco in reference namelist : Sigma-stretching parameters 
     
    20462046#endif 
    20472047 
     2048!!gm   I don't like that HERE we are supposed to set the reference coordinate (i.e. _0 arrays) 
     2049!!gm   and only that !!!!! 
     2050!!gm   THIS should be removed from here ! 
    20482051      gdept_n(:,:,:) = gdept_0(:,:,:) 
    20492052      gdepw_n(:,:,:) = gdepw_0(:,:,:) 
     
    20562059      e3uw_n (:,:,:) = e3uw_0 (:,:,:) 
    20572060      e3vw_n (:,:,:) = e3vw_0 (:,:,:) 
     2061!! gm end 
    20582062!! 
    20592063      ! HYBRID :  
     
    20632067               IF( scobot(ji,jj) >= gdept_n(ji,jj,jk) )   mbathy(ji,jj) = MAX( 2, jk ) 
    20642068            END DO 
    2065             IF( scobot(ji,jj) == 0._wp               )   mbathy(ji,jj) = 0 
     2069            IF( scobot(ji,jj) == 0._wp                )   mbathy(ji,jj) = 0 
    20662070         END DO 
    20672071      END DO 
     
    21142118         END DO 
    21152119      ENDIF 
    2116  
    2117 !================================================================================ 
    2118 ! check the coordinate makes sense 
    2119 !================================================================================ 
     2120      ! 
     2121      !================================================================================ 
     2122      ! check the coordinate makes sense 
     2123      !================================================================================ 
    21202124      DO ji = 1, jpi 
    21212125         DO jj = 1, jpj 
    2122  
     2126            ! 
    21232127            IF( hbatt(ji,jj) > 0._wp) THEN 
    21242128               DO jk = 1, mbathy(ji,jj) 
     
    21472151                 ENDIF 
    21482152               END DO 
    2149  
     2153               ! 
    21502154               DO jk = 1, mbathy(ji,jj)-1 
    21512155                 ! and check it never exceeds the total depth 
     
    21572161                 ENDIF 
    21582162               END DO 
    2159  
    21602163            ENDIF 
    2161  
    21622164         END DO 
    21632165      END DO 
     
    21692171   END SUBROUTINE zgr_sco 
    21702172 
    2171 !!====================================================================== 
     2173 
    21722174   SUBROUTINE s_sh94() 
    2173  
    21742175      !!---------------------------------------------------------------------- 
    21752176      !!                  ***  ROUTINE s_sh94  *** 
     
    21822183      !! Reference : Song and Haidvogel 1994.  
    21832184      !!---------------------------------------------------------------------- 
    2184       ! 
    21852185      INTEGER  ::   ji, jj, jk           ! dummy loop argument 
    21862186      REAL(wp) ::   zcoeft, zcoefw   ! temporary scalars 
     
    21882188      REAL(wp), POINTER, DIMENSION(:,:,:) :: z_gsigw3, z_gsigt3, z_gsi3w3 
    21892189      REAL(wp), POINTER, DIMENSION(:,:,:) :: z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3            
    2190  
    2191       CALL wrk_alloc( jpi, jpj, jpk, z_gsigw3, z_gsigt3, z_gsi3w3                                      ) 
    2192       CALL wrk_alloc( jpi, jpj, jpk, z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 ) 
     2190      !!---------------------------------------------------------------------- 
     2191 
     2192      CALL wrk_alloc( jpi,jpj,jpk,  z_gsigw3, z_gsigt3, z_gsi3w3                                      ) 
     2193      CALL wrk_alloc( jpi,jpj,jpk,  z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 ) 
    21932194 
    21942195      z_gsigw3  = 0._wp   ;   z_gsigt3  = 0._wp   ;   z_gsi3w3  = 0._wp 
     
    21962197      z_esigtu3 = 0._wp   ;   z_esigtv3 = 0._wp   ;   z_esigtf3 = 0._wp 
    21972198      z_esigwu3 = 0._wp   ;   z_esigwv3 = 0._wp 
    2198  
     2199      ! 
    21992200      DO ji = 1, jpi 
    22002201         DO jj = 1, jpj 
    2201  
     2202            ! 
    22022203            IF( hbatt(ji,jj) > rn_hc ) THEN    !deep water, stretched sigma 
    22032204               DO jk = 1, jpk 
     
    22622263        END DO 
    22632264      END DO 
    2264  
    2265       CALL wrk_dealloc( jpi, jpj, jpk, z_gsigw3, z_gsigt3, z_gsi3w3                                      ) 
    2266       CALL wrk_dealloc( jpi, jpj, jpk, z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 ) 
    2267  
     2265      ! 
     2266      CALL wrk_dealloc( jpi,jpj,jpk,  z_gsigw3, z_gsigt3, z_gsi3w3                                      ) 
     2267      CALL wrk_dealloc( jpi,jpj,jpk,  z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 ) 
     2268      ! 
    22682269   END SUBROUTINE s_sh94 
    22692270 
     2271 
    22702272   SUBROUTINE s_sf12 
    2271  
    22722273      !!---------------------------------------------------------------------- 
    22732274      !!                  ***  ROUTINE s_sf12 ***  
     
    22852286      !! Reference : Siddorn and Furner 2012 (submitted Ocean modelling). 
    22862287      !!---------------------------------------------------------------------- 
    2287       ! 
    22882288      INTEGER  ::   ji, jj, jk           ! dummy loop argument 
    22892289      REAL(wp) ::   zsmth               ! smoothing around critical depth 
     
    22922292      REAL(wp), POINTER, DIMENSION(:,:,:) :: z_gsigw3, z_gsigt3, z_gsi3w3 
    22932293      REAL(wp), POINTER, DIMENSION(:,:,:) :: z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3            
    2294  
     2294      !!---------------------------------------------------------------------- 
    22952295      ! 
    22962296      CALL wrk_alloc( jpi, jpj, jpk, z_gsigw3, z_gsigt3, z_gsi3w3                                      ) 
     
    23852385             e3f_0(ji,jj,jk)=(scosrf(ji,jj)+hbatf(ji,jj))*z_esigtf3(ji,jj,jk) 
    23862386             ! 
    2387              e3w_0(ji,jj,jk)=hbatt(ji,jj)*z_esigw3(ji,jj,jk) 
     2387             e3w_0 (ji,jj,jk)=hbatt(ji,jj)*z_esigw3(ji,jj,jk) 
    23882388             e3uw_0(ji,jj,jk)=hbatu(ji,jj)*z_esigwu3(ji,jj,jk) 
    23892389             e3vw_0(ji,jj,jk)=hbatv(ji,jj)*z_esigwv3(ji,jj,jk) 
     
    23982398      CALL lbc_lnk(e3uw_0,'T',1.) ; CALL lbc_lnk(e3vw_0,'T',1.) 
    23992399      ! 
    2400       !                                               ! ============= 
    2401  
    2402       CALL wrk_dealloc( jpi, jpj, jpk, z_gsigw3, z_gsigt3, z_gsi3w3                                      ) 
    2403       CALL wrk_dealloc( jpi, jpj, jpk, z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 ) 
    2404  
     2400      CALL wrk_dealloc( jpi,jpj,jpk,   z_gsigw3, z_gsigt3, z_gsi3w3                                      ) 
     2401      CALL wrk_dealloc( jpi,jpj,jpk,   z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 ) 
     2402      ! 
    24052403   END SUBROUTINE s_sf12 
    24062404 
     2405 
    24072406   SUBROUTINE s_tanh() 
    2408  
    24092407      !!---------------------------------------------------------------------- 
    24102408      !!                  ***  ROUTINE s_tanh***  
Note: See TracChangeset for help on using the changeset viewer.