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 7646 for trunk/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90 – NEMO

Ignore:
Timestamp:
2017-02-06T10:25:03+01:00 (7 years ago)
Author:
timgraham
Message:

Merge of dev_merge_2016 into trunk. UPDATE TO ARCHFILES NEEDED for XIOS2.
LIM_SRC_s/limrhg.F90 to follow in next commit due to change of kind (I'm unable to do it in this commit).
Merged using the following steps:

1) svn merge --reintegrate svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk .
2) Resolve minor conflicts in sette.sh and namelist_cfg for ORCA2LIM3 (due to a change in trunk after branch was created)
3) svn commit
4) svn switch svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk
5) svn merge svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/2016/dev_merge_2016 .
6) At this stage I checked out a clean copy of the branch to compare against what is about to be committed to the trunk.
6) svn commit #Commit code to the trunk

In this commit I have also reverted a change to Fcheck_archfile.sh which was causing problems on the Paris machine.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90

    r6505 r7646  
    204204      !! 
    205205      !!     ln_teos10 : polynomial TEOS-10 equation of state is used for rho(t,s,z). 
    206       !!         Check value: rho = 1028.21993233072 kg/m^3 for z=3000 dbar, ct=3 Celcius, sa=35.5 g/kg 
     206      !!         Check value: rho = 1028.21993233072 kg/m^3 for z=3000 dbar, ct=3 Celsius, sa=35.5 g/kg 
    207207      !! 
    208208      !!     ln_eos80 : polynomial EOS-80 equation of state is used for rho(t,s,z). 
    209       !!         Check value: rho = 1028.35011066567 kg/m^3 for z=3000 dbar, pt=3 Celcius, sp=35.5 psu 
     209      !!         Check value: rho = 1028.35011066567 kg/m^3 for z=3000 dbar, pt=3 Celsius, sp=35.5 psu 
    210210      !! 
    211211      !!     ln_seos : simplified equation of state 
     
    221221      !!                TEOS-10 Manual, 2010 
    222222      !!---------------------------------------------------------------------- 
    223       REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) ::   pts   ! 1 : potential temperature  [Celcius] 
     223      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) ::   pts   ! 1 : potential temperature  [Celsius] 
    224224      !                                                               ! 2 : salinity               [psu] 
    225225      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(  out) ::   prd   ! in situ density            [-] 
     
    316316      !! 
    317317      !!---------------------------------------------------------------------- 
    318       REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) ::   pts    ! 1 : potential temperature  [Celcius] 
     318      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) ::   pts    ! 1 : potential temperature  [Celsius] 
    319319      !                                                                ! 2 : salinity               [psu] 
    320320      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(  out) ::   prd    ! in situ density            [-] 
     
    481481      !! 
    482482      !!---------------------------------------------------------------------- 
    483       REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in   ) ::   pts   ! 1 : potential temperature  [Celcius] 
     483      REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in   ) ::   pts   ! 1 : potential temperature  [Celsius] 
    484484      !                                                           ! 2 : salinity               [psu] 
    485485      REAL(wp), DIMENSION(jpi,jpj)     , INTENT(in   ) ::   pdep  ! depth                      [m] 
     
    907907      !! 
    908908      !!---------------------------------------------------------------------- 
    909       REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) ::  pts   ! pot. temperature and salinity   [Celcius,psu] 
    910       REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) ::  pab   ! thermal/haline expansion coef.  [Celcius-1,psu-1] 
     909      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) ::  pts   ! pot. temperature and salinity   [Celsius,psu] 
     910      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) ::  pab   ! thermal/haline expansion coef.  [Celsius-1,psu-1] 
    911911      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(  out) ::  pn2   ! Brunt-Vaisala frequency squared [1/s^2] 
    912912      ! 
     
    944944      !!                 ***  ROUTINE eos_pt_from_ct  *** 
    945945      !! 
    946       !! ** Purpose :   Compute pot.temp. from cons. temp. [Celcius] 
     946      !! ** Purpose :   Compute pot.temp. from cons. temp. [Celsius] 
    947947      !! 
    948948      !! ** Method  :   rational approximation (5/3th order) of TEOS-10 algorithm 
     
    952952      !!                Rational approximation to TEOS10 algorithm (rms error on WOA13 values: 4.0e-5 degC) 
    953953      !!---------------------------------------------------------------------- 
    954       REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) ::   ctmp   ! Cons. Temp [Celcius] 
    955       REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) ::   psal   ! salinity   [psu] 
     954      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) ::   ctmp   ! Cons. Temp   [Celsius] 
     955      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) ::   psal   ! salinity     [psu] 
    956956      ! Leave result array automatic rather than making explicitly allocated 
    957       REAL(wp), DIMENSION(jpi,jpj) ::   ptmp   ! potential temperature [Celcius] 
     957      REAL(wp), DIMENSION(jpi,jpj) ::   ptmp   ! potential temperature [Celsius] 
    958958      ! 
    959959      INTEGER  ::   ji, jj               ! dummy loop indices 
     
    10031003      !!                 ***  ROUTINE eos_fzp  *** 
    10041004      !! 
    1005       !! ** Purpose :   Compute the freezing point temperature [Celcius] 
    1006       !! 
    1007       !! ** Method  :   UNESCO freezing point (ptf) in Celcius is given by 
     1005      !! ** Purpose :   Compute the freezing point temperature [Celsius] 
     1006      !! 
     1007      !! ** Method  :   UNESCO freezing point (ptf) in Celsius is given by 
    10081008      !!       ptf(t,z) = (-.0575+1.710523e-3*sqrt(abs(s))-2.154996e-4*s)*s - 7.53e-4*z 
    10091009      !!       checkvalue: tf=-2.588567 Celsius for s=40psu, z=500m 
     
    10131013      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   )           ::   psal   ! salinity   [psu] 
    10141014      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ), OPTIONAL ::   pdep   ! depth      [m] 
    1015       REAL(wp), DIMENSION(jpi,jpj), INTENT(out  )           ::   ptf    ! freezing temperature [Celcius] 
     1015      REAL(wp), DIMENSION(jpi,jpj), INTENT(out  )           ::   ptf    ! freezing temperature [Celsius] 
    10161016      ! 
    10171017      INTEGER  ::   ji, jj          ! dummy loop indices 
     
    10561056      !!                 ***  ROUTINE eos_fzp  *** 
    10571057      !! 
    1058       !! ** Purpose :   Compute the freezing point temperature [Celcius] 
    1059       !! 
    1060       !! ** Method  :   UNESCO freezing point (ptf) in Celcius is given by 
     1058      !! ** Purpose :   Compute the freezing point temperature [Celsius] 
     1059      !! 
     1060      !! ** Method  :   UNESCO freezing point (ptf) in Celsius is given by 
    10611061      !!       ptf(t,z) = (-.0575+1.710523e-3*sqrt(abs(s))-2.154996e-4*s)*s - 7.53e-4*z 
    10621062      !!       checkvalue: tf=-2.588567 Celsius for s=40psu, z=500m 
     
    10661066      REAL(wp), INTENT(in )           ::   psal         ! salinity   [psu] 
    10671067      REAL(wp), INTENT(in ), OPTIONAL ::   pdep         ! depth      [m] 
    1068       REAL(wp), INTENT(out)           ::   ptf          ! freezing temperature [Celcius] 
     1068      REAL(wp), INTENT(out)           ::   ptf          ! freezing temperature [Celsius] 
    10691069      ! 
    10701070      REAL(wp) :: zs   ! local scalars 
Note: See TracChangeset for help on using the changeset viewer.