New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 6489 for trunk/NEMOGCM/NEMO/OPA_SRC – NEMO

Ignore:
Timestamp:
2016-04-20T14:26:29+02:00 (8 years ago)
Author:
timgraham
Message:

#1687 fixed. Namelist changes required.

Location:
trunk/NEMOGCM/NEMO/OPA_SRC
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r6416 r6489  
    10061006      IF( srcv(jpr_toce)%laction ) THEN                      ! received by sas in case of opa <-> sas coupling 
    10071007         sst_m(:,:) = frcv(jpr_toce)%z3(:,:,1) 
    1008          IF( srcv(jpr_soce)%laction .AND. ln_useCT ) THEN    ! make sure that sst_m is the potential temperature 
     1008         IF( srcv(jpr_soce)%laction .AND. l_useCT ) THEN    ! make sure that sst_m is the potential temperature 
    10091009            sst_m(:,:) = eos_pt_from_ct( sst_m(:,:), sss_m(:,:) ) 
    10101010         ENDIF 
     
    17631763          
    17641764         IF ( nn_components == jp_iam_opa ) THEN 
    1765             ztmp1(:,:) = tsn(:,:,1,jp_tem)   ! send temperature as it is (potential or conservative) -> use of ln_useCT on the received part 
     1765            ztmp1(:,:) = tsn(:,:,1,jp_tem)   ! send temperature as it is (potential or conservative) -> use of l_useCT on the received part 
    17661766         ELSE 
    17671767            ! we must send the surface potential temperature  
    1768             IF( ln_useCT )  THEN    ;   ztmp1(:,:) = eos_pt_from_ct( tsn(:,:,1,jp_tem), tsn(:,:,1,jp_sal) ) 
     1768            IF( l_useCT )  THEN    ;   ztmp1(:,:) = eos_pt_from_ct( tsn(:,:,1,jp_tem), tsn(:,:,1,jp_sal) ) 
    17691769            ELSE                    ;   ztmp1(:,:) = tsn(:,:,1,jp_tem) 
    17701770            ENDIF 
  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90

    r6140 r6489  
    7070         ssu_m(:,:) = ub(:,:,1) 
    7171         ssv_m(:,:) = vb(:,:,1) 
    72          IF( ln_useCT )  THEN    ;   sst_m(:,:) = eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
     72         IF( l_useCT )  THEN    ;   sst_m(:,:) = eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
    7373         ELSE                    ;   sst_m(:,:) = zts(:,:,jp_tem) 
    7474         ENDIF 
     
    9292            ssu_m(:,:) = zcoef * ub(:,:,1) 
    9393            ssv_m(:,:) = zcoef * vb(:,:,1) 
    94             IF( ln_useCT )  THEN    ;   sst_m(:,:) = zcoef * eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
     94            IF( l_useCT )  THEN    ;   sst_m(:,:) = zcoef * eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
    9595            ELSE                    ;   sst_m(:,:) = zcoef * zts(:,:,jp_tem) 
    9696            ENDIF 
     
    120120         ssu_m(:,:) = ssu_m(:,:) + ub(:,:,1) 
    121121         ssv_m(:,:) = ssv_m(:,:) + vb(:,:,1) 
    122          IF( ln_useCT )  THEN    ;   sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
     122         IF( l_useCT )  THEN    ;   sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
    123123         ELSE                    ;   sst_m(:,:) = sst_m(:,:) + zts(:,:,jp_tem) 
    124124         ENDIF 
     
    241241         ssu_m(:,:) = ub(:,:,1) 
    242242         ssv_m(:,:) = vb(:,:,1) 
    243          IF( ln_useCT )  THEN    ;   sst_m(:,:) = eos_pt_from_ct( tsn(:,:,1,jp_tem), tsn(:,:,1,jp_sal) ) 
     243         IF( l_useCT )  THEN    ;   sst_m(:,:) = eos_pt_from_ct( tsn(:,:,1,jp_tem), tsn(:,:,1,jp_sal) ) 
    244244         ELSE                    ;   sst_m(:,:) = tsn(:,:,1,jp_tem) 
    245245         ENDIF 
  • trunk/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90

    r6140 r6489  
    7575 
    7676   !                               !!** Namelist nameos ** 
    77    INTEGER , PUBLIC ::   nn_eos     ! = 0/1/2 type of eq. of state and Brunt-Vaisala frequ. 
    78    LOGICAL , PUBLIC ::   ln_useCT   ! determine if eos_pt_from_ct is used to compute sst_m 
     77   LOGICAL , PUBLIC ::   ln_TEOS10   ! determine if eos_pt_from_ct is used to compute sst_m 
     78   LOGICAL , PUBLIC ::   ln_EOS80   ! determine if eos_pt_from_ct is used to compute sst_m 
     79   LOGICAL , PUBLIC ::   ln_SEOS   ! determine if eos_pt_from_ct is used to compute sst_m 
     80 
     81   ! Parameters 
     82   LOGICAL , PUBLIC    ::   l_useCT         ! =T in ln_TEOS10=T (i.e. use eos_pt_from_ct to compute sst_m), =F otherwise 
     83   INTEGER , PUBLIC    ::   neos            ! Identifier for equation of state used 
     84 
     85   INTEGER , PARAMETER ::   np_teos10 = -1  ! parameter for using TEOS10 
     86   INTEGER , PARAMETER ::   np_eos80  =  0  ! parameter for using EOS80 
     87   INTEGER , PARAMETER ::   np_seos   = 1   ! parameter for using Simplified Equation of state 
    7988 
    8089   !                               !!!  simplified eos coefficients (default value: Vallis 2006) 
     
    184193      !! ** Purpose :   Compute the in situ density (ratio rho/rau0) from 
    185194      !!       potential temperature and salinity using an equation of state 
    186       !!       defined through the namelist parameter nn_eos. 
     195      !!       selected in the nameos namelist 
    187196      !! 
    188197      !! ** Method  :   prd(t,s,z) = ( rho(t,s,z) - rau0 ) / rau0 
     
    194203      !!                rau0   reference density            kg/m^3 
    195204      !! 
    196       !!     nn_eos = -1 : polynomial TEOS-10 equation of state is used for rho(t,s,z). 
     205      !!     ln_teos10 : polynomial TEOS-10 equation of state is used for rho(t,s,z). 
    197206      !!         Check value: rho = 1028.21993233072 kg/m^3 for z=3000 dbar, ct=3 Celcius, sa=35.5 g/kg 
    198207      !! 
    199       !!     nn_eos =  0 : polynomial EOS-80 equation of state is used for rho(t,s,z). 
     208      !!     ln_eos80 : polynomial EOS-80 equation of state is used for rho(t,s,z). 
    200209      !!         Check value: rho = 1028.35011066567 kg/m^3 for z=3000 dbar, pt=3 Celcius, sp=35.5 psu 
    201210      !! 
    202       !!     nn_eos =  1 : simplified equation of state 
     211      !!     ln_seos : simplified equation of state 
    203212      !!              prd(t,s,z) = ( -a0*(1+lambda/2*(T-T0)+mu*z+nu*(S-S0))*(T-T0) + b0*(S-S0) ) / rau0 
    204213      !!              linear case function of T only: rn_alpha<>0, other coefficients = 0 
     
    224233      IF( nn_timing == 1 )   CALL timing_start('eos-insitu') 
    225234      ! 
    226       SELECT CASE( nn_eos ) 
    227       ! 
    228       CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
     235      SELECT CASE( neos ) 
     236      ! 
     237      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    229238         ! 
    230239         DO jk = 1, jpkm1 
     
    266275         END DO 
    267276         ! 
    268       CASE( 1 )                !==  simplified EOS  ==! 
     277      CASE( np_seos )                !==  simplified EOS  ==! 
    269278         ! 
    270279         DO jk = 1, jpkm1 
     
    300309      !! ** Purpose :   Compute the in situ density (ratio rho/rau0) and the 
    301310      !!      potential volumic mass (Kg/m3) from potential temperature and 
    302       !!      salinity fields using an equation of state defined through the 
    303       !!     namelist parameter nn_eos. 
     311      !!      salinity fields using an equation of state selected in the 
     312      !!     namelist. 
    304313      !! 
    305314      !! ** Action  : - prd  , the in situ density (no units) 
     
    322331      IF( nn_timing == 1 )   CALL timing_start('eos-pot') 
    323332      ! 
    324       SELECT CASE ( nn_eos ) 
    325       ! 
    326       CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
     333      SELECT CASE ( neos ) 
     334      ! 
     335      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    327336         ! 
    328337         ! Stochastic equation of state 
     
    430439         ENDIF 
    431440          
    432       CASE( 1 )                !==  simplified EOS  ==! 
     441      CASE( np_seos )                !==  simplified EOS  ==! 
    433442         ! 
    434443         DO jk = 1, jpkm1 
     
    467476      !! ** Purpose :   Compute the in situ density (ratio rho/rau0) from 
    468477      !!      potential temperature and salinity using an equation of state 
    469       !!      defined through the namelist parameter nn_eos. * 2D field case 
     478      !!      selected in the nameos namelist. * 2D field case 
    470479      !! 
    471480      !! ** Action  : - prd , the in situ density (no units) (unmasked) 
     
    486495      prd(:,:) = 0._wp 
    487496      ! 
    488       SELECT CASE( nn_eos ) 
    489       ! 
    490       CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
     497      SELECT CASE( neos ) 
     498      ! 
     499      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    491500         ! 
    492501         DO jj = 1, jpjm1 
     
    527536         CALL lbc_lnk( prd, 'T', 1. )                    ! Lateral boundary conditions 
    528537         ! 
    529       CASE( 1 )                !==  simplified EOS  ==! 
     538      CASE( np_seos )                !==  simplified EOS  ==! 
    530539         ! 
    531540         DO jj = 1, jpjm1 
     
    576585      IF( nn_timing == 1 )   CALL timing_start('rab_3d') 
    577586      ! 
    578       SELECT CASE ( nn_eos ) 
    579       ! 
    580       CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
     587      SELECT CASE ( neos ) 
     588      ! 
     589      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    581590         ! 
    582591         DO jk = 1, jpkm1 
     
    635644         END DO 
    636645         ! 
    637       CASE( 1 )                  !==  simplified EOS  ==! 
     646      CASE( np_seos )                  !==  simplified EOS  ==! 
    638647         ! 
    639648         DO jk = 1, jpkm1 
     
    657666      CASE DEFAULT 
    658667         IF(lwp) WRITE(numout,cform_err) 
    659          IF(lwp) WRITE(numout,*) '          bad flag value for nn_eos = ', nn_eos 
     668         IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    660669         nstop = nstop + 1 
    661670         ! 
     
    690699      pab(:,:,:) = 0._wp 
    691700      ! 
    692       SELECT CASE ( nn_eos ) 
    693       ! 
    694       CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
     701      SELECT CASE ( neos ) 
     702      ! 
     703      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    695704         ! 
    696705         DO jj = 1, jpjm1 
     
    750759         CALL lbc_lnk( pab(:,:,jp_sal), 'T', 1. )                     
    751760         ! 
    752       CASE( 1 )                  !==  simplified EOS  ==! 
     761      CASE( np_seos )                  !==  simplified EOS  ==! 
    753762         ! 
    754763         DO jj = 1, jpjm1 
     
    773782      CASE DEFAULT 
    774783         IF(lwp) WRITE(numout,cform_err) 
    775          IF(lwp) WRITE(numout,*) '          bad flag value for nn_eos = ', nn_eos 
     784         IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    776785         nstop = nstop + 1 
    777786         ! 
     
    806815      pab(:) = 0._wp 
    807816      ! 
    808       SELECT CASE ( nn_eos ) 
    809       ! 
    810       CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
     817      SELECT CASE ( neos ) 
     818      ! 
     819      CASE( np_teos10, np_eos80 )      !==  polynomial TEOS-10 / EOS-80 ==! 
    811820         ! 
    812821         ! 
     
    859868         ! 
    860869         ! 
    861       CASE( 1 )                  !==  simplified EOS  ==! 
     870      CASE( np_seos )                  !==  simplified EOS  ==! 
    862871         ! 
    863872         zt    = pts(jp_tem) - 10._wp   ! pot. temperature anomaly (t-T0) 
    864873         zs    = pts(jp_sal) - 35._wp   ! abs. salinity anomaly (s-S0) 
    865          zh    = pdep                    ! depth at the partial step level 
     874         zh    = pdep                   ! depth at the partial step level 
    866875         ! 
    867876         zn  = rn_a0 * ( 1._wp + rn_lambda1*zt + rn_mu1*zh ) + rn_nu*zs 
     
    873882      CASE DEFAULT 
    874883         IF(lwp) WRITE(numout,cform_err) 
    875          IF(lwp) WRITE(numout,*) '          bad flag value for nn_eos = ', nn_eos 
     884         IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    876885         nstop = nstop + 1 
    877886         ! 
     
    10091018      !!---------------------------------------------------------------------- 
    10101019      ! 
    1011       SELECT CASE ( nn_eos ) 
    1012       ! 
    1013       CASE ( -1, 1 )                !==  CT,SA (TEOS-10 formulation) ==! 
     1020      SELECT CASE ( neos ) 
     1021      ! 
     1022      CASE ( np_teos10, np_eos80 )                !==  CT,SA (TEOS-10 formulation) ==! 
    10141023         ! 
    10151024         DO jj = 1, jpj 
     
    10241033         IF( PRESENT( pdep ) )   ptf(:,:) = ptf(:,:) - 7.53e-4 * pdep(:,:) 
    10251034         ! 
    1026       CASE ( 0 )                     !==  PT,SP (UNESCO formulation)  ==! 
     1035      CASE ( np_seos )                     !==  PT,SP (UNESCO formulation)  ==! 
    10271036         ! 
    10281037         ptf(:,:) = ( - 0.0575_wp + 1.710523e-3_wp * SQRT( psal(:,:) )   & 
     
    10331042      CASE DEFAULT 
    10341043         IF(lwp) WRITE(numout,cform_err) 
    1035          IF(lwp) WRITE(numout,*) '          bad flag value for nn_eos = ', nn_eos 
     1044         IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    10361045         nstop = nstop + 1 
    10371046         ! 
     
    10591068      !!---------------------------------------------------------------------- 
    10601069      ! 
    1061       SELECT CASE ( nn_eos ) 
    1062       ! 
    1063       CASE ( -1, 1 )                !==  CT,SA (TEOS-10 formulation) ==! 
     1070      SELECT CASE ( neos ) 
     1071      ! 
     1072      CASE ( np_teos10, np_eos80 )                !==  CT,SA (TEOS-10 formulation) ==! 
    10641073         ! 
    10651074         zs  = SQRT( ABS( psal ) * r1_S0 )           ! square root salinity 
     
    10701079         IF( PRESENT( pdep ) )   ptf = ptf - 7.53e-4 * pdep 
    10711080         ! 
    1072       CASE ( 0 )                     !==  PT,SP (UNESCO formulation)  ==! 
     1081      CASE ( np_seos )                     !==  PT,SP (UNESCO formulation)  ==! 
    10731082         ! 
    10741083         ptf = ( - 0.0575_wp + 1.710523e-3_wp * SQRT( psal )   & 
     
    10791088      CASE DEFAULT 
    10801089         IF(lwp) WRITE(numout,cform_err) 
    1081          IF(lwp) WRITE(numout,*) '          bad flag value for nn_eos = ', nn_eos 
     1090         IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    10821091         nstop = nstop + 1 
    10831092         ! 
     
    11181127      IF( nn_timing == 1 )   CALL timing_start('eos_pen') 
    11191128      ! 
    1120       SELECT CASE ( nn_eos ) 
    1121       ! 
    1122       CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
     1129      SELECT CASE ( neos ) 
     1130      ! 
     1131      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    11231132         ! 
    11241133         DO jk = 1, jpkm1 
     
    11831192         END DO 
    11841193         ! 
    1185       CASE( 1 )                !==  Vallis (2006) simplified EOS  ==! 
     1194      CASE( np_seos )                !==  Vallis (2006) simplified EOS  ==! 
    11861195         ! 
    11871196         DO jk = 1, jpkm1 
     
    12051214      CASE DEFAULT 
    12061215         IF(lwp) WRITE(numout,cform_err) 
    1207          IF(lwp) WRITE(numout,*) '          bad flag value for nn_eos = ', nn_eos 
     1216         IF(lwp) WRITE(numout,*) '          bad flag value for neos = ', neos 
    12081217         nstop = nstop + 1 
    12091218         ! 
     
    12241233      !!---------------------------------------------------------------------- 
    12251234      INTEGER  ::   ios   ! local integer 
    1226       !! 
    1227       NAMELIST/nameos/ nn_eos, ln_useCT, rn_a0, rn_b0, rn_lambda1, rn_mu1,   & 
     1235      INTEGER  ::   ioptio   ! local integer 
     1236      !! 
     1237      NAMELIST/nameos/ ln_TEOS10, ln_EOS80, ln_SEOS, rn_a0, rn_b0, rn_lambda1, rn_mu1,   & 
    12281238         &                                             rn_lambda2, rn_mu2, rn_nu 
    12291239      !!---------------------------------------------------------------------- 
     
    12451255         WRITE(numout,*) 'eos_init : equation of state' 
    12461256         WRITE(numout,*) '~~~~~~~~' 
    1247          WRITE(numout,*) '          Namelist nameos : set eos parameters' 
    1248          WRITE(numout,*) '             flag for eq. of state and N^2  nn_eos   = ', nn_eos 
    1249          IF( ln_useCT )   THEN 
    1250             WRITE(numout,*) '             model uses Conservative Temperature' 
    1251             WRITE(numout,*) '             Important: model must be initialized with CT and SA fields' 
    1252          ELSE 
    1253             WRITE(numout,*) '             model does not use Conservative Temperature' 
    1254          ENDIF 
     1257         WRITE(numout,*) '   Namelist nameos : Chosen the Equation Of Seawater (EOS)' 
     1258         WRITE(numout,*) '      TEOS-10 : rho=F(Conservative Temperature, Absolute  Salinity, depth)   ln_TEOS10 = ', ln_TEOS10 
     1259         WRITE(numout,*) '      EOS-80  : rho=F(Potential    Temperature, Practical Salinity, depth)   ln_EOS80  = ', ln_EOS80 
     1260         WRITE(numout,*) '      S-EOS   : rho=F(Conservative Temperature, Absolute  Salinity, depth)   ln_SEOS   = ', ln_SEOS 
    12551261      ENDIF 
    1256       ! 
    1257       SELECT CASE( nn_eos )         ! check option 
    1258       ! 
    1259       CASE( -1 )                       !==  polynomial TEOS-10  ==! 
     1262 
     1263      ! Check options for equation of state & set neos based on logical flags 
     1264      ioptio = 0 
     1265      IF( ln_TEOS10 ) THEN   ;   ioptio = ioptio+1   ;   neos = np_teos10   ;   ENDIF 
     1266      IF( ln_EOS80  ) THEN   ;   ioptio = ioptio+1   ;   neos = np_eos80    ;   ENDIF 
     1267      IF( ln_SEOS   ) THEN   ;   ioptio = ioptio+1   ;   neos = np_seos     ;   ENDIF 
     1268      IF( ioptio /= 1 )   CALL ctl_stop("Exactly one equation of state option must be selected") 
     1269      ! 
     1270      SELECT CASE( neos )         ! check option 
     1271      ! 
     1272      CASE( np_teos10 )                       !==  polynomial TEOS-10  ==! 
    12601273         IF(lwp) WRITE(numout,*) 
    12611274         IF(lwp) WRITE(numout,*) '          use of TEOS-10 equation of state (cons. temp. and abs. salinity)' 
     1275         ! 
     1276         l_useCT = .TRUE.                          ! model temperature is Conservative temperature  
    12621277         ! 
    12631278         rdeltaS = 32._wp 
     
    14461461         BPE002 = 1.7269476440e-04_wp 
    14471462         ! 
    1448       CASE( 0 )                        !==  polynomial EOS-80 formulation  ==! 
     1463      CASE( np_eos80 )                        !==  polynomial EOS-80 formulation  ==! 
    14491464         ! 
    14501465         IF(lwp) WRITE(numout,*) 
    14511466         IF(lwp) WRITE(numout,*) '          use of EOS-80 equation of state (pot. temp. and pract. salinity)' 
    14521467         ! 
     1468         l_useCT = .FALSE.                         ! model temperature is Potential temperature 
    14531469         rdeltaS = 20._wp 
    14541470         r1_S0  = 1._wp/40._wp 
     
    16361652         BPE002 = 5.3661089288e-04_wp 
    16371653         ! 
    1638       CASE( 1 )                        !==  Simplified EOS     ==! 
     1654      CASE( np_seos )                        !==  Simplified EOS     ==! 
    16391655         IF(lwp) THEN 
    16401656            WRITE(numout,*) 
     
    16511667            WRITE(numout,*) '               Caution: rn_beta0=0 incompatible with ddm parameterization ' 
    16521668         ENDIF 
    1653          ! 
    1654       CASE DEFAULT                     !==  ERROR in nn_eos  ==! 
    1655          WRITE(ctmp1,*) '          bad flag value for nn_eos = ', nn_eos 
     1669         l_useCT = .TRUE.          ! Use conservative temperature 
     1670         ! 
     1671      CASE DEFAULT                     !==  ERROR in neos  ==! 
     1672         WRITE(ctmp1,*) '          bad flag value for neos = ', neos, '. You should never see this error' 
    16561673         CALL ctl_stop( ctmp1 ) 
    16571674         ! 
     
    16621679      r1_rcp      = 1._wp / rcp 
    16631680      r1_rau0_rcp = 1._wp / rau0_rcp  
     1681      ! 
     1682      IF(lwp) THEN 
     1683         IF( l_useCT )   THEN 
     1684            WRITE(numout,*) '             model uses Conservative Temperature' 
     1685            WRITE(numout,*) '             Important: model must be initialized with CT and SA fields' 
     1686         ELSE 
     1687            WRITE(numout,*) '             model does not use Conservative Temperature' 
     1688         ENDIF 
     1689      ENDIF 
    16641690      ! 
    16651691      IF(lwp) WRITE(numout,*) 
  • trunk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90

    r6140 r6489  
    3131   USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)   
    3232 
    33    USE eosbn2, ONLY : nn_eos 
     33   USE eosbn2, ONLY : neos 
    3434 
    3535   IMPLICIT NONE 
     
    175175      !  Compute Ekman depth from wind stress forcing. 
    176176      ! ------------------------------------------------------- 
    177       zflageos = ( 0.5 + SIGN( 0.5, nn_eos - 1. ) ) * rau0 
     177      zflageos = ( 0.5 + SIGN( 0.5, neos - 1. ) ) * rau0 
    178178      DO jj = 2, jpjm1 
    179179            DO ji = fs_2, fs_jpim1 
Note: See TracChangeset for help on using the changeset viewer.